Developmental Gene Expression Research Articles

Nucleosomal asymmetry: a novel mechanism to regulate nucleosome function.

Nucleosomes constitute the fundamental building blocks of chromatin. They are comprised of DNA wrapped around a histone octamer formed of two copies each of the four core histones H2A, H2B, H3, and H4. Nucleosomal histones undergo a plethora of posttranslational modifications that regulate gene expression and other chromatin-templated processes by altering chromatin structure or by recruiting effector proteins. Given their symmetric arrangement, the sister histones within a nucleosome have commonly been considered to be equivalent and to carry the same modifications. However, it is now clear that nucleosomes can exhibit asymmetry, combining differentially modified sister histones or different variants of the same histone within a single nucleosome. Enabled by the development of novel tools that allow generating asymmetrically modified nucleosomes, recent biochemical and cell-based studies have begun to shed light on the origins and functional consequences of nucleosomal asymmetry. These studies indicate that nucleosomal asymmetry represents a novel regulatory mechanism in the establishment and functional readout of chromatin states. Asymmetry expands the combinatorial space available for setting up complex sets of histone marks at individual nucleosomes, regulating multivalent interactions with histone modifiers and readers. The resulting functional consequences of asymmetry regulate transcription, poising of developmental gene expression by bivalent chromatin, and the mechanisms by which oncohistones deregulate chromatin states in cancer. Here, we review recent progress and current challenges in uncovering the mechanisms and biological functions of nucleosomal asymmetry.

Biologically informed NeuralODEs for genome-wide regulatory dynamics

BackgroundGene regulatory network (GRN) models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such gene regulatory ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the underlying GRN governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impede either scalability, explainability, or both.ResultsWe developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that overcomes limitations of other methods by flexibly incorporating prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of GRN ODEs. We tested the accuracy of PHOENIX in a series of in silico experiments, benchmarking it against several currently used tools. We demonstrated PHOENIX’s flexibility by modeling regulation of oscillating expression profiles obtained from synchronized yeast cells. We also assessed the scalability of PHOENIX by modeling genome-scale GRNs for breast cancer samples ordered in pseudotime and for B cells treated with Rituximab.ConclusionsPHOENIX uses a combination of user-defined prior knowledge and functional forms from systems biology to encode biological “first principles” as soft constraints on the GRN allowing us to predict subsequent gene expression patterns in a biologically explainable manner.

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

Insulin-like growth factor-1 effects on kidney development in preterm piglets.

Preterm birth disrupts fetal kidney development, potentially leading to postnatal acute kidney injury. Preterm infants are deficient in insulin-like growth factor 1 (IGF-1), a growth factor that stimulates organ development. By utilizing a preterm pig model, this study investigated whether IGF-1 supplementation enhances preterm kidney maturation. Cesarean-delivered preterm pigs were treated systemically IGF-1 or vehicle control for 5, 9 or 19 days after birth. Blood, urine, and kidney tissue were collected for biochemical, histological and gene expression analyses. Age-matched term-born pigs were sacrificed at similar postnatal ages and served as thereference group. Compared with term pigs, preterm pigs exhibited impaired kidney maturation, as indicated by analyses of renal morphology, histopathology, and inflammatory and injury markers. Supplementation with IGF-1 reduced signs of kidney immaturity, particularly in the first week of life, as indicated by improved morphology, upregulated expression of key developmental genes, reduced severity and incidence of microscopic lesions, and decreased levels of inflammatory and injury markers. No association was seen between the symptoms of necrotizing enterocolitis and kidney defects. Preterm birth in pigs impairs kidney maturation and exogenous IGF-1 treatment partially reverses this impairment. Early IGF-1 supplementation could support the development of preterm kidneys. Preterm birth may disrupt kidney development in newborns, potentially leading to morphological changes, injury, and inflammation. Preterm pigs have previously been used as models for preterm infants, but not for kidney development. IGF-1 supplementation promotes kidney maturation and alleviates renal impairments in the first week of life in preterm pigs. IGF-1 may hold potential as a supportive therapy for preterm infants sensitive to acute kidney injury.

Evolutionary dynamics of chromatin structure and duplicate gene expression in diploid and allopolyploid cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n=4X=52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this inter-genomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Using a hybrid neural network architecture for DNA sequence representation: A study on [formula omitted]-methylcytosine sites

N4-methylcytosine (4mC) is a modified form of cytosine found in DNA, contributing to epigenetic regulation. It exists in various genomes, including the Rosaceae family encompassing significant fruit crops like apples, cherries, and roses. Previous investigations have examined the distribution and functional implications of 4mC sites within the Rosaceae genome, focusing on their potential roles in gene expression regulation, environmental adaptation, and evolution. This research aims to improve the accuracy of predicting 4mC sites within the genome of Fragaria vesca, a Rosaceae plant species. Building upon the original 4mc-w2vec method, which combines word embedding processing and a convolutional neural network (CNN), we have incorporated additional feature encoding techniques and leveraged pre-trained natural language processing (NLP) models with different deep learning architectures including different forms of CNN, recurrent neural networks (RNN) and long short-term memory (LSTM). Our assessments have shown that the best model is derived from a CNN model using fastText encoding. This model demonstrates enhanced performance, achieving a sensitivity of 0.909, specificity of 0.77, and accuracy of 0.879 on an independent dataset. Furthermore, our model surpasses previously published works on the same dataset, thus showcasing its superior predictive capabilities.

Spatial transcriptomics in embryonic mouse diaphragm muscle reveals regional gradients and subdomains of developmental gene expression

Spatial transcriptomics in embryonic mouse diaphragm muscle reveals regional gradients and subdomains of developmental gene expression

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects.

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

The quantitative genetics of gene expression in Mimulus guttatus.

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05). cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTLs make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules," which has important implications for the evolution of gene expression and how it is studied by geneticists.

FEVER: an interactive web-based resource for evolutionary transcriptomics across fishes.

Teleost fish represent one of the largest and most diverse clades of vertebrates, which makes them great models in various research areas such as ecology and evolution. Recent sequencing endeavors provided high-quality genomes for species covering the main fish evolutionary lineages, opening up large-scale comparative genomics studies. However, transcriptomic data across fish species and organs are heterogenous and have not been integrated with newly sequenced genomes making gene expression quantification and comparative analyses particularly challenging. Thus, resources integrating genomic and transcriptomic data across fish species and organs are still lacking. Here, we present FEVER, a web-based resource allowing evolutionary transcriptomics across species and tissues. First, based on query genes FEVER reconstructs gene trees providing orthologous and paralogous relationships as well as their evolutionary dynamics across 13 species covering the major fish lineages, and 4 model species as evolutionary outgroups. Second, it provides unbiased gene expression across 11 tissues using up-to-date fish genomes. Finally, genomic and transcriptomic data are combined together allowing the exploration of gene expression evolution following speciation and duplication events. FEVER is freely accessible at https://fever.sk8.inrae.fr/.

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense).

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

Evolution of phenotypic variance provides insights into the genetic basis of adaptation.

Most traits are polygenic, and the contributing loci can be identified by GWAS. The genetic basis of adaptation is, however, difficult to characterize. Here, we propose to study the genetic basis of trait evolution by monitoring the evolution of their phenotypic variance during adaptation to a new environment in well-defined laboratory conditions. Extensive computer simulations show that the evolution of phenotypic variance in a replicated experimental evolution setting can distinguish between oligogenic and polygenic adaptive architectures. We compared gene expression variance in male Drosophila simulans before and after 100 generations of adaptation to a novel hot environment. The variance change in gene expression was indistinguishable for genes with and without a significant change in mean expression after 100 generations of evolution. We suggest that the majority of adaptive gene expression evolution can be explained by a polygenic architecture. We propose that tracking the evolution of phenotypic variance across generations can provide an approach to characterize the adaptive architecture.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis.

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Arabidopsis BECLIN1-induced autophagy mediates reprogramming in tapetal programmed cell death by altering the gross cellular homeostasis

In flowering plants, the tapetum degeneration in post-meiotic anther occurs through developmental programmed cell death (dPCD), which is one of the most critical and sensitive steps for the proper development of male gametophytes and fertility. Yet the pathways of dPCD, its regulation, and its interaction with autophagy remain elusive. Here, we report that high-level expression of Arabidopsis autophagy-related gene BECLIN1 (BECN1 or AtATG6) in the tobacco tapetum prior to their dPCD resulted in developmental defects. BECN1 induces severe autophagy and multiple cytoplasm-to-vacuole pathways, which alters tapetal cell reactive oxygen species (ROS)-homeostasis that represses the tapetal dPCD. The transcriptome analysis reveals that BECN1- expression caused major changes in the pathway, resulting in altered cellular homeostasis in the tapetal cell. Moreover, BECN1-mediated autophagy reprograms the execution of tapetal PCD by altering the expression of the key developmental PCD marker genes: SCPL48, CEP1, DMP4, BFN1, MC9, EXI1, and Bcl-2 member BAG5, and BAG6. This study demonstrates that BECN1-mediated autophagy is inhibitory to the dPCD of the tapetum, but the severity of autophagy leads to autophagic death in the later stages. The delayed and altered mode of tapetal degeneration resulted in male sterility.

H2AK119ub1 differentially fine-tunes gene expression by modulating canonical PRC1- and H1-dependent chromatin compaction

Polycomb repressive complex 1 (PRC1) is a key transcriptional regulator in development via modulating chromatin structure and catalyzing histone H2A ubiquitination at Lys119 (H2AK119ub1). H2AK119ub1 isone of the most abundant histone modifications in mammalian cells. However, the function of H2AK119ub1 in polycomb-mediated gene silencing remains debated. In this study, we reveal that H2AK119ub1 has two distinct roles in gene expression, through differentially modulating chromatin compaction mediated by canonical PRC1 and the linker histone H1. Interestingly, we find that H2AK119ub1 plays a positive role in transcription through interfering with the binding of canonical PRC1 to nucleosomes and therefore counteracting chromatin condensation. Conversely, we demonstrate that H2AK119ub1 facilitates H1-dependent chromatin condensation and enhances the silencing of developmental genes in mouse embryonic stem cells, suggesting that H1 may be one of several possible pathways for H2AK119ub1 in repressing transcription. These results provide insights and molecular mechanisms by which H2AK119ub1 differentially fine-tunes developmental gene expression.

ING4 and ING5 are essential for histone H3 lysine 14 acetylation and epicardial cell lineage development.

Inhibitor of growth 4 and 5 (ING4, ING5) are structurally similar chromatin-binding proteins in the KAT6A, KAT6B and KAT7 histone acetyltransferase protein complexes. Heterozygous mutations in the KAT6A or KAT6B gene cause human disorders with cardiac defects, but the contribution of their chromatin-adaptor proteins to development is unknown. We found that Ing5-/- mice had isolated cardiac ventricular septal defects. Ing4-/-Ing5-/- embryos failed to undergo chorioallantoic fusion and arrested in development at embryonic day 8.5, displaying loss of histone H3 lysine 14 acetylation, reduction in H3 lysine 23 acetylation levels and reduced developmental gene expression. Embryonic day 12.5 Ing4+/-Ing5-/- hearts showed a paucity of epicardial cells and epicardium-derived cells, failure of myocardium compaction, and coronary vasculature defects, accompanied by reduced expression of epicardium genes. Cell adhesion gene expression and proepicardium outgrowth were defective in the ING4- and ING5-deficient state. Our findings suggest that ING4 and ING5 are essential for heart development and promote epicardium and epicardium-derived cell fates and imply mutation of the human ING5 gene as a possible cause of isolated ventricular septal defects.

The evolution of seminal fluid gene expression and postmating reproductive isolation in Drosophila.

Seminal fluid protein (Sfp) genes show, in general, a higher rate of sequence divergence than genes from other categories, which is often attributed to forms of postcopulatory sexual selection or sexual conflict. Recently, the relaxation of selective constraints has been proposed as an alternative explanation for the rapid sequence evolution of Sfps and other genes with sex-limited expression. The expression of Sfp genes is a likely target of selection, but the evolution of differences in their expression levels is less understood. Here, we explore both polymorphism and divergence in Sfp gene expression between Drosophila melanogaster and Drosophila simulans, how selection might have influenced their expression, and whether changes in expression might trigger the evolution of reproductive isolating barriers. In our analysis, Sfp genes showed higher divergence, but not higher polymorphism, in expression than the average male reproductive glands gene. Sfp genes with reproductive-tissue-specific expression were enriched for both directional and stabilizing selection, while relaxed selection was the predominant mode of evolution among Sfp genes with any other nonreproductive tissue-specific or nontissue-specific expression. The knockdown of single genes known to affect intraspecific sperm competition, and with patterns of expression divergence and polymorphism suggestive of directional selection, was not enough to break down postmating reproductive isolation barriers between species. Our results identify the expression of male-specific Sfp genes as an enriched target of selection and suggest a complex molecular relationship between postcopulatory sexual selection on a single gene's expression and its effect on the onset of speciation.

Hybridization and gene expression: Beyond differentially expressed genes.

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.