Developmental Gene Expression Patterns Research Articles

CBIO-14. SINGLE CELL RNA SEQUENCING ANALYSIS OF HUMAN GLIOBLASTOMA STEM-LIKE CELL CULTURES AND XENOGRAFT TUMORS

Abstract Single cell RNA-seq (scRNA-seq) studies for glioma have yielded critical insight into intratumoral heterogeneity and developmental gene expression patterns for primary gliomas. One key conclusion from these studies is that each tumor represents a complex, yet maligned, neuro-developmental ecosystem, harboring diverse cell types, which presumably contribute to tumor growth and homeostasis in specific ways (e.g., vascular mimicry, immune evasion, recreating NSC niches, neural injury responses, etc.). Here, to better understand experimental models of human glioblastoma (GB), we performed single cell RNA-seq analysis of human GB stem-like cells (GSCs) of distinct tumor subtypes (mesenchymal and proneural) during their in vitro culture in serum-free conditions and also during tumor formation in immunocompromised mice. This analysis revealed surprising differences between in vitro and in vivo grown GSCs. Among our results, we find that in vivo mesenchymal GSCs are capable of transitioning to proneural-like states, while proneural GSCs are capable of transitioning to mesenchymal states. We characterize cycling cells based on expression of and G2/M and S phase makers, estimate RNA velocity, and examine different developmental trajectories arising in vitro and in vivo. We also compare and discuss different analysis pipelines for scRNA-seq data.

Single-cell transcriptomics reveals hidden information in trypanosomatids

Dixenic parasites often encounter environmental extremes during the transition from vector to host. Preadapted transmission stages overcome these challenges to promote parasites' survival and ensure life cycle progression. Recently, Vigneron et al. and Briggs et al. used single-cell transcriptomics to investigate developmental stage specific gene expression patterns during parasite differentiation.

Transcriptome Reveals Allele Contribution to Heterosis in Maize.

Heterosis, which has greatly increased maize yields, is associated with gene expression patterns during key developmental stages that enhance hybrid phenotypes relative to parental phenotypes. Before heterosis can be more effectively used for crop improvement, hybrid maize developmental gene expression patterns must be better understood. Here, six maize hybrids, including the popular hybrid Zhengdan958 (ZC) from China, were studied. Maize hybrids created in-house were generated using an incomplete diallel cross (NCII)-based strategy from four elite inbred parental lines. Differential gene expression (DEG) profiles corresponding to three developmental stages revealed that hybrid partial expression patterns exhibited complementarity of expression of certain parental genes, with parental allelic expression patterns varying both qualitatively and quantitatively in hybrids. Single-parent expression (SPE) and parent-specific expression (PSE) types of qualitative variation were most prevalent, 43.73 and 41.07% of variation, respectively. Meanwhile, negative super-dominance (NSD) and positive super-dominance (PSD) types of quantitative variation were most prevalent, 31.06 and 24.30% of variation, respectively. During the early reproductive growth stage, the gene expression pattern differed markedly from other developmental stage patterns, with allelic expression patterns during seed development skewed toward low-value parental alleles in hybrid seeds exhibiting significant quantitative variation-associated superiority. Comparisons of qualitative gene expression variation rates between ZC and other hybrids revealed proportions of SPE-DEGs (41.36%) in ZC seed DEGs that significantly exceeded the average proportion of SPE-DEGs found in seeds of other hybrids (28.36%). Importantly, quantitative gene expression variation rate comparisons between ZC and hybrids, except for transgressive expression, revealed that the ZC rate exceeded the average rate for other hybrids, highlighting the importance of partial gene expression in heterosis. Moreover, enriched ZC DEGs exhibiting distinct tissue-specific expression patterns belonged to four biological pathways, including photosynthesis, plant hormone signal transduction, biology metabolism and biosynthesis. These results provide valuable technical insights for creating hybrids exhibiting strong heterosis.

Developmental expression patterns of toolkit genes in male accessory gland of Drosophila parallels those of mammalian prostate.

ABSTRACTConservation of genetic toolkits in disparate phyla may help reveal commonalities in organ designs transcending their extreme anatomical disparities. A male accessory sexual organ in mammals, the prostate, for instance, is anatomically disparate from its analogous, phylogenetically distant counterpart – the male accessory gland (MAG) – in insects like Drosophila. It has not been ascertained if the anatomically disparate Drosophila MAG shares developmental parallels with those of the mammalian prostate. Here we show that the development of Drosophila mesoderm-derived MAG entails recruitment of similar genetic toolkits of tubular organs like that seen in endoderm-derived mammalian prostate. For instance, like mammalian prostate, Drosophila MAG morphogenesis is marked by recruitment of fibroblast growth factor receptor (FGFR) – a signalling pathway often seen recruited for tubulogenesis – starting early during its adepithelial genesis. A specialisation of the individual domains of the developing MAG tube, on the other hand, is marked by the expression of a posterior Hox gene transcription factor, Abd-B, while Hh-Dpp signalling marks its growth. Drosophila MAG, therefore, reveals the developmental design of a unitary bud-derived tube that appears to have been co-opted for the development of male accessory sexual organs across distant phylogeny and embryonic lineages.This article has an associated First Person interview with the first author of the paper.

First person – Jaya Kumari

ABSTRACT First Person is a series of interviews with the first authors of a selection of papers published in Biology Open, helping early-career researchers promote themselves alongside their papers. Jaya Kumari is first author on ‘ Developmental expression patterns of toolkit genes in male accessory gland of Drosophila parallels those of mammalian prostate’, published in BiO. Jaya conducted the research described in this article while a graduate student in Pradip Sinha's lab at the Department of Biological Sciences and Bioengineering, Indian Institute of Technology Kanpur, Kanpur, India. Jaya is currently an early-career postdoc in the lab of Pradip Sinha, investigating development and carcinogenesis in Drosophila organs with squamous epithelial cells.

Molecular architecture of the developing mouse brain.

The mammalian brain develops through a complex interplay of spatial cues generated by diffusible morphogens, cell-cell interactions and intrinsic genetic programs that result in probably more than a thousand distinct cell types. A complete understanding of this process requires a systematic characterization of cell states over the entire spatiotemporal range of brain development. The ability of single-cell RNA sequencing and spatial transcriptomics to reveal the molecular heterogeneity of complex tissues has therefore been particularly powerful in the nervous system. Previous studies have explored development in specific brain regions1-8, the whole adult brain9 and even entire embryos10. Here we report a comprehensive single-cell transcriptomic atlas of the embryonic mouse brain between gastrulation and birth. We identified almost eight hundred cellular states that describe a developmental program for the functional elements of the brain and its enclosing membranes, including the early neuroepithelium, region-specific secondary organizers, and both neurogenic and gliogenic progenitors. We also used in situ mRNA sequencing to map the spatial expression patterns of key developmental genes. Integrating the in situ data with our single-cell clusters revealed the precise spatial organization of neural progenitors during the patterning of the nervous system.

Nuclear Organization during Hepatogenesis in Zebrafish Requires Uhrf1.

Acquisition of cellular fate during development is initiated and maintained by well-coordinated patterns of gene expression that are dictated by the epigenetic landscape and genome organization in the nucleus. While the epigenetic marks that mediate developmental gene expression patterns during organogenesis have been well studied, less is known about how epigenetic marks influence nuclear organization during development. This study examines the relationship between nuclear structure, chromatin accessibility, DNA methylation, and gene expression during hepatic outgrowth in zebrafish larvae. We investigate the relationship between these features using mutants that lack DNA methylation. Hepatocyte nuclear morphology was established coincident with hepatocyte differentiation at 80 h post-fertilization (hpf), and nuclear shape and size continued to change until the conclusion of outgrowth and morphogenesis at 120 hpf. Integrating ATAC-Seq analysis with DNA methylation profiling of zebrafish livers at 120 hpf showed that closed and highly methylated chromatin occupies most transposable elements and that open chromatin correlated with gene expression. DNA hypomethylation, due to mutation of genes encoding ubiquitin-like, containing PHD and RING Finger Domains 1 (uhrf1) and DNA methyltransferase (dnmt1), did not block hepatocyte differentiation, but had dramatic effects on nuclear organization. Hepatocytes in uhrf1 mutants have large, deformed nuclei with multiple nucleoli, downregulation of nucleolar genes, and a complete lack of the nuclear lamina. Loss of lamin B2 staining was phenocopied by dnmt1 mutation. Together, these data show that hepatocyte nuclear morphogenesis coincides with organ morphogenesis and outgrowth, and that DNA methylation directs chromatin organization, and, in turn, hepatocyte nuclear shape and size during liver development.

Enhancing Gonadotrope Gene Expression Through Regulatory lncRNAs.

The world of long non-coding RNAs (lncRNAs) has opened up massive new prospects in understanding the regulation of gene expression. Not only are there seemingly almost infinite numbers of lncRNAs in the mammalian cell, but they have highly diverse mechanisms of action. In the nucleus, some are chromatin-associated, transcribed from transcriptional enhancers (eRNAs) and/or direct changes in the epigenetic landscape with profound effects on gene expression. The pituitary gonadotrope is responsible for activation of reproduction through production and secretion of appropriate levels of the gonadotropic hormones. As such, it exemplifies a cell whose function is defined through changes in developmental and temporal patterns of gene expression, including those that are hormonally induced. Roles for diverse distal regulatory elements and eRNAs in gonadotrope biology have only just begun to emerge. Here, we will present an overview of the different kinds of lncRNAs that alter gene expression, and what is known about their roles in regulating some of the key gonadotrope genes. We will also review various screens that have detected differentially expressed pituitary lncRNAs associated with changes in reproductive state and those whose expression is found to play a role in gonadotrope-derived nonfunctioning pituitary adenomas. We hope to shed light on this exciting new field, emphasize the open questions, and encourage research to illuminate the roles of lncRNAs in various endocrine systems.

Genome-wide analysis of the Universal stress protein A gene family in Vitis and expression in response to abiotic stress

Universal Stress Protein A (USPA) plays critical roles in the regulation of growth, development and response to abiotic stress in plants. To date, most research related to the role of USPA in plants has been carried out in herbaceous models such as Arabidopsis, rice and soybean. Here, we used bioinformatics approaches to identify 21 USPA genes in the genome of Vitis vinifera L. Phylogenetic analysis revealed that VvUSPAs could be divided into eight clades. Based on predicted chromosomal locations, we identified 16 pairs of syntenic, orthologous genes between A. thaliana and V. vinifera. Further promoter cis-elements analysis, together with identification of potential microRNA (miRNA) binding sites, suggested that at least some of the VvUSPAs participate in response to phytohormones and abiotic stress. To add support for this, we analyzed the developmental and stress-responsive expression patterns of the homologous USPA genes in the drought-resistant wild Vitis yeshanensis accession ‘Yanshan-1’ and the drought-sensitive Vitis riparia accession ‘He'an’. Most of the USPA genes were upregulated in different degrees in the two genotypes after drought stress and exposure to ethephon (ETH), abscisic acid (ABA) and methyl jasmonate (MeJA). Individual USPA genes showed various tissue-specific expression patterns. Heterologous expression of five selected genes (VvUSPA2, VvUSPA3, VvUSPA11, VvUSPA13 and VvUSPA16) in Escherichia coli (E. coli) enhanced resistance to drought stress. Our study provides a model for mapping gene function in response to abiotic stress and identified three candidate genes, VvUSPA3, VvUSPA11 and VvUSPA16, as regulators of drought response in V. vinifera.

Recapitulating Developmental Condensation and Constructing Self-organised Cartilaginous Tissue for Cartilage Regeneration.

To develop a novel chondrocyte condensation culture strategy recapitulating developmental condensation and construct self-organised cartilaginous tissue for cartilage regeneration. Cell-condensation aggregate (CCA) was generated using the condensation culture method by sequential cell seeding. The chondrification capacities and biocompatibilities of CCA were assessed by comparison with the cell-scaffold complex (CSC), which was constructed by cell-scaffold coculture. Preclinical studies including implantation into nude mice subcutaneously and cartilage defect repair in rabbits were performed. CCA constructed by condensation culture exhibited a morphology of self-organised cartilaginous tissue. Meanwhile, the condensation culture inhibited or abolished expression of HOX genes including HOXC4 and HOXD8, which was partially consistent with developmental HOX gene expression patterns and associated with enhanced regeneration capacities. Compared with CSC, CCA showed a higher capacity for chondrification and regeneration of rabbit cartilage defects. The therapeutic assessments indicate that CCA is an efficient therapeutic tool for cartilage regeneration, providing a new strategy for tissue engineering by mimicking developmental events.

Dynamic regulation of anterior-posterior patterning genes in living Drosophila embryos

Expression of the gap and pair-rule genes plays an essential role in body segmentation during Drosophila embryogenesis.1-5 However, it remains unclear how precise expression patterns of these key developmental genes arise from stochastic transcriptional activation at the single-cell level. Here, I employed genome-editing and live-imaging approaches to comprehensively visualize regulation of the gap and pair-rule genes at the endogenous loci. Quantitative image analysis revealed that the total duration of active transcription (transcription period) is a major determinant of spatial patterning of gene expression in early embryos. The length of the transcription period is determined by the continuity of bursting activities in individual nuclei, with the core expression domain producing more bursts than boundary regions. Each gene exhibits a distinct rate of nascent RNA production during transcriptional bursting, which contributes to gene-to-gene variability in the total output. I also provide evidence for "enhancer interference," wherein a distal weak enhancer interferes with transcriptional activation by a strong proximal enhancer to downregulate the length of the transcription period without changing the transcription rate. Analysis of the endogenous hunchback (hb) locus revealed that the removal of the distal shadow enhancer induces strong ectopic transcriptional activation, which suppresses refinement of the initial broad expression domain into narrower stripe patterns at the anterior part of embryos. This study provides key insights into the link between transcriptional bursting, enhancer-promoter interaction, and spatiotemporal patterning of gene expression during animal development.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species.

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5′-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

Decoding Human Megakaryocyte Development.

Despite our growing understanding of embryonic immune development, rare early megakaryocytes (MKs) remain relatively understudied. Here we used single-cell RNA sequencing of human MKs from embryonic yolk sac (YS) and fetal liver (FL) to characterize the transcriptome, cellular heterogeneity, and developmental trajectories of early megakaryopoiesis. In the YS and FL, we found heterogeneous MK subpopulations with distinct developmental routes and patterns of gene expression that could reflect early functional specialization. Intriguingly, we identified a subpopulation of CD42b+CD14+ MKs invivo that exhibit high expression of genes associated with immune responses and can also be derived from human embryonic stem cells (hESCs) invitro. Furthermore, we identified THBS1 as an early marker for MK-biased embryonic endothelial cells. Overall, we provide important insights and invaluable resources for dissection of the molecular and cellular programs underlying early human megakaryopoiesis.

PCB126 Exposure Revealed Alterations in m6A RNA Modifications in Transcripts Associated With AHR Activation.

Chemical modifications of proteins, DNA, and RNA moieties play critical roles in regulating gene expression. Emerging evidence suggests the RNA modifications (epitranscriptomics) have substantive roles in basic biological processes. One of the most common modifications in mRNA and noncoding RNAs is N6-methyladenosine (m6A). In a subset of mRNAs, m6A sites are preferentially enriched near stop codons, in 3' UTRs, and within exons, suggesting an important role in the regulation of mRNA processing and function including alternative splicing and gene expression. Very little is known about the effect of environmental chemical exposure on m6A modifications. As many of the commonly occurring environmental contaminants alter gene expression profiles and have detrimental effects on physiological processes, it is important to understand the effects of exposure on this important layer of gene regulation. Hence, the objective of this study was to characterize the acute effects of developmental exposure to PCB126, an environmentally relevant dioxin-like PCB, on m6A methylation patterns. We exposed zebrafish embryos to PCB126 for 6 h starting from 72 h post fertilization and profiled m6A RNA using methylated RNA immunoprecipitation followed by sequencing (MeRIP-seq). Our analysis revealed 117 and 217 m6A peaks in the DMSO and PCB126 samples (false discovery rate5%), respectively. The majority of the peaks were preferentially located around the 3' UTR and stop codons. Statistical analysis revealed 15 m6A marked transcripts to be differentially methylated by PCB126 exposure. These include transcripts that are known to be activated by AHR agonists (eg, ahrra, tiparp, nfe2l2b) as well as others that are important for normal development (vgf, cebpd, sned1). These results suggest that environmental chemicals such as dioxin-like PCBs could affect developmental gene expression patterns by altering m6A levels. Further studies are necessary to understand the functional consequences of exposure-associated alterations in m6A levels.

The Dynamics, Causes, and Impacts of Mammalian Evolutionary Rates Revealed by the Analyses of Capybara Draft Genome Sequences.

Capybara (Hydrochoerus hydrochaeri) is the largest species among the extant rodents. The draft genome of capybara was sequenced with the estimated genome size of 2.6 Gb. Although capybara is about 60 times larger than guinea pig, comparative analyses revealed that the neutral evolutionary rates of the two species were not substantially different. However, analyses of 39 mammalian genomes revealed very heterogeneous evolutionary rates. The highest evolutionary rate, 8.5 times higher than the human rate, was found in the Cricetidae-Muridae common ancestor after the divergence of Spalacidae. Muridae, the family with the highest number of species among mammals, emerged after the rate acceleration. Factors responsible for the evolutionary rate heterogeneity were investigated through correlations between the evolutionary rate and longevity, gestation length, litter frequency, litter size, body weight, generation interval, age at maturity, and taxonomic order. The regression analysis of these factors showed that the model with three factors (taxonomic order, generation interval, and litter size) had the highest predictive power (R2 = 0.74). These three factors determine the number of meiosis per unit time. We also conducted transcriptome analysis and found that the evolutionary rate dynamics affects the evolution of gene expression patterns.

Developmental expression patterns and correlation analysis of TLR4 and its downstream genes in the intestinal and immune tissues of Meishan pigs

ABSTRACT To understand the developmental expression patterns of key genes in the Toll-like receptor 4 ( TLR4 ) pathway and their regulatory characteristics in the immune response in pigs, we examined TLR4 and its downstream genes expression levels in the intestinal and immune tissues of Meishan pigs. The genes were expressed in all examined tissues at the different developmental stages. TLR4 expression was higher in spleen and lower in other tissues. Spleen and lymph TLR4 expression was significantly lower in [...]

Selection of reference genes for the study of relative gene expression during development of Litopenaeus vannamei larvae

Pacific whiteleg shrimp (Litopenaeus vannamei) production is one of the most economically important aquaculture industries around the world. Gene expression studies during larval development are a key tool to elucidate biological aspects during metamorphosis and the effects of pathological conditions; the adequate quantification of gene expression can provide molecular markers and basic knowledge to optimize culture conditions, including the development of strategies for the prevention and control of diseases. In this context, the selection and validation of reference genes is a critical process to avoid bias due to inaccurate measurements and the consequent misinterpretation of biological processes. In this study, nine candidate reference genes (Clathrin, Cyt-c, SubF0, EF1α, β-Actin, GAPDH, TBP, AK and PK) were selected to test their expression stability during larval development of L. vannamei using geNorm, NormFinder, BestKeeper and the integrated tool RefFinder algorithms. Based on our analysis the most stable gene was SubF0, followed by GAPDH and EF1α. Relative expression of developmental genes (Twist, Mef-2 and Ubx) was quantified using the three most stable reference genes (individually and combined), taking special attention to the expected expression and biological function of these genes. The expression of Mef-2 was the most sensitive to the reference gene. However, the combination of the three most stable reference genes provided consistent results according to the expected expression patterns of developmental genes.

Immunostaining of the Embryonic and Larval Drosophila Brain.

Immunostaining is used to visualize the spatiotemporal expression pattern of developmental control genes that regulate the genesis and specification of the embryonic and larval brain of Drosophila. It is also used to visualize the effects of targeted misexpression or inactivation of disease-related genes. Immunostaining uses specific antibodies to mark expressed proteins and allows their localization to be traced. This method reveals insights into gene regulation, cell type specification, neuron and glial differentiation, axonal and synaptic scaffolding and posttranslational protein modifications underlying the patterning and specification of the maturing brain. Depending on the targeted protein, it is possible to visualize a multitude of regions of the Drosophila brain, such as small groups of neurons or glia, defined subcomponents of the brain's axon scaffold, or pre- and postsynaptic structures of neurons. Thus, antibody probes that recognize defined tissues, cells, or subcellular structures like axons or synaptic terminals can be used as markers to identify and analyze phenotypes in embryos and larvae. Several antibodies, combined with different labels can be used concurrently to examine protein colocalization. This protocol spans over 3-4days.

Accurate genome-wide predictions of spatio-temporal gene expression during embryonic development.

Comprehensive information on the timing and location of gene expression is fundamental to our understanding of embryonic development and tissue formation. While high-throughput in situ hybridization projects provide invaluable information about developmental gene expression patterns for model organisms like Drosophila, the output of these experiments is primarily qualitative, and a high proportion of protein coding genes and most non-coding genes lack any annotation. Accurate data-centric predictions of spatio-temporal gene expression will therefore complement current in situ hybridization efforts. Here, we applied a machine learning approach by training models on all public gene expression and chromatin data, even from whole-organism experiments, to provide genome-wide, quantitative spatio-temporal predictions for all genes. We developed structured in silico nano-dissection, a computational approach that predicts gene expression in >200 tissue-developmental stages. The algorithm integrates expression signals from a compendium of 6,378 genome-wide expression and chromatin profiling experiments in a cell lineage-aware fashion. We systematically evaluated our performance via cross-validation and experimentally confirmed 22 new predictions for four different embryonic tissues. The model also predicts complex, multi-tissue expression and developmental regulation with high accuracy. We further show the potential of applying these genome-wide predictions to extract tissue specificity signals from non-tissue-dissected experiments, and to prioritize tissues and stages for disease modeling. This resource, together with the exploratory tools are freely available at our webserver http://find.princeton.edu, which provides a valuable tool for a range of applications, from predicting spatio-temporal expression patterns to recognizing tissue signatures from differential gene expression profiles.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species.

Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis. By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.