Ye Zhang

The human brain is a tissue of vast complexity in terms of the cell types it comprises. Conventional approaches to classifying cell types in the human brain at single cell resolution have been limited to exploring relatively few markers and therefore have provided a limited molecular characterization of any given cell type. We used single cell RNA sequencing on 466 cells to capture the cellular complexity of the adult and fetal human brain at a whole transcriptome level. Healthy adult temporal lobe tissue was obtained during surgical procedures where otherwise normal tissue was removed to gain access to deeper hippocampal pathology in patients with medical refractory seizures. We were able to classify individual cells into all of the major neuronal, glial, and vascular cell types in the brain. We were able to divide neurons into individual communities and show that these communities preserve the categorization of interneuron subtypes that is typically observed with the use of classic interneuron markers. We then used single cell RNA sequencing on fetal human cortical neurons to identify genes that are differentially expressed between fetal and adult neurons and those genes that display an expression gradient that reflects the transition between replicating and quiescent fetal neuronal populations. Finally, we observed the expression of major histocompatibility complex type I genes in a subset of adult neurons, but not fetal neurons. The work presented here demonstrates the applicability of single cell RNA sequencing on the study of the adult human brain and constitutes a first step toward a comprehensive cellular atlas of the human brain.

Glioblastoma (GBM) is the most common primary brain cancer in adults and is notoriously difficult to treat because of its diffuse nature. We performed single-cell RNA sequencing (RNA-seq) on 3,589 cells in a cohort of four patients. We obtained cells from the tumor core as well as surrounding peripheral tissue. Our analysis revealed cellular variation in the tumor's genome and transcriptome. We were also able to identify infiltrating neoplastic cells in regions peripheral to the core lesions. Despite the existence of significant heterogeneity among neoplastic cells, we found that infiltrating GBM cells share a consistent gene signature between patients, suggesting a common mechanism of infiltration. Additionally, in investigating the immunological response to the tumors, we found transcriptionally distinct myeloid cell populations residing in the tumor core and the surrounding peritumoral space. Our data provide a detailed dissection of GBM cell types, revealing an abundance of information about tumor formation and migration.

Our previous reports demonstrated an alarming increase in resistance to adamantanes among influenza A(H3N2) viruses isolated in 2001-2005. To continue monitoring drug resistance, we conducted a comprehensive analysis of influenza A(H3N2) and A(H1N1) viruses isolated globally in 2005-2006. The results obtained by pyrosequencing indicate that 96.4% (n=761) of A(H3N2) viruses circulating in the United States were adamantane resistant. Drug resistance has reached 100% among isolates from some Asian countries. Analysis of correlation between the appearance of drug resistance and the evolutionary pathway of the hemagglutinin (HA) gene suggests at least 2 separate introductions of resistance into circulating populations that gave rise to identifiable subclades. It also indicates that resistant A(H3N2) viruses may have emerged in Asia in late 2001. Among A(H1N1) viruses isolated worldwide, resistance reached 15.5% in 2005-2006; in the United States alone, it was 4.0%. Phylogenetic analysis of the HA and M genes indicates that the acquisition of resistance in A(H1N1) viruses can be linked to a specific genetic group and was not a result of reassortment between A(H3N2) and A(H1N1) viruses. The results of the study highlight the necessity of close monitoring of resistance to existing antivirals as wells as the need for new therapeutics.

BACKGROUND: Eucalyptus is a highly diverse genus of the Myrtaceae family and widely planted in the world for timber and pulp production. Tissue culture induced callus has become a common tool for Eucalyptus breeding, however, our knowledge about the genes related to the callus maturation and shoot regeneration is still poor. RESULTS: We set up an experiment to monitor the callus induction and callus development of two Eucalyptus species - E. camaldulensis (high embryogenic potential) and E. grandis x urophylla (low embryogenic potential). Then, we performed transcriptome sequencing for primary callus, mature callus, shoot regeneration stage callus and senescence callus. We identified 707 upregulated and 694 downregulated genes during the maturation process of the two Eucalyptus species and most of them were involved in the signaling pathways like plant hormone and MAPK. Next, we identified 135 and 142 genes that might play important roles during the callus development of E. camaldulensis and E. grandis x urophylla, respectively. Further, we found 15 DEGs shared by these two Eucalyptus species during the callus development, including Eucgr.D00640 (stem-specific protein TSJT1), Eucgr.B00171 (BTB/POZ and TAZ domain-containing protein 1), Eucgr.C00948 (zinc finger CCCH domain-containing protein 20), Eucgr.K01667 (stomatal closure-related actinbinding protein 3), Eucgr.C00663 (glutaredoxin-C10) and Eucgr.C00419 (UPF0481 protein At3g47200). Interestingly, the expression patterns of these genes displayed "N" shape in the samples. Further, we found 51 genes that were dysregulated during the callus development of E. camaldulensis but without changes in E. grandis x urophylla, such as Eucgr.B02127 (GRF1-interacting factor 1), Eucgr.C00947 (transcription factor MYB36), Eucgr.B02752 (laccase-7), Eucgr.B03985 (transcription factor MYB108), Eucgr.D00536 (GDSL esterase/lipase At5g45920) and Eucgr.B02347 (scarecrow-like protein 34). These 51 genes might be associated with the high propagation ability of Eucalyptus and 22 might be induced after the dedifferentiation. Last, we performed WGCNA to identify the co-expressed genes during the callus development of Eucalyptus and qRT-PCR experiment to validate the gene expression patterns. CONCLUSIONS: This is the first time to globally study the gene profiles during the callus development of Eucalyptus. The results will improve our understanding of gene regulation and molecular mechanisms in the callus maturation and shoot regeneration.

Plants produce a vast array of specialized metabolites, many of which are used as pharmaceuticals, flavors, fragrances, and other high-value fine chemicals. However, most of these compounds occur in non-model plants for which genomic sequence information is not yet available. The production of a large amount of nucleotide sequence data using next-generation technologies is now relatively fast and cost-effective, especially when using the latest Roche-454 and Illumina sequencers with enhanced base-calling accuracy. To investigate specialized metabolite biosynthesis in non-model plants we have established a data-mining framework, employing next-generation sequencing and computational algorithms, to construct and analyze the transcriptomes of 75 non-model plants that produce compounds of interest for biotechnological applications. After sequence assembly an extensive annotation approach was applied to assign functional information to over 800,000 putative transcripts. The annotation is based on direct searches against public databases, including RefSeq and InterPro. Gene Ontology (GO), Enzyme Commission (EC) annotations and associated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps are also collected. As a proof-of-concept, the selection of biosynthetic gene candidates associated with six specialized metabolic pathways is described. A web-based BLAST server has been established to allow public access to assembled transcriptome databases for all 75 plant species of the PhytoMetaSyn Project (www.phytometasyn.ca).