Exon size and sequence conservation improves identification of splice-altering nucleotides

Splicing accuracy varies across human introns, tissues, age and disease

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigate splicing accuracy using RNA-sequencing data from >14k control samples and 40 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that splicing inaccuracies occur at different rates across introns and tissues and are affected by the abundance of core components of the spliceosome assembly and its regulators. We find that age is positively correlated with a global decline in splicing fidelity, mostly affecting genes implicated in neurodegenerative diseases. We find support for the latter by observing a genome-wide increase in splicing inaccuracies in samples affected with Alzheimer's disease as compared to neurologically normal individuals. In this work, we provide an in-depth characterisation of splicing accuracy, with implications for our understanding of the role of inaccuracies in ageing and neurodegenerative disorders.

Evolutionarily Developed Alternatively Spliced Exons Containing Translation Initiation Sites

Alternative splicing is essential for the generation of various protein isoforms that are involved in cell differentiation and tissue development. In addition to internal coding exons, alternative splicing affects the exons with translation initiation codons; however, little is known about these exons. Here, we performed a systematic classification of human alternative exons using coding information. The analysis showed that more than 5% of cassette exons contain translation initiation codons (alternatively skipped exons harboring a 5' untranslated region and coding region, 5UC-ASEs) although their skipping causes the deletion of translation initiation sites essential for protein synthesis. The splicing of 5UC-ASEs is under the repressive control of MATR3, a DNA/RNA-binding protein associated with neurodegeneration, and is distinctly regulated particularly in the human brain, muscle, and testis. Interestingly, MATR3 represses its own translation by skipping a 5UC-ASE in MATR3 to autoregulate its expression level. 5UC-ASEs are larger than other types of alternative exons. Furthermore, evolutionary analysis revealed that 5UC-ASEs have already appeared in cartilaginous fishes, have increased in amphibians, and are concentrated in the genes involved in transcription in mammals. Taken together, our analysis identified a unique set of alternative exons, 5UC-ASEs, that have evolutionarily acquired a repression mechanism for gene expression through association with MATR3.

Challenges in identifying mRNA transcript starts and ends from long-read sequencing data

Long-read sequencing (LRS) technologies have the potential to revolutionize scientific discoveries in RNA biology through the comprehensive identification and quantification of full-length mRNA isoforms. Despite great promise, challenges remain in the widespread implementation of LRS technologies for RNA-based applications, including concerns about low coverage, high sequencing error, and robust computational pipelines. Although much focus has been placed on defining mRNA exon composition and structure with LRS data, less careful characterization has been done of the ability to assess the terminal ends of isoforms, specifically, transcription start and end sites. Such characterization is crucial for completely delineating full mRNA molecules and regulatory consequences. However, there are substantial inconsistencies in both start and end coordinates of LRS reads spanning a gene, such that LRS reads often fail to accurately recapitulate annotated or empirically derived terminal ends of mRNA molecules. Here, we describe the specific challenges of identifying and quantifying mRNA terminal ends with LRS technologies and how these issues influence biological interpretations of LRS data. We then review recent experimental and computational advances designed to alleviate these problems, with ideal use cases for each approach. Finally, we outline anticipated developments and necessary improvements for the characterization of terminal ends from LRS data.

Association of SPP1 and NCAPG genes with milk production traits in Chinese Holstein cows: polymorphism and functional validation analysis

Milk production traits play an important role in dairy cattle breeding, and single nucleotide polymorphisms can be used as effective molecular markers for milk production trait marker-assisted breeding in dairy cattle. Based on the results of the preliminary GWAS, candidate genes SPP1 and NCAPG associated with milk production traits were screened. In this study, the aim was to screen and characterize the SNPs of SPP1 and NCAPG genes about milk production traits. Two SNPs and one haplotype block of the SPP1 gene and four SNPs and one haplotype block of the NCAPG gene were obtained by amplification, sequencing and association analysis, and all six SNPs were located in the exon region. Association analysis showed that all six SNPs were significantly associated with milk protein percentage. Linkage disequilibrium analysis showed that 2 SNPs of SPP1 (g. 36,700,265 C > T and g. 36,693,596 C > A) constituted a haplotype that correlated with milk protein percentage, and the dominant haplotype was H2H2, which was CCTT. 4 SNPs of NCAPG (g. 37,342,705 C > A, g. 37,343,379 G > T, g. 37,374,314 C > A and g. 37,377,857 G > A) constituted a haplotype associated with milk protein percentage, 305-days milk protein yield and 305 days milk yield. Tissue expression profiling results revealed that SPP1 and NCAPG had the highest expression in mammary tissue. Interference with SPP1 and NCAPG inhibited the proliferation of Bovine mammary epithelial cells. (BMECs), down-regulated the expression of PCNA, CDK2 and CCND1, up-regulated the expression of BAX and BAD, and promoted apoptosis. Reduced triglyceride synthesis in BMECs, down-regulated the expression of DGAT1, DGAT2, LPIN1, and AGPAT6.SPP1 and NCAPG are involved in the synthesis of milk proteins, and interfering with SPP1 and NCAPG decreased the secretion of β-casein, κ-casein, and αs1-casein, as well as up-regulated the CSN2 and CSN3 expression. The above results indicate that the SNP loci of SPP1 and NCAPG can be used as potential molecular markers to improve milk production traits in dairy cows, laying the foundation for marker-assisted selection. It also proves that SPP1 and NCAPG can be used as candidate key genes for milk production traits in dairy cows, providing new insights into the physiological mechanisms of lactation regulation in dairy cows.

Branch site recognition by the spliceosome

The spliceosome is a eukaryotic multimegadalton RNA-protein complex that removes introns from transcripts. The spliceosome ensures the selection of each exon-intron boundary through multiple recognition events. Initially, the 5' splice site (5' SS) and branch site (BS) are bound by the U1 small nuclear ribonucleoprotein (snRNP) and the U2 snRNP, respectively, while the 3' SS is mostly determined by proximity to the branch site. A large number of splicing factors recognize the splice sites and recruit the snRNPs before the stable binding of the snRNPs occurs by base-pairing the snRNA to the transcript. Fidelity of this process is crucial, as mutations in splicing factors and U2 snRNP components are associated with many diseases. In recent years, major advances have been made in understanding how splice sites are selected in Saccharomyces cerevisiae and humans. Here, I review and discuss the current understanding of the recognition of splice sites by the spliceosome with a focus on recognition and binding of the branch site by the U2 snRNP in humans.

Association of Gene Variant Type and Location with Breast Cancer Risk in the General Population

Importance

Pathogenic variants (PVs) in ATM, BRCA1, BRCA2, CHEK2 , and PALB2 are associated with increased breast cancer risk. However, it is unknown whether breast cancer risk differs by PV type or location in carriers ascertained from the general population.

Objective

To evaluate breast cancer risks associated with PV type and location in ATM, BRCA1, BRCA2, CHEK2 , and PALB2 .

Design

Age adjusted case-control association analysis for all participants, subsets of PV carriers, and women with no breast cancer family history in population-based and clinical testing cohorts.

Setting

Twelve US population-based studies within the Cancer Risk Estimates Related to Susceptibility (CARRIERS) Consortium, and breast cancer cases from the UK-Biobank and an Ambry Genetics clinical testing cohort.

Participants

32,247 women with and 32,544 age-matched women without a breast cancer diagnosis from CARRIERS; 237 and 1351 women with BRCA2 PVs and breast cancer from the UKBB and Ambry Genetics, respectively.

Exposures

PVs in ATM, BRCA1, BRCA2, CHEK2, and PALB2.

Main Outcomes and Measures

PVs were grouped by type and location within genes and assessed for risks of breast cancer (odds ratios (OR), 95% confidence intervals (CI), and p-values) using logistic regression. Mean ages at diagnosis were compared using linear regression.

Results

Compared to women carrying BRCA2 exon 11 protein truncating variants (PTVs) in the CARRIERS population-based study, women with BRCA2 ex13-27 PTVs (OR=2.7, 95%CI 1.1-7.9) and ex1-10 PTVs (OR=1.6, 95%CI 0.8-3.5) had higher breast cancer risks, lower rates of ER-negative breast cancer (ex13-27 OR=0.5, 95%CI 0.2-0.9; ex1-10 OR=0.5, 95%CI 0.1-1.0), and earlier age of breast cancer diagnosis (ex13-27 5.5 years, p<0.001; ex1-10 2.4 years, p=0.17). These associations with ER-negative breast cancer and age replicated in a high-risk clinical cohort and the population-based UK Biobank cohort. No differences in risk or age at diagnosis by gene region were observed for PTVs in other predisposition genes.

Conclusions and Relevance

Population-based and clinical high-risk cohorts establish that PTVs in exon 11 of BRCA2 are associated with reduced risk of breast cancer, later age at diagnosis, and greater risk of ER-negative disease. These differential risks may improve individualized risk prediction and clinical management for women carrying BRCA2 PTVs.

Key Points

Question: Does ATM , BRCA1 , BRCA2 , CHEK2 and PALB2 pathogenic variant type and location influence breast cancer risk in population-based studies? Findings: Breast cancer risk and estrogen receptor status differ based on the type and location of pathogenic variants in BRCA2 . Women carrying protein truncating variants in exon 11 have a lower breast cancer risk in the population-based cohorts, older age at diagnosis and higher rates of estrogen receptor negative breast cancer than women with exon 1-10 or exon 13-27 truncation variants in population-based and clinical testing cohorts. Meaning: Incorporating pathogenic variant type and location in cancer risk models may improve individualized risk prediction.

Splice site recognition - deciphering Exon-Intron transitions for genetic insights using Enhanced integrated Block-Level gated LSTM model

Bioinformatics is a contemporary interdisciplinary area focused on analyzing the growing number of genome sequences. Gene variants are differences in DNA sequences among individuals within a population. Splice site recognition is a crucial step in the process of gene expression, where the coding sequences of genes are joined together to form mature messenger RNA (mRNA). These genetic variants that disrupt genes are believed to be the primary reason for neuro-developmental disorders like ASD (Autism Spectrum Disorder) is a neuro-developmental disorder that is diagnosed in individuals, families, and society and occurs as the developmental delay in one among the hundred genes that are associated with these disorders. Missense variants, premature stop codons, or deletions alter both the quality and quantity of encoded proteins. Predicting genes within exons and introns presents main challenges, such as dealing with sequencing errors, short reads, incomplete genes, overlapping, and more. Although many traditional techniques have been utilized in creating an exon prediction system, the primary challenge lies in accurately identifying the length and spliced strand location classification of exons in conjunction with introns. From now on, the suggested approach utilizes a Deep Learning algorithm to analyze intricate and extensive genomic datasets. M-LSTM is utilized to categorize three binary combinations (EI as 1, IE as 2, and none as 3) using spliced DNA strands. The M-LSTM system is able to sequence extensive datasets, ensuring that long information can be stored without any impact on the current input or output. This enables it to recognize and address long-term connections and problems with rapidly increasing gradients. The proposed model is compared internally with Naïve Bayes and Random Forest to assess its efficacy. Additionally, the proposed model's performance is forecasted by utilizing probabilistic parameters like recall, F1-score, precision, and accuracy to assess the effectiveness of the proposed system.

The Human Developing Cerebral Cortex Is Characterized by an Elevated De Novo Expression of Long Noncoding RNAs in Excitatory Neurons

The outstanding human cognitive capacities are computed in the cerebral cortex, a mammalian-specific brain region and the place of massive biological innovation. Long noncoding RNAs have emerged as gene regulatory elements with higher evolutionary turnover than mRNAs. The many long noncoding RNAs identified in neural tissues make them candidates for molecular sources of cerebral cortex evolution and disease. Here, we characterized the genomic and cellular shifts that occurred during the evolution of the long noncoding RNA repertoire expressed in the developing cerebral cortex and explored putative roles for these long noncoding RNAs in the evolution of the human brain. Using transcriptomics and comparative genomics, we comprehensively annotated the cortical transcriptomes of humans, rhesus macaques, mice, and chickens and classified human cortical long noncoding RNAs into evolutionary groups as a function of their predicted minimal ages. Long noncoding RNA evolutionary groups showed differences in expression levels, splicing efficiencies, transposable element contents, genomic distributions, and transcription factor binding to their promoters. Furthermore, older long noncoding RNAs showed preferential expression in germinative zones, outer radial glial cells, and cortical inhibitory (GABAergic) neurons. In comparison, younger long noncoding RNAs showed preferential expression in cortical excitatory (glutamatergic) neurons, were enriched in primate and human-specific gene co-expression modules, and were dysregulated in neurodevelopmental disorders. These results suggest different evolutionary routes for older and younger cortical long noncoding RNAs, highlighting old long noncoding RNAs as a possible source of molecular evolution of conserved developmental programs; conversely, we propose that the de novo expression of primate- and human-specific young long noncoding RNAs is a putative source of molecular evolution and dysfunction of cortical excitatory neurons, warranting further investigation.

Copper-binding proteins and exonic splicing enhancers and silencers

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.

Deciphering the role of alternative splicing in neoplastic diseases for immune-oncological therapies

Alternative splicing (AS) is an important molecular biological mechanism regulated by complex mechanisms involving a plethora of cis and trans-acting elements. Furthermore, AS is tissue specific and altered in various pathologies, including infectious, inflammatory, and neoplastic diseases. Recently developed immuno-oncological therapies include monoclonal antibodies (mAbs) and chimeric antigen receptor (CAR) T cells targeting, among others, immune checkpoint (ICP) molecules. Despite therapeutic successes have been demonstrated, only a limited number of patients showed long-term benefit from these therapies with tumor entity-related differential response rates were observed. Interestingly, splice variants of common immunotherapeutic targets generated by AS are able to completely escape and/or reduce the efficacy of mAb- and/or CAR-based tumor immunotherapies. Therefore, the analyses of splicing patterns of targeted molecules in tumor specimens prior to therapy might help correct stratification, thereby increasing therapy success by antibody panel selection and antibody dosages. In addition, the expression of certain splicing factors has been linked with the patients' outcome, thereby highlighting their putative prognostic potential. Outstanding questions are addressed to translate the findings into clinical application. This review article provides an overview of the role of AS in (tumor) diseases, its molecular mechanisms, clinical relevance, and therapy response.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Evolution and phylogenetic diversity of the aquaporin gene family in arachnids

Water flux across cells predominantly occurs through the pore formed by the aquaporin channels. Since water balance is one of the most important challenges to terrestrial animals, aquaporin evolution and diversity is known to play roles in animal terrestrialisation. Arachnids (Arthropoda: Chelicerata: Arachnida) are the second most diverse group and represent the pioneer land colonists in animals; however, there remains no thorough investigation on aquaporin evolution and diversity in this evolutionarily important lineage. Here we reported a phylogenetic study of aquaporin evolution and diversity using genomic data from 116 arachnid species covering almost all (15/16) extant orders. A previously unrecognised subfamily related to aquaporin-4 (i.e. Aqp4-like subfamily) via phylogenetic analysis was identified, suggesting certain underestimate of the arachnid aquaporin diversity in earlier studies probably due to limited taxonomic sampling. Further analysis indicates that this subfamily emerged deep within the life tree of arthropods. Gene tree of another Aqp4-like subfamily (PripL) shows an unexpected basal split between acariform mites (Acariformes) and other arachnids. A closer inspection demonstrated that the PripL evolved quickly and has been under differential selection pressure in acariform mites. Evidence is provided that the evolutionarily ancient Glp subfamily (i.e. aquaglyceroporin) is significantly expanded in terrestrial arachnids compared with their marine relatives. Finally, in spite of the phylogenetic diversity, there exists conservation of some exons in size, functional domain, and intron-insertion phase: an 81-bp and a 218-bp exon, respectively, in apq4-like and glp genes across Eumetazoa lineages including arachnids and human beings. Both exons encode the carboxyl-terminal NPA motif, implying the coding and splicing pressure during hundreds of million years of animal evolution. Hypotheses were tested to explore the possible link between these findings and arachnid terrestrialisation.

Novornabreak: Local Assembly for Novel Splice Junction and Fusion Transcript Detection from RNA-Seq Data

We present novoRNABreak, a unified framework for cancer specific novel splice junction and fusion transcript detection in RNA-seq data obtained from human cancer samples. novoRNABreak is based on a local assembly model, which offers a tradeoff between the alignment-based and de novo whole transcriptome assembly (WTA) methods. This approach is accurate and sensitive in assembling novel junctions that are difficult to directly align or have multiple alignments. Additionally, it is more efficient due to the strategy that focuses on junctions rather than full length transcripts. The performance of novoRNABreak is demonstrated by a comprehensive set of experiments using synthetic data generated based on genome reference, as well as real RNA-seq data from breast cancer and prostate cancer samples. The results show that our tool has a better performance by fully utilizing unmapped reads and precisely identifying the junctions where short reads or small exons have multiple alignments. novoRNABreak is a fully-fledged program available on GitHub (https://github.com/KChen-lab/novoRNABreak).

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Hidden codes in mRNA: Control of gene expression by m6A

Information in mRNA has largely been thought to be confined to its nucleotide sequence. However, the advent of mapping techniques to detect modified nucleotides has revealed that mRNA contains additional information in the form of chemical modifications. The most abundant modified nucleotide is N⁶-methyladenosine (m⁶A), a methyl modification of adenosine. Although early studies viewed m⁶A as a dynamic and tissue-specific modification, it is now clear that the mRNAs that contain m⁶A and the location of m⁶A in those transcripts are largely universal and are influenced by gene architecture, i.e., the size and location of exons and introns. m⁶A can affect nuclear processes such as splicing and epigenetic regulation, but the major effect of m⁶A on mRNAs is to promote degradation in the cytoplasm. m⁶A marks a functionally related cohort of mRNAs linked to certain biological processes, including cell differentiation and cell fate determination. m⁶A is also enriched in other cohorts of mRNAs and can therefore affect their respective cellular processes and pathways. Future work will focus on understanding how the m⁶A pathway is regulated to achieve control of m⁶A-containing mRNAs.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Splice site proximity influences alternative exon definition

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

The Dynamics of NO3- and NH4+ Uptake in Duckweed Are Coordinated with the Expression of Major Nitrogen Assimilation Genes

Duckweed plants play important roles in aquatic ecosystems worldwide. They rapidly accumulate biomass and have potential uses in bioremediation of water polluted by fertilizer runoff or other chemicals. Here we studied the assimilation of two major sources of inorganic nitrogen, nitrate (NO3- ) and ammonium (NH4+), in six duckweed species: Spirodela polyrhiza, Landoltia punctata, Lemna aequinoctialis, Lemna turionifera, Lemna minor, and Wolffia globosa. All six duckweed species preferred NH4+ over NO3- and started using NO3- only when NH4+ was depleted. Using the available genome sequence, we analyzed the molecular structure and expression of eight key nitrogen assimilation genes in S. polyrhiza. The expression of genes encoding nitrate reductase and nitrite reductase increased about 10-fold when NO3- was supplied and decreased when NH4+ was supplied. NO3- and NH4+ induced the glutamine synthetase (GS) genes GS1;2 and the GS2 by 2- to 5-fold, respectively, but repressed GS1;1 and GS1;3. NH4+ and NO3- upregulated the genes encoding ferredoxin- and NADH-dependent glutamate synthases (Fd-GOGAT and NADH-GOGAT). A survey of nitrogen assimilation gene promoters suggested complex regulation, with major roles for NRE-like and GAATC/GATTC cis-elements, TATA-based enhancers, GA/CTn repeats, and G-quadruplex structures. These results will inform efforts to improve bioremediation and nitrogen use efficiency.

HIV-1: To Splice or Not to Splice, That Is the Question

The transcription of the HIV-1 provirus results in only one type of transcript-full length genomic RNA. To make the mRNA transcripts for the accessory proteins Tat and Rev, the genomic RNA must completely splice. The mRNA transcripts for Vif, Vpr, and Env must undergo splicing but not completely. Genomic RNA (which also functions as mRNA for the Gag and Gag/Pro/Pol precursor polyproteins) must not splice at all. HIV-1 can tolerate a surprising range in the relative abundance of individual transcript types, and a surprising amount of aberrant and even odd splicing; however, it must not over-splice, which results in the loss of full-length genomic RNA and has a dramatic fitness cost. Cells typically do not tolerate unspliced/incompletely spliced transcripts, so HIV-1 must circumvent this cell policing mechanism to allow some splicing while suppressing most. Splicing is controlled by RNA secondary structure, cis-acting regulatory sequences which bind splicing factors, and the viral protein Rev. There is still much work to be done to clarify the combinatorial effects of these splicing regulators. These control mechanisms represent attractive targets to induce over-splicing as an antiviral strategy. Finally, splicing has been implicated in latency, but to date there is little supporting evidence for such a mechanism. In this review we apply what is known of cellular splicing to understand splicing in HIV-1, and present data from our newer and more sensitive deep sequencing assays quantifying the different HIV-1 transcript types.