Aberrant splicing in maize rough endosperm3 reveals a conserved role for U12 splicing in eukaryotic multicellular development

Taxonomy of introns and the evolution of minor introns

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, major and non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing.

The temporal protein signature analyses of developing human deciduous molar tooth germ

The tooth serves as an exemplary model for developmental studies, encompassing epithelial-mesenchymal transition and cell differentiation. The essential factors and pathways identified in tooth development will help understand the natural development process and the malformations of mineralized tissues such as skeleton. The time-dependent proteomic changes were investigated through the proteomics of healthy human molars during embryonic stages, ranging from the cap-to-early bell stage. A comprehensive analysis revealed 713 differentially expressed proteins (DEPs) exhibiting five distinct temporal expression patterns. Through the application of weighted gene co-expression network analysis (WGCNA), 24 potential driver proteins of tooth development were screened, including CHID1, RAP1GDS1, HAPLN3, AKAP12, WLS, GSS, DDAH1, CLSTN1, AFM, RBP1, AGO1, SET, HMGB2, HMGB1, ANP32A, SPON1, FREM1, C8B, PRPS2, FCHO2, PPP1R12A, GPALPP1, U2AF2, and RCC2. Then, the proteomics and transcriptomics expression patterns of these proteins were further compared, complemented by single-cell RNA-sequencing (scRNA-seq). In summary, this study not only offers a wealth of information regarding the molecular intricacies of human embryonic epithelial and mesenchymal cell differentiation but also serves as an invaluable resource for future mechanistic inquiries into tooth development.

Where the minor things are: a pan-eukaryotic survey suggests neutral processes may explain much of minor intron evolution

Spliceosomal introns are gene segments removed from RNA transcripts by ribonucleoprotein machineries called spliceosomes. In some eukaryotes a second 'minor' spliceosome is responsible for processing a tiny minority of introns. Despite its seemingly modest role, minor splicing has persisted for roughly 1.5 billion years of eukaryotic evolution. Identifying minor introns in over 3000 eukaryotic genomes, we report diverse evolutionary histories including surprisingly high numbers in some fungi and green algae, repeated loss, as well as general biases in their positional and genic distributions. We estimate that ancestral minor intron densities were comparable to those of vertebrates, suggesting a trend of long-term stasis. Finally, three findings suggest a major role for neutral processes in minor intron evolution. First, highly similar patterns of minor and major intron evolution contrast with both functionalist and deleterious model predictions. Second, observed functional biases among minor intron-containing genes are largely explained by these genes' greater ages. Third, no association of intron splicing with cell proliferation in a minor intron-rich fungus suggests that regulatory roles are lineage-specific and thus cannot offer a general explanation for minor splicing's persistence. These data constitute the most comprehensive view of minor introns and their evolutionary history to date, and provide a foundation for future studies of these remarkable genetic elements.

Genetic analysis of human RNA binding motif protein 48 (RBM48) reveals an essential role in U12-type intron splicing

U12-type or minor introns are found in most multicellular eukaryotes and constitute ∼0.5% of all introns in species with a minor spliceosome. Although the biological significance for the evolutionary conservation of U12-type introns is debated, mutations disrupting U12 splicing cause developmental defects in both plants and animals. In human hematopoietic stem cells, U12 splicing defects disrupt proper differentiation of myeloid lineages and are associated with myelodysplastic syndrome, predisposing individuals to acute myeloid leukemia. Mutants in the maize ortholog of RNA binding motif protein 48 (RBM48) have aberrant U12-type intron splicing. Human RBM48 was recently purified biochemically as part of the minor spliceosome and shown to recognize the 5' end of the U6atac snRNA. In this report, we use CRISPR/Cas9-mediated ablation of RBM48 in human K-562 cells to show the genetic function of RBM48. RNA-seq analysis comparing wild-type and mutant K-562 genotypes found that 48% of minor intron-containing genes have significant U12-type intron retention in RBM48 mutants. Comparing these results to maize rbm48 mutants defined a subset of minor intron-containing genes disrupted in both species. Mutations in the majority of these orthologous minor intron-containing genes have been reported to cause developmental defects in both plants and animals. Our results provide genetic evidence that the primary defect of human RBM48 mutants is aberrant U12-type intron splicing, while a comparison of human and maize RNA-seq data identifies candidate genes likely to mediate mutant phenotypes of U12-type splicing defects.

Zrsr2 Is Essential for the Embryonic Development and Splicing of Minor Introns in RNA and Protein Processing Genes in Zebrafish

ZRSR2 (zinc finger CCCH-type, RNA binding motif and serine/arginine rich 2) is an essential splicing factor involved in 3' splice-site recognition as a component of both the major and minor spliceosomes that mediate the splicing of U2-type (major) and U12-type (minor) introns, respectively. Studies of ZRSR2-depleted cell lines and ZRSR2-mutated patient samples revealed its essential role in the U12-dependent minor spliceosome. However, the role of ZRSR2 during embryonic development is not clear, as its function is compensated for by Zrsr1 in mice. Here, we utilized the zebrafish model to investigate the role of zrsr2 during embryonic development. Using CRISPR/Cas9 technology, we generated a zrsr2-knockout zebrafish line, termed zrsr2hg129/hg129 (p.Trp167Argfs*9) and examined embryo development in the homozygous mutant embryos. zrsr2hg129/hg129 embryos displayed multiple developmental defects starting at 4 days post fertilization (dpf) and died after 8 dpf, suggesting that proper Zrsr2 function is required during embryonic development. The global transcriptome analysis of 3 dpf zrsr2hg129/hg129 embryos revealed that the loss of Zrsr2 results in the downregulation of essential metabolic pathways and the aberrant retention of minor introns in about one-third of all minor intron-containing genes in zebrafish. Overall, our study has demonstrated that the role of Zrsr2 as a component of the minor spliceosome is conserved and critical for proper embryonic development in zebrafish.

Ubiquitination and Ubiquitin-Like Modifications as Mediators of Alternative Pre-mRNA Splicing in Arabidopsis thaliana

Alternative splicing (AS) is a common post-transcriptional regulatory process in eukaryotes. AS has an irreplaceable role during plant development and in response to environmental stress as it evokes differential expression of downstream genes or splicing factors (e.g., serine/arginine-rich proteins). Numerous studies have reported that loss of AS capacity leads to defects in plant growth and development, and induction of stress-sensitive phenotypes. A role for post-translational modification (PTM) of AS components has emerged in recent years. These modifications are capable of regulating the activity, stability, localization, interaction, and folding of spliceosomal proteins in human cells and yeast, indicating that PTMs represent another layer of AS regulation. In this review, we summarize the recent reports concerning ubiquitin and ubiquitin-like modification of spliceosome components and analyze the relationship between spliceosome and the ubiquitin/26S proteasome pathway in plants. Based on the totality of the evidence presented, we further speculate on the roles of protein ubiquitination mediated AS in plant development and environmental response.

Zrsr2 and functional U12-dependent spliceosome are necessary for follicular development

ZRSR2 is a splicing factor involved in recognition of 3'-intron splice sites that is frequently mutated in myeloid malignancies and several tumors; however, the role of mutations of Zrsr2 in other tissues has not been analyzed. To explore the biological role of ZRSR2, we generated three Zrsr2 mutant mouse lines. All Zrsr2 mutant lines exhibited blood cell anomalies, and in two lines, oogenesis was blocked at the secondary follicle stage. RNA-seq of Zrsr2 mu secondary follicles showed aberrations in gene expression and showed altered alternative splicing (AS) events involving enrichment of U12-type intron retention (IR), supporting the functional Zrsr2 action in minor spliceosomes. IR events were preferentially associated with centriole replication, protein phosphorylation, and DNA damage checkpoint. Notably, we found alterations in AS events of 50 meiotic genes. These results indicate that ZRSR2 mutations alter splicing mainly in U12-type introns, which may affect peripheral blood cells, and impede oogenesis and female fertility.

Alternative splicing during fruit development among fleshy fruits

BACKGROUND

Alternative splicing (AS) is an important mechanism of posttranscriptional modification and dynamically regulates multiple physiological processes in plants, including fruit ripening. However, little is known about alternative splicing during fruit development in fleshy fruits.

RESULTS

We studied the alternative splicing at the immature and ripe stages during fruit development in cucumber, melon, papaya and peach. We found that 14.96-17.48% of multiexon genes exhibited alternative splicing. Intron retention was not always the most frequent event, indicating that the alternative splicing pattern during different developmental process differs. Alternative splicing was significantly more prevalent at the ripe stage than at the immature stage in cucumber and melon, while the opposite trend was shown in papaya and peach, implying that developmental stages adopt different alternative splicing strategies for their specific functions. Some genes involved in fruit ripening underwent stage-specific alternative splicing, indicating that alternative splicing regulates fruits ripening. Conserved alternative splicing events did not appear to be stage-specific. Clustering fruit developmental stages across the four species based on alternative splicing profiles resulted in species-specific clustering, suggesting that diversification of alternative splicing contributes to lineage-specific evolution in fleshy fruits.

CONCLUSIONS

We obtained high quality transcriptomes and alternative splicing events during fruit development across the four species. Dynamics and nonconserved alternative splicing were discovered. The candidate stage-specific AS genes involved in fruit ripening will provide valuable insight into the roles of alternative splicing during the developmental processes of fleshy fruits.

Post-transcriptional regulation of seed dormancy and germination: Current understanding and future directions

Seed dormancy is a developmental checkpoint that prevents mature seeds from germinating under conditions that are otherwise favorable for germination. Temperature and light are the most relevant environmental factors that regulate seed dormancy and germination. These environmental cues can trigger molecular and physiological responses including hormone signaling, particularly that of abscisic acid and gibberellin. The balance between the content and sensitivity of these hormones is the key to the regulation of seed dormancy. Temperature and light tightly regulate the transcription of thousands of genes, as well as other aspects of gene expression such as mRNA splicing, translation, and stability. Chromatin remodeling determines specific transcriptional outputs, and alternative splicing leads to different outcomes and produces transcripts that encode proteins with altered or lost functions. Proper regulation of chromatin remodeling and alternative splicing may be highly relevant to seed germination. Moreover, microRNAs are also critical for the control of gene expression in seeds. This review aims to discuss recent updates on post-transcriptional regulation during seed maturation, dormancy, germination, and post-germination events. We propose future prospects for understanding how different post-transcriptional processes in crop seeds can contribute to the design of genotypes with better performance and higher productivity.

Leaf senescence: progression, regulation, and application

Leaf senescence, the last stage of leaf development, is a type of postmitotic senescence and is characterized by the functional transition from nutrient assimilation to nutrient remobilization which is essential for plants' fitness. The initiation and progression of leaf senescence are regulated by a variety of internal and external factors such as age, phytohormones, and environmental stresses. Significant breakthroughs in dissecting the molecular mechanisms underpinning leaf senescence have benefited from the identification of senescence-altered mutants through forward genetic screening and functional assessment of hundreds of senescence-associated genes (SAGs) via reverse genetic research in model plant Arabidopsis thaliana as well as in crop plants. Leaf senescence involves highly complex genetic programs that are tightly tuned by multiple layers of regulation, including chromatin and transcription regulation, post-transcriptional, translational and post-translational regulation. Due to the significant impact of leaf senescence on photosynthesis, nutrient remobilization, stress responses, and productivity, much effort has been made in devising strategies based on known senescence regulatory mechanisms to manipulate the initiation and progression of leaf senescence, aiming for higher yield, better quality, or improved horticultural performance in crop plants. This review aims to provide an overview of leaf senescence and discuss recent advances in multi-dimensional regulation of leaf senescence from genetic and molecular network perspectives. We also put forward the key issues that need to be addressed, including the nature of leaf age, functional stay-green trait, coordination between different regulatory pathways, source-sink relationship and nutrient remobilization, as well as translational researches on leaf senescence.

Accelerating biological insight for understudied genes

The rapid expansion of genome sequence data is increasing the discovery of protein-coding genes across all domains of life. Annotating these genes with reliable functional information is necessary to understand evolution, to define the full biochemical space accessed by nature, and to identify target genes for biotechnology improvements. The vast majority of proteins are annotated based on sequence conservation with no specific biological, biochemical, genetic, or cellular function identified. Recent technical advances throughout the biological sciences enable experimental research on these understudied protein-coding genes in a broader collection of species. However, scientists have incentives and biases to continue focusing on well documented genes within their preferred model organism. This perspective suggests a research model that seeks to break historic silos of research bias by enabling interdisciplinary teams to accelerate biological functional annotation. We propose an initiative to develop coordinated projects of collaborating evolutionary biologists, cell biologists, geneticists, and biochemists that will focus on subsets of target genes in multiple model organisms. Concurrent analysis in multiple organisms takes advantage of evolutionary divergence and selection, which causes individual species to be better suited as experimental models for specific genes. Most importantly, multisystem approaches would encourage transdisciplinary critical thinking and hypothesis testing that is inherently slow in current biological research.

Minor intron retention drives clonal hematopoietic disorders and diverse cancer predisposition

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

Maize kernel development

Maize (Zea mays) is a leading cereal crop in the world. The maize kernel is the storage organ and the harvest portion of this crop and is closely related to its yield and quality. The development of maize kernel is initiated by the double fertilization event, leading to the formation of a diploid embryo and a triploid endosperm. The embryo and endosperm are then undergone independent developmental programs, resulting in a mature maize kernel which is comprised of a persistent endosperm, a large embryo, and a maternal pericarp. Due to the well-characterized morphogenesis and powerful genetics, maize kernel has long been an excellent model for the study of cereal kernel development. In recent years, with the release of the maize reference genome and the development of new genomic technologies, there has been an explosive expansion of new knowledge for maize kernel development. In this review, we overviewed recent progress in the study of maize kernel development, with an emphasis on genetic mapping of kernel traits, transcriptome analysis during kernel development, functional gene cloning of kernel mutants, and genetic engineering of kernel traits.

Engineering 6-phosphogluconate dehydrogenase improves grain yield in heat-stressed maize

Endosperm starch synthesis is a primary determinant of grain yield and is sensitive to high-temperature stress. The maize chloroplast-localized 6-phosphogluconate dehydrogenase (6PGDH), PGD3, is critical for endosperm starch accumulation. Maize also has two cytosolic isozymes, PGD1 and PGD2, that are not required for kernel development. We found that cytosolic PGD1 and PGD2 isozymes have heat-stable activity, while amyloplast-localized PGD3 activity is labile under heat stress conditions. We targeted heat-stable 6PGDH to endosperm amyloplasts by fusing the Waxy1 chloroplast targeting the peptide coding sequence to the Pgd1 and Pgd2 open reading frames (ORFs). These WPGD1 and WPGD2 fusion proteins import into isolated chloroplasts, demonstrating a functional targeting sequence. Transgenic maize plants expressing WPGD1 and WPGD2 with an endosperm-specific promoter increased 6PGDH activity with enhanced heat stability in vitro. WPGD1 and WPGD2 transgenes complement the pgd3-defective kernel phenotype, indicating the fusion proteins are targeted to the amyloplast. In the field, the WPGD1 and WPGD2 transgenes can mitigate grain yield losses in high-nighttime-temperature conditions by increasing kernel number. These results provide insight into the subcellular distribution of metabolic activities in the endosperm and suggest the amyloplast pentose phosphate pathway is a heat-sensitive step in maize kernel metabolism that contributes to yield loss during heat stress.

Engineering 6-phosphogluconate dehydrogenase improves grain yield in heat-stressed maize

Endosperm starch synthesis is a primary determinant of grain yield and is sensitive to high-temperature stress. The maize chloroplast-localized 6-phosphogluconate dehydrogenase (6PGDH), PGD3, is critical for endosperm starch accumulation. Maize also has two cytosolic isozymes, PGD1 and PGD2, that are not required for kernel development. We found that cytosolic PGD1 and PGD2 isozymes have heat-stable activity, while amyloplast-localized PGD3 activity is labile under heat stress conditions. We targeted heat-stable 6PGDH to endosperm amyloplasts by fusing the Waxy1 chloroplast targeting the peptide coding sequence to the Pgd1 and Pgd2 open reading frames (ORFs). These WPGD1 and WPGD2 fusion proteins import into isolated chloroplasts, demonstrating a functional targeting sequence. Transgenic maize plants expressing WPGD1 and WPGD2 with an endosperm-specific promoter increased 6PGDH activity with enhanced heat stability in vitro. WPGD1 and WPGD2 transgenes complement the pgd3-defective kernel phenotype, indicating the fusion proteins are targeted to the amyloplast. In the field, the WPGD1 and WPGD2 transgenes can mitigate grain yield losses in high-nighttime-temperature conditions by increasing kernel number. These results provide insight into the subcellular distribution of metabolic activities in the endosperm and suggest the amyloplast pentose phosphate pathway is a heat-sensitive step in maize kernel metabolism that contributes to yield loss during heat stress.

Splicing impact of deep exonic missense variants in CAPN3 explored systematically by minigene functional assay

Improving the accuracy of variant interpretation during diagnostic sequencing is a major goal for genomic medicine. To explore an often-overlooked splicing effect of missense variants, we developed the functional assay ("minigene") for the majority of exons of CAPN3, the gene responsible for limb girdle muscular dystrophy. By systematically screening 21 missense variants distributed along the gene, we found that eight clinically relevant missense variants located at a certain distance from the exon-intron borders (deep exonic missense variants) disrupted normal splicing of CAPN3 exons. Several recent machine learning-based computational tools failed to predict splicing impact for the majority of these deep exonic missense variants, highlighting the importance of including variants of this type in the training sets during the future algorithm development. Overall, 24 variants in CAPN3 gene were explored, leading to the change in the American College of Medical Genetics and Genomics classification of seven of them when results of the "minigene" functional assay were considered. Our findings reveal previously unknown splicing impact of several clinically important variants in CAPN3 and draw attention to the existence of deep exonic variants with a disruptive effect on gene splicing that could be overlooked by the current approaches in clinical genetics.

Minor Splicing Factors Zrsr1 and Zrsr2 Are Essential for Early Embryo Development and 2-Cell-Like Conversion

Minor splicing plays an important role in vertebrate development. Zrsr1 and Zrsr2 paralog genes have essential roles in alternative splicing, mainly participating in the recognition of minor (U12) introns. To further explore their roles during early embryo development, we produced Zrsr1mu and Zrsr2mu mutant mice, containing truncating mutations within the second zinc finger domain. Both homozygous mutant mice were viable with a normal lifespan. When we crossed a homozygous Zrsr2mu/mu female with Zrsr1mu/mu male, the double heterozygotes were non-viable, giving rise to embryos that stopped developing mainly between the 2- and 4-cell stages, just after zygotic gene activation. RNA-seq analysis of Zrsr1/2mu 2-cell embryos showed altered gene and isoform expression of thousands of genes enriched in gene ontology terms and biological pathways related to ribosome, RNA transport, spliceosome, and essential zygotic gene activation steps. Alternative splicing was analyzed, showing a significant increase in intron retention in both U2 and U12 intron-containing genes related to cell cycle and mitotic nuclear division. Remarkably, both Zrsr1 and Zrsr2 were required for the conversion of mouse-induced pluripotent stem cells into 2C-like cells. According to our results, Zrsr1 or Zrsr2 are necessary for ZGA and both are indispensable for the conversion of induced pluripotent stem cells into 2C-like cells.

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.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

Minor intron splicing revisited: identification of new minor intron-containing genes and tissue-dependent retention and alternative splicing of minor introns

BACKGROUND

Mutations in minor spliceosome components such as U12 snRNA (cerebellar ataxia) and U4atac snRNA (microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1)) result in tissue-specific symptoms. Given that the minor spliceosome is ubiquitously expressed, we hypothesized that these restricted phenotypes might be caused by the tissue-specific regulation of the minor spliceosome targets, i.e. minor intron-containing genes (MIGs). The current model of inefficient splicing is thought to apply to the regulation of the ~ 500 MIGs identified in the U12DB. However this database was created more than 10 years ago. Therefore, we first wanted to revisit the classification of minor introns in light of the most recent reference genome. We then sought to address specificity of MIG expression, minor intron retention, and alternative splicing (AS) across mouse and human tissues.

RESULTS

We employed position-weight matrices to obtain a comprehensive updated list of minor introns, consisting of 722 mouse and 770 human minor introns. These can be found in the Minor Intron DataBase (MIDB). Besides identification of 99% of the minor introns found in the U12DB, we also discovered ~ 150 new MIGs. We then analyzed the RNAseq data from eleven different mouse tissues, which revealed tissue-specific MIG expression and minor intron retention. Additionally, many minor introns were efficiently spliced compared to their flanking major introns. Finally, we identified several novel AS events across minor introns in both mouse and human, which were also tissue-dependent. Bioinformatics analysis revealed that several of the AS events could result in the production of novel tissue-specific proteins. Moreover, like the major introns, we found that these AS events were more prevalent in long minor introns, while retention was favoured in shorter introns.

CONCLUSION

Here we show that minor intron splicing and AS across minor introns is a highly organised process that might be regulated in coordination with the major spliceosome in a tissue-specific manner. We have provided a framework to further study the impact of the minor spliceosome and the regulation of MIG expression. These findings may shed light on the mechanism underlying tissue-specific phenotypes in diseases associated with minor spliceosome inactivation. MIDB can be accessed at https://midb.pnb.uconn.edu .

Sex-Dimorphic Behavioral Alterations and Altered Neurogenesis in U12 Intron Splicing-Defective Zrsr1 Mutant Mice

Mutant mice with respect to the splicing factor Zrsr1 present altered spermatogenesis and infertility. To investigate whether Zrsr1 is involved in the homeostatic control that the hypothalamus exerts over reproductive functions, we first analyzed both differential gene and isoform expression and alternative splicing alterations in Zrsr1 mutant (Zrsr1^mu) hypothalamus; second, we analyzed the spontaneous and social behavior of Zrsr1^mu mice; and third, we analyzed adult cell proliferation and survival in the Zrsr1^mu hypothalamus. The Zrsr1^mu hypothalamus showed altered expression of genes and isoforms related to the glutathione metabolic process, synaptonemal complex assembly, mRNA transport, and altered splicing events involving the enrichment of U12-type intron retention (IR). Furthermore, increased IR in U12-containing genes related with the prolactin, progesterone, and gonadotropin-releasing hormone (GnRH) reproductive signaling pathway was observed. This was associated with a hyperactive phenotype in both males and females, with an anxious phenotype in females, and with increased social interaction in males, instead of the classical aggressive behavior. In addition, Zrsr1^mu females but not males exhibited reduced cell proliferation in both the hypothalamus and the subventricular zone. Overall, these results suggest that Zrsr1 expression and function are relevant to organization of the hypothalamic cell network controlling behavior.

Dek42 encodes an RNA-binding protein that affects alternative pre-mRNA splicing and maize kernel development

RNA-binding proteins (RBPs) play an important role in post-transcriptional gene regulation. However, the functions of RBPs in plants remain poorly understood. Maize kernel mutant dek42 has small defective kernels and lethal seedlings. Dek42 was cloned by Mutator tag isolation and further confirmed by an independent mutant allele and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 materials. Dek42 encodes an RRM_RBM48 type RNA-binding protein that localizes to the nucleus. Dek42 is constitutively expressed in various maize tissues. The dek42 mutation caused a significant reduction in the accumulation of DEK42 protein in mutant kernels. RNA-seq analysis showed that the dek42 mutation significantly disturbed the expression of thousands of genes during maize kernel development. Sequence analysis also showed that the dek42 mutation significantly changed alternative splicing in expressed genes, which were especially enriched for the U12-type intron-retained type. Yeast two-hybrid screening identified SF3a1 as a DEK42-interacting protein. DEK42 also interacts with the spliceosome component U1-70K. These results suggested that DEK42 participates in the regulation of pre-messenger RNA splicing through its interaction with other spliceosome components. This study showed the function of a newly identified RBP and provided insights into alternative splicing regulation during maize kernel development.

RNA Binding Motif Protein 48 Is Required for U12 Splicing and Maize Endosperm Differentiation

The last eukaryotic common ancestor had two classes of introns that are still found in most eukaryotic lineages. Common U2-type and rare U12-type introns are spliced by the major and minor spliceosomes, respectively. Relatively few splicing factors have been shown to be specific to the minor spliceosome. We found that the maize (Zea mays) RNA binding motif protein 48 (RBM48) is a U12 splicing factor that functions to promote cell differentiation and repress cell proliferation. RBM48 is coselected with the U12 splicing factor, zinc finger CCCH-type, RNA binding motif, and Ser/Arg rich 2/Rough endosperm 3 (RGH3). Protein-protein interactions between RBM48, RGH3, and U2 Auxiliary Factor (U2AF) subunits suggest major and minor spliceosome factors required for intron recognition form complexes with RBM48. Human RBM48 interacts with armadillo repeat containing 7 (ARMC7). Maize RBM48 and ARMC7 have a conserved protein-protein interaction. These data predict that RBM48 is likely to function in U12 splicing throughout eukaryotes and that U12 splicing promotes endosperm cell differentiation in maize.

IntEREst: intron-exon retention estimator

BACKGROUND

In-depth study of the intron retention levels of transcripts provide insights on the mechanisms regulating pre-mRNA splicing efficiency. Additionally, detailed analysis of retained introns can link these introns to post-transcriptional regulation or identify aberrant splicing events in human diseases.

RESULTS

We present IntEREst, Intron-Exon Retention Estimator, an R package that supports rigorous analysis of non-annotated intron retention events (in addition to the ones annotated by RefSeq or similar databases), and support intra-sample in addition to inter-sample comparisons. It accepts binary sequence alignment/map (.bam) files as input and determines genome-wide estimates of intron retention or exon-exon junction levels. Moreover, it includes functions for comparing subsets of user-defined introns (e.g. U12-type vs U2-type) and its plotting functions allow visualization of the distribution of the retention levels of the introns. Statistical methods are adapted from the DESeq2, edgeR and DEXSeq R packages to extract the significantly more or less retained introns. Analyses can be performed either sequentially (on single core) or in parallel (on multiple cores). We used IntEREst to investigate the U12- and U2-type intron retention in human and plant RNAseq dataset with defects in the U12-dependent spliceosome due to mutations in the ZRSR2 component of this spliceosome. Additionally, we compared the retained introns discovered by IntEREst with that of other methods and studies.

CONCLUSION

IntEREst is an R package for Intron retention and exon-exon junction levels analysis of RNA-seq data. Both the human and plant analyses show that the U12-type introns are retained at higher level compared to the U2-type introns already in the control samples, but the retention is exacerbated in patient or plant samples carrying a mutated ZRSR2 gene. Intron retention events caused by ZRSR2 mutation that we discovered using IntEREst (DESeq2 based function) show considerable overlap with the retained introns discovered by other methods (e.g. IRFinder and edgeR based function of IntEREst). Our results indicate that increase in both the number of biological replicates and the depth of sequencing library promote the discovery of retained introns, but the effect of library size gradually decreases with more than 35 million reads mapped to the introns.

The emerging role of minor intron splicing in neurological disorders

Pre-mRNA splicing is an essential step in eukaryotic gene expression. Mutations in cis-acting sequence elements within pre-mRNA molecules or trans-acting factors involved in pre-mRNA processing have both been linked to splicing dysfunction that give rise to a large number of human diseases. These mutations typically affect the major splicing pathway, which excises more than 99% of all introns in humans. However, approximately 700-800 human introns feature divergent intron consensus sequences at their 5' and 3' ends and are recognized by a separate pre-mRNA processing machinery denoted as the minor spliceosome. This spliceosome has been studied less than its major counterpart, but has received increasing attention during the last few years as a novel pathomechanistic player on the stage in neurodevelopmental and neurodegenerative diseases. Here, we review the current knowledge on minor spliceosome function and discuss its potential pathomechanistic role and impact in neurodegeneration.

Minor spliceosome and disease

The U12-dependent (minor) spliceosome excises a rare group of introns that are characterized by a highly conserved 5' splice site and branch point sequence. Several new congenital or somatic diseases have recently been associated with mutations in components of the minor spliceosome. A common theme in these diseases is the detection of elevated levels of transcripts containing U12-type introns, of which a subset is associated with other splicing defects. Here we review the present understanding of minor spliceosome diseases, particularly those associated with the specific components of the minor spliceosome. We also present a model for interpreting the molecular-level consequences of the different diseases.

Isobaric tags for relative and absolute quantification-based proteomic analysis of testis biopsies in rhesus monkeys treated with transient scrotal hyperthermia

This study aimed to examine the cellular and molecular events that occur in rhesus monkey testes after scrotal hyperthermia. Eight male adult rhesus monkeys were subjected to scrotal hyperthermia at 43°C for 30 min daily for 6 consecutive days. Sperm concentration, reproductive hormones, and testis histology were examined before hyperthermia (day 0), and at 8, 15, 30, 45, 60, 75, and 90 days after the initiation of hyperthermia. iTRAQ-based proteomic analysis was conducted on testicular tissues collected on days 0, 8, and 60 to identify differentially expressed proteins at the early and recovery stages of testicular damage. The sperm concentration was significantly decreased at days 30 and 45 after treatment (p < 0.01) and recovered to baseline at day 60. When compared with day 0, 101 and 24 differentially expressed proteins were identified at days 8 and 60 after heat treatment, respectively. The molecular functions of the differentially expressed proteins at day 8 were mainly nucleic acid binding, unfolded protein binding, nucleotide binding, and nucleoside phosphate binding. Spliceosome was enriched as the most significant pathway at day 8. CIRBP, PSIP1, Sam68, and Decorin were validated and found to be consistent with the proteomic data, indicating the reliability of the proteomic profiles identified in this study. In summary, we suggest that the proteins identified in this study may play important roles in heat-induced spermatogenic impairment. Some of these proteins, such as CIRBP, PSIP1, Sam68, and Decorin, may be early molecular targets responsible for spermatogenesis suppression induced by heat treatment.

A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize

Identifying and characterizing alternative splicing (AS) enables our understanding of the biological role of transcript isoform diversity. This study describes the use of publicly available RNA-Seq data to identify and characterize the global diversity of AS isoforms in maize using the inbred lines B73 and Mo17, and a related species, sorghum. Identification and characterization of AS within maize tissues revealed that genes expressed in seed exhibit the largest differential AS relative to other tissues examined. Additionally, differences in AS between the two genotypes B73 and Mo17 are greatest within genes expressed in seed. We demonstrate that changes in the level of alternatively spliced transcripts (intron retention and exon skipping) do not solely reflect differences in total transcript abundance, and we present evidence that intron retention may act to fine-tune gene expression across seed development stages. Furthermore, we have identified temperature sensitive AS in maize and demonstrate that drought-induced changes in AS involve distinct sets of genes in reproductive and vegetative tissues. Examining our identified AS isoforms within B73 × Mo17 recombinant inbred lines (RILs) identified splicing QTL (sQTL). The 43.3% of cis-sQTL regulated junctions are actually identified as alternatively spliced junctions in our analysis, while 10 Mb windows on each side of 48.2% of trans-sQTLs overlap with splicing related genes. Using sorghum as an out-group enabled direct examination of loss or conservation of AS between homeologous genes representing the two subgenomes of maize. We identify several instances where AS isoforms that are conserved between one maize homeolog and its sorghum ortholog are absent from the second maize homeolog, suggesting that these AS isoforms may have been lost after the maize whole genome duplication event. This comprehensive analysis provides new insights into the complexity of AS in maize.