Rapid, transcript-specific changes in splicing in response to environmental stress

The spliceosome impacts morphogenesis in the human fungal pathogen Candida albicans

At human body temperature, the fungal pathogen Candida albicans can transition from yeast to filamentous morphologies in response to host-relevant cues. Additionally, elevated temperatures encountered during febrile episodes can independently induce C. albicans filamentation. However, the underlying genetic pathways governing this developmental transition in response to elevated temperatures remain largely unexplored. Here, we conducted a functional genomic screen to unravel the genetic mechanisms orchestrating C. albicans filamentation specifically in response to elevated temperature, implicating 45% of genes associated with the spliceosome or pre-mRNA splicing in this process. Employing RNA-Seq to elucidate the relationship between mRNA splicing and filamentation, we identified greater levels of intron retention in filaments compared to yeast, which correlated with reduced expression of the affected genes. Intriguingly, homozygous deletion of a gene encoding a spliceosome component important for filamentation (PRP19) caused even greater levels of intron retention compared with wild type and displayed globally dysregulated gene expression. This suggests that intron retention is a mechanism for fine-tuning gene expression during filamentation, with perturbations of the spliceosome exacerbating this process and blocking filamentation. Overall, this study unveils a novel biological process governing C. albicans filamentation, providing new insights into the complex regulation of this key virulence trait.IMPORTANCEFungal pathogens such as Candida albicans can cause serious infections with high mortality rates in immunocompromised individuals. When C. albicans is grown at temperatures encountered during human febrile episodes, yeast cells undergo a transition to filamentous cells, and this process is key to its virulence. Here, we expanded our understanding of how C. albicans undergoes filamentation in response to elevated temperature and identified many genes involved in mRNA splicing that positively regulate filamentation. Through transcriptome analyses, we found that intron retention is a mechanism for fine-tuning gene expression in filaments, and perturbation of the spliceosome exacerbates intron retention and alters gene expression substantially, causing a block in filamentation. This work adds to the growing body of knowledge on the role of introns in fungi and provides new insights into the cellular processes that regulate a key virulence trait in C. albicans.

Splicing factor Prp18p promotes genome-wide fidelity of consensus 3'-splice sites

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is challenging considering the low complexity of the 3'SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3'SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3'SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3'SS. The role of Prp18p in 3'SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3'SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3'SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3'SS sequences.

Single-cell analysis reveals region-heterogeneous responses in rhesus monkey spinal cord with complete injury

Spinal cord injury (SCI) leads to severe sensory and motor dysfunction below the lesion. However, the cellular dynamic responses and heterogeneity across different regions below the lesion remain to be elusive. Here, we used single-cell transcriptomics to investigate the region-related cellular responses in female rhesus monkeys with complete thoracic SCI from acute to chronic phases. We found that distal lumbar tissue cells were severely impacted, leading to degenerative microenvironments characterized by disease-associated microglia and oligodendrocytes activation alongside increased inhibitory interneurons proportion following SCI. By implanting scaffold into the injury sites, we could improve the injury microenvironment through glial cells and fibroblast regulation while remodeling spared lumbar tissues via reduced inhibitory neurons proportion and improved phagocytosis and myelination. Our findings offer crucial pathological insights into the spared distal tissues and proximal tissues after SCI, emphasizing the importance of scaffold-based treatment approaches targeting heterogeneous microenvironments.

Moderate maternal nutrient reduction in pregnancy alters fatty acid oxidation and RNA splicing in the nonhuman primate fetal liver

Fetal liver tissue collected from a nonhuman primate (NHP) baboon model of maternal nutrient reduction (MNR) at four gestational time points (90, 120, 140, and 165 days gestation [dG], term in the baboon is ∼185 dG) was used to quantify MNR effects on the fetal liver transcriptome. 28 transcripts demonstrated different expression patterns between MNR and control livers during the second half of gestation, a developmental period when the fetus undergoes rapid weight gain and fat accumulation. Differentially expressed transcripts were enriched for fatty acid oxidation and RNA splicing-related pathways. Increased RNA splicing activity in MNR was reflected in greater abundances of transcript splice variant isoforms in the MNR group. It can be hypothesized that the increase in splice variants is deployed in an effort to adapt to the poor in utero environment and ensure near-normal development and energy metabolism. This study is the first to study developmental programming across four critical gestational stages during primate fetal liver development and reveals a potentially novel cellular response mechanism mediating fetal programming in response to MNR.

SUMOylation regulates pre-mRNA splicing to overcome DNA damage in fungi

Pathogenic fungi are subject to DNA damage stress derived from host immune responses during infection. Small ubiquitin-like modifier (SUMO) modification and precursor (pre)-mRNA splicing are both involved in DNA damage response (DDR). However, the mechanisms of how SUMOylation and splicing coordinated in DDR remain largely unknown. Combining with biochemical analysis, RNA-Seq method, and biological analysis, we report that SUMO pathway participates in DDR and virulence in Fusarium graminearum, a causal agent of Fusarium head blight of cereal crops world-wide. Interestingly, a key transcription factor FgSR is SUMOylated upon DNA damage stress. SUMOylation regulates FgSR nuclear-cytoplasmic partitioning and its phosphorylation by FgMec1, and promotes its interaction with chromatin remodeling complex SWI/SNF for activating the expression of DDR-related genes. Moreover, the SWI/SNF complex was found to further recruit splicing-related NineTeen Complex, subsequently modulates pre-mRNA splicing during DDR. Our findings reveal a novel function of SUMOylation in DDR by regulating a transcription factor to orchestrate gene expression and pre-mRNA splicing to overcome DNA damage during the infection of F. graminearum, which advances the understanding of the delicate regulation of DDR by SUMOylation in pathogenic fungi, and extends the knowledge of cooperation of SUMOylation and pre-mRNA splicing in DDR in eukaryotes.

Decline of regenerative potential of old muscle stem cells: contribution to muscle aging

Muscle stem cells (MuSCs) are required for life-long muscle regeneration. In general, aging has been linked to a decline in the numbers and the regenerative potential of MuSCs. Muscle regeneration depends on the proper functioning of MuSCs, which is itself dependent on intricate interactions with its niche components. Aging is associated with both cell-intrinsic and niche-mediated changes, which can be the result of transcriptional, posttranscriptional, or posttranslational alterations in MuSCs or in the components of their niche. The interplay between cell intrinsic alterations in MuSCs and changes in the stem cell niche environment during aging and its impact on the number and the function of MuSCs is an important emerging area of research. In this review, we discuss whether the decline in the regenerative potential of MuSCs with age is the cause or the consequence of aging skeletal muscle. Understanding the effect of aging on MuSCs and the individual components of their niche is critical to develop effective therapeutic approaches to diminish or reverse the age-related defects in muscle regeneration.

Spliceosomal mutations decouple 3' splice site fidelity from cellular fitness

The fidelity of splice site selection is thought to be critical for proper gene expression and cellular fitness. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is a daunting task considering the low complexity of the 3'SS consensus sequence YAG. Here we show that inactivating the near-essential splicing factor Prp18p results in a global activation of alternative 3'SS, many of which harbor sequences that highly diverge from the YAG consensus, including some highly unusual non-AG 3'SS. We show that the role of Prp18p in 3'SS fidelity is promoted by physical interactions with the essential splicing factors Slu7p and Prp8p and synergized by the proofreading activity of the Prp22p helicase. Strikingly, structure-guided point mutations that disrupt Prp18p-Slu7p and Prp18p-Prp8p interactions mimic the loss of 3'SS fidelity without any impact on cellular growth, suggesting that accumulation of incorrectly spliced transcripts does not have a major deleterious effect on cellular viability. These results show that spliceosomes exhibit remarkably relaxed fidelity in the absence of Prp18p, and that new 3'SS sampling can be achieved genome-wide without a major negative impact on cellular fitness, a feature that could be used during evolution to explore new productive alternative splice sites.

Mosaic translation hypothesis: chimeric polypeptides produced via multiple ribosomal frameshifting as a basis for adaptability

How many different proteins can be produced from a single spliced transcript? Genome annotation projects overlook the coding potential of reading frames other than that of the reference open reading frames (refORFs). Recently, alternative open reading frames (altORFs) and their translational products, alternative proteins, have been shown to carry out important functions in various organisms. AltORFs overlapping refORFs or other altORFs in a different reading frame may be involved in one fundamental mechanism so far overlooked. A few years ago, it was proposed that altORFs may act as building blocks for chimeric (mosaic) polypeptides, which are produced via multiple ribosomal frameshifting events from a single mature transcript. We adopt terminology from that earlier discussion and call this mechanism mosaic translation. This way of extracting and combining genetic information may significantly increase proteome diversity. Thus, we hypothesize that this mechanism may have contributed to the flexibility and adaptability of organisms to a variety of environmental conditions. Specialized ribosomes acting as sensors probably played a central role in this process. Importantly, mosaic translation may be the main source of protein diversity in genomes that lack alternative splicing. The idea of mosaic translation is a testable hypothesis, although its direct demonstration is challenging. Should mosaic translation occur, we would currently highly underestimate the complexity of translation mechanisms and thus the proteome.

Phylogenomics of a Saccharomyces cerevisiae cocoa strain reveals adaptation to a West African fermented food population

Various yeast strains have been proposed as candidate starter cultures for cocoa fermentation, especially strains of Saccharomyces cerevisiae. In the current study, the genome of the cocoa strain S. cerevisiae IMDO 050523 was unraveled based on a combination of long- and short-read sequencing. It consisted of 16 nuclear chromosomes and a mitochondrial chromosome, which were organized in 20 contigs, with only two small gaps. A phylogenomic analysis of this genome together with another 105 S cerevisiae genomes, among which 20 from cocoa strains showed a geographical distribution of the latter, including S. cerevisiae IMDO 050523. Its genome clustered together with that of a West African fermented food population, indicating a wider adaptation to West African food niches than cocoa. Furthermore, S. cerevisiae IMDO 050523 contained genetic signatures involved in sucrose hydrolysis, pectin degradation, osmotolerance, and conserved amino acid changes in key ester-producing enzymes that could point toward specific niche adaptations.

Transcript-specific determinants of pre-mRNA splicing revealed through in vivo kinetic analyses of the 1st and 2nd chemical steps

We generate high-precision measurements of the in vivo rates of both chemical steps of pre-mRNA splicing across the genome-wide complement of substrates in yeast by coupling metabolic labeling, multiplexed primer-extension sequencing, and kinetic modeling. We demonstrate that the rates of intron removal vary widely, splice-site sequences are primary determinants of 1st step but have little apparent impact on 2nd step rates, and the 2nd step is generally faster than the 1st step. Ribosomal protein genes (RPGs) are spliced faster than non-RPGs at each step, and RPGs share evolutionarily conserved properties that may contribute to their faster splicing. A genetic variant defective in the 1st step of the pathway reveals a genome-wide defect in the 1st step but an unexpected, transcript-specific change in the 2nd step. Our work demonstrates that extended co-transcriptional association is an important determinant of splicing rate, a conclusion at odds with recent claims of ultra-fast splicing.

Splicing inactivation generates hybrid mRNA-snoRNA transcripts targeted by cytoplasmic RNA decay

Many small nucleolar RNAs (snoRNA)s are processed from introns of host genes, but the importance of splicing for proper biogenesis and the fate of the snoRNAs is not well understood. Here, we show that inactivation of splicing factors or mutation of splicing signals leads to the accumulation of partially processed hybrid messenger RNA-snoRNA (hmsnoRNA) transcripts. hmsnoRNAs are processed to the mature 3' ends of the snoRNAs by the nuclear exosome and bound by small nucleolar ribonucleoproteins. hmsnoRNAs are unaffected by translation-coupled RNA quality-control pathways, but they are degraded by the major cytoplasmic exonuclease Xrn1p, due to their messenger RNA (mRNA)-like 5' extensions. These results show that completion of splicing is required to promote complete and accurate processing of intron-encoded snoRNAs and that splicing defects lead to degradation of hybrid mRNA-snoRNA species by cytoplasmic decay, underscoring the importance of splicing for the biogenesis of intron-encoded snoRNAs.

Intron-mediated induction of phenotypic heterogeneity

Intragenic regions that are removed during maturation of the RNA transcript—introns—are universally present in the nuclear genomes of eukaryotes¹. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein genes that differ primarily in their introns^2,3. Although studies have shed light on the role of ribosomal protein introns under stress and starvation^4–6, understanding the contribution of introns to ribosome regulation remains challenging. Here, by combining isogrowth profiling⁷ with single-cell protein measurements⁸, we show that introns can mediate inducible phenotypic heterogeneity that confers a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, which is mediated by an intron in the 5′ untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low levels of Rps22B protein result in prolonged survival under sustained starvation, whereas high levels of Rps22B enable cells to grow faster after transient starvation. Furthermore, yeasts growing at high concentrations of sugar, similar to those in ripe grapes, exhibit bimodal expression of Rps22B when approaching the stationary phase. Differential intron-mediated regulation of ribosomal protein genes thus provides a way to diversify the population when starvation threatens in natural environments. Our findings reveal a role for introns in inducing phenotypic heterogeneity in changing environments, and suggest that duplicated ribosomal protein genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness⁹.

Experiments in yeast show that introns have a role in inducing phenotypic heterogeneity and that intron-mediated regulation of ribosomal proteins confers a fitness advantage by enabling yeast populations to diversify under nutrient-scarce conditions.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

Heat increases full-length SMN splicing: promise for splice-augmenting therapies for SMA

Spinal muscular atrophy (SMA) is a debilitating neurodegenerative pediatric disease characterized by low levels of the survival motor protein (SMN). Humans have two SMN genes that produce identical SMN proteins, but they differ at a key nucleotide in exon 7 that induces differential mRNA splicing. SMN1 primarily produces full-length SMN protein, but due to the spliceosome's inability to efficiently recognize exon 7, SMN2 transcripts are often truncated. SMA occurs primarily through mutations or deletions in the SMN1 gene; therefore, current therapies use antisense oligonucleotides (ASOs) to target exon 7 inclusion in SMN2 mRNA and promote full-length SMN protein production. Here, we explore additional methods that can target SMN splicing and therapeutically increase full-length SMN protein. We demonstrate that in vitro heat treatment of cells increases exon 7 inclusion and relative abundance of full-length SMN2 mRNA and protein, a response that is modulated through the upregulation of the positive splicing factor TRA2 beta. We also observe that HSP90, but not HSP40 or HSP70, in the heat shock response is essential for SMN2 exon 7 splicing under hyperthermic conditions. Finally, we show that pulsatile heat treatments for one hour in vitro and in vivo are effective in increasing full-length SMN2 levels. These findings suggest that timed interval treatments could be a therapeutic alternative for SMA patients who do not respond to current ASO-based therapies or require a unique combination regimen.

Multiple functions of SWI/SNF chromatin remodeling complex in plant-pathogen interactions

The SWI/SNF chromatin remodeling complex utilizes the energy of ATP hydrolysis to facilitate chromatin access and plays essential roles in DNA-based events. Studies in animals, plants and fungi have uncovered sophisticated regulatory mechanisms of this complex that govern development and various stress responses. In this review, we summarize the composition of SWI/SNF complex in eukaryotes and discuss multiple functions of the SWI/SNF complex in regulating gene transcription, mRNA splicing, and DNA damage response. Our review further highlights the importance of SWI/SNF complex in regulating plant immunity responses and fungal pathogenesis. Finally, the potentials in exploiting chromatin remodeling for management of crop disease are presented.

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Global Transcriptomics Uncovers Distinct Contributions From Splicing Regulatory Proteins to the Macrophage Innate Immune Response

Pathogen sensing via pattern recognition receptors triggers massive reprogramming of macrophage gene expression. While the signaling cascades and transcription factors that activate these responses are well-known, the role of post-transcriptional RNA processing in modulating innate immune gene expression remains understudied. Given their crucial role in regulating pre-mRNA splicing and other RNA processing steps, we hypothesized that members of the SR/hnRNP protein families regulate innate immune gene expression in distinct ways. We analyzed steady state gene expression and alternatively spliced isoform production in ten SR/hnRNP knockdown RAW 264.7 macrophage-like cell lines following infection with the bacterial pathogen Salmonella enterica serovar Typhimurium (Salmonella). We identified thousands of transcripts whose abundance is increased or decreased by SR/hnRNP knockdown in macrophages. Notably, we observed that SR and hnRNP proteins influence expression of different genes in uninfected versus Salmonella-infected macrophages, suggesting functionalization of these proteins upon pathogen sensing. Likewise, we found that knockdown of SR/hnRNPs promoted differential isoform usage (DIU) for thousands of macrophage transcripts and that these alternative splicing changes were distinct in uninfected and Salmonella-infected macrophages. Finally, having observed a surprising degree of similarity between the differentially expressed genes (DEGs) and DIUs in hnRNP K and U knockdown macrophages, we found that hnRNP K and U knockdown macrophages are both more restrictive to Vesicular Stomatitis Virus (VSV), while hnRNP K knockdown macrophages are more permissive to Salmonella Typhimurium. Based on these findings, we conclude that many innate immune genes evolved to rely on one or more SR/hnRNPs to ensure the proper magnitude of their induction, supporting a model wherein pre-mRNA splicing is critical for regulating innate immune gene expression and controlling infection outcomes in macrophages ex vivo.

Reciprocal hemizygosity analysis reveals that the Saccharomyces cerevisiae CGI121 gene affects lag time duration in synthetic grape must

It is standard practice to ferment white wines at low temperatures (10–18°C). However, low temperatures increase fermentation duration and risk of problem ferments, leading to significant costs. The lag duration at fermentation initiation is heavily impacted by temperature; therefore, identification of Saccharomyces cerevisiae genes influencing fermentation kinetics is of interest for winemaking. We selected 28 S. cerevisiae BY4743 single deletants, from a prior list of open reading frames (ORFs) mapped to quantitative trait loci (QTLs) on Chr. VII and XIII, influencing the duration of fermentative lag time. Five BY4743 deletants, Δapt1, Δcgi121, Δclb6, Δrps17a, and Δvma21, differed significantly in their fermentative lag duration compared to BY4743 in synthetic grape must (SGM) at 15 °C, over 72 h. Fermentation at 12.5°C for 528 h confirmed the longer lag times of BY4743 Δcgi121, Δrps17a, and Δvma21. These three candidates ORFs were deleted in S. cerevisiae RM11-1a and S288C to perform single reciprocal hemizygosity analysis (RHA). RHA hybrids and single deletants of RM11-1a and S288C were fermented at 12.5°C in SGM and lag time measurements confirmed that the S288C allele of CGI121 on Chr. XIII, encoding a component of the EKC/KEOPS complex, increased fermentative lag phase duration. Nucleotide sequences of RM11-1a and S288C CGI121 alleles differed by only one synonymous nucleotide, suggesting that intron splicing, codon bias, or positional effects might be responsible for the impact on lag phase duration. This research demonstrates a new role of CGI121 and highlights the applicability of QTL analysis for investigating complex phenotypic traits in yeast.

Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

Does Saccharomyces cerevisiae Require Specific Post-Translational Silencing against Leaky Translation of Hac1up?

HAC1 encodes a key transcription factor that transmits the unfolded protein response (UPR) from the endoplasmic reticulum (ER) to the nucleus and regulates downstream UPR genes in Saccharomyces cerevisiae. In response to the accumulation of unfolded proteins in the ER, Ire1p oligomers splice HAC1 pre-mRNA (HAC1^u) via a non-conventional process and allow the spliced HAC1 (HAC1ⁱ) to be translated efficiently. However, leaky splicing and translation of HAC1^u may occur in non-UPR cells to induce undesirable UPR. To control accidental UPR activation, multiple fail-safe mechanisms have been proposed to prevent leaky HAC1 splicing and translation and to facilitate rapid degradation of translated Hac1^up and Hac1ⁱp. Among proposed regulatory mechanisms is a degron sequence encoded at the 5′ end of the HAC1 intron that silences Hac1^up expression. To investigate the necessity of an intron-encoded degron sequence that specifically targets Hac1^up for degradation, we employed publicly available transcriptomic data to quantify leaky HAC1 splicing and translation in UPR-induced and non-UPR cells. As expected, we found that HAC1^u is only efficiently spliced into HAC1ⁱ and efficiently translated into Hac1ⁱp in UPR-induced cells. However, our analysis of ribosome profiling data confirmed frequent occurrence of leaky translation of HAC1^u regardless of UPR induction, demonstrating the inability of translation fail-safe to completely inhibit Hac1^up production. Additionally, among 32 yeast HAC1 surveyed, the degron sequence is highly conserved by Saccharomyces yeast but is poorly conserved by all other yeast species. Nevertheless, the degron sequence is the most conserved HAC1 intron segment in yeasts. These results suggest that the degron sequence may indeed play an important role in mitigating the accumulation of Hac1^up to prevent accidental UPR activation in the Saccharomyces yeast.

RNA-binding protein syncrip regulates starvation-induced hyperactivity in adult Drosophila

How to respond to starvation determines fitness. One prominent behavioral response is increased locomotor activities upon starvation, also known as Starvation-Induced Hyperactivity (SIH). SIH is paradoxical as it promotes food seeking but also increases energy expenditure. Despite its importance in fitness, the genetic contributions to SIH as a behavioral trait remains unexplored. Here, we examined SIH in the Drosophila melanogaster Genetic Reference Panel (DGRP) and performed genome-wide association studies. We identified 23 significant loci, corresponding to 14 genes, significantly associated with SIH in adult Drosophila. Gene enrichment analyses indicated that genes encoding ion channels and mRNA binding proteins (RBPs) were most enriched in SIH. We are especially interested in RBPs because they provide a potential mechanism to quickly change protein expression in response to environmental challenges. Using RNA interference, we validated the role of syp in regulating SIH. syp encodes Syncrip (Syp), an RBP. While ubiquitous knockdown of syp led to semi-lethality in adult flies, adult flies with neuron-specific syp knockdown were viable and exhibited decreased SIH. Using the Temporal and Regional Gene Expression Targeting (TARGET) system, we further confirmed the role of Syp in adult neurons in regulating SIH. To determine how syp is regulated by starvation, we performed RNA-seq using the heads of flies maintained under either food or starvation conditions. RNA-seq analyses revealed that syp was alternatively spliced under starvation while its expression level was unchanged. We further generated an alternatively-spliced-exon-specific knockout (KO) line and found that KO flies showed reduced SIH. Together, this study demonstrates a significant genetic contribution to SIH as a behavioral trait, identifies syp as a SIH gene, and highlights the significance of RBPs and post-transcriptional processes in the brain in regulating behavioral responses to starvation.

Animals living in the wild often face periods of starvation. How to physiologically and behaviorally respond to starvation is essential for survival. One behavioral response is Starvation-Induced Hyperactivity (SIH). We used the Drosophila melanogaster Genetic Reference Panel, derived from a wild population, to study the genetic basis of SIH. Our results show that there is a significant genetic contribution to SIH in this population, and that genes encoding RNA binding proteins (RBPs) are especially important. Using RNA interference and the TARGET system, we confirmed the role of an RBP Syp in adult neurons in SIH. Using RNA-seq and Western blotting, we found that syp was alternatively spliced under starvation while its expression level was unchanged. Further studies from syp exon-specific knockout flies showed that alternative splicing involving two exons in syp was important for SIH. Together, this study identifies syp as a SIH gene and highlights an essential role of post-transcriptional modification in regulating this behavior.

Alternative mRNA Processing of Innate Response Pathways in Respiratory Syncytial Virus (RSV) Infection

The innate immune response (IIR) involves rapid genomic expression of protective interferons (IFNs) and inflammatory cytokines triggered by intracellular viral replication. Although the transcriptional control of the innate pathway is known in substantial detail, little is understood about the complexity of alternative splicing (AS) and alternative polyadenylation (APA) of mRNAs underlying the cellular IIR. In this study, we applied single-molecule, real-time (SMRT) sequencing with mRNA quantitation using short-read mRNA sequencing to characterize changes in mRNA processing in the epithelial response to respiratory syncytial virus (RSV) replication. Mock or RSV-infected human small-airway epithelial cells (hSAECs) were profiled using SMRT sequencing and the curated transcriptome analyzed by structural and quality annotation of novel transcript isoforms (SQANTI). We identified 113,082 unique isoforms; 28,561 represented full splice matches, and 45% of genes expressed six or greater AS mRNA isoforms. Identification of differentially expressed AS isoforms was accomplished by mapping a short-read RNA sequencing expression matrix to the curated transcriptome, and 905 transcripts underwent differential polyadenylation site analysis enriched in protein secretion, translation, and mRNA degradation. We focused on 355 genes showing differential isoform utilization (DIU), indicating where a new AS isoform becomes a major fraction of mRNA isoforms expressed. In pathway and network enrichment analyses, we observed that DIU transcripts are substantially enriched in cell cycle control and IIR pathways. Interestingly, the RelA/IRF7 innate regulators showed substantial DIU where major transcripts included distinct isoforms with exon occlusion, intron inclusion, and alternative transcription start site utilization. We validated the presence of RelA and IRF7 AS isoforms as well as their induction by RSV using eight isoform-specific RT-PCR assays. These isoforms were identified in both immortalized and primary small-airway epithelial cells. We concluded that the cell cycle and IIR are differentially spliced in response to RSV. These data indicate that substantial post-transcriptional complexity regulates the antiviral response.

CtIP -mediated alternative mRNA splicing finetunes the DNA damage response

In order to survive to the exposure of DNA damaging agents, cells activate a complex response that coordinates the cellular metabolism, cell cycle progression and DNA repair. Among many other events, recent evidence has described global changes in mRNA splicing in cells treated with genotoxic agents. Here, we explore further this DNA damage-dependent alternative splicing. Indeed, we show that both the splicing factor SF3B2 and the repair protein CtIP contribute to the global pattern of splicing both in cells treated or not to DNA damaging agents. Additionally, we focus on a specific DNA damage- and CtIP-dependent alternative splicing event of the helicase PIF1 and explore its relevance for the survival of cells upon exposure to ionizing radiation. Indeed, we described how the nuclear, active form of PIF1 is substituted by a splicing variant, named vPIF1, in a fashion that requires both the presence of DNA damage and CtIP. Interestingly, timely expression of vPIF1 is required for optimal survival to exposure to DNA damaging agents, but early expression of this isoform delays early events of the DNA damage response. On the contrary, expression of the full length PIF1 facilitates those early events, but increases the sensitivity to DNA damaging agents if the expression is maintained long-term.

Identification of Genomewide Alternative Splicing Events in Sequential, Isogenic Clinical Isolates of Candida albicans Reveals a Novel Mechanism of Drug Resistance and Tolerance to Cellular Stresses

The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

Alternative splicing (AS)—a process by which a single gene gives rise to different protein isoforms in eukaryotes—has been implicated in many basic cellular processes, but little is known about its role in drug resistance and fungal pathogenesis. The most common human fungal pathogen, Candida albicans, has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Here, we report AS regulating drug resistance in C. albicans. Comparative RNA-sequencing of two different sets of sequential, isogenic azole-sensitive and -resistant isolates of C. albicans revealed differential expression of splice isoforms of 14 genes. One of these was the superoxide dismutase gene SOD3, which contains a single intron. The sod3Δ/Δ mutant was susceptible to the antifungals amphotericin B (AMB) and menadione (MND). While AMB susceptibility was rescued by overexpression of both the spliced and unspliced SOD3 isoforms, only the spliced isoform could overcome MND susceptibility, demonstrating the functional relevance of this splicing in developing drug resistance. Furthermore, unlike AMB, MND inhibits SOD3 splicing and acts as a splicing inhibitor. Consistent with these observations, MND exposure resulted in increased levels of unspliced SOD3 isoform that are unable to scavenge reactive oxygen species (ROS), resulting in increased drug susceptibility. Collectively, these observations suggest that AS is a novel mechanism for stress adaptation and overcoming drug susceptibility in C. albicans.

IMPORTANCE The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

tappAS: a comprehensive computational framework for the analysis of the functional impact of differential splicing

Recent advances in long-read sequencing solve inaccuracies in alternative transcript identification of full-length transcripts in short-read RNA-Seq data, which encourages the development of methods for isoform-centered functional analysis. Here, we present tappAS, the first framework to enable a comprehensive Functional Iso-Transcriptomics (FIT) analysis, which is effective at revealing the functional impact of context-specific post-transcriptional regulation. tappAS uses isoform-resolved annotation of coding and non-coding functional domains, motifs, and sites, in combination with novel analysis methods to interrogate different aspects of the functional readout of transcript variants and isoform regulation. tappAS software and documentation are available at https://app.tappas.org.

Full-length transcriptome reconstruction reveals a large diversity of RNA and protein isoforms in rat hippocampus

Gene annotation is a critical resource in genomics research. Many computational approaches have been developed to assemble transcriptomes based on high-throughput short-read sequencing, however, only with limited accuracy. Here, we combine next-generation and third-generation sequencing to reconstruct a full-length transcriptome in the rat hippocampus, which is further validated using independent 5´ and 3´-end profiling approaches. In total, we detect 28,268 full-length transcripts (FLTs), covering 6,380 RefSeq genes and 849 unannotated loci. Based on these FLTs, we discover co-occurring alternative RNA processing events. Integrating with polysome profiling and ribosome footprinting data, we predict isoform-specific translational status and reconstruct an open reading frame (ORF)-eome. Notably, a high proportion of the predicted ORFs are validated by mass spectrometry-based proteomics. Moreover, we identify isoforms with subcellular localization pattern in neurons. Collectively, our data advance our knowledge of RNA and protein isoform diversity in the rat brain and provide a rich resource for functional studies.

Nutritional Stress in Head and Neck Cancer Originating Cell Lines: The Sensitivity of the NRF2-NQO1 Axis

Nutritional stress disturbs the cellular redox-status, which is characterized by the increased generation of reactive oxygen species (ROS). The NRF2-NQO1 axis represents a protective mechanism against ROS. Its strength is cell type-specific. FaDu, Cal 27 and Detroit 562 cells differ with respect to basal NQO1 activity. These cells were grown for 48 hours in nutritional conditions (NC): (a) Low glucose-NC2, (b) no glucose, no glutamine-NC3, (c) no glucose with glutamine-NC4. After determining the viability, proliferation and ROS generation, NC2 and NC3 were chosen for further exploration. These conditions were also applied to IMR-90 fibroblasts. The transcripts/transcript variants of NRF2 and NQO1 were quantified and transcript variants were characterized. The proteins (NRF2, NQO1 and TP53) were analyzed by a western blot in both cellular fractions. Under NC2, the NRF2-NQO1 axis did not appear activated in the cancer cell lines. Under NC3, the NRF2-NQO1axis appeared slightly activated in Detroit 562. There are opposite trends with respect to TP53 nuclear signal when comparing Cal 27 and Detroit 562 to FaDu, under NC2 and NC3. The strong activation of the NRF2-NQO1 axis in IMR-90 resulted in an increased expression of catalytically deficient NQO1, due to NQO1*2/*2 polymorphism (rs1800566). The presented results call for a comprehensive exploration of the stress response in complex biological systems.

RNA processing errors triggered by cadmium and integrator complex disruption are signals for environmental stress

BACKGROUND

Adaptive responses to stress are essential for cell and organismal survival. In metazoans, little is known about the impact of environmental stress on RNA homeostasis.

RESULTS

By studying the regulation of a cadmium-induced gene named numr-1 in Caenorhabditis elegans, we discovered that disruption of RNA processing acts as a signal for environmental stress. We find that NUMR-1 contains motifs common to RNA splicing factors and influences RNA splicing in vivo. A genome-wide screen reveals that numr-1 is strongly and specifically induced by silencing of genes that function in basal RNA metabolism including subunits of the metazoan integrator complex. Human integrator processes snRNAs for functioning with splicing factors, and we find that silencing of C. elegans integrator subunits disrupts snRNA processing, causes aberrant pre-mRNA splicing, and induces the heat shock response. Cadmium, which also strongly induces numr-1, has similar effects on RNA and the heat shock response. Lastly, we find that heat shock factor-1 is required for full numr-1 induction by cadmium.

CONCLUSION

Our results are consistent with a model in which disruption of integrator processing of RNA acts as a molecular damage signal initiating an adaptive stress response mediated by heat shock factor-1. When numr-1 is induced via this pathway in C. elegans, its function in RNA metabolism may allow it to mitigate further damage and thereby promote tolerance to cadmium.

Translation Control of HAC1 by Regulation of Splicing in Saccharomyces cerevisiae

Hac1p is a key transcription factor regulating the unfolded protein response (UPR) induced by abnormal accumulation of unfolded/misfolded proteins in the endoplasmic reticulum (ER) in Saccharomyces cerevisiae. The accumulation of unfolded/misfolded proteins is sensed by protein Ire1p, which then undergoes trans-autophosphorylation and oligomerization into discrete foci on the ER membrane. HAC1 pre-mRNA, which is exported to the cytoplasm but is blocked from translation by its intron sequence looping back to its 5’UTR to form base-pair interaction, is transported to the Ire1p foci to be spliced, guided by a cis-acting bipartite element at its 3’UTR (3’BE). Spliced HAC1 mRNA can be efficiently translated. The resulting Hac1p enters the nucleus and activates, together with coactivators, a large number of genes encoding proteins such as protein chaperones to restore and maintain ER homeostasis and secretary protein quality control. This review details the translation regulation of Hac1p production, mediated by the nonconventional splicing, in the broad context of translation control and summarizes the evolution and diversification of the UPR signaling pathway among fungal, metazoan and plant lineages.

Pharmacological and physiological activation of AMPK improves the spliceopathy in DM1 mouse muscles

Myotonic dystrophy type 1 (DM1) is a debilitating multisystemic disorder caused by a triplet repeat expansion in the 3' untranslated region of dystrophia myotonica protein kinase mRNAs. Mutant mRNAs accumulate in the nucleus of affected cells and misregulate RNA-binding proteins, thereby promoting characteristic missplicing events. However, little is known about the signaling pathways that may be affected in DM1. Here, we investigated the status of activated protein kinase (AMPK) signaling in DM1 skeletal muscle and found that the AMPK pathway is markedly repressed in a DM1 mouse model (human skeletal actin-long repeat, HSALR) and patient-derived DM1 myoblasts. Chronic pharmacological activation of AMPK signaling in DM1 mice with 5-aminoimidazole-4-carboxamide-1-β-D-ribofuranoside (AICAR) has multiple beneficial effects on the DM1 phenotype. Indeed, a 6-week AICAR treatment of DM1 mice promoted expression of a slower, more oxidative phenotype, improved muscle histology and corrected several events associated with RNA toxicity. Importantly, AICAR also had a dose-dependent positive effect on the spliceopathy in patient-derived DM1 myoblasts. In separate experiments, we also show that chronic treatment of DM1 mice with resveratrol as well as voluntary wheel running also rescued missplicing events in muscle. Collectively, our findings demonstrate the therapeutic potential of chronic AMPK stimulation both physiologically and pharmacologically for DM1 patients.

RNA-Seq approach for accurate characterization of splicing efficiency of yeast introns

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N_EE) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N_EI5 and N_EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N_Total = N_EE + p * N_EI5 + (1-p) * N_EI3, with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N_EI5 is expected to be smaller than N_EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N_EE/N_Total). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated.

Differential expression of alternative transcripts of soluble guanylyl cyclase, GYCY1a3 and GUCY1b3 genes, in the malignant and benign breast tumors

Extensive alterations in splicing is one of the molecular indicator for human cancers. Soluble guanylyl cyclase (sGC), an obligatory heterodimer, is composed of α1 and β1 subunits. Each subunit is encoded by a separate gene, GUCY1a3 and GUCY1b3, correspondingly. sGC activity has been regulated by an alternative splicing and it has an important effect on the breast cancer. sGC alternative splicing has been evaluated in the 55 malignant, 25 benign and 30 normal breast tissues using qRT-PCR and RT-PCR. The differences between groups were analyzed by Mann-Whitney U. The expression of six different splice forms have been detected, three for α1 and three for β1 sGC. Expressions of Tr1, Tr2 β1 sGC and Tr7, Tr6 α1 sGC mRNA in the malignant breast tumors were significantly lower than those of benign and normal breast tissues. However, the expression of Tr3 α1 sGC mRNA was significantly higher than that of benign and normal tissues. Present data have provided some evidences for an alteration in the expression of α1 and β1 sGC alternative splicing forms which may contribute to the loss of sGC functions in the breast cancer. The observed information might be discussed by the cGMP status.

Effect of ethanol on lipid metabolism

Hepatic lipid metabolism is a series of complex processes that control influx and efflux of not only hepatic lipid pools, but also organismal pools. Lipid homeostasis is usually tightly controlled by expression, substrate supply, oxidation and secretion that keep hepatic lipid pools relatively constant. However, perturbations of any of these processes can lead to lipid accumulation in the liver. Although it is thought that these responses are hepatic arms of the 'thrifty genome', they are maladaptive in the context of chronic fatty liver diseases. Ethanol is likely unique among toxins, in that it perturbs almost all aspects of hepatic lipid metabolism. This complex response is due in part to the large metabolic demand placed on the organ by alcohol metabolism, but also appears to involve more nuanced changes in expression and substrate supply. The net effect is that steatosis is a rapid response to alcohol abuse. Although transient steatosis is largely an inert pathology, the chronicity of alcohol-related liver disease seems to require steatosis. Better and more specific understanding of the mechanisms by which alcohol causes steatosis may therefore translate into targeted therapies to treat alcohol-related liver disease and/or prevent its progression.

Excised linear introns regulate growth in yeast

Spliceosomal introns are ubiquitous non-coding RNAs that are typically destined for rapid debranching and degradation. Here we describe 34 excised introns in Saccharomyces cerevisiae that-despite being rapidly degraded in log-phase growth-accumulate as linear RNAs under either saturated-growth conditions or other stresses that cause prolonged inhibition of TORC1, which is a key integrator of growth signalling. Introns that become stabilized remain associated with components of the spliceosome and differ from other spliceosomal introns in having a short distance between their lariat branch point and 3' splice site, which is necessary and sufficient for their stabilization. Deletion of these unusual introns is disadvantageous in saturated conditions and causes aberrantly high growth rates in yeast that are chronically challenged with the TORC1 inhibitor rapamycin. The reintroduction of native or engineered stable introns suppresses this aberrant rapamycin response. Thus, excised introns function within the TOR growth-signalling network of S. cerevisiae and, more generally, excised spliceosomal introns can have biological functions.

Non-coding transcript variants of protein-coding genes - what are they good for?

The total number of protein-coding genes in the human genome is not significantly higher than those in much simpler eukaryotes, despite a general increase in genome size proportionate to the organismal complexity. The large non-coding transcriptome and extensive differential splicing, are increasingly being accepted as the factors contributing to the complex mammalian physiology and architecture. Recent studies reveal additional layers of functional complexity: some long non-coding RNAs have been re-defined as micropeptide or microprotein encoding transcripts, and in turn some protein-coding RNAs are bifunctional and display also non-coding functions. Moreover, several protein-coding genes express long non-coding RNA splice-forms and generate circular RNAs in addition to their canonical mRNA transcripts, revoking the strict definition of a gene as coding or non-coding. In this mini review, we discuss the current understanding of these hybrid genes and their possible roles and relevance.

Alternatively spliced variants in Atlantic cod (Gadus morhua) support response to variable salinity environment

Analysis of gill transcriptome of the Atlantic cod from the Baltic Sea demonstrated that alternatively spliced (AS) variants may be actively involved in the process of adaptation to altered salinity. Some AS variants of different genes, like phospholipase A2 group IVC (PLA2G4C), appeared only in fish exposed to altered salinity, while other isoforms of the same genes were present in all experimental groups. Novel sequence arrangements represent 89% of all AS in the Baltic cod compared to the Atlantic population. Profiles of modified pathways suggest that regulation by AS can afford specific changes of genes expressed in response to the environment. The AS variants appear to be involved in the response to stress by modifications of signalling in apoptosis pathways, an innate immunological response and pro-inflammatory process. Present results support the hypothesis that developing new AS variants could support genome complexity and reinforce the ability to fast adapt to local environments.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Dual RNA Processing Roles of Pat1b via Cytoplasmic Lsm1-7 and Nuclear Lsm2-8 Complexes

Pat1 RNA-binding proteins, enriched in processing bodies (P bodies), are key players in cytoplasmic 5' to 3' mRNA decay, activating decapping of mRNA in complex with the Lsm1-7 heptamer. Using co-immunoprecipitation and immunofluorescence approaches coupled with RNAi, we provide evidence for a nuclear complex of Pat1b with the Lsm2-8 heptamer, which binds to the spliceosomal U6 small nuclear RNA (snRNA). Furthermore, we establish the set of interactions connecting Pat1b/Lsm2-8/U6 snRNA/SART3 and additional U4/U6.U5 tri-small nuclear ribonucleoprotein particle (tri-snRNP) components in Cajal bodies, the site of snRNP biogenesis. RNA sequencing following Pat1b depletion revealed the preferential upregulation of mRNAs normally found in P bodies and enriched in 3' UTR AU-rich elements. Changes in >180 alternative splicing events were also observed, characterized by skipping of regulated exons with weak donor sites. Our data demonstrate the dual role of a decapping enhancer in pre-mRNA processing as well as in mRNA decay via distinct nuclear and cytoplasmic Lsm complexes.

Thyroid transcriptome analysis reveals different adaptive responses to cold environmental conditions between two chicken breeds

Selection for cold tolerance in chickens is important for improving production performance and animal welfare. The identification of chicken breeds with higher cold tolerance and production performance will help to target candidates for the selection. The thyroid gland plays important roles in thermal adaptation, and its function is influenced by breed differences and transcriptional plasticity, both of which remain largely unknown in the chicken thyroid transcriptome. In this study, we subjected Bashang Long-tail (BS) and Rhode Island Red (RIR) chickens to either cold or warm environments for 21 weeks and investigated egg production performance, body weight changes, serum thyroid hormone concentrations, and thyroid gland transcriptome profiles. RIR chickens had higher egg production than BS chickens under warm conditions, but BS chickens produced more eggs than RIRs under cold conditions. Furthermore, BS chickens showed stable body weight gain under cold conditions while RIRs did not. These results suggested that BS breed is a preferable candidate for cold-tolerance selection and that the cold adaptability of RIRs should be improved in the future. BS chickens had higher serum thyroid hormone concentrations than RIRs under both environments. RNA-Seq generated 344.3 million paired-end reads from 16 sequencing libraries, and about 90% of the processed reads were concordantly mapped to the chicken reference genome. Differential expression analysis identified 46-1,211 genes in the respective comparisons. With regard to breed differences in the thyroid transcriptome, BS chickens showed higher cell replication and development, and immune response-related activity, while RIR chickens showed higher carbohydrate and protein metabolism activity. The cold environment reduced breed differences in the thyroid transcriptome compared with the warm environment. Transcriptional plasticity analysis revealed different adaptive responses in BS and RIR chickens to cope with the cold, and showed higher responsiveness in BS compared with RIR chickens, suggesting greater adaptability of the thyroid in BS chickens. Moreover, 10,053 differential splicing events were revealed among the groups, with RNA splicing and processing, gene expression, transport, and metabolism being the main affected biological processes, identifying a valuable alternative splicing repertoire for the chicken thyroid. A short isoform of TPO (encoding thyroid peroxidase) containing multiple open reading frames was generated in both breeds by skipping exons 4 and 5 in the cold environment. These findings provide novel clues for future studies of the molecular mechanisms underlying cold adaptation and/or acclimation in chickens.

Introns provide a platform for intergenic regulatory feedback of RPL22 paralogs in yeast

Ribosomal protein genes (RPGs) in Saccharomyces cerevisiae are a remarkable regulatory group that may serve as a model for understanding genetic redundancy in evolutionary adaptations. Most RPGs exist as pairs of highly conserved functional paralogs with divergent untranslated regions and introns. We examined the roles of introns in strains with various combinations of intron and gene deletions in RPL22, RPL2, RPL16, RPL37, RPL17, RPS0, and RPS18 paralog pairs. We found that introns inhibited the expression of their genes in the RPL22 pair, with the RPL22B intron conferring a much stronger effect. While the WT RPL22A/RPL22B mRNA ratio was 93/7, the rpl22aΔi/RPL22B and RPL22A/rpl22bΔi ratios were >99/<1 and 60/40, respectively. The intron in RPL2A stimulated the expression of its own gene, but the removal of the other introns had little effect on expression of the corresponding gene pair. Rpl22 protein abundances corresponded to changes in mRNAs.

Using splicing reporters containing endogenous intron sequences, we demonstrated that these effects were due to the inhibition of splicing by Rpl22 proteins but not by their RNA-binding mutant versions. Indeed, only WT Rpl22A/Rpl22B proteins (but not the mutants) interacted in a yeast three-hybrid system with an RPL22B intronic region between bp 165 and 236. Transcriptome analysis showed that both the total level of Rpl22 and the A/B ratio were important for maintaining the WT phenotype. The data presented here support the contention that the Rpl22B protein has a paralog-specific role.

The RPL22 singleton of Kluyveromyces lactis, which did not undergo whole genome duplication, also responded to Rpl22-mediated inhibition in K. lactis cells. Vice versa, the overproduction of the K. lactis protein reduced the expression of RPL22A/B in S. cerevisiae. The extraribosomal function of of the K. lactis Rpl22 suggests that the loop regulating RPL22 paralogs of S. cerevisiae evolved from autoregulation.

Intron specificity in pre-mRNA splicing

The occurrence of spliceosomal introns in eukaryotic genomes is highly diverse and ranges from few introns in an organism to multiple introns per gene. Introns vary with respect to their lengths, strengths of splicing signals, and position in resident genes. Higher intronic density and diversity in genetically complex organisms relies on increased efficiency and accuracy of spliceosomes for pre-mRNA splicing. Since intron diversity is critical for functions in RNA stability, regulation of gene expression and alternative splicing, RNA-binding proteins, spliceosomal regulatory factors and post-translational modifications of splicing factors ought to make the splicing process intron-specific. We recently reported function and regulation of a ubiquitin fold harboring splicing regulator, Sde2, which following activation by ubiquitin-specific proteases facilitates excision of selected introns from a subset of multi-intronic genes in Schizosaccharomyces pombe (Thakran et al. EMBO J, https://doi.org/10.15252/embj.201796751 , 2017). By reviewing our findings with understandings of intron functions and regulated splicing processes, we propose possible functions and mechanism of intron-specific pre-mRNA splicing and suggest that this process is crucial to highlight importance of introns in eukaryotic genomes.

Transcription rate strongly affects splicing fidelity and cotranscriptionality in budding yeast

The functional consequences of alternative splicing on altering the transcription rate have been the subject of intensive study in mammalian cells but less is known about effects of splicing on changing the transcription rate in yeast. We present several lines of evidence showing that slow RNA polymerase II elongation increases both cotranscriptional splicing and splicing efficiency and that faster elongation reduces cotranscriptional splicing and splicing efficiency in budding yeast, suggesting that splicing is more efficient when cotranscriptional. Moreover, we demonstrate that altering the RNA polymerase II elongation rate in either direction compromises splicing fidelity, and we reveal that splicing fidelity depends largely on intron length together with secondary structure and splice site score. These effects are notably stronger for the highly expressed ribosomal protein coding transcripts. We propose that transcription by RNA polymerase II is tuned to optimize the efficiency and accuracy of ribosomal protein gene expression, while allowing flexibility in splice site choice with the nonribosomal protein transcripts.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5'ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5'ss in dtd1+ intron 1 and of an upstream alternative 3'ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5'ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3'ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

ARSDA: A New Approach for Storing, Transmitting and Analyzing Transcriptomic Data

Two major stumbling blocks exist in high-throughput sequencing (HTS) data analysis. The first is the sheer file size, typically in gigabytes when uncompressed, causing problems in storage, transmission, and analysis. However, these files do not need to be so large, and can be reduced without loss of information. Each HTS file, either in compressed .SRA or plain text .fastq format, contains numerous identical reads stored as separate entries. For example, among 44,603,541 forward reads in the SRR4011234.sra file (from a Bacillus subtilis transcriptomic study) deposited at NCBI’s SRA database, one read has 497,027 identical copies. Instead of storing them as separate entries, one can and should store them as a single entry with the SeqID_NumCopy format (which I dub as FASTA+ format). The second is the proper allocation of reads that map equally well to paralogous genes. I illustrate in detail a new method for such allocation. I have developed ARSDA software that implement these new approaches. A number of HTS files for model species are in the process of being processed and deposited at http://coevol.rdc.uottawa.ca to demonstrate that this approach not only saves a huge amount of storage space and transmission bandwidth, but also dramatically reduces time in downstream data analysis. Instead of matching the 497,027 identical reads separately against the B. subtilis genome, one only needs to match it once. ARSDA includes functions to take advantage of HTS data in the new sequence format for downstream data analysis such as gene expression characterization. I contrasted gene expression results between ARSDA and Cufflinks so readers can better appreciate the strength of ARSDA. ARSDA is freely available for Windows, Linux. and Macintosh computers at http://dambe.bio.uottawa.ca/ARSDA/ARSDA.aspx.

The conserved AU dinucleotide at the 5' end of nascent U1 snRNA is optimized for the interaction with nuclear cap-binding-complex

Splicing is initiated by a productive interaction between the pre-mRNA and the U1 snRNP, in which a short RNA duplex is established between the 5' splice site of a pre-mRNA and the 5' end of the U1 snRNA. A long-standing puzzle has been why the AU dincucleotide at the 5'-end of the U1 snRNA is highly conserved, despite the absence of an apparent role in the formation of the duplex. To explore this conundrum, we varied this AU dinucleotide into all possible permutations and analyzed the resulting molecular consequences. This led to the unexpected findings that the AU dinucleotide dictates the optimal binding of cap-binding complex (CBC) to the 5' end of the nascent U1 snRNA, which ultimately influences the utilization of U1 snRNP in splicing. Our data also provide a structural interpretation as to why the AU dinucleotide is conserved during evolution.

The chromatin remodeling complex Swi/Snf regulates splicing of meiotic transcripts in Saccharomyces cerevisiae

Despite its relatively streamlined genome, there are important examples of regulated RNA splicing in Saccharomyces cerevisiae, such as splicing of meiotic transcripts. Like other eukaryotes, S. cerevisiae undergoes a dramatic reprogramming of gene expression during meiosis, including regulated splicing of a number of crucial meiosis-specific RNAs. Splicing of a subset of these is dependent upon the splicing activator Mer1. Here we show a crucial role for the chromatin remodeler Swi/Snf in regulation of splicing of meiotic genes and find that the complex affects meiotic splicing in two ways. First, we show that Swi/Snf regulates nutrient-dependent downregulation of ribosomal protein encoding RNAs, leading to the redistribution of spliceosomes from this abundant class of intron-containing RNAs (the ribosomal protein genes) to Mer1-regulated transcripts. We also demonstrate that Mer1 expression is dependent on Snf2, its acetylation state and histone H3 lysine 9 acetylation at the MER1 locus. Hence, Snf2 exerts systems level control of meiotic gene expression through two temporally distinct mechanisms, demonstrating that it is a key regulator of meiotic splicing in S. cerevisiae. We also reveal an evolutionarily conserved mechanism whereby the cell redirects its energy from maintaining its translational capacity to the process of meiosis.