Characterizing the interplay between gene nucleotide composition bias and splicing

Reduced protein-coding transcript diversity in severe dengue emphasises the role of alternative splicing

Dengue fever, a neglected tropical arboviral disease, has emerged as a global health concern in the past decade. Necessitating a nuanced comprehension of the intricate dynamics of host-virus interactions influencing disease severity, we analysed transcriptomic patterns using bulk RNA-seq from 112 age- and gender-matched NS1 antigen-confirmed hospital-admitted dengue patients with varying severity. Severe cases exhibited reduced platelet count, increased lymphocytosis, and neutropenia, indicating a dysregulated immune response. Using bulk RNA-seq, our analysis revealed a minimal overlap between the differentially expressed gene and transcript isoform, with a distinct expression pattern across the disease severity. Severe patients showed enrichment in retained intron and nonsense-mediated decay transcript biotypes, suggesting altered splicing efficiency. Furthermore, an up-regulated programmed cell death, a haemolytic response, and an impaired interferon and antiviral response at the transcript level were observed. We also identified the potential involvement of the RBM39 gene among others in the innate immune response during dengue viral pathogenesis, warranting further investigation. These findings provide valuable insights into potential therapeutic targets, underscoring the importance of exploring transcriptomic landscapes between different disease sub-phenotypes in infectious diseases.

Genetic Diversity in Bronchial Asthma Susceptibility: Exploring the Role of Vitamin D Receptor Gene Polymorphisms in Varied Geographic Contexts

Bronchial asthma (BA) exhibits varying prevalence across global populations, prompting a comprehensive investigation into genetic and environmental determinants. Vitamin D is a potent immunomodulator capable of suppressing inflammatory signals in several cell types involved in the asthmatic response; it exerts effects on the immune system by binding to the nuclear vitamin D receptor (VDR). VDR gene genetic variations are affecting serum vitamin D levels with a possible role in the BA risk. The current study aimed to examine the complex interaction of various factors (genetic background, serum vitamin D levels, and geographic location) to identify differences in the influence of these factors on the susceptibility to asthma between populations at different latitudes. Focusing on Eastern European cohorts from Latvia and Lithuania and comparing them with published data on East Asian populations, we explore the impact of VDR gene polymorphisms on BA susceptibility. Genotyping four key VDR SNPs and assessing their association with 25-hydroxyvitamin D levels, our study unveils significant associations of the studied loci with the risk of asthma-both risk-reducing and increasing effects, differently distributed between Baltic and East Asian populations. The functional effects of in silico VDR gene genetic variations are also identified and discussed.

SplicingLore: a web resource for studying the regulation of cassette exons by human splicing factors

One challenge faced by scientists from the alternative RNA splicing field is to decode the cooperative or antagonistic effects of splicing factors (SFs) to understand and eventually predict splicing outcomes on a genome-wide scale. In this manuscript, we introduce SplicingLore, an open-access database and web resource that help to fill this gap in a straightforward manner. The database contains a collection of RNA-sequencing-derived lists of alternative exons regulated by a total of 75 different SFs. All datasets were processed in a standardized manner, ensuring valid comparisons and correlation analyses. The user can easily retrieve a factor-specific set of differentially included exons from the database or provide a list of exons and search which SF(s) control(s) their inclusion. Our simple workflow is fast and easy to run, and it ensures a reliable calculation of correlation scores between the tested datasets. As a proof of concept, we predicted and experimentally validated a novel functional cooperation between the RNA helicases DDX17 and DDX5 and the heterogeneous nuclear ribonucleoprotein C (HNRNPC) protein. SplicingLore is available at https://splicinglore.ens-lyon.fr/. Database URL:  https://splicinglore.ens-lyon.fr/.

Glucose dissociates DDX21 dimers to regulate mRNA splicing and tissue differentiation

Glucose is a universal bioenergy source; however, its role in controlling protein interactions is unappreciated, as are its actions during differentiation-associated intracellular glucose elevation. Azido-glucose click chemistry identified glucose binding to a variety of RNA binding proteins (RBPs), including the DDX21 RNA helicase, which was found to be essential for epidermal differentiation. Glucose bound the ATP-binding domain of DDX21, altering protein conformation, inhibiting helicase activity, and dissociating DDX21 dimers. Glucose elevation during differentiation was associated with DDX21 re-localization from the nucleolus to the nucleoplasm where DDX21 assembled into larger protein complexes containing RNA splicing factors. DDX21 localized to specific SCUGSDGC motif in mRNA introns in a glucose-dependent manner and promoted the splicing of key pro-differentiation genes, including GRHL3, KLF4, OVOL1, and RBPJ. These findings uncover a biochemical mechanism of action for glucose in modulating the dimerization and function of an RNA helicase essential for tissue differentiation.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

False gene and chromosome losses in genome assemblies caused by GC content variation and repeats

BACKGROUND

Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements.

RESULTS

Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies.

CONCLUSIONS

Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

RNA helicase-dependent gene looping impacts messenger RNA processing

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

The endogenous HBZ interactome in ATL leukemic cells reveals an unprecedented complexity of host interacting partners involved in RNA splicing

Adult T-cell leukemia/lymphoma (ATL) is a T-cell lymphoproliferative neoplasm caused by the human T-cell leukemia virus type 1 (HTLV-1). Two viral proteins, Tax-1 and HBZ play important roles in HTLV-1 infectivity and in HTLV-1-associated pathologies by altering key pathways of cell homeostasis. However, the molecular mechanisms through which the two viral proteins, particularly HBZ, induce and/or sustain the oncogenic process are still largely elusive. Previous results suggested that HBZ interaction with nuclear factors may alter cell cycle and cell proliferation. To have a more complete picture of the HBZ interactions, we investigated in detail the endogenous HBZ interactome in leukemic cells by immunoprecipitating the HBZ-interacting complexes of ATL-2 leukemic cells, followed by tandem mass spectrometry analyses. RNA seq analysis was performed to decipher the differential gene expression and splicing modifications related to HTLV-1. Here we compared ATL-2 with MOLT-4, a non HTLV-1 derived leukemic T cell line and further compared with HBZ-induced modifications in an isogenic system composed by Jurkat T cells and stably HBZ transfected Jurkat derivatives. The endogenous HBZ interactome of ATL-2 cells identified 249 interactors covering three main clusters corresponding to protein families mainly involved in mRNA splicing, nonsense-mediated RNA decay (NMD) and JAK-STAT signaling pathway. Here we analyzed in detail the cluster involved in RNA splicing. RNAseq analysis showed that HBZ specifically altered the transcription of many genes, including crucial oncogenes, by affecting different splicing events. Consistently, the two RNA helicases, members of the RNA splicing family, DDX5 and its paralog DDX17, recently shown to be involved in alternative splicing of cellular genes after NF-κB activation by HTLV-1 Tax-1, interacted and partially co-localized with HBZ. For the first time, a complete picture of the endogenous HBZ interactome was elucidated. The wide interaction of HBZ with molecules involved in RNA splicing and the subsequent transcriptome alteration strongly suggests an unprecedented complex role of the viral oncogene in the establishment of the leukemic state.

EWS splicing regulation contributes to balancing Foxp1 isoforms required for neuronal differentiation

Alternative splicing is a key regulatory process underlying the amplification of genomic information and the expansion of proteomic diversity, particularly in brain. Here, we identify the Ewing sarcoma protein (EWS) as a new player of alternative splicing regulation during neuronal differentiation. Knockdown of EWS in neuronal progenitor cells leads to premature differentiation. Transcriptome profiling of EWS-depleted cells revealed global changes in splicing regulation. Bioinformatic analyses and biochemical experiments demonstrated that EWS regulates alternative exons in a position-dependent fashion. Notably, several EWS-regulated splicing events are physiologically modulated during neuronal differentiation and EWS depletion in neuronal precursors anticipates the splicing-pattern of mature neurons. Among other targets, we found that EWS controls the alternative splicing of the forkhead family transcription factor FOXP1, a pivotal transcriptional regulator of neuronal differentiation, possibly contributing to the switch of gene expression underlying the neuronal differentiation program.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

The oncogenic kinase NEK2 regulates an RBFOX2-dependent pro-mesenchymal splicing program in triple-negative breast cancer cells

BACKGROUND

Triple-negative breast cancer (TNBC) is the most heterogeneous and malignant subtype of breast cancer (BC). TNBC is defined by the absence of expression of estrogen, progesterone and HER2 receptors and lacks efficacious targeted therapies. NEK2 is an oncogenic kinase that is significantly upregulated in TNBC, thereby representing a promising therapeutic target. NEK2 localizes in the nucleus and promotes oncogenic splice variants in different cancer cells. Notably, alternative splicing (AS) dysregulation has recently emerged as a featuring trait of TNBC that contributes to its aggressive phenotype.

METHODS

To investigate whether NEK2 modulates TNBC transcriptome we performed RNA-sequencing analyses in a representative TNBC cell line (MDA-MB-231) and results were validated in multiple TNBC cell lines. Bioinformatics and functional analyses were carried out to elucidate the mechanism of splicing regulation by NEK2. Data from The Cancer Genome Atlas were mined to evaluate the potential of NEK2-sensitive exons as markers to identify the TNBC subtype and to assess their prognostic value.

RESULTS

Transcriptome analysis revealed a widespread impact of NEK2 on the transcriptome of TNBC cells, with 1830 AS events that are susceptible to its expression. NEK2 regulates the inclusion of cassette exons in splice variants that discriminate TNBC from other BC and that correlate with poor prognosis, suggesting that this kinase contributes to the TNBC-specific splicing program. NEK2 elicits its effects by modulating the expression of the splicing factor RBFOX2, a well-known regulator of epithelial to mesenchymal transition (EMT). Accordingly, NEK2 splicing-regulated genes are enriched in functional terms related to cell adhesion and contractile cytoskeleton and NEK2 depletion in mesenchymal TNBC cells induces phenotypic and molecular traits typical of epithelial cells. Remarkably, depletion of select NEK2-sensitive splice-variants that are prognostic in TNBC patients is sufficient to interfere with TNBC cell morphology and motility, suggesting that NEK2 orchestrates a pro-mesenchymal splicing program that modulates migratory and invasive properties of TNBC cells.

CONCLUSIONS

Our study uncovers an extensive splicing program modulated by NEK2 involving splice variants that confer an invasive phenotype to TNBCs and that might represent, together with NEK2 itself, valuable therapeutic targets for this disease.

A Dynamic Splicing Program Ensures Proper Synaptic Connections in the Developing Cerebellum

Tight coordination of gene expression in the developing cerebellum is crucial for establishment of neuronal circuits governing motor and cognitive function. However, transcriptional changes alone do not explain all of the switches underlying neuronal differentiation. Here we unveiled a widespread and highly dynamic splicing program that affects synaptic genes in cerebellar neurons. The motifs enriched in modulated exons implicated the splicing factor Sam68 as a regulator of this program. Sam68 controls splicing of exons with weak branchpoints by directly binding near the 3' splice site and competing with U2AF recruitment. Ablation of Sam68 disrupts splicing regulation of synaptic genes associated with neurodevelopmental diseases and impairs synaptic connections and firing of Purkinje cells, resulting in motor coordination defects, ataxia, and abnormal social behavior. These findings uncover an unexpectedly dynamic splicing regulatory network that shapes the synapse in early life and establishes motor and cognitive circuitry in the developing cerebellum.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Identification of long regulatory elements in the genome of Plasmodium falciparum and other eukaryotes

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

An Intergenic rs9275596 Polymorphism on Chr. 6p21 Is Associated with Multiple Sclerosis in Latvians

Background and objectives: Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system, leading to demyelination of neurons and potentially debilitating physical and mental symptoms. The disease is more prevalent in women than in men. The major histocompatibility complex (MHC) region has been identified as a major genetic determinant for autoimmune diseases, and its role in some neurological disorders including MS was evaluated. An intergenic single-nucleotide polymorphism (SNP), rs9275596, located between the HLA-DQB1 and HLA-DQA2 genes, is in significant association with various autoimmune diseases according to genome-wide association studies (GWASs). A cumulative effect of this SNP with other polymorphisms from this region was revealed. The aim of the study was to verify the data on rs9275596 association in multiple sclerosis in a case/control study of the Latvian population and to evaluate eventual functional significance of allele substitutions. Materials and Methods: rs9275596 (chr6:32713854; GRCh38.p12) was genotyped in 273 MS patients and 208 controls on main and sex-specific associations. Eventual functional significance of allele substitutions was evaluated in silico using publicly available tools. Results: The rs9275596 rare alleles were identified as a disease susceptibility factor in association with the MS main group and in affected females (p < 0.001 and p < 0.01, respectively). Risk factor genotypes with rare alleles included were associated with the MS common cohort (p < 0.002) and female cohort (odds ratio, OR = 2.24) and were identified as disease susceptible in males (OR = 2.41). It was shown that structural changes of rs9275596 affect the secondary structure of DNA. Functional significance of allele substitutions was evaluated on the eventual sequence affinity to transcription factors (TFs) and splicing signals similarity. A possible impact of the particular polymorphisms on the transcription and splicing efficiency is discussed. Conclusions: Our results suggest susceptibility of rs9275596 to multiple sclerosis in Latvians.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.