Evolutionary innovation in conserved regulatory elements across the mammalian tree of life.

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

CpG island turnover events predict evolutionary changes in enhancer activity.

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Trp53 ablation fails to prevent microcephaly in mouse pallium with impaired minor intron splicing.

Minor spliceosome inhibition due to mutations in RNU4ATAC are linked to primary microcephaly. Ablation of Rnu11, which encodes a minor spliceosome snRNA, inhibits the minor spliceosome in the developing mouse pallium, causing microcephaly. There, cell cycle defects and p53-mediated apoptosis in response to DNA damage resulted in loss of radial glial cells (RGCs), underpinning microcephaly. Here, we ablated Trp53 to block cell death in Rnu11 cKO mice. We report that Trp53 ablation failed to prevent microcephaly in these double knockout (dKO) mice. We show that the transcriptome of the dKO pallium was more similar to the control compared with the Rnu11 cKO. We find aberrant minor intron splicing in minor intron-containing genes involved in cell cycle regulation, resulting in more severely impaired mitotic progression and cell cycle lengthening of RGCs in the dKO that was detected earlier than in the Rnu11 cKO. Furthermore, we discover a potential role of p53 in causing DNA damage in the developing pallium, as detection of γH2aX+ was delayed in the dKO. Thus, we postulate that microcephaly in minor spliceosome-related diseases is primarily caused by cell cycle defects.

Role of NDE1 in the Development and Evolution of the Gyrified Cortex.

An expanded cortex is a hallmark of human neurodevelopment and endows increased cognitive capabilities. Recent work has shown that the cell cycle-related gene NDE1 is essential for proper cortical development. Patients who have mutations in NDE1 exhibit congenital microcephaly as a primary phenotype. At the cellular level, NDE1 is essential for interkinetic nuclear migration and mitosis of radial glial cells, which translates to an indispensable role in neurodevelopment. The nuclear migration function of NDE1 is well conserved across Opisthokonta. In mammals, multiple isoforms containing alternate terminal exons, which influence the functionality of NDE1, have been reported. It has been noted that the pattern of terminal exon usage mirrors patterns of cortical complexity in mammals. To provide context to these findings, here, we provide a comprehensive review of the literature regarding NDE1, its molecular biology and physiological relevance at the cellular and organismal levels. In particular, we outline the potential roles of NDE1 in progenitor cell behavior and explore the spectrum of NDE1 pathogenic variants. Moreover, we assessed the evolutionary conservation of NDE1 and interrogated whether the usage of alternative terminal exons is characteristic of species with gyrencephalic cortices. We found that gyrencephalic species are more likely to express transcripts that use the human-associated terminal exon, whereas lissencephalic species tend to express transcripts that use the mouse-associated terminal exon. Among gyrencephalic species, the human-associated terminal exon was preferentially expressed by those with a high order of gyrification. These findings underscore phylogenetic relationships between the preferential usage of NDE1 terminal exon and high-order gyrification, which provide insight into cortical evolution underlying high-order brain functions.

An Integrated Model of Minor Intron Emergence and Conservation.

Minor introns constitute <0.5% of the introns in the human genome and have remained an enigma since their discovery. These introns are removed by a distinct splicing complex, the minor spliceosome. Both are ancient, tracing back to the last eukaryotic common ancestor (LECA), which is reflected by minor intron enrichment in specific gene families, such as the mitogen activated-protein kinase kinases, voltage-gated sodium and calcium ion channels, and E2F transcription factors. Most minor introns occur as single introns in genes with predominantly major introns. Due to this organization, minor intron-containing gene (MIG) expression requires the coordinated action of two spliceosomes, which increases the probability of missplicing. Thus, one would expect loss of minor introns via purifying selection. This has resulted in complete minor intron loss in at least nine eukaryotic lineages. However, minor introns are highly conserved in land plants and metazoans, where their importance is underscored by embryonic lethality when the minor spliceosome is inactivated. Conditional inactivation of the minor spliceosome has shown that rapidly dividing progenitor cells are highly sensitive to minor spliceosome loss. Indeed, we found that MIGs were significantly enriched in a screen for genes essential for survival in 341 cycling cell lines. Here, we propose that minor introns inserted randomly into genes in LECA or earlier and were subsequently conserved in genes crucial for cycling cell survival. We hypothesize that the essentiality of MIGs allowed minor introns to endure through the unicellularity of early eukaryotic evolution. Moreover, we identified 59 MIGs that emerged after LECA, and that many of these are essential for cycling cell survival, reinforcing our essentiality model for MIG conservation. This suggests that minor intron emergence is dynamic across eukaryotic evolution, and that minor introns should not be viewed as molecular fossils. We also posit that minor intron splicing was co-opted in multicellular evolution as a regulatory switch for en masse control of MIG expression and the biological processes they regulate. Specifically, this mode of regulation could control cell proliferation and thus body size, an idea supported by domestication syndrome, wherein MIGs are enriched in common candidate animal domestication genes.

Minor spliceosome inactivation causes microcephaly, owing to cell cycle defects and death of self-amplifying radial glial cells.

Mutation in minor spliceosome components is linked to the developmental disorder microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1). Here, we inactivated the minor spliceosome in the developing mouse cortex (pallium) by ablating Rnu11, which encodes the crucial minor spliceosome small nuclear RNA (snRNA) U11. Rnu11 conditional knockout mice were born with microcephaly, which was caused by the death of self-amplifying radial glial cells (RGCs), while intermediate progenitor cells and neurons were produced. RNA sequencing suggested that this cell death was mediated by upregulation of p53 (Trp53 - Mouse Genome Informatics) and DNA damage, which were both observed specifically in U11-null RGCs. Moreover, U11 loss caused elevated minor intron retention in genes regulating the cell cycle, which was consistent with fewer RGCs in S-phase and cytokinesis, alongside prolonged metaphase in RGCs. In all, we found that self-amplifying RGCs are the cell type most sensitive to loss of minor splicing. Together, these findings provide a potential explanation of how disruption of minor splicing might cause microcephaly in MOPD1.

Network-based bioinformatics analysis of spatio-temporal RNA-Seq data reveals transcriptional programs underpinning normal and aberrant retinal development.

BACKGROUND: The retina as a model system with extensive information on genes involved in development/maintenance is of great value for investigations employing deep sequencing to capture transcriptome change over time. This in turn could enable us to find patterns in gene expression across time to reveal transition in biological processes. METHODS: We developed a bioinformatics pipeline to categorize genes based on their differential expression and their alternative splicing status across time by binning genes based on their transcriptional kinetics. Genes within same bins were then leveraged to query gene annotation databases to discover molecular programs employed by the developing retina. RESULTS: Using our pipeline on RNA-Seq data obtained from fractionated (nucleus/cytoplasm) developing retina at embryonic day (E) 16 and postnatal day (P) 0, we captured high-resolution as in the difference between the cytoplasm and the nucleus at the same developmental time. We found de novo transcription of genes whose transcripts were exclusively found in the nuclear transcriptome at P0. Further analysis showed that these genes enriched for functions that are known to be executed during postnatal development, thus showing that the P0 nuclear transcriptome is temporally ahead of that of its cytoplasm. We extended our strategy to perform temporal analysis comparing P0 data to either P21-Nrl-wildtype (WT) or P21-Nrl-knockout (KO) retinae, which predicted that the KO retina would have compromised vasculature. Indeed, histological manifestation of vasodilation has been reported at a later time point (P60). CONCLUSIONS: Thus, our approach was predictive of a phenotype before it presented histologically. Our strategy can be extended to investigating the development and/or disease progression of other tissue types.

Minor splicing snRNAs are enriched in the developing mouse CNS and are crucial for survival of differentiating retinal neurons.

In eukaryotes, gene expression requires splicing, which starts with the identification of exon-intron boundaries by the small, nuclear RNA (snRNAs) of the spliceosome, aided by associated proteins. In the mammalian genome, <1% of introns lack canonical exon-intron boundary sequences and cannot be spliced by the canonical splicing machinery. These introns are spliced by the minor spliceosome, consisting of unique snRNAs (U11, U12, U4atac, and U6atac). The importance of the minor spliceosome is underscored by the disease microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1), which is caused by mutation in U4atac. Thus, it is important to understand the expression and function of the minor spliceosome and its targets in mammalian development, for which we used the mouse as our model. Here, we report enrichment of the minor snRNAs in the developing head/central nervous system (CNS) between E9.5 and E12.5, along with enrichment of these snRNAs in differentiating retinal neurons. Moreover, dynamic expression kinetics of minor intron-containing genes (MIGs) was observed across retinal development. DAVID analysis of MIGs that were cotranscriptionally upregulated embryonically revealed enrichment for RNA metabolism and cell cycle regulation. In contrast, MIGs that were cotranscriptionally upregulated postnatally revealed enrichment for protein localization/transport, vesicle-mediated transport, and calcium transport. Finally, we used U12 morpholino to inactivate the minor spliceosome in the postnatal retina, which resulted in apoptosis of differentiating retinal neurons. Taken together, our data suggest that the minor spliceosome may have distinct functions in embryonic versus postnatal development. Importantly, we show that the minor spliceosome is crucial for the survival of terminally differentiating retinal neurons.

Replication-dependent histone genes are actively transcribed in differentiating and aging retinal neurons.

In the mammalian genome, each histone family contains multiple replication-dependent paralogs, which are found in clusters where their transcription is thought to be coupled to the cell cycle. Here, we wanted to interrogate the transcriptional regulation of these paralogs during retinal development and aging. We employed deep sequencing, quantitative PCR, in situ hybridization (ISH), and microarray analysis, which revealed that replication-dependent histone genes were not only transcribed in progenitor cells but also in differentiating neurons. Specifically, by ISH analysis we found that different histone genes were actively transcribed in a subset of neurons between postnatal day 7 and 14. Interestingly, within a histone family, not all paralogs were transcribed at the same level during retinal development. For example, expression of Hist1h1b was higher embryonically, while that of Hist1h1c was higher postnatally. Finally, expression of replication-dependent histone genes was also observed in the aging retina. Moreover, transcription of replication-dependent histones was independent of rapamycin-mediated mTOR pathway inactivation. Overall, our data suggest the existence of variant nucleosomes produced by the differential expression of the replication-dependent histone genes across retinal development. Also, the expression of a subset of replication-dependent histone isotypes in senescent neurons warrants re-examining these genes as "replication-dependent." Thus, our findings underscore the importance of understanding the transcriptional regulation of replication-dependent histone genes in the maintenance and functioning of neurons.