The first intron of EIJ1 confers a specific response to wounding and herbivore stresses

Plants are constantly exposed to different kinds of biotic stress, such as herbivore attack and wounding. To deal with these stresses, plants have evolved sophisticated defence mechanisms to protect themselves. Previously, we found that EIJ1 (EDS1-interacting J protein 1) plays a negative regulatory role in plant disease resistance in Arabidopsis thaliana. Follow-up studies revealed that EIJ1 specifically responds to wounding and herbivore stresses. The expression of EIJ1 was specifically induced by wounding or herbivore stress, as demonstrated by similar results in EIJ1 protein assay. Interestingly, GUS staining found that the promoter of EIJ1 is not involved in the induction of expression under wounding stress. Instead, we identified the first intron of EIJ1 as a key factor in response to wounding stress. Deleting the first intron of EIJ1 resulted in a loss of response to wounding stress in plants. Our results broaden the role of EIJ1 in plant resistance to biotic stress and provide new insights into plant responses to biotic stress.

Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening", is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200-bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene as well as the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

Cytochrome c oxidase subunit I gene in Thalassiosira nordenskioeldii strains inhabiting in cold and warm sea waters

From the biota beneath the sea ice in Lake Saroma, which is adjacent to Sea of Okhotsk, a diatom culture of Saroma 16 was isolated. Strutted processes and a labiate process in Saroma 16 were characteristic of those in Thalassiosira nordenskioeldii. Similarity search analysis showed that the 826-bp rbcL-3P region sequence of this strain was 100% identical to multiple sequences registered as T. nordenskioeldii in a public database. The 4305-bp PCR-amplified mitochondrial cytochrome c oxidase subunit I (COI) gene (COI)-5P region of Saroma 16 included a 1060-bp open reading frame (ORF), which was interrupted by 934-bp and 2311-bp introns that included frame-shifted ORFs encoding reverse-transcriptase (RTase)-like proteins. Previous reports showed that a strain of the same species, CNS00052, originating from the East China Sea included no introns in the COI, whereas North Atlantic Ocean strains of the same species, such as CCMP992, CCMP993, and CCMP997, included a 2.3-kb intron in the same position as Saroma 16.

Functional Characterization of Two Novel Intron 4 SERPING1 Gene Splice Site Pathogenic Variants in Families with Hereditary Angioedema

Variants that affect splice sites comprise 14.3% of all pathogenic variants in the SERPING1 gene; more than half of them are located outside the canonical sites. To make a clinical decision concerning patients with such variants, it is essential to know the exact way in which the effect of the variant would be realized. The optimal approach to determine the consequences is considered to be mRNA analysis. In the current study, we present the results of functional analysis of two previously non-described variants in the SERPING1 gene (NM_000062.3) affecting intron 4: c.686-1G>A and c.685+4dup, which were detected in members of two Russian families with autosomal dominant inheritance of angioedema type 1. Analysis of the patients' mRNA (extracted from whole blood) showed that the SERPING1(NM_000062.3):c.685+4dup variant leads to the loss of the donor splice site and the activation of the cryptic site in exon 4: r.710_745del (p.Gly217_Pro228del), while the SERPING1(NM_000062.3):c.686-1G>A variant leads to the skipping of exon 5: r.746_949del (p.Asp229_Ser296del).

Impact on splicing in Saccharomyces cerevisiae of random 50-base sequences inserted into an intron

Intron splicing is a key regulatory step in gene expression in eukaryotes. Three sequence elements required for splicing-5' and 3' splice sites and a branchpoint-are especially well-characterized in Saccharomyces cerevisiae, but our understanding of additional intron features that impact splicing in this organism is incomplete, due largely to its small number of introns. To overcome this limitation, we constructed a library in S. cerevisiae of random 50-nt (N50) elements individually inserted into the intron of a reporter gene and quantified canonical splicing and the use of cryptic splice sites by sequencing analysis. More than 70% of approximately 140,000 N50 elements reduced splicing by at least 20%. N50 features, including higher GC content, presence of GU repeats, and stronger predicted secondary structure of its pre-mRNA, correlated with reduced splicing efficiency. A likely basis for the reduced splicing of such a large proportion of variants is the formation of RNA structures that pair N50 bases-such as the GU repeats-with other bases specifically within the reporter pre-mRNA analyzed. However, multiple models were unable to explain more than a small fraction of the variance in splicing efficiency across the library, suggesting that complex nonlinear interactions in RNA structures are not accurately captured by RNA structure prediction methods. Our results imply that the specific context of a pre-mRNA may determine the bases allowable in an intron to prevent secondary structures that reduce splicing. This large data set can serve as a resource for further exploration of splicing mechanisms.

α-tubulin regulation by 5' introns in Saccharomyces cerevisiae

Across eukaryotic genomes, multiple α- and β-tubulin genes require regulation to ensure sufficient production of tubulin heterodimers. Features within these gene families that regulate expression remain underexplored. Here, we investigate the role of the 5' intron in regulating α-tubulin expression in Saccharomyces cerevisiae. We find that the intron in the α-tubulin, TUB1, promotes α-tubulin expression and cell fitness during microtubule stress. The role of the TUB1 intron depends on proximity to the TUB1 promoter and sequence features that are distinct from the intron in the alternative α-tubulin isotype, TUB3. These results lead us to perform a screen to identify genes that act with the TUB1 intron. We identified several genes involved in chromatin remodeling, α/β-tubulin heterodimer assembly, and the spindle assembly checkpoint. We propose a model where the TUB1 intron promotes expression from the chromosomal locus and that this may represent a conserved mechanism for tubulin regulation under conditions that require high levels of tubulin production.

Transcriptomic and genomic identification of spliceosomal genes from Euglena gracilis

Diverse splicing types in nuclear and chloroplast genes of protist Euglena gracilis have been recognized for decades. However, the splicing machinery responsible for processing nuclear precursor messenger RNA introns, including trans-splicing of the 5' terminal outron and spliced leader (SL) RNA, remains elusive. Here, we identify 166 spliceosomal protein genes and two snRNA genes from E. gracilis by performing bioinformatics analysis from a combination of next-generation and full-length transcriptomic RNA sequencing (RNAseq) data as well as draft genomic data. With the spliceosomal proteins we identified in hand, the insensitivity of E. gracilis to some splicing modulators is revealed at the sequence level. The prevalence of SL RNA-mediated trans-splicing is estimated to be more than 70% from our full-length RNAseq data. Finally, the splicing proteomes between E. gracilis and its three evolutionary cousins within the same Excavata group are compared. In conclusion, our study characterizes the spliceosomal components in E. gracilis and provides the molecular basis for further exploration of underlying splicing mechanisms.

Improved recovery and annotation of genes in metagenomes through the prediction of fungal introns

Metagenomics provides a tool to assess the functional potential of environmental and host-associated microbiomes based on the analysis of environmental DNA: assembly, gene prediction and annotation. While gene prediction is straightforward for most bacterial and archaeal taxa, it has limited applicability in the majority of eukaryotic organisms, including fungi that contain introns in gene coding sequences. As a consequence, eukaryotic genes are underrepresented in metagenomics datasets and our understanding of the contribution of fungi and other eukaryotes to microbiome functioning is limited. Here, we developed a machine intelligence-based algorithm that predicts fungal introns in environmental DNA with reasonable precision and used it to improve the annotation of environmental metagenomes. Intron removal increased the number of predicted genes by up to 9.1% and improved the annotation of several others. The proportion of newly predicted genes increased with the share of eukaryotic genes in the metagenome and-within fungal taxa-increased with the number of introns per gene. Our approach provides a tool named SVMmycointron for improved metagenome annotation, especially of microbiomes with a high proportion of eukaryotes. The scripts described in the paper are made publicly available and can be readily utilized by microbiome researchers analysing metagenomics data.

Intron Regions as Genetic Markers for Population Genetic Investigations of Opisthorchis viverrini sensu lato and Clonorchis sinensis

Opisthorchiasis and clonorchiasis are prevalent in Southeast and Far-East Asia, which are caused by the group 1 carcinogenic liver flukes Opisthorchis viverrini sensu lato and Clonorchis sinensis infection. There have been comprehensive investigations of systematics and genetic variation of these liver flukes. Previous studies have shown that O. viverrini is a species complex, called "O. viverrini sensu lato". More comprehensive investigations of molecular systematics and population genetics of each of the species that make up the species complex are required. Thus, other polymorphic genetic markers need to be developed. Therefore, this study aimed to characterize the intron regions of taurocyamine kinase gene (TK) to examine the genetic variation and population genetics of O. viverrini and C. sinensis collected from different geographical isolates and from a range of animal hosts. We screened seven intron regions embedded in TK. Of these, we selected an intron 5 of domain 1 (TkD1Int5) region to investigate the genetic variation and population genetics of theses liver flukes. The high nucleotide and haplotype diversity of TkD1Int5 was detected in O. viverrine. Heterozygosity with several insertion/deletion (indel) regions were detected in TkD1Int5 of the O. viverrine samples, whereas only an indel nucleotide was detected in one C. sinensis sample. Several O. viverrine samples contained three different haplotypes within a particular heterozygous sample. There were no genetic differences between C. sinensis isolated from various animal host. Heterozygous patterns specifically detected in humans was observed in C. sinensis. Thus, TkD1Int5 is a high polymorphic genetic marker, which could be an alternative marker for further population genetic investigations of these carcinogenic liver flukes and other related species from a wide geographical distribution and variety of animal hosts.

SAFB associates with nascent RNAs and can promote gene expression in mouse embryonic stem cells

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Whole F8 gene sequencing combined with splicing functional analyses led to a substantial increase of the molecular diagnosis yield for non-severe haemophilia A

INTRODUCTION

Conventional genetic investigation fails to identify the F8 causal variant in 2.5%-10% of haemophilia A (HA) patients with non-severe phenotypes. In these cases, F8 deep intronic variants could be causal.

AIM

To identify pathogenic F8 deep intronic variants in genetically unresolved families with non-severe HA analysed in the haematology laboratory of the Hospices Civils de Lyon.

METHODS

The whole F8 was analysed by next generation sequencing. The pathogenic impact of candidate variants identified was assessed using both in silico analysis (MaxEntScan and spliceAI) and functional analysis (RNA or minigene assay).

RESULTS

Sequencing was performed in 49/55 families included for which a DNA sample from a male propositus was available. In total, 33 candidate variants from 43 propositi were identified. These variants corresponded to 31 single nucleotide substitutions, one 173-bp deletion, and an 869-bp tandem triplication. No candidate variant was found in six propositi. The most frequent variants found were the association of [c.2113+1154G>C and c.5374-304C>T], identified in five propositi, and the c.2114-6529C>G identified in nine propositi. Four variants had been previously described as HA-causing. Splicing functional assay found a deleterious impact for 11 substitutions (c.671-94G>A, c.788-312A>G, c.2113+1154G>C, c.2114-6529C>G, c.5999-820A>T, c.5999-786C>A, c.5999-669G>T, c.5999-669G>A, c.5999-669G>C, c.6900+4104A>C, and c.6901-2992A>G). The HA-causing variant was identified in 33/49 (67%) cases. In total, F8 deep intronic variants caused 8.8% of the non-severe HA among the 1643 families analysed in our laboratory.

CONCLUSION

The results emphasise the value of whole F8 gene sequencing combined with splicing functional analyses to improve the diagnosis yield for non-severe HA.

Identification of a novel mitochondria-localized LKB1 variant required for the regulation of the oxidative stress response

The tumor suppressor Liver Kinase B1 (LKB1) is a multifunctional serine/threonine protein kinase that regulates cell metabolism, polarity, and growth and is associated with Peutz-Jeghers Syndrome and cancer predisposition. The LKB1 gene comprises 10 exons and 9 introns. Three spliced LKB1 variants have been documented, and they reside mainly in the cytoplasm, although two possess a nuclear-localization sequence (NLS) and are able to shuttle into the nucleus. Here, we report the identification of a fourth and novel LKB1 isoform that is, interestingly, targeted to the mitochondria. We show that this mitochondria-localized LKB1 (mLKB1) is generated from alternative splicing in the 5' region of the transcript and translated from an alternative initiation codon encoded by a previously unknown exon 1b (131 bp) hidden within the long intron 1 of LKB1 gene. We found by replacing the N-terminal NLS of the canonical LKB1 isoform, the N-terminus of the alternatively spliced mLKB1 variant encodes a mitochondrial transit peptide that allows it to localize to the mitochondria. We further demonstrate that mLKB1 colocalizes histologically with mitochondria-resident ATP Synthase and NAD-dependent deacetylase sirtuin-3, mitochondrial (SIRT3) and that its expression is rapidly and transiently upregulated by oxidative stress. We conclude that this novel LKB1 isoform, mLKB1, plays a critical role in regulating mitochondrial metabolic activity and oxidative stress response.

Mitofusin-2 gene polymorphisms and metabolic dysfunction associated fatty liver disease: a case-control study in a Chinese population

OBJECTIVES

Mitofusion-2 (Mfn2) may have a role in mitochondrial oxidative stress and insulin resistance that can promote the development of metabolic dysfunction associated fatty liver disease (MAFLD). This retrospective and case control study aimed to explore the relationships between common Mfn2 single nucleotide polymorphisms (SNPs) and MAFLD in a northern Han Chinese population.

METHODS

Six Mfn2 SNPs (rs2336384, rs873458, rs873457, rs4846085, rs2878677, and rs2236057) were genotyped using the ligase detection reaction in 466 MAFLD patients and 423 healthy controls. Genotype and allele frequencies were calculated, along with haplotype analysis and pairwise linkage disequilibrium.

RESULTS

The genotype distribution of rs2336384, rs2878677, and rs2236057 among the MAFLD patients showed a significantly different pattern from that of healthy controls. The data showed that an increased risk of MAFLD was significantly correlated with patients carrying the GG genotype of rs2336384, CC genotype of rs873457, TT genotype of rs4846085, TT genotype of rs2878677, and the AA genotype of rs2236057. Moreover, The GGCTTA haplotype was found to be adversely linked with MAFLD by haplotype analysis.

CONCLUSION

The current findings suggest a strong link between certain Mfn2 gene polymorphisms and MAFLD.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Role of the proline-rich disordered domain of DROSHA in intronic microRNA processing

DROSHA serves as a gatekeeper of the microRNA (miRNA) pathway by processing primary transcripts (pri-miRNAs). While the functions of structured domains of DROSHA have been well documented, the contribution of N-terminal proline-rich disordered domain (PRD) remains elusive. Here we show that the PRD promotes the processing of miRNA hairpins located within introns. We identified a DROSHA isoform (p140) lacking the PRD, which is produced by proteolytic cleavage. Small RNA sequencing revealed that p140 is significantly impaired in the maturation of intronic miRNAs. Consistently, our minigene constructs demonstrated that PRD enhances the processing of intronic hairpins, but not those in exons. Splice site mutations did not affect the PRD's enhancing effect on intronic constructs, suggesting that the PRD acts independently of splicing reaction by interacting with sequences residing within introns. The N-terminal regions from zebrafish and Xenopus DROSHA can replace the human counterpart, indicating functional conservation despite poor sequence alignment. Moreover, we found that rapidly evolving intronic miRNAs are generally more dependent on PRD than conserved ones, suggesting a role of PRD in miRNA evolution. Our study reveals a new layer of miRNA regulation mediated by a low-complexity disordered domain that senses the genomic contexts of miRNA loci.

New Insights into the Evolution and Gene Structure of the Mitochondrial Carrier Family Unveiled by Analyzing the Frequent and Conserved Intron Positions

Mitochondrial carriers (MCs) belong to a eukaryotic protein family of transporters that in higher organisms is called the solute carrier family 25 (SLC25). All MCs have characteristic triplicated sequence repeats forming a 3-fold symmetrical structure of a six-transmembrane α-helix bundle with a centrally located substrate-binding site. Biochemical characterization has shown that MCs altogether transport a wide variety of substrates but can be divided into subfamilies, each transporting a few specific substrates. We have investigated the intron positions in the human MC genes and their orthologs of highly diversified organisms. The results demonstrate that several intron positions are present in numerous MC sequences at the same specific points, of which some are 3-fold symmetry related. Many of these frequent intron positions are also conserved in subfamilies or in groups of subfamilies transporting similar substrates. The analyses of the frequent and conserved intron positions in MCs suggest phylogenetic relationships not only between close but also distant homologs as well as a possible involvement of the intron positions in the evolution of the substrate specificity diversification of the MC family members.

The evolution of pre-mRNA splicing and its machinery revealed by reduced extremophilic red algae

The Cyanidiales are a group of mostly thermophilic and acidophilic red algae that thrive near volcanic vents. Despite their phylogenetic relationship, the reduced genomes of Cyanidioschyzon merolae and Galdieria sulphuraria are strikingly different with respect to pre-mRNA splicing, a ubiquitous eukaryotic feature. Introns are rare and spliceosomal machinery is extremely reduced in C. merolae, in contrast to G. sulphuraria. Previous studies also revealed divergent spliceosomes in the mesophilic red alga Porphyridium purpureum and the red algal derived plastid of Guillardia theta (Cryptophyta), along with unusually high levels of unspliced transcripts. To further examine the evolution of splicing in red algae, we compared C. merolae and G. sulphuraria, investigating splicing levels, intron position, intron sequence features, and the composition of the spliceosome. In addition to identifying 11 additional introns in C. merolae, our transcriptomic analysis also revealed typical eukaryotic splicing in G. sulphuraria, whereas most transcripts in C. merolae remain unspliced. The distribution of intron positions within their host genes was examined to provide insight into patterns of intron loss in red algae. We observed increasing variability of 5' splice sites and branch donor regions with increasing intron richness. We also found these relationships to be connected to reductions in and losses of corresponding parts of the spliceosome. Our findings highlight patterns of intron and spliceosome evolution in related red algae under the pressures of genome reduction.

A novel approach to analyze the association characteristics between post-spliced introns and their corresponding mRNA

Studies have shown that post-spliced introns promote cell survival when nutrients are scarce, and intron loss/gain can influence many stages of mRNA metabolism. However, few approaches are currently available to study the correlation between intron sequences and their corresponding mature mRNA sequences. Here, based on the results of the improved Smith-Waterman local alignment-based algorithm method (SW method) and binding free energy weighted local alignment algorithm method (BFE method), the optimal matched segments between introns and their corresponding mature mRNAs in Caenorhabditis elegans (C.elegans) and their relative matching frequency (RF) distributions were obtained. The results showed that although the distributions of relative matching frequencies on mRNAs obtained by the BFE method were similar to the SW method, the interaction intensity in 5'and 3'untranslated regions (UTRs) regions was weaker than the SW method. The RF distributions in the exon-exon junction regions were comparable, the effects of long and short introns on mRNA and on the five functional sites with BFE method were similar to the SW method. However, the interaction intensity in 5'and 3'UTR regions with BFE method was weaker than with SW method. Although the matching rate and length distribution shape of the optimal matched fragment were consistent with the SW method, an increase in length was observed. The matching rates and the length of the optimal matched fragments were mainly in the range of 60%-80% and 20-30bp, respectively. Although we found that there were still matching preferences in the 5'and 3'UTR regions of the mRNAs with BFE, the matching intensities were significantly lower than the matching intensities between introns and their corresponding mRNAs with SW method. Overall, our findings suggest that the interaction between introns and mRNAs results from synergism among different types of sequences during the evolutionary process.

Interplay of heavy chain introns influences efficient transcript splicing and affects product quality of recombinant biotherapeutic antibodies from CHO cells

Introns are included in genes encoding therapeutic proteins for their well-documented function of boosting expression. However, mis-splicing of introns in recombinant immunoglobulin (IgG) heavy chain (HC) transcripts can produce amino acid sequence product variants. These variants can affect product quality; therefore, purification process optimization may be needed to remove them, or if they cannot be removed, then in-depth characterization must be carried out to understand their effects on biological activity. In this study, HC transgene engineering approaches were investigated and were successful in significantly reducing the previously identified IgG HC splice variants to <0.5%. Subsequently, a comprehensive evaluation was conducted to understand the influence of the different introns in the HC genes on the expression of recombinant biotherapeutic antibodies. The data revealed an unexpected cooperation between specific introns for efficient splicing, where intron retention led to significant reductions in IgG expression of up to 75% for some intron combinations. Furthermore, it was shown that HC introns could be fully removed without significantly affecting productivity. This work paves the way for future biotherapeutic antibody transgene design with regard to inclusion of HC introns. By removing unnecessary introns, transgene mRNA transcript will no longer be mis-spliced, thereby eliminating HC splice variants and improving antibody product quality.

Whole F9 gene sequencing identified deep intronic variations in genetically unresolved hemophilia B patients

BACKGROUND

The disease-causative variant remains unidentified in approximately 0.5% to 2% of hemophilia B patients using conventional genetic investigations, and F9 deep intronic variations could be responsible for these phenotypes.

OBJECTIVES

This study aimed to characterize deep intronic variants in hemophilia B patients for whom genetic investigations failed.

METHODS

We performed whole F9 sequencing in 17 genetically unsolved hemophilia B patients. The pathogenic impact of the candidate variants identified was studied using both in silico analysis (MaxEntScan and spliceAI) and minigene assay.

RESULTS

In total, 9 candidate variants were identified in 15 patients; 7 were deep intronic substitutions and 2 corresponded to insertions of mobile elements. The most frequent variants found were c.278-1806A>C and the association of c.278-1244A>G and c.392-864T>G, identified in 4 and 6 unrelated individuals, respectively. In silico analysis predicted splicing impact for 4 substitutions (c.278-1806A>C, c.392-864T>G, c.724-2385G>T, c.723+4297T>A). Minigene assay showed a deleterious splicing impact for these 4 substitutions and also for the c.278-1786_278-1785insLINE. In the end, 5 variants were classified as likely pathogenic using the American College of Medical Genetics and Genomics guidelines, and 4 as of unknown significance. Thus, the hemophilia B-causing variant was identified in 13/17 (76%) families.

CONCLUSION

We elucidated the causing defect in three-quarters of the families included in this study, and we reported new F9 deep intronic variants that can cause hemophilia B.

Exon Elongation Added Intrinsically Disordered Regions to the Encoded Proteins and Facilitated the Emergence of the Last Eukaryotic Common Ancestor

Most prokaryotic proteins consist of a single structural domain (SD) with little intrinsically disordered regions (IDRs) that by themselves do not adopt stable structures, whereas the typical eukaryotic protein comprises multiple SDs and IDRs. How eukaryotic proteins evolved to differ from prokaryotic proteins has not been fully elucidated. Here, we found that the longer the internal exons are, the more frequently they encode IDRs in eight eukaryotes including vertebrates, invertebrates, a fungus, and plants. Based on this observation, we propose the "small bang" model from the proteomic viewpoint: the protoeukaryotic genes had no introns and mostly encoded one SD each, but a majority of them were subsequently divided into multiple exons (step 1). Many exons unconstrained by SDs elongated to encode IDRs (step 2). The elongated exons encoding IDRs frequently facilitated the acquisition of multiple SDs to make the last common ancestor of eukaryotes (step 3). One prediction of the model is that long internal exons are mostly unconstrained exons. Analytical results of the eight eukaryotes are consistent with this prediction. In support of the model, we identified cases of internal exons that elongated after the rat-mouse divergence and discovered that the expanded sections are mostly in unconstrained exons and preferentially encode IDRs. The model also predicts that SDs followed by long internal exons tend to have other SDs downstream. This prediction was also verified in all the eukaryotic species analyzed. Our model accounts for the dichotomy between prokaryotic and eukaryotic proteins and proposes a selective advantage conferred by IDRs.

A common intronic single nucleotide variant modifies PKD1 expression level

Autosomal dominant polycystic kidney disease (ADPKD), caused by mutations in PKD1 and PKD2 (PKD1/2), has unexplained phenotypic variability likely affected by environmental and other genetic factors. Approximately 10% of individuals with ADPKD phenotype have no causal mutation detected, possibly due to unrecognized risk variants of PKD1/2. This study was designed to identify risk variants of PKD genes through population genetic analyses. We used Wright's F-statistics (Fst) to evaluate common single nucleotide variants (SNVs) potentially favored by positive natural selection in PKD1 from 1000 Genomes Project (1KG) and genotyped 388 subjects from the Rogosin Institute ADPKD Data Repository. The variants with >90th percentile Fst scores underwent further investigation by in silico analysis and molecular genetics analyses. We identified a deep intronic SNV, rs3874648G> A, located in a conserved binding site of the splicing regulator Tra2-β in PKD1 intron 30. Reverse-transcription PCR (RT-PCR) of peripheral blood leukocytes (PBL) from an ADPKD patient homozygous for rs3874648-A identified an atypical PKD1 splice form. Functional analyses demonstrated that rs3874648-A allele increased Tra2-β binding affinity and activated a cryptic acceptor splice-site, causing a frameshift that introduced a premature stop codon in mRNA, thereby decreasing PKD1 full-length transcript level. PKD1 transcript levels were lower in PBL from rs3874648-G/A carriers than in rs3874648-G/G homozygotes in a small cohort of normal individuals and patients with PKD2 inactivating mutations. Our findings indicate that rs3874648G > A is a PKD1 expression modifier attenuating PKD1 expression through Tra2-β, while the derived G allele advantageously maintains PKD1 expression and is predominant in all subpopulations.

The evolutionary dynamics that retain long neutral genomic sequences in face of indel deletion bias: a model and its application to human introns

Insertions and deletions (indels) of short DNA segments are common evolutionary events. Numerous studies showed that deletions occur more often than insertions in both prokaryotes and eukaryotes. It raises the question why neutral sequences are not eradicated from the genome. We suggest that this is due to a phenomenon we term border-induced selection. Accordingly, a neutral sequence is bordered between conserved regions. Deletions occurring near the borders occasionally protrude to the conserved region and are thereby subject to strong purifying selection. Thus, for short neutral sequences, an insertion bias is expected. Here, we develop a set of increasingly complex models of indel dynamics that incorporate border-induced selection. Furthermore, we show that short conserved sequences within the neutrally evolving sequence help explain: (i) the presence of very long sequences; (ii) the high variance of sequence lengths; and (iii) the possible emergence of multimodality in sequence length distributions. Finally, we fitted our models to the human intron length distribution, as introns are thought to be mostly neutral and bordered by conserved exons. We show that when accounting for the occurrence of short conserved sequences within introns, we reproduce the main features, including the presence of long introns and the multimodality of intron distribution.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Confirmation of the first intronic sequence in the bivalvian mitochondrial genome of Macoma balthica (Linnaeus, 1758)

In 2020, the first male-type mitochondrial genome from the clam Macoma balthica was published. Apart from the unusual doubly uniparental inheritance of mtDNA, scientists observed a unique (over 4k bp long) extension in the middle of the cox2 gene. We have attempted to replicate these data by NGS DNA sequencing and explore further the expression of the long cox2 gene. In our study, we report an even longer cox2 gene (over 5.5 kbp) with no stop codon separating conserved cox2 domains, as well as, based on the rtPCR, a lower relative gene expression pattern of the middle part of the gene (5' = 1; mid = 0.46; 3' = 0.89). Lastly, we sequenced the cox2 gene transcript proving the excision of the intronic sequence.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

Dynamic plastid and mitochondrial genomes in Chaetopeltidales (Chlorophyceae) and characterization of a new chlorophyte taxon

PREMISE

Chaetopeltidales is a poorly characterized order in the Chlorophyceae, with only two plastid and no mitochondrial genomes published. Here we describe a new taxon in Chaetopeltidales, Gormaniella terricola gen. et sp. nov. and characterize both of its organellar genomes.

METHODS

Gormaniella terricola was inadvertently isolated from a surface-sterilized hornwort thallus. Light microscopy was used to characterize its vegetative morphology. Organellar genomes were assembled, annotated, and analyzed using a variety of software packages.

RESULTS

The mitochondrial genome (66,927 bp) represents the first complete mitochondrial genome published for Chaetopeltidales. The chloroplast genome, measuring 428,981 bp, is one of the largest plastid genomes published to date and shares this large size and an incredible number of short, dispersed repeats with the other sequenced chloroplast genomes in Chaetopeltidales. Despite these shared features, the chloroplast genomes of Chaetopeltidales appear to be highly rearranged when compared to one another, with numerous inversions, translocations, and duplications, suggesting a particularly dynamic chloroplast genome. Both the chloroplast and mitochondrial genomes of G. terricola contain a number of mobile group I and group II introns, which appear to have invaded separately. Three of the introns within the mitochondrial genome encode homing endonucleases that are phylogenetically nested within those found in fungi, rather than algae, suggesting a possible case of horizontal gene transfer.

CONCLUSIONS

These results help to shed light on a poorly understood group of algae and their unusual organellar genomes, raising additional questions about the unique patterns of genome evolution within Chaetopeltidales.

Splicing modulators elicit global translational repression by condensate-prone proteins translated from introns

Chemical splicing modulators that bind to the spliceosome have provided an attractive avenue for cancer treatment. Splicing modulators induce accumulation and subsequent translation of a subset of intron-retained mRNAs. However, the biological effect of proteins containing translated intron sequences remains unclear. Here, we identify a number of truncated proteins generated upon treatment with the splicing modulator spliceostatin A (SSA) via genome-wide ribosome profiling and bio-orthogonal noncanonical amino acid tagging (BONCAT) mass spectrometry. A subset of these truncated proteins has intrinsically disordered regions, forms insoluble cellular condensates, and triggers the proteotoxic stress response through c-Jun N-terminal kinase (JNK) phosphorylation, thereby inhibiting the mTORC1 pathway. In turn, this reduces global translation. These findings indicate that creating an overburden of condensate-prone proteins derived from introns represses translation and prevents further production of harmful truncated proteins. This mechanism appears to contribute to the antiproliferative and proapoptotic activity of splicing modulators.

RHD Genotypes in a Chinese Cohort of Pregnant Women

RHD variants in D¯ Chinese pregnant women arose difficulties in management during pregnancy. Therefore, this study aims to precisely manage D¯ pregnant women by evaluating the spectrum of RHD mutations in D¯ pregnant women and getting insight into the possible rare alleles of RHD. A total of 76 D¯ pregnant women were analyzed by performing polymerase chain reactions with sequence-specific primers (PCR-SSP), the 10 RHD exons Sanger sequencing, RHD zygosity detection, and mRNA sequencing (mRNA-seq). About 40% of alleles are variations of RHD, including RHD 1227A homozygous, RHD-CE(2-9)-D, et al. Therefore, we developed a molecular diagnostic strategy for Chinese women, and most D¯ pregnant women can be diagnosed with this simple decision tree. After RHD genotyping for D¯ pregnancy women, we eliminated at least 15% unnecessary ante- and postpartum injections of Rh immunoglobulin (RhIG). As the first pedigree study and the first functional analysis under physiological conditions, mRNA-seq revealed that c.336-1G>A mutation mainly led to the inclusion of the intron 2, which indirectly explained the D¯ phenotype in this family. We also developed a robust protocol for determining fetal RhD status from maternal plasma. All 31 fetuses were predicted as RhD positive and confirmed the RhD status after birth.

An overview of RNA splicing and functioning of splicing factors in land plant chloroplasts

RNA splicing refers to a process by which introns of a pre-mRNA are excised and the exons at both ends are joined together. Chloroplast introns are inherently self-splicing ribozymes, but over time, they have lost self-splicing ability due to the degeneration of intronic elements. Thus, the splicing of chloroplast introns relies heavily on nuclear-encoded splicing factors, which belong to diverse protein families. Different splicing factors and their shared intron targets are supposed to form ribonucleoprotein particles (RNPs) to facilitate intron splicing. As characterized in a previous review, around 14 chloroplast intron splicing factors were identified until 2010. However, only a few genetic and biochemical evidence has shown that these splicing factors are required for the splicing of one or several introns. The roles of splicing factors are generally believed to facilitate intron folding; however, the precise role of each protein in RNA splicing remains ambiguous. This may be because the precise binding site of most of these splicing factors remains unexplored. In the last decade, several new splicing factors have been identified. Also, several splicing factors were found to bind to specific sequences within introns, which enhanced the understanding of splicing factors. Here, we summarize recent progress on the splicing factors in land plant chloroplasts and discuss their possible roles in chloroplast RNA splicing based on previous studies.

First Identification of Rare Exonic and Deep Intronic Splice-Altering Variants in Patients With Beta-Sarcoglycanopathy

BACKGROUND

The precise genetic diagnosis of a sarcoglycanopathy or dystrophinopathy is sometimes extremely challenging, as pathogenic non-coding variants and/or complex structural variants do exist in DMD or sarcoglycan genes. This study aimed to determine the genetic diagnosis of three patients from two unrelated families with a suspected sarcoglycanopathy or dystrophinopathy based on their clinical, radiological, and pathological features, for whom routine genomic detection approaches failed to yield a definite genetic diagnosis.

METHODS

Muscle-derived reverse transcription-polymerase chain reaction analysis and/or TA cloning of DMD, SGCA, SGCB, SGCD, and SGCG mRNA were performed to identify aberrant transcripts. Genomic Sanger sequencing around the aberrant transcripts was performed to detect possible splice-altering variants. Bioinformatic and segregation studies of the detected genomic variants were performed in both families.

RESULTS

In patients F1-II1 and F1-II2, we identified two novel pathogenic compound heterozygous variants in SGCB. One is a deep intronic splice-altering variant (DISV), c.243 + 1558C > T in intron 2 causing the activation of an 87-base pair (bp) pseudoexon, and the other one is a non-canonical splicing site variant, c.243 + 6T > A leading to the partial intron inclusion of 10-bp sequence. A novel DISV, c.243 + 1576C > G causing a 106-bp pseudoexon activation, and a nonsense variant in SGCB were identified in compound heterozygous state in patient F2-II1. Unexpectedly, the predicted nonsense variant, c.334C > T in exon 3, created a new donor splice site in exon 3 that was stronger than the natural one, resulting in a 97-bp deletion of exon 3 (r.333_429del).

CONCLUSION

This is the first identification of rare exonic and DISVs in the SGCB gene.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

Information-theoretic approach in allometric scaling relations of DNA and proteins

Allometric scaling relations can be observed in between molecular parameters. Hence, we looked for presence of such relation among sizes (i.e., lengths) of proteins and genes. Protein lengths exist in the literature as the number of amino acids. They can also be derived from the mRNA lengths. Here, we looked for allometric scaling relation by using such data and simultaneously, the data was compared with the sizes of genes and proteins that were obtained from our modified information-theoretic approach. Results implied presence of scaling relation in the calculated results. This was expected due to the implemented modification in the information-theoretic calculation. Relation in the literature-based data was lacking high goodness of fit value. It could be due to physical factors and selective pressures, which ended up in deviations of the literature-sourced values from those in the model. Genome size is correlated with cell size. Intracellular volume, which is related to the DNA size, would require certain number of proteins, the sizes of which can therefore be correlated with the protein sizes. Cell sizes, genome sizes, and average protein and gene sizes, along with the number of proteins, namely the expression levels of the genes, are the physical factors, and the molecular factors influence those physical factors. The selective pressures on those can act through the connection between those physical factors and limit the dynamic ranges. Biological measures could be prone to such forces and are likely to deviate from expected models, regardless of the validity of assumptions, unless those are also implemented in the models. Yet, present discrepancies could be pointing at the need for model improvement, data imperfection, invalid assumptions, etc. Still, current work highlights possible use of information-theoretic approach in allometric scaling relations' studies.

Evidence for the Association between the Intronic Haplotypes of Ionotropic Glutamate Receptors and First-Episode Schizophrenia

The heritable component of schizophrenia (SCH) as a polygenic trait is represented by numerous variants from a heterogeneous group of genes each contributing a relatively small effect. Various SNPs have already been found and analyzed in genes encoding the NMDAR subunits. However, less is known about genetic variations of genes encoding the AMPA and kainate receptor subunits. We analyzed sixteen iGluR genes in full length to determine the sequence variability of iGluR genes. Our aim was to describe the rate of genetic variability, its distribution, and the co-occurrence of variants and to identify new candidate risk variants or haplotypes. The cumulative effect of genetic risk was then estimated using a simple scoring model. GRIN2A-B, GRIN3A-B, and GRIK4 genes showed significantly increased genetic variation in SCH patients. The fixation index statistic revealed eight intronic haplotypes and an additional four intronic SNPs within the sequences of iGluR genes associated with SCH (p < 0.05). The haplotypes were used in the proposed simple scoring model and moreover as a test for genetic predisposition to schizophrenia. The positive likelihood ratio for the scoring model test reached 7.11. We also observed 41 protein-altering variants (38 missense variants, four frameshifts, and one nonsense variant) that were not significantly associated with SCH. Our data suggest that some intronic regulatory regions of iGluR genes and their common variability are among the components from which the genetic predisposition to SCH is composed.

Evidence for exon shuffling is sensitive to model choice

The exon shuffling theory posits that intronic recombination creates new domain combinations, facilitating the evolution of novel protein function. This theory predicts that introns will be preferentially situated near domain boundaries. Many studies have sought evidence for exon shuffling by testing the correspondence between introns and domain boundaries against chance intron positioning. Here, we present an empirical investigation of how the choice of null model influences significance. Although genome-wide studies have used a uniform null model, exclusively, more realistic null models have been proposed for single gene studies. We extended these models for genome-wide analyses and applied them to 21 metazoan and fungal genomes. Our results show that compared with the other two models, the uniform model does not recapitulate genuine exon lengths, dramatically underestimates the probability of chance agreement, and overestimates the significance of intron-domain correspondence by as much as 100 orders of magnitude. Model choice had much greater impact on the assessment of exon shuffling in fungal genomes than in metazoa, leading to different evolutionary conclusions in seven of the 16 fungal genomes tested. Genome-wide studies that use this overly permissive null model may exaggerate the importance of exon shuffling as a general mechanism of multidomain evolution.

The Gene Rearrangement, Loss, Transfer, and Deep Intronic Variation in Mitochondrial Genomes of Conidiobolus

The genus Conidiobolus s.s. was newly delimited from Conidiobolus s.l. In order to gain insight into its mitochondrial genetic background, this study sequenced six mitochondrial genomes of the genus Conidiobolus s.s. These mitogenomes were all composed of circular DNA molecules, ranging from 29,253 to 48,417 bp in size and from 26.61 to 27.90% in GC content. The order and direction for 14 core protein-coding genes (PCGs) were identical, except for the atp8 gene lost in Conidiobolus chlamydosporus, Conidiobolus polyspermus, and Conidiobolus polytocus, and rearranged in the other Conidiobolus s.s. species. Besides, the atp8 gene split the cox1 gene in Conidiobolus taihushanensis. Phylogenomic analysis based on the 14 core PCGs confirmed that all Conidiobolus s.s. species formed a monophyly in the Entomophthoromycotina lineage. The number and length of introns were the main factors contributing to mitogenomic size, and deep variations and potential transfer were detected in introns. In addition, gene transfer occurred between the mitochondrial and nuclear genomes. This study promoted the understanding of the evolution and phylogeny of the Conidiobolus s.s. genus.

Expression Elements Derived From Plant Sequences Provide Effective Gene Expression Regulation and New Opportunities for Plant Biotechnology Traits

Plant biotechnology traits provide a means to increase crop yields, manage weeds and pests, and sustainably contribute to addressing the needs of a growing population. One of the key challenges in developing new traits for plant biotechnology is the availability of expression elements for efficacious and predictable transgene regulation. Recent advances in genomics, transcriptomics, and computational tools have enabled the generation of new expression elements in a variety of model organisms. In this study, new expression element sequences were computationally generated for use in crops, starting from native Arabidopsis and maize sequences. These elements include promoters, 5' untranslated regions (5' UTRs), introns, and 3' UTRs. The expression elements were demonstrated to drive effective transgene expression in stably transformed soybean plants across multiple tissues types and developmental stages. The expressed transcripts were characterized to demonstrate the molecular function of these expression elements. The data show that the promoters precisely initiate transcripts, the introns are effectively spliced, and the 3' UTRs enable predictable processing of transcript 3' ends. Overall, our results indicate that these new expression elements can recapitulate key functional properties of natural sequences and provide opportunities for optimizing the expression of genes in future plant biotechnology traits.

Alu RNA and their roles in human disease states

Alu RNA are implicated in the poor prognosis of several human disease states. These RNA are transcription products of primate specific transposable elements called Alu elements. These elements are extremely abundant, comprising over 10% of the human genome, and 100 to 1000 cytoplasmic copies of Alu RNA per cell. Alu RNA do not have a single universal functional role aside from selfish self-propagation. Despite this, Alu RNA have been found to operate in a diverse set of translational and transcriptional mechanisms. This review will focus on the current knowledge of Alu RNA involved in human disease states and known mechanisms of action. Examples of Alu RNA that are transcribed in a variety of contexts such as introns, mature mRNA, and non-coding transcripts will be discussed. Past and present challenges in studying Alu RNA, and the future directions of Alu RNA in basic and clinical research will also be examined.

Eukaryotic Genomes Show Strong Evolutionary Conservation of k-mer Composition and Correlation Contributions between Introns and Intergenic Regions

Several strongly conserved DNA sequence patterns in and between introns and intergenic regions (IIRs) consisting of short tandem repeats (STRs) with repeat lengths <3 bp have already been described in the kingdom of Animalia. In this work, we expanded the search and analysis of conserved DNA sequence patterns to a wider range of eukaryotic genomes. Our aims were to confirm the conservation of these patterns, to support the hypothesis on their functional constraints and/or the identification of unknown patterns. We pairwise compared genomic DNA sequences of genes, exons, CDS, introns and intergenic regions of 34 Embryophyta (land plants), 30 Protista and 29 Fungi using established k-mer-based (alignment-free) comparison methods. Additionally, the results were compared with values derived for Animalia in former studies. We confirmed strong correlations between the sequence structures of IIRs spanning over the entire domain of Eukaryotes. We found that the high correlations within introns, intergenic regions and between the two are a result of conserved abundancies of STRs with repeat units ≤2 bp (e.g., (AT)n). For some sequence patterns and their inverse complementary sequences, we found a violation of equal distribution on complementary DNA strands in a subset of genomes. Looking at mismatches within the identified STR patterns, we found specific preferences for certain nucleotides stable over all four phylogenetic kingdoms. We conclude that all of these conserved patterns between IIRs indicate a shared function of these sequence structures related to STRs.

Systematic Identification and Functional Validation of New snoRNAs in Human Muscle Progenitors

Small non-coding RNAs (sncRNAs) represent an important class of regulatory RNAs involved in the regulation of transcription, RNA splicing or translation. Among these sncRNAs, small nucleolar RNAs (snoRNAs) mostly originate from intron splicing in humans and are central to posttranscriptional regulation of gene expression. However, the characterization of the complete repertoire of sncRNAs in a given cellular context and the functional annotation of the human transcriptome are far from complete. Here, we report the large-scale identification of sncRNAs in the size range of 50 to 200 nucleotides without a priori on their biogenesis, structure and genomic origin in the context of normal human muscle cells. We provided a complete set of experimental validation of novel candidate snoRNAs by evaluating the prerequisites for their biogenesis and functionality, leading to their validation as genuine snoRNAs. Interestingly, we also found intergenic snoRNAs, which we showed are in fact integrated into candidate introns of unannotated transcripts or degraded by the Nonsense Mediated Decay pathway. Hence, intergenic snoRNAs represent a new type of landmark for the identification of new transcripts that have gone undetected because of low abundance or degradation after the release of the snoRNA.

Noncoding sequence variants define a novel regulatory element in the first intron of the N-acetylglutamate synthase gene

N-acetylglutamate synthase deficiency is an autosomal recessive urea cycle disorder caused either by decreased expression of the NAGS gene or defective NAGS enzyme resulting in decreased production of N-acetylglutamate (NAG), an allosteric activator of carbamylphosphate synthetase 1 (CPS1). NAGSD is the only urea cycle disorder that can be effectively treated with a single drug, N-carbamylglutamate (NCG), a stable NAG analog, which activates CPS1 to restore ureagenesis. We describe three patients with NAGSD due to four novel noncoding sequence variants in the NAGS regulatory regions. All three patients had hyperammonemia that resolved upon treatment with NCG. Sequence variants NM_153006.2:c.427-222G>A and NM_153006.2:c.427-218A>C reside in the 547 bp-long first intron of NAGS and define a novel NAGS regulatory element that binds retinoic X receptor α. Sequence variants NC_000017.10:g.42078967A>T (NM_153006.2:c.-3065A>T) and NC_000017.10:g.42078934C>T (NM_153006.2:c.-3098C>T) reside in the NAGS enhancer, within known HNF1 and predicted glucocorticoid receptor binding sites, respectively. Reporter gene assays in HepG2 and HuH-7 cells demonstrated that all four substitutions could result in reduced expression of NAGS. These findings show that analyzing noncoding regions of NAGS and other urea cycle genes can reveal molecular causes of disease and identify novel regulators of ureagenesis.

Contribution of Noncanonical Splice Variants to TTN Truncating Variant Cardiomyopathy

BACKGROUND

Heterozygous TTN truncating variants cause 10% to 20% of idiopathic dilated cardiomyopathy (DCM). Although variants which disrupt canonical splice signals (ie, invariant dinucleotide of the splice donor site, invariant dinucleotide of the splice acceptor site) at exon-intron junctions are readily recognized as TTN truncating variants, the effects of other nearby sequence variations on splicing and their contribution to disease is uncertain.

METHODS

Rare variants of unknown significance located in the splice regions of highly expressed TTN exons from 203 DCM cases, 3329 normal subjects, and clinical variant databases were identified. The effects of these variants on splicing were assessed using an in vitro splice assay.

RESULTS

Splice-altering variants of unknown significance were enriched in DCM cases over controls and present in 2% of DCM patients (P=0.002). Application of this method to clinical variant databases demonstrated 20% of similar variants of unknown significance in TTN splice regions affect splicing. Noncanonical splice-altering variants were most frequently located at position +5 of the donor site (P=4.4×10⁷) and position -3 of the acceptor site (P=0.002). SpliceAI, an emerging in silico prediction tool, had a high positive predictive value (86%-95%) but poor sensitivity (15%-50%) for the detection of splice-altering variants. Alternate exons spliced out of most TTN transcripts frequently lacked the consensus base at +5 donor and -3 acceptor positions.

CONCLUSIONS

Noncanonical splice-altering variants in TTN explain 1-2% of DCM and offer a 10-20% increase in the diagnostic power of TTN sequencing in this disease. These data suggest rules that may improve efforts to detect splice-altering variants in other genes and may explain the low percent splicing observed for many alternate TTN exons.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

Proximity to the Promoter and Terminator Regions Regulates the Transcription Enhancement Potential of an Intron

An evolutionarily conserved feature of introns is their ability to enhance expression of genes that harbor them. Introns have been shown to regulate gene expression at the transcription and post-transcription level. The general perception is that a promoter-proximal intron is most efficient in enhancing gene expression and the effect diminishes with the increase in distance from the promoter. Here we show that the intron regains its positive influence on gene expression when in proximity to the terminator. We inserted ACT1 intron into different positions within IMD4 and INO1 genes. Transcription Run-On (TRO) analysis revealed that the transcription of both IMD4 and INO1 was maximal in constructs with a promoter-proximal intron and decreased with the increase in distance of the intron from the promoter. However, activation was partially restored when the intron was placed close to the terminator. We previously demonstrated that the promoter-proximal intron stimulates transcription by affecting promoter directionality through gene looping-mediated recruitment of termination factors in the vicinity of the promoter region. Here we show that the terminator-proximal intron also enhances promoter directionality and results in compact gene architecture with the promoter and terminator regions in close physical proximity. Furthermore, we show that both the promoter and terminator-proximal introns facilitate assembly or stabilization of the preinitiation complex (PIC) on the promoter. On the basis of these findings, we propose that proximity to both the promoter and the terminator regions affects the transcription regulatory potential of an intron, and the terminator-proximal intron enhances transcription by affecting both the assembly of preinitiation complex and promoter directionality.

Efficient expression and function of a receptor-like kinase in wheat powdery mildew defence require an intron-located MYB binding site

The LRK10-like receptor kinases (LRK10L-RLKs) are ubiquitously present in higher plants, but knowledge of their expression and function is still limited. Here, we report expression and functional analysis of TtdLRK10L-1, a typical LRK10L-RLK in durum wheat (Triticum turgidum L. ssp. durum). The introns of TtdLRK10L-1 contained multiple kinds of predicted cis-elements. To investigate the potential effect of these cis-elements on TtdLRK10L-1 expression and function, two types of transgenic wheat lines were prepared, which expressed a GFP-tagged TtdLRK10L-1 protein (TtdLRK10L-1:GFP) from the cDNA or genomic DNA (gDNA) sequence of TtdLRK10L-1 under the native promoter. TtdLRK10L-1:GFP expression was up-regulated by the powdery mildew pathogen Blumeria graminis f. sp. tritici (Bgt) in both types of transgenic plants, with the scale of the elevation being much stronger in the gDNA lines. Both types of transgenic plants exhibited enhanced resistance to Bgt infection relative to wild type control. Notably, the Bgt defence activated in the gDNA lines was significantly stronger than that in the cDNA lines. Further analysis revealed that a putative MYB transcription factor binding site (MYB-BS, CAGTTA) located in TtdLRK10L-1 intron I was critical for the efficient expression and function of TtdLRK10L-1 in Bgt defence. This MYB-BS could also increase the activity of a superpromoter widely used in ectopic gene expression studies in plants. Together, our results deepen the understanding of the expression and functional characteristics of LRK10L-RLKs. TtdLRK10L-1 is likely useful for further dissecting the molecular processes underlying wheat defence against Bgt and for developing Bgt resistant wheat crops.

Gene Architecture Facilitates Intron-Mediated Enhancement of Transcription

Introns impact several vital aspects of eukaryotic organisms like proteomic plasticity, genomic stability, stress response and gene expression. A role for introns in the regulation of gene expression at the level of transcription has been known for more than thirty years. The molecular basis underlying the phenomenon, however, is still not entirely clear. An important clue came from studies performed in budding yeast that indicate that the presence of an intron within a gene results in formation of a multi-looped gene architecture. When looping is defective, these interactions are abolished, and there is no enhancement of transcription despite normal splicing. In this review, we highlight several potential mechanisms through which looping interactions may enhance transcription. The promoter-5′ splice site interaction can facilitate initiation of transcription, the terminator-3′ splice site interaction can enable efficient termination of transcription, while the promoter-terminator interaction can enhance promoter directionality and expedite reinitiation of transcription. Like yeast, mammalian genes also exhibit an intragenic interaction of the promoter with the gene body, especially exons. Such promoter-exon interactions may be responsible for splicing-dependent transcriptional regulation. Thus, the splicing-facilitated changes in gene architecture may play a critical role in regulation of transcription in yeast as well as in higher eukaryotes.

Atypical acute myeloid leukemia-specific transcripts generate shared and immunogenic MHC class-I-associated epitopes

Acute myeloid leukemia (AML) has not benefited from innovative immunotherapies, mainly because of the lack of actionable immune targets. Using an original proteogenomic approach, we analyzed the major histocompatibility complex class I (MHC class I)-associated immunopeptidome of 19 primary AML samples and identified 58 tumor-specific antigens (TSAs). These TSAs bore no mutations and derived mainly (86%) from supposedly non-coding genomic regions. Two AML-specific aberrations were instrumental in the biogenesis of TSAs, intron retention, and epigenetic changes. Indeed, 48% of TSAs resulted from intron retention and translation, and their RNA expression correlated with mutations of epigenetic modifiers (e.g., DNMT3A). AML TSA-coding transcripts were highly shared among patients and were expressed in both blasts and leukemic stem cells. In AML patients, the predicted number of TSAs correlated with spontaneous expansion of cognate T cell receptor clonotypes, accumulation of activated cytotoxic T cells, immunoediting, and improved survival. These TSAs represent attractive targets for AML immunotherapy.

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.