A Bioinformatics-Based Alternative mRNA Splicing Code that May Explain Some Disease Mutations Is Conserved in Animals

Statistical analysis of synonymous and stop codons in pseudo-random and real sequences as a function of GC content

Knowledge of the frequencies of synonymous triplets in protein-coding and non-coding DNA stretches can be used in gene finding. These frequencies depend on the GC content of the genome or parts of it. An example of interest is provided by stop codons. This is relevant for the definition of Open Reading Frames. A generic case is provided by pseudo-random sequences, especially when they code for complex proteins or when they are non-coding and not subject to selection pressure. Here, we calculate, for such sequences and for all 25 known genetic codes, the frequency of each amino acid and stop codon based on their set of codons and as a function of GC content. The amino acids can be classified into five groups according to the GC content where their expected frequency reaches its maximum. We determine the overall Shannon information based on groups of synonymous codons and show that it becomes maximum at a percent GC of 43.3% (for the standard code). This is in line with the observation that in most fungi, plants, and animals, this genomic parameter is in the range from 35 to 50%. By analysing natural sequences, we show that there is a clear bias for triplets corresponding to stop codons near the 5'- and 3'-splice sites in the introns of various clades.

10 Years of Toxicogenomics section in Frontiers in Genetics: Past discoveries and Future Perspectives

The Frontiers Media family has over 200 journals, which are each headed by usually one Field Chief Editor, and several thousand specialty sections, which are each headed by one or more Specialty Chief Editors. The year 2021 was the 10th anniversary of the founding of the Frontiers in Genetics journal and the Frontiers in Toxicogenomics specialty section of this journal. In 2021, we also announce one of the newest of the Frontiers journals-Frontiers in Toxicology which is part of the Frontiers Media family of journals but independent of Frontiers in Genetics. Dr. Ruden is the founding, and currently sole, Specialty Chief Editor of Frontiers in Toxicogenomics and one of 9 Specialty Chief Editors of Frontiers in Toxicology. As of 2021, Frontiers in Toxicogenomics has published over 138 articles and has over 370 Editors including 90 Associate Editors and 280 Review Editors. The Frontiers in Genetics impact factor was initially approximately 2.5 when it was first listed in PubMed in 2015 and has risen steadily to its current value of 4.8, which is typical for the majority of the over 200 Frontiers journals that have established impact factors. In this overview of the first decade of Frontiers in Toxicogenomics, we discuss the top 5 articles with the highest Scopus citations, which were all written in the first few years of the journal. The article with the highest number of citations, with 353 Scopus over 600 Google Scholar citations, and the highest average number of citations (67) that steadily increased from 10 citations in 2013 to 119 citations in 2021, was written in 2012 by Dr. Ruden's laboratory and titled, "Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift." The five most influential authors who published in the journal in the past 10 years based on Scopus citations of a particular paper are Dr. Ruden's laboratory, with 353 Scopus citations for the SnpSift paper mentioned above; Drs. Brock Christensen and Carmen J. Marsit, with 86 Scopus citations for their review, "Epigenomics in environmental health"; Dr. Michael Aschner and colleagues, with 61 Scopus citations for their paper "Genetic factors and manganese-induced neurotoxicity"; and Dr. Sandra C. dos Santos and colleagues, with 59 Scopus citations for their paper, "Yeast toxicogenomics: genome-wide responses to chemical stresses with impact in environmental health, pharmacology, and biotechnology." While the top 5 papers were published in the early years of the journal, we will also discuss a more recent article published in 2018 on a comparison of RNA-seq and microarray methods by Dr. Michael Liguori's laboratory, "Comparison of RNA-Seq and Microarray Gene Expression Platforms for the Toxicogenomic Evaluation of Liver From Short-Term Rat Toxicity Studies," that far exceeds the number of downloads and views of all the other articles published in the first 10 years of the journal and will likely be a top cited paper in the second decade highlights of this journal. Finally, we discuss where the Frontiers in Toxicogenomics specialty journal and the Frontiers in Toxicology journal will go to advance the field of toxicogenomics, and more generally, toxicology, in the future.

Delivering the Promise of Gene Therapy with Nanomedicines in Treating Central Nervous System Diseases

Central Nervous System (CNS) diseases, such as Alzheimer's diseases (AD), Parkinson's Diseases (PD), brain tumors, Huntington's disease (HD), and stroke, still remain difficult to treat by the conventional molecular drugs. In recent years, various gene therapies have come into the spotlight as versatile therapeutics providing the potential to prevent and treat these diseases. Despite the significant progress that has undoubtedly been achieved in terms of the design and modification of genetic modulators with desired potency and minimized unwanted immune responses, the efficient and safe in vivo delivery of gene therapies still poses major translational challenges. Various non-viral nanomedicines have been recently explored to circumvent this limitation. In this review, an overview of gene therapies for CNS diseases is provided and describes recent advances in the development of nanomedicines, including their unique characteristics, chemical modifications, bioconjugations, and the specific applications that those nanomedicines are harnessed to deliver gene therapies.

CoCoView - A Codon Conservation Viewer via Sequence Logos

Sequence logos are a simple way to display a set of aligned sequences, and they are useful to identify conserved patterns. Since their introduction, several tools have been developed for generating these representations at the single residue level (amino acids or nucleotides). We have developed a tool to build sequence logos of protein-coding sequences at the codon level, allowing more accurate analysis of coding-sequences as they represent synonymous and non-synonymous changes instead of showing only changes that imply on amino acid substitutions. We built CoCoView on top of the Logomaker Python API. It creates codon sequence logos from a multiple sequence alignment of protein-coding sequences. Some properties of the data and the generated logos can be controlled by the end-users, such as data redundancy, plot type and alphabet color.

• Split aligned sequences into codon positions;

• For each position compute codon frequency and information content;

• Use the computed information to plot the graphic.

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.

A new murine Rpl5 (uL18) mutation provides a unique model of variably penetrant Diamond-Blackfan anemia

Ribosome dysfunction is implicated in multiple abnormal developmental and disease states in humans. Heterozygous germline mutations in genes encoding ribosomal proteins are found in most individuals with Diamond-Blackfan anemia (DBA), whereas somatic mutations have been implicated in a variety of cancers and other disorders. Ribosomal protein-deficient animal models show variable phenotypes and penetrance, similar to human patients with DBA. In this study, we characterized a novel ENU mouse mutant (Skax23m1Jus) with growth and skeletal defects, cardiac malformations, and increased mortality. After genetic mapping and whole-exome sequencing, we identified an intronic Rpl5 mutation, which segregated with all affected mice. This mutation was associated with decreased ribosome generation, consistent with Rpl5 haploinsufficiency. Rpl5Skax23-Jus/+ animals had a profound delay in erythroid maturation and increased mortality at embryonic day (E) 12.5, which improved by E14.5. Surviving mutant animals had macrocytic anemia at birth, as well as evidence of ventricular septal defect (VSD). Surviving adult and aged mice exhibited no hematopoietic defect or VSD. We propose that this novel Rpl5Skax23-Jus/+ mutant mouse will be useful in studying the factors influencing the variable penetrance that is observed in DBA.

GHRL gene-based genotyping of ovine and caprine breeds reveals highly polymorphic intronic sequences in Awassi sheep with several RNA motifs

BACKGROUND

The current study was conducted to identify the genetic polymorphism of ghrelin (GHRL) gene of sheep and goats, as well as to determine whether these polymorphisms were associated with the evolutionary genetic differences in the involved species. This study was performed on 233 sheep and 91 goats. Two genetic loci of 113 bp and 262 bp partially spanning over exon 2/intron 2 and intron 4/exon 5 of GHRL gene respectively were amplified and genotyped using polymerase chain reaction-single-strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods.

RESULTS

The SSCP banding pattern of 262-bp locus indicated the presence of four diplotypes (BC, BB, AC, and AB) in Awassi sheep, three diplotypes (BC, BB, and AB) in Karadi sheep, and only two diplotypes (BC and BB) in all goats' samples. The current study detected several novel SNPs in the ovine-caprine populations as well as two SNPs that are observed only in sheep, including intron4:119 C>A and intron4:123 T>G. The phylogenetic analysis revealed that the observed diplotypes resided within ovine sequences and were closely related to caprine counterparts. Computational analyses indicated the presence of various intronic RNA motifs. However, all these motifs were gathered in Awassi breed.

CONCLUSION

It is stated that the intron 4 is highly diverse amongst goats and sheep as well as within sheep with a particular emphasis on Awassi. This genetic peculiarity may in turn suggest a high polymorphic pattern of this breed in comparison with other related counterparts.

Genome-wide analyses supported by RNA-Seq reveal non-canonical splice sites in plant genomes

BACKGROUND

Most eukaryotic genes comprise exons and introns thus requiring the precise removal of introns from pre-mRNAs to enable protein biosynthesis. U2 and U12 spliceosomes catalyze this step by recognizing motifs on the transcript in order to remove the introns. A process which is dependent on precise definition of exon-intron borders by splice sites, which are consequently highly conserved across species. Only very few combinations of terminal dinucleotides are frequently observed at intron ends, dominated by the canonical GT-AG splice sites on the DNA level.

RESULTS

Here we investigate the occurrence of diverse combinations of dinucleotides at predicted splice sites. Analyzing 121 plant genome sequences based on their annotation revealed strong splice site conservation across species, annotation errors, and true biological divergence from canonical splice sites. The frequency of non-canonical splice sites clearly correlates with their divergence from canonical ones indicating either an accumulation of probably neutral mutations, or evolution towards canonical splice sites. Strong conservation across multiple species and non-random accumulation of substitutions in splice sites indicate a functional relevance of non-canonical splice sites. The average composition of splice sites across all investigated species is 98.7% for GT-AG, 1.2% for GC-AG, 0.06% for AT-AC, and 0.09% for minor non-canonical splice sites. RNA-Seq data sets of 35 species were incorporated to validate non-canonical splice site predictions through gaps in sequencing reads alignments and to demonstrate the expression of affected genes.

CONCLUSION

We conclude that bona fide non-canonical splice sites are present and appear to be functionally relevant in most plant genomes, although at low abundance.

RBM4a modulates the impact of PRDM16 on development of brown adipocytes through an alternative splicing mechanism

Brown adipocytes (BAs) exhibit an energy-expending signature that is important in balancing metabolic homeostasis. In this study, results of transcriptome analyses revealed the reprogrammed splicing profile of the PR domain containing 16 (PRDM16) gene, a key transcription factor involved in brown adipogenesis, throughout development of wild-type brown adipose tissues (BATs). Moreover, discriminative splicing patterns of PRDM16 transcripts were noted in embryonic and postnatal RBM4a^-/- BATs. Overexpression of RBM4a enhanced the relative levels of PRDM16^{-ex 16} transcripts by simultaneously interacting with exonic and intronic CU elements, which encoded the PRDM16_S isoform containing a distinct C-terminus. The presence of the overexpressed PRDM16_S isoform showed a stronger effect than the overexpressed PRDM16_L isoform on enhancing transcriptional activity of the RBM4a and the PGC-1α promoter. Overexpression of the PRDM16_S isoform exerted more-prominent effects on enhancing the BAT-related gene program and energy expenditure compared to those of PRDM16_L-overexpressing cells. Our studies demonstrated that RBM4a-regulated alternative splicing constituted another regulatory mechanism for strengthening the influence of PRDM16 on the development of brown adipocytes.