Regulation by 3'-Untranslated Regions

Molecular Mimicry of Transposable Elements in Plants

Transposable elements (TEs) are mobile DNA elements that are particularly abundant in the plant genomes. They have long been considered as junk DNA; however, a growing body of evidence suggests that TE insertions promote genetic diversity that is essential for the adaptive evolution of a species. Thus far, studies have mainly investigated the cis-acting regulatory roles of TEs generated by their insertions nearby or within the host genes. However, the trans-acting effects of TE-derived RNA and DNA remained obscure to date. TEs contain various regulatory elements within their sequences that can accommodate the binding of specific RNAs and proteins. Recently, it was suggested that some of these cellular regulators are shared between TEs and the host genes, and the competition for the common host factors underlies the fine-tuned developmental reprogramming. In this review, we will highlight and discuss the latest discoveries on the biological functions of plant TEs, with a particular focus on their competitive binding with specific developmental regulators.

RNA processing in innate immunity: regulation by RNA-binding proteins

RNA processing is an important but often overlooked process in determining protein expression. Alternative polyadenylation and regulated mRNA stability control the amount and duration of protein expression, while alternative splicing also controls protein identity and function. Much work in innate immunity has focused on the activation of transcription factors and the downstream consequences in gene expression. However, there is increasing evidence indicating that regulation of RNA processing by RNA-binding proteins (RBPs) contributes significantly to tuning the innate immune response. Herein we review work highlighting the impact of RNA processing in innate immunity and describe the RBPs and mechanisms driving this regulation. We conclude with a discussion of unanswered questions to motivate continued research in this important and understudied field.

The phenotypic spectrum of the Cornelia de Lange-like "Alazami-Yuan syndrome": A case report of the 7th diagnosed individual and review of the literature

We present a 17-year-old female with a childhood clinical diagnosis of Cornelia de Lange Syndrome (CdLS), however later genetic testing identified compound heterozygous variants in TAF6, consistent with AYS. This case report adds to the phenotypic spectrum observed in AYS, and draws connections to transcriptional pathways between CdLS and AYS.

Morphogenesis, starvation, and light responses in a mushroom-forming fungus revealed by long-read sequencing and extensive expression profiling

Mushroom-forming fungi (Agaricomycetes) are emerging as pivotal players in several fields of science and industry. Genomic data for Agaricomycetes are accumulating rapidly; however, this is not paralleled by improvements of gene annotations, which leave gene function notoriously poorly understood. We set out to improve our functional understanding of the model mushroom Coprinopsis cinerea by integrating a new, chromosome-level assembly, high-quality gene predictions, and functional information derived from broad gene-expression profiling data. The new annotation includes 5' and 3' untranslated regions (UTRs), polyadenylation sites (PASs), upstream open reading frames (uORFs), splicing isoforms, and microexons, as well as core gene sets corresponding to carbon starvation, light response, and hyphal differentiation. As a result, the genome of C. cinerea has now become the most comprehensively annotated genome among mushroom-forming fungi, which will contribute to multiple rapidly expanding fields, including research on their life history, light and stress responses, as well as multicellular development.

Profiling the regulatory landscape of sialylation through miRNA targeting of CMP- sialic acid synthetase

Cell surface sialic acid is an important glycan modification that contributes to both normal and pathological physiology. The enzyme cytidine monophosphate N-acetylneuraminic acid synthetase (CMAS) biosynthesizes the activated sugar donor cytidine monophosphate (CMP) sialic acid, which is required for all sialylation. CMAS levels impact sialylation with corresponding biological effects. The mechanisms that regulate CMAS are relatively uncharacterized. Herein, we use a high throughput genetically encoded fluorescence assay (miRFluR) to comprehensively profile the posttranscriptional regulation of CMAS by miRNA. These small non-coding RNAs have been found to impact glycosylation. Mapping the interactions of the human miRNAome with the 3'-untranslated region of CMAS, we identified miRNA whose impact on CMAS expression was either downregulatory or upregulatory. This follows previous work from our laboratory and others showing that miRNA regulation is bidirectional. Validation of the high-throughput results confirmed our findings. We also identified the direct binding sites for two upregulatory and two downregulatory miRNAs. Functional enrichment analysis for miRNAs upregulating CMAS revealed associations with pancreatic cancer, where sialic acid metabolism and the α-2,6-sialyltransferase ST6GAL1 have been found to be important. We found that miRNA associated with the enriched signature enhanced pancreatic cell-surface α-2,6-sialylation via CMAS expression in the absence of effects on ST6GAL1. We also find overlap between the miRNA regulation of CMAS and that of previously analyzed sialyltransferases. Overall, our work points to the importance of miRNA in regulating sialylation levels in disease and add further evidence to the bidirectional nature of miRNA regulation.

Pumilio differentially binds to mRNA 3' UTR isoforms to regulate localization of synaptic proteins

In neuronal cells, the regulation of RNA is crucial for the spatiotemporal control of gene expression, but how the correct localization, levels, and function of synaptic proteins are achieved is not well understood. In this study, we globally investigate the role of alternative 3' UTRs in regulating RNA localization in the synaptic regions of the Drosophila brain. We identify direct mRNA targets of the translational repressor Pumilio, finding that mRNAs bound by Pumilio encode proteins enriched in synaptosomes. Pumilio differentially binds to RNA isoforms of the same gene, favoring long, neuronal 3' UTRs. These longer 3' UTRs tend to remain in the neuronal soma, whereas shorter UTR isoforms localize to the synapse. In cultured pumilio mutant neurons, axon outgrowth defects are accompanied by mRNA isoform mislocalization, and proteins encoded by these Pumilio target mRNAs display excessive abundance at synaptic boutons. Our study identifies an important mechanism for the spatiotemporal regulation of protein function in neurons.

3' untranslated region somatic variants connect alternative polyadenylation dysregulation in human cancers

Somatic variants in the cancer genome influence gene expression through diverse mechanisms depending on their specific locations. However, a systematic evaluation of the effects of somatic variants located in 3' untranslated regions (3' UTRs) on alternative polyadenylation (APA) of mRNA remains lacking. In this study, we analyze 10,199 tumor samples across 32 cancer types and identify 1,333 somatic single nucleotide variants (SNVs) associated with abnormal 3' UTR APA. Mechanistically, these 3' UTR SNVs can alter cis-regulatory elements, such as the poly(A) signal and UGUA motif, leading to changes in APA. Minigene assays confirm that 3' UTR SNVs in multiple genes, including RPS23 and CHTOP, induce aberrant APA. Among affected genes, 62 exhibit differential stability between tandem 3' UTR isoforms, including HSPA4 and UCK2, validated by experimental assays. Finally, we establish that SNV-related abnormal APA usage serves as an additional layer of expression regulation for tumor-suppressor gene HMGN2 in breast cancer. Collectively, this study reveals 3' UTR APA as a critical mechanism mediating the functional impact of somatic noncoding variants in human cancers.

Axonal RNA localization is essential for long-term memory

Localization of mRNAs to neuronal terminals, coupled to local translation, has emerged as a prevalent mechanism controlling the synaptic proteome. However, the physiological regulation and function of this process in the context of mature in vivo memory circuits has remained unclear. Here, we combined synaptosome RNA profiling with whole brain high-resolution imaging to uncover mRNAs with different localization patterns in the axons of Drosophila Mushroom Body memory neurons, some exhibiting regionalized, input-dependent, recruitment along axons. By integrating transcriptome-wide binding approaches and functional assays, we show that the conserved Imp RNA binding protein controls the transport of mRNAs to Mushroom Body axons and characterize a mutant in which this transport is selectively impaired. Using this unique mutant, we demonstrate that axonal mRNA localization is required for long-term, but not short-term, behavioral memory. This work uncovers circuit-dependent mRNA targeting in vivo and demonstrates the importance of local RNA regulation in memory consolidation.

Exploiting the SunTag system to study the developmental regulation of mRNA translation

The ability to quantitatively study mRNA translation using SunTag imaging is transforming our understanding of the translation process. Here, we expand the SunTag method to study new aspects of translation regulation in Drosophila. Repression of the maternal hunchback (hb) mRNA in the posterior of the Drosophila embryo is a textbook example of translational control. Using SunTag imaging to quantify translation of maternal SunTag-hb mRNAs, we show that repression in the posterior is leaky, as ∼5% of SunTag-hb mRNAs are translated. In the anterior of the embryo, the maternal and zygotic SunTag-hb mRNAs show similar translation efficiency despite having different untranslated regions (UTRs). We demonstrate that the SunTag-hb mRNA can be used as a reporter to study ribosome pausing at single-mRNA resolution, by exploiting the conserved xbp1 mRNA and A60 pausing sequences. Finally, we adapt the detector component of the SunTag system to visualise and quantify translation of the short gastrulation (sog) mRNA, encoding an essential secreted extracellular BMP regulator, at the endoplasmic reticulum in fixed and live embryos. Together, these tools will facilitate the future dissection of translation regulatory mechanisms during development.

An atlas of RNA-dependent proteins in cell division reveals the riboregulation of mitotic protein-protein interactions

Ribonucleoprotein complexes are dynamic assemblies of RNA with RNA-binding proteins, which modulate the fate of RNA. Inversely, RNA riboregulates the interactions and functions of the associated proteins. Dysregulation of ribonucleoprotein functions is linked to diseases such as cancer and neurological disorders. In dividing cells, RNA and RNA-binding proteins are present in mitotic structures, but their impact on cell division remains unclear. By applying the proteome-wide R-DeeP strategy to cells synchronized in mitosis versus interphase integrated with the RBP2GO knowledge, we provided an atlas of RNA-dependent proteins in cell division, accessible at R-DeeP3.dkfz.de. We uncovered AURKA, KIFC1 and TPX2 as unconventional RNA-binding proteins. KIFC1 was identified as a new substrate of AURKA, and new TPX2-interacting protein. Their pair-wise interactions were RNA dependent. In addition, RNA stimulated AURKA kinase activity and stabilized its conformation. In this work, we highlighted riboregulation of major mitotic factors as an additional complexity level of cell division.

The LSmAD Domain of Ataxin-2 Modulates the Structure and RNA Binding of Its Preceding LSm Domain

Ataxin-2 (Atx2), an RNA-binding protein, plays a pivotal role in the regulation of RNA, intracellular metabolism, and translation within the cellular environment. Although both the Sm-like (LSm) and LSm-associated (LSmAD) domains are considered to associated with RNA binding, there is still a lack of experimental evidence supporting their functions. To address this, we designed and constructed several recombinants containing the RNA-binding domain (RBD) of Atx2. By employing biophysical and biochemical techniques, such as EMSA and SHAPE chemical detection, we identified that LSm is responsible for RNA binding, whereas LSmAD alone does not bind RNA. NMR and small-angle X-ray scattering (SAXS) analyses have revealed that the LSmAD domain exhibits limited structural integrity and poor folding capability. The EMSA data confirmed that both LSm and LSm-LSmAD bind RNA, whereas LSmAD alone cannot, suggesting that LSmAD may serve as an auxiliary role to the LSm domain. SHAPE chemical probing further demonstrates that LSm binds to the AU-rich, GU-rich, or CU-rich sequence, but not to the CA-rich sequence. These findings indicate that Atx2 can interact with the U-rich sequences in the 3'-UTR, implicating its role in poly(A) tailing and the regulation of mRNA translation and degradation.

Identification of cis-sQTL demonstrates genetic associations and functional implications of inflammatory processes in Nelore cattle muscle tissue

This study aimed to identify splicing quantitative trait loci (cis-sQTL) in Nelore cattle muscle tissue and explore the involvement of spliced genes (sGenes) in immune system-related biological processes. Genotypic data from 80 intact male Nelore cattle were obtained using SNP-Chip technology, while RNA-Seq analysis was performed to measure gene expression levels, enabling the integration of genomic and transcriptomic datasets. The normalized expression levels of spliced transcripts were associated with single nucleotide polymorphisms (SNPs) through an analysis of variance using an additive linear model with the MatrixEQTL package. A permutation analysis then assessed the significance of the best SNPs for each spliced transcript. Functional enrichment analysis was performed on the sGenes to investigate their roles in the immune system. In total, 3,187 variants were linked to 3,202 spliced transcripts, with 83 sGenes involved in immune system processes. Of these, 31 sGenes were enriched for five transcription factors. Most cis-sQTL effects were found in intronic regions, with 27 sQTL variants associated with disease susceptibility and resistance in cattle. Key sGenes identified, such as GSDMA, NLRP6, CASP6, GZMA, CASP4, CASP1, TREM2, NLRP1, and NAIP, were related to inflammasome formation and pyroptosis. Additionally, genes like PIDD1, OPTN, NFKBIB, STAT1, TNIP3, and TREM2 were involved in regulating the NF-kB pathway. These findings lay the groundwork for breeding disease-resistant cattle and enhance our understanding of genetic mechanisms in immune responses.

Transcriptome-wide mapping of N3-methylcytidine modification at single-base resolution

3-Methylcytidine (m3C), a prevalent modification of transfer RNAs (tRNAs), was recently identified in eukaryotic messenger RNAs (mRNAs). However, its precise distribution and formation mechanisms in mRNAs remain elusive. Here, we develop a novel approach, m3C immunoprecipitation and sequencing (m3C-IP-seq), utilizing antibody enrichment to profile the m3C methylome at single-nucleotide resolution. m3C-IP-seq captures 12 cytoplasmic tRNA isoacceptors and 2 mitochondrial tRNA isoacceptors containing m3C modifications. Moreover, m3C-IP-seq permits the comprehensive profiling of m3C sites in mRNAs and long noncoding RNAs, with their presence reliant on a nuclear isoform of METTL8. A significant proportion of m3C sites is concentrated in the 3' untranslated region (3' UTR) of mRNAs and is associated with mRNA degradation. Additionally, m3C methylation is dynamic and responds to hypoxia. Collectively, our data demonstrate the widespread presence of m3C modification in the human transcriptome and provide a resource for functional studies of m3C-mediated RNA metabolism.

Molecular underpinnings of plasticity and supergene-mediated polymorphism in fire ant queens

Characterizing molecular underpinnings of plastic traits and balanced polymorphisms represent 2 important goals of evolutionary biology. Fire ant gynes (pre-reproductive queens) provide an ideal system to study potential links between these phenomena because they exhibit both supergene-mediated polymorphism and nutritional plasticity in weight and colony-founding behaviour. Gynes with the inversion supergene haplotype are lightweight and depend on existing workers to initiate reproduction. Gynes with only the ancestral, non-inverted gene arrangement accumulate more nutrient reserves as adults and, in a distinct colony-founding behaviour, initiate reproduction without help from workers. However, when such gynes overwinter in the natal nest they develop an environmentally induced lightweight phenotype and colony-founding behaviour, similar to gynes with the inversion haplotype that have not overwintered. To evaluate the extent of shared mechanisms between plasticity and balanced polymorphism in fire ant gyne traits, we assessed whether genes with expression variation linked to overwintering plasticity may be affected by the evolutionary divergence between supergene haplotypes. To do so, we first compared transcriptional profiles of brains and ovaries from overwintered and non-overwintered gynes to identify plasticity-associated genes. These genes were enriched for metabolic and behavioural functions. Next, we compared plasticity-associated genes to those differentially expressed by supergene genotype, revealing a significant overlap of the 2 sets in ovarian tissues. We also identified sequence substitutions between supergene variants of multiple plasticity-associated genes, consistent with a scenario in which an ancestrally plastic phenotype responsive to an environmental condition became increasingly genetically regulated.

Multilevel Gene Expression Changes in Lineages Containing Adaptive Copy Number Variants

Copy number variants (CNVs) are an important class of genetic variation that can mediate rapid adaptive evolution. Whereas, CNVs can increase the relative fitness of the organism, they can also incur a cost due to the associated increased gene expression and repetitive DNA. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that gene expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- chemostats. We find CNV-amplified genes correlate with higher mRNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of gene expression by the abundance of the preceding step we were able to identify widespread differences in the efficiency of each level of gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA-binding sites for Ssd1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which genome evolution shapes gene expression.

Targeting the CLK2/SRSF9 splicing axis in prostate cancer leads to decreased ARV7 expression

In advanced prostate cancer (PC), in particular after acquisition of resistance to androgen receptor (AR) signaling inhibitors (ARSI), upregulation of AR splice variants compromises endocrine therapy efficiency. Androgen receptor splice variant-7 (ARV7) is clinically the most relevant and has a distinct 3' untranslated region (3'UTR) compared to the AR full-length variant, suggesting a unique post-transcriptional regulation. Here, we set out to evaluate the applicability of the ARV7 3'UTR as a therapy target. A common single nucleotide polymorphism, rs5918762, was found to affect the splicing rate and thus the expression of ARV7 in cellular models and patient specimens. Serine/arginine-rich splicing factor 9 (SRSF9) was found to bind to and increase the inclusion of the cryptic exon 3 of ARV7 during the splicing process in the alternative C allele of rs5918762. The dual specificity protein kinase CLK2 interferes with the activity of SRSF9 by regulating its expression. Inhibition of the Cdc2-like kinase (CLK) family by the small molecules cirtuvivint or lorecivivint results in the decreased expression of ARV7. Both inhibitors show potent anti-proliferative effects in enzalutamide-treated or -naive PC models. Thus, targeting aberrant alternative splicing at the 3'UTR of ARV7 by disturbing the CLK2/SRSF9 axis might be a valuable therapeutic approach in late stage, ARSI-resistant PC.

Alternative polyadenylation in cancer: Molecular mechanisms and clinical application

Alternative polyadenylation (APA) serves as a crucial mechanism for the posttranscriptional regulation of gene expression and influences gene expression by generating diverse mRNA isoforms. This process is regulated by a diverse array of RNA-binding proteins (RBPs), which selectively bind to specific sequences or structures within the pre-mRNA molecule. Dysregulation of APA and its associated RBPs has been implicated in numerous diseases, including cardiovascular diseases, nervous system disease, and cancer. For instance, aberrant APA events have been observed in several types of tumors, contributing to tumor heterogeneity and affecting key cellular pathways involved in cell proliferation, invasion, metastasis, and response to therapy. This review critically evaluates the current understanding of APA mechanisms and the multifaceted roles of RBPs in orchestrating this intricate process. We highlight recent advancements in high-throughput sequencing and bioinformatics tools that have enhanced our ability to study APA on a genome-wide scale. Moreover, we explored the pathological consequences of APA dysregulation, emphasizing its role in oncogenesis. By elucidating the intricate relationships between APA and RBPs, this review aims to underscore the potential of targeting the APA machinery and RBPs for therapeutic intervention. Understanding these molecular processes holds promise for developing novel diagnostic markers and treatment strategies for a range of human cancers.

Characterization and molecular targeting of CFIm25 (NUDT21/CPSF5) mRNA using miRNAs

Changes in protein levels of the mammalian cleavage factor, CFIm25, play a role in regulating pathological processes including neural dysfunction, fibrosis, and tumorigenesis. However, despite these effects, little is known about how CFIm25 (NUDT21) expression is regulated at the RNA level. A potential regulator of NUDT21 mRNA are small non-coding microRNAs (miRNAs). In general, miRNAs bind to the 3'untranslated regions (3'UTRs) and can target the bound mRNA for degradation or inhibit translation thus affecting the levels of protein in cells. Interestingly, a mechanism known as alternative polyadenylation (APA) enables mRNAs to escape miRNA regulation by generating mRNAs with 3'UTRs of different sizes. As many miRNA target sites are located within the 3'UTR, shortening the 3'UTR allows mRNAs to evade miRNAs targeting this region. The differences in the lengths and the sequence composition of the 3'UTRs may also impact the mRNA's translatability and subcellular localization. APA has been reported to regulate over 70% of protein coding genes, thus increasing the transcript repertoire. Several proteins, including mammalian cleavage factor, CFIm25 (NUDT21), have been shown to regulate APA. In this study we wanted to determine whether CFIm25 (NUDT21), itself a regulator of APA, undergoes APA to evade miRNA regulation. We used the blood cancer mantle cell lymphoma (MCL) cells as a model and showed that in these cells, NUDT21 is relatively stable with a long half-life. In addition, the NUDT21 pre-mRNA undergoes alternative APA within the same terminal exon. The three different sized NUDT21 mRNAs have different 3'UTR lengths and they each use a different canonical polyadenylation signal, AAUAAA, for 3'end cleavage and polyadenylation. Use of miRNA mimics and inhibitors showed that miR-23a, miR-222, and miR-323a play a significant role in regulating NUDT21 expression. Hence, these results suggest that NUDT21 mRNA is stable and the different 3'UTRs generated through APA of NUDT21 play an important role in evading miRNA regulation and offers insights into how levels of CFIm25 (NUDT21) may be fine-tuned as needed under different physiological and pathological conditions.

mRNA-LM: full-length integrated SLM for mRNA analysis

The success of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) messenger RNA (mRNA) vaccine has led to increased interest in the design and use of mRNA for vaccines and therapeutics. Still, selecting the most appropriate mRNA sequence for a protein remains a challenge. Several recent studies have shown that the specific mRNA sequence can have a significant impact on the translation efficiency, half-life, degradation rates, and other issues that play a major role in determining vaccine efficiency. To enable the selection of the most appropriate sequence, we developed mRNA-LM, an integrated small language model for modeling the entire mRNA sequence. mRNA-LM uses the contrastive language-image pretraining integration technology to combine three separate language models for the different mRNA segments. We trained mRNA-LM on millions of diverse mRNA sequences from several different species. The unsupervised model was able to learn meaningful biology related to evolution and host-pathogen interactions. Fine-tuning of mRNA-LM allowed us to use it in several mRNA property prediction tasks. As we show, using the full-length integrated model led to accurate predictions, improving on prior methods proposed for this task.

Impact of rare non-coding variants on human diseases through alternative polyadenylation outliers

Although rare non-coding variants (RVs) play crucial roles in complex traits and diseases, understanding their mechanisms and identifying disease-associated RVs continue to be major challenges. Here we constructed a comprehensive atlas of alternative polyadenylation (APA) outliers (aOutliers), including 1334 3' UTR and 200 intronic aOutliers, from 15,201 samples across 49 human tissues. These aOutliers exhibit unique characteristics from transcription or splicing outliers, with a pronounced RV enrichment. Mechanistically, aOutlier-RVs alter poly(A) signals and splicing sites, and perturbation indeed triggers APA events. Furthermore, we developed a Bayesian-based APA RV prediction model, which successfully pinpointed a specific set of 1799 RVs impacting 278 genes with significantly large disease effect sizes. Notably, we observed a convergence effect between rare and common cancer variants, exemplified by regulation in the DDX18 gene. Together, this study introduced an APA-enhanced framework for genome annotation, underscoring APA's role in uncovering functional RVs linked to complex traits and diseases.

Heterogeneous nuclear ribonucleoprotein C promotes non-small cell lung cancer progression by enhancing XB130 mRNA stability and translation

BACKGROUND

XB130, a classical adaptor protein, exerts a critical role in diverse cellular processes. Aberrant expression of XB130 is closely associated with tumorigenesis and aggressiveness. However, the mechanisms governing its expression regulation remain poorly understood. Heterogeneous nuclear ribonucleoprotein C (hnRNPC), as an RNA-binding protein, is known to modulate multiple aspects of RNA metabolism and has been implicated in the pathogenesis of various cancers. We have previously discovered that hnRNPC is one of the candidate proteins that interact with the 3' untranslated region (3'UTR) of XB130 in non-small cell lung cancer (NSCLC). Therefore, this study aims to comprehensively elucidate how hnRNPC regulates the expression of XB130 in NSCLC.

MATERIALS AND METHODS

We evaluated the expression of hnRNPC in cancer and assessed the correlation between hnRNPC expression and prognosis in cancer patients using public databases. Subsequently, several stable cell lines were constructed. The proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT) of these cells were detected through Real-time cellular analysis, adherent colony formation, wound healing assay, invasion assay, and Western blotting. The specific regulatory manner between hnRNPC and XB130 was investigated by Real-time quantitative PCR, Western blotting, RNA pull‑down assay, dual‑luciferase reporter assay, RNA immunoprecipitation, and Co-Immunoprecipitation.

RESULTS

We identified that hnRNPC expression is significantly elevated in NSCLC and correlates with poor prognosis in patients with lung adenocarcinoma. HnRNPC overexpression in NSCLC cells increased the expression of XB130, subsequently activating the PI3K/Akt signaling pathway and ultimately promoting cell proliferation and EMT. Additionally, overexpressing XB130 in hnRNPC-silenced cells partially restored cell proliferation and EMT. Mechanistically, hnRNPC specifically bound to the 3'UTR segments of XB130 mRNA, enhancing mRNA stability by inhibiting the recruitment of nucleases 5'-3' exoribonuclease 1 (XRN1) and DIS3-like 3'-5' exoribonuclease 2 (DIS3L2). Furthermore, hnRNPC simultaneously interacted with the eukaryotic initiation factor 4E (eIF4E), a component of the eIF4F complex, facilitating the circularization of XB130 mRNA and thereby increasing its translation efficiency.

CONCLUSIONS

HnRNPC overexpression promotes NSCLC progression by enhancing XB130 mRNA stability and translation, suggesting that hnRNPC might be a potential therapeutic and prognostic target for NSCLC.

Screening the human miRNA interactome reveals coordinated up-regulation in melanoma, adding bidirectional regulation to miRNA networks

Cellular protein expression is coordinated posttranscriptionally by an intricate regulatory network. The current presumption is that microRNAs (miRNAs) work by repression of functionally related targets within a system. In recent work, up-regulation of protein expression via direct interactions of messenger RNA with miRNA has been found in dividing cells, providing an additional mechanism of regulation. Herein, we demonstrate coordinated up-regulation of functionally coupled proteins by miRNA. We focused on CD98hc, the heavy chain of the amino acid transporter LAT-1, and α-2,3-sialyltransferases ST3GAL1 and ST3GAL2, which are critical for CD98hc stability in melanoma. Profiling miRNA regulation using our high-throughput miRFluR assay, we identified miRNA that up-regulated the expression of both CD98hc and either ST3GAL1 or ST3GAL2. These co-up-regulating miRNAs were enriched in melanoma datasets associated with transformation and progression. Our findings add co-up-regulation by miRNA into miRNA regulatory networks and add a bidirectional twist to the impact miRNAs have on protein regulation and glycosylation.

Detecting gene expression in Caenorhabditis elegans

Reliable methods for detecting and analyzing gene expression are necessary tools for understanding development and investigating biological responses to genetic and environmental perturbation. With its fully sequenced genome, invariant cell lineage, transparent body, wiring diagram, detailed anatomy, and wide array of genetic tools, Caenorhabditis elegans is an exceptionally useful model organism for linking gene expression to cellular phenotypes. The development of new techniques in recent years has greatly expanded our ability to detect gene expression at high resolution. Here, we provide an overview of gene expression methods for C. elegans, including techniques for detecting transcripts and proteins in situ, bulk RNA sequencing of whole worms and specific tissues and cells, single-cell RNA sequencing, and high-throughput proteomics. We discuss important considerations for choosing among these techniques and provide an overview of publicly available online resources for gene expression data.

Pharmacogenetic and pharmacokinetic factors for dexmedetomidine-associated hemodynamic instability in pediatric patients

Purpose

The incidence of hemodynamic instability associated with dexmedetomidine (DEX) sedation has been reported to exceed 50%, with substantial inter-individual variability in response. Genetic factors have been suggested to contribute significantly to such variation. The aim of this study was to identify the clinical, pharmacokinetic, and genetic factors associated with DEX-induced hemodynamic instability in pediatric anesthesia patients.

Methods

A cohort of 270 pediatric patients scheduled for elective interventional surgery received an intranasal dose of 3 mcg·kg-1 of dexmedetomidine, and subsequent propofol induction was conducted when patients had a UMSS of 2-4. The primary endpoint was hemodynamic instability-defined as a composite of hypotension and/or bradycardia, which is characterized by a 20% reduction from age-specific baseline values. Plasma concentrations of dexmedetomidine were determined, and single-nucleotide polymorphisms (SNPs) were genotyped. A validated population pharmacokinetic model was used to estimate pharmacokinetic parameters. LASSO regression was used to identify significant factors, and a Cox's proportional hazards model-derived nomogram for hemodynamic instability was developed.

Results

Hemodynamic instability was observed in 52 out of 270 patients (209 events), resulting in a cumulative incidence of 16.30% at 90 min, as estimated by Kaplan-Meier estimation, and it was associated with a median time to event of 35 min. The interval time between DEX initiation and propofol induction was 16 min (IQR: 12-22 min). The cumulative incidence was 8.2% within 22 min after DEX initiation. The identified significant risk factors for DEX-associated hemodynamic instability included weight, DEX clearance, concomitant propofol use, and the following gene variants UGT2B10 rs1841042 (hazard ratio (HR):1.41, 95% confidence interval (CI): 1.12-1.79), CYP2A6 rs8192733 (HR:0.28, 95%CI:0.09-0.88), ADRA2B rs3813662 (HR:1.39,95%CI:1.02-1.89), CACNA2D2 rs2236957 (HR:1.46, 95%CI:1.09-1.96), NR1I2 rs3814057 (HR:0.64, 95%CI:0.43-0.95), and CACNB2 rs10764319 (HR:1.40,95%CI:1.05-1.87). The areas under the curve for the training and test cohorts were 0.881 and 0.762, respectively. The calibration curve indicated excellent agreement.

Conclusion

The predictive nomogram, which incorporates genetic variants (UGT2B10, CYP2A6, ADRA2B, CACNA2D2, NR1I2, and CACNB2) along with clinical factors such as weight, DEX clearance, and propofol use, may help prevent DEX-associated hemodynamic instability. Delayed hemodynamic instability is likely to occur after 35-min DEX initiation in patients with lower DEX clearance after propofol induction.

Mechanistic Evaluation of miRNAs and Their Targeted Genes in the Pathogenesis and Therapeutics of Parkinson's Disease

MicroRNA (miRNA) are usually 18-25 nucleotides long non-coding RNA targeting post-transcriptional regulation of genes involved in various biological processes. The function of miRNA is essential for maintaining a homeostatic cellular condition, regulating autophagy, cellular motility, and inflammation. Dysregulation of miRNA is responsible for multiple disorders, including neurodegeneration, which has emerged as a severe problem in recent times and has verified itself as a life-threatening condition that can be understood by the continuous destruction of neurons affecting various cognitive and motor functions. Parkinson's disease (PD) is the second most common, permanently debilitating neurodegenerative disorder after Alzheimer's, mainly characterized by uncontrolled tremor, stiffness, bradykinesia or akinesia (slowness in movement), and post-traumatic stress disorder. PD is mainly caused by the demolition of the primary dopamine neurotransmitter secretory cells and dopaminergic or dopamine secretory neurons in the substantia nigra pars compacta of the midbrain, which are majorly responsible for motor functions. In this study, a systematic evaluation of research articles from year 2017 to 2022 was performed on multiple search engines, and lists of miRNA being dysregulated in PD in different body components were generated. This study highlighted miR-7, miR-124, miR-29 family, and miR-425, showing altered expression levels during PD's progression, further regulating the expression of multiple genes responsible for PD.

Assessment of mRNA Decay and Calculation of Codon Occurrence to mRNA Stability Correlation Coefficients after 5-EU Metabolic Labeling

mRNA translation and decay are tightly connected. This chapter describes a method to assess the influence of each codon identity on mRNA stability in cultured cells. The technique involves metabolic labeling of the nascent mRNAs by addition of the nucleoside analog 5-ethynyluridine (5-EU), purification of the RNA at different time-points after chase of the 5-EU, then biotinylation with Click chemistry, pull-down, and sequencing. The transcripts' half-lives are calculated from the expression level of each mRNA at the different time-points. Finally, the method describes the calculation of the Codon occurrence to mRNA Stability correlation Coefficient, or CSC, as a correlation between the codon occurrence in a transcript and the transcript half-life, for each codon.

Recent progress in mRNA cancer vaccines

The research and development of messenger RNA (mRNA) cancer vaccines have gradually overcome numerous challenges through the application of personalized cancer antigens, structural optimization of mRNA, and the development of alternative RNA-based vectors and efficient targeted delivery vectors. Clinical trials are currently underway for various cancer vaccines that encode tumor-associated antigens (TAAs), tumor-specific antigens (TSAs), or immunomodulators. In this paper, we summarize the optimization of mRNA and the emergence of RNA-based expression vectors in cancer vaccines. We begin by reviewing the advancement and utilization of state-of-the-art targeted lipid nanoparticles (LNPs), followed by presenting the primary classifications and clinical applications of mRNA cancer vaccines. Collectively, mRNA vaccines are emerging as a central focus in cancer immunotherapy, offering the potential to address multiple challenges in cancer treatment, either as standalone therapies or in combination with current cancer treatments.

Crosstalk between metabolic and epigenetic modifications during cell carcinogenesis

Genetic mutations arising from various internal and external factors drive cells to become cancerous. Cancerous cells undergo numerous changes, including metabolic reprogramming and epigenetic modifications, to support their abnormal proliferation. This metabolic reprogramming leads to the altered expression of many metabolic enzymes and the accumulation of metabolites. Recent studies have shown that these enzymes and metabolites can serve as substrates or cofactors for chromatin-modifying enzymes, thereby participating in epigenetic modifications and promoting carcinogenesis. Additionally, epigenetic modifications play a role in the metabolic reprogramming and immune evasion of cancer cells, influencing cancer progression. This review focuses on the origins of cancer, particularly the metabolic reprogramming of cancer cells and changes in epigenetic modifications. We discuss how metabolites in cancer cells contribute to epigenetic remodeling, including lactylation, acetylation, succinylation, and crotonylation. Finally, we review the impact of epigenetic modifications on tumor immunity and the latest advancements in cancer therapies targeting these modifications.

Selective translation of nuclear mitochondrial respiratory proteins reprograms succinate metabolism in AML development and chemoresistance

Mitochondrial adaptations dynamically reprogram cellular bioenergetics and metabolism and confer key properties for human cancers. However, the selective regulation of these mitochondrial responses remains largely elusive. Here, inspired by a genetic screening in acute myeloid leukemia (AML), we identify RAS effector RREB1 as a translational regulator and uncover a unique translation control system for nuclear-encoded mitochondrial proteins in human cancers. RREB1 deletion reduces mitochondrial activities and succinate metabolism, thereby damaging leukemia stem cell (LSC) function and AML development. Replenishing complex II subunit SDHD rectifies these deficiencies. Notably, inhibition of complex II re-sensitizes AML cells to venetoclax treatment. Mechanistically, a short RREB1 variant binds to a conserved motif in the 3' UTRs and cooperates with elongation factor eEF1A1 to enhance protein translation of nuclear-encoded mitochondrial mRNAs. Overall, our findings reveal a unique translation control mechanism for mitochondrial adaptations in AML pathogenesis and provide a potential strategy for targeting this vulnerability of LSCs.

Landscape of alternative splicing and polyadenylation during growth and development of muscles in pigs

Alternative polyadenylation (APA) is emerging as a post-transcriptional regulatory mechanism, similar as that of alternative splicing (AS), and plays a prominent role in regulating gene expression and increasing the complexity of the transcriptome and proteome. We use polyadenylation selected long-read isoform sequencing to obtain full-length transcript sequences in porcine muscles at five developmental stages. We identify numerous novel transcripts unannotated in the existing pig genome, including transcripts mapping to known and unknown gene loci, and widespread transcript diversity in porcine muscles. The top 100 most isoformic genes are mainly enriched in Gene Ontology terms related to muscle growth and development. It is revealed that intron retention/exon inclusion and the usage of distal polyadenylation site (PAS) are associated with ageing through analyzing changes of AS and PAS during muscle development. We also identify developmental changes in major transcripts and major PASs. Furthermore, genes/transcripts important for muscle development are identified. The results confirm the importance of AS and APA in pig muscles, substantially increasing transcriptional diversity and showing an important mechanism underlying gene regulation in muscles.

Dissecting non-B DNA structural motifs in untranslated regions of eukaryotic genomes

The untranslated regions (UTRs) of genes significantly impact various biological processes, including transcription, posttranscriptional control, mRNA stability, localization, and translation efficiency. In functional areas of genomes, non-B DNA structures such as cruciform, curved, triplex, G-quadruplex, and Z-DNA structures are common and have an impact on cellular physiology. Although the role of these structures in cis-regulatory regions such as promoters is well established in eukaryotic genomes, their prevalence within UTRs across different eukaryotic classes has not been extensively documented. Our study investigated the prevalence of various non-B DNA motifs within the 5' and 3' UTRs across diverse eukaryotic species. Our comparative analysis encompassed the 5'-UTRs and 3'UTRs of 360 species representing diverse eukaryotic domains of life, including Arthropoda (Diptera, Hemiptera, and Hymenoptera), Chordata (Artiodactyla, Carnivora, Galliformes, Passeriformes, Primates, Rodentia, Squamata, Testudines), Magnoliophyta (Brassicales), Fabales (Poales), and Nematoda (Rhabditida), on the basis of datasets derived from the UTRdb. We observed that species belonging to taxonomic orders such as Rhabditida, Diptera, Brassicales, and Hemiptera present a prevalence of curved DNA motifs in their UTRs, whereas orders such as Testudines, Galliformes, and Rodentia present a preponderance of G-quadruplexes in both UTRs. The distribution of motifs is conserved across different taxonomic classes, although species-specific variations in motif preferences were also observed. Our research unequivocally illuminates the prevalence and potential functional implications of non-B DNA motifs, offering invaluable insights into the evolutionary and biological significance of these structures.

Generation of Transcript Length Variants and Reprogramming of mRNA Splicing During Atherosclerosis Progression in ApoE-Deficient Mice

Background. Variant 3'UTRs provide mRNAs with different binding sites for miRNAs or RNA-binding proteins (RBPs) allowing the establishment of new regulatory environments. Regulation of 3'UTR length impacts on the control of gene expression by regulating accessibility of miRNAs or RBPs to homologous sequences in mRNAs. Objective. Studying the dynamics of mRNA length variations in atherosclerosis (ATS) progression and reversion in ApoE-deficient mice exposed to a high-fat diet and treated with an αCD40-specific siRNA or with a sequence-scrambled siRNA as control. Methods. We gathered microarray mRNA expression data from the aortas of mice after 2 or 16 weeks of treatments, and used these data in a Bioinformatics analysis. Results. Here, we report the lengthening of the 5'UTR/3'UTRs and the shortening of the CDS in downregulated mRNAs during ATS progression. Furthermore, treatment with the αCD40-specific siRNA resulted in the partial reversion of the 3'UTR lengthening. Exon analysis showed that these length variations were actually due to changes in the number of exons embedded in mRNAs, and the further examination of transcripts co-expressed at weeks 2 and 16 in mice treated with the control siRNA revealed a process of mRNA isoform switching in which transcript variants differed in the patterns of alternative splicing or activated latent/cryptic splice sites. Conclusion. We document length variations in the 5'UTR/3'UTR and CDS of mRNAs downregulated during atherosclerosis progression and suggest a role for mRNA splicing reprogramming and transcript isoform switching in the generation of disease-related mRNA sequence diversity and variability.

APAtizer: a tool for alternative polyadenylation analysis of RNA-Seq data

SUMMARY

APAtizer is a tool designed to analyze alternative polyadenylation events on RNA-sequencing data. The tool handles different file formats, including BAM, htseq, and DaPars bedGraph files. It provides a user-friendly interface that allows users to generate informative visualizations, including Volcano plots, heatmaps, and gene lists. These outputs allow the user to retrieve useful biological insights such as the occurrence of polyadenylation events when comparing two biological conditions. In addition, it can perform differential gene expression, gene ontology analysis, visualization of Venn diagram intersections, and correlation analysis.

AVAILABILITY AND IMPLEMENTATION

Source code and example case studies for APAtizer are available at https://github.com/GeneRegulationi3S/APAtizer/.

Decoding biology with massively parallel reporter assays and machine learning

Massively parallel reporter assays (MPRAs) are powerful tools for quantifying the impacts of sequence variation on gene expression. Reading out molecular phenotypes with sequencing enables interrogating the impact of sequence variation beyond genome scale. Machine learning models integrate and codify information learned from MPRAs and enable generalization by predicting sequences outside the training data set. Models can provide a quantitative understanding of cis-regulatory codes controlling gene expression, enable variant stratification, and guide the design of synthetic regulatory elements for applications from synthetic biology to mRNA and gene therapy. This review focuses on cis-regulatory MPRAs, particularly those that interrogate cotranscriptional and post-transcriptional processes: alternative splicing, cleavage and polyadenylation, translation, and mRNA decay.

A genome-wide association study reveals molecular mechanism underlying powdery mildew resistance in cucumber

BACKGROUND

Powdery mildew is a disease with one of the most substantial impacts on cucumber production globally. The most efficient approach for controlling powdery mildew is the development of genetic resistance; however, few genes associated with inherent variations in cucumber powdery mildew resistance have been identified as of yet.

RESULTS

In this study, we re-sequence 299 cucumber accessions, which are divided into four geographical groups. A genome-wide association study identifies 50 sites significantly associated with natural variations in powdery mildew resistance. Linkage disequilibrium analysis further divides these 50 sites into 32 linkage disequilibrium blocks containing 41 putative genes. Virus-induced gene silencing and gene expression analysis implicate CsGy5G015960, which encodes a phosphate transporter, as the candidate gene regulating powdery mildew resistance. On the basis of the resequencing data, we generate five CsGy5G015960 haplotypes, identifying Hap.1 as the haplotype most likely associated with powdery mildew resistance. In addition, we determine that a 29-bp InDel in the 3' untranslated region of CsGy5G015960 is responsible for mRNA stability. Overexpression of CsGy5G015960Hap.1 in the susceptible line enhances powdery mildew resistance and phosphorus accumulation. Further comparative RNA-seq analysis demonstrates that CsGy5G015960Hap.1 may regulate cucumber powdery mildew resistance by maintaining a higher H2O2 level through the depletion of multiple class III peroxidases.

CONCLUSIONS

Here we identify a candidate powdery mildew-resistant gene in cucumber using GWAS. The identified gene may be a promising target for molecular breeding and genetic engineering in cucumber to enhance powdery mildew resistance.

Determination of the Frequency of BCL-2 Polymorphisms (c.-717C>A and c.*2364G>A) and LIF Polymorphism (c.*1414T>G) in Patients with Congenital Anomalies of the Kidney and Urinary Tract

Introduction

Congenital anomalies of the kidney and urinary tract (CAKUT) are characterized by several malformations. Its prevalence is 0.3-0.6% in live births. The B-cell lymphoma (BCL-2) gene regulates apoptosis, and the Leukemia Inhibitory Factor (LIF) gene plays a role in many biological processes, such as blastocyst growth and uterine preparation for implantation. In this study, two single nucleotide polymorphisms (SNPs) of the BCL-2 gene (rs2279115 and rs4987856) and one SNP of the LIF gene (rs929271) were investigated in CAKUT patients for the first time.

Methods

Hundred and twenty-nine CAKUT patients and 105 controls were enrolled in this study. We used polymerase chain reaction-restriction fragment length polymorphism for rs2279115 and rs929271 and SNaPshot for rs4987856. The χ2 test was used to compare discrete variables, and the independent sample t test was used to compare continuous variables.

Results

The allele frequencies for the rs2279115 and rs4987856 polymorphisms of BCL-2 and the rs929271 polymorphism of LIF were not significantly different between the patient and control groups (p = 0.162, p = 0.053, p = 0.635, respectively). However, the co-segregation analysis revealed a significant difference in the distribution of allele frequencies between the patient and control groups for two genetic variations: LIF rs929271 SNP and BCL-2 rs4987856 SNP (p = 0.034). The relative odds ratio was 2.444 (95% Confidence Interval (CI) 1.054-5.671).

Conclusion

This study, which is the first time in the literature, showed that changes in BCL-2 and LIF genes are associated with CAKUT disease.

An atlas of RNA-dependent proteins in cell division reveals the riboregulation of mitotic protein-protein interactions

Ribonucleoprotein complexes are dynamic assemblies of RNA with RNA-binding proteins (RBPs), which can modulate the fate of the RNA molecules from transcription to degradation. Vice versa, RNA can regulate the interactions and functions of the associated proteins. Dysregulation of RBPs is linked to diseases such as cancer and neurological disorders. RNA and RBPs are present in mitotic structures like the centrosomes and spindle microtubules, but their influence on mitotic spindle integrity remains unknown. Thus, we applied the R-DeeP strategy for the proteome-wide identification of RNA-dependent proteins and complexes to cells synchronized in mitosis versus interphase. The resulting atlas of RNA-dependent proteins in cell division can be accessed through the R-DeeP 3.0 database (R-DeeP3.dkfz.de). It revealed key mitotic factors as RNA-dependent such as AURKA, KIFC1 and TPX2 that were linked to RNA despite their lack of canonical RNA-binding domains. KIFC1 was identified as a new interaction partner and phosphorylation substrate of AURKA at S349 and T359. In addition, KIFC1 interacted with both, AURKA and TPX2, in an RNA-dependent manner. Our data suggest a riboregulation of mitotic protein-protein interactions during spindle assembly, offering new perspectives on the control of cell division processes by RNA-protein complexes.

The Role of Ryanodine Receptor 2 Polymorphisms in Oral Squamous Cell Carcinoma Susceptibility and Clinicopathological Features

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide, and oral squamous cell carcinoma (OSCC) is one of the most common types. There is strong evidence that ryanodine receptor 2 (RYR2) plays an important role in different types of cancer according to previous studies. Its expression is associated with survival in patients with HNSCC, but it is unknown whether altered RYR2 expression contributes to tumorigenesis. Therefore, we examined how RYR2 polymorphisms affect OSCC susceptibility and clinicopathological characteristics. Five single nucleotide polymorphisms (SNPs) of RYR2, rs12594, rs16835904, rs2779359, rs3765097, and rs3820216, were analyzed in 562 cases of OSCC and 332 healthy controls using real-time PCR. We demonstrated that RYR2 SNP rs12594 was significantly different between the case and control groups, but this difference was not significant after adjusting for personal habits. In contrast, we found that different genotypes of SNP rs2779359 were significantly associated with the characteristics of clinical stage and tumor size in OSCC patients, according to the odds ratios and the adjusted odds ratios; specifically, patients with the T genotype had 1.477-fold (95% CI, 1.043 to 2.091; p = 0.028) and 1.533-fold (95% CI, 1.087-2.162; p = 0.015) increases in clinical stage and tumor size, respectively, compared with patients with the C allele. The results of our study, in which RYR2 SNPs associated with OSCC progression and development were examined for the first time, suggest that clinicopathological characteristics may alter OSCC susceptibility. Finally, RYR2 SNP rs2779359 not only plays a role in both the prognosis and diagnosis of oral cancer but is also likely an important predictive factor for recurrence, response to treatment, and medication toxicity.

HybriSeq: Probe-based Device-free Single-cell RNA Profiling

We have developed the HybriSeq method for single-cell RNA profiling, which utilizes in situ hybridization of multiple probes for targeted transcripts, followed by split-pool barcoding and sequencing analysis of the probes. We have shown that HybriSeq can achieve high sensitivity for RNA detection with multiple probes and profile differential splicing. The utility of HybriSeq is demonstrated in characterizing cell-to-cell heterogeneities of a panel of 95 cell-cycle-related genes and the detection of misannotated transcripts.

Variation within the non-coding genome influences genetic and epigenetic regulation of the human leukocyte antigen genes

Variation within the non-coding genome may influence the regulation and expression of important genes involved in immune control such as the human leukocyte antigen (HLA) system. Class I and Class II HLA molecules are essential for peptide presentation which is required for T lymphocyte activation. Single nucleotide polymorphisms within non-coding regions of HLA Class I and Class II genes may influence the expression of these genes by affecting the binding of transcription factors and chromatin modeling molecules. Furthermore, an interplay between genetic and epigenetic factors may also influence HLA expression. Epigenetic factors such as DNA methylation and non-coding RNA, regulate gene expression without changing the DNA sequence. However, genetic variation may promote or allow genes to escape regulation by epigenetic factors, resulting in altered expression. The HLA system is central to most diseases, therefore, understanding the role of genetics and epigenetics on HLA regulation will tremendously impact healthcare. The knowledge gained from these studies may lead to novel and cost-effective diagnostic approaches and therapeutic interventions. This review discusses the role of non-coding variants on HLA regulation. Furthermore, we discuss the interplay between genetic and epigenetic factors on the regulation of HLA by evaluating literature based on polymorphisms within DNA methylation and miRNA regulatory sites within class I and Class II HLA genes. We also provide insight into the importance of the HLA non-coding genome on disease, discuss ethnic-specific differences across the HLA region and provide guidelines for future HLA studies.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

In-locus gene silencing in plants using genome editing

Gene silencing is crucial in crop breeding for desired trait development. RNA interference (RNAi) has been used widely but is limited by ectopic expression of transgenes and genetic instability. Introducing an upstream start codon (uATG) into the 5'untranslated region (5'UTR) of a target gene may 'silence' the target gene by inhibiting protein translation from the primary start codon (pATG). Here, we report an efficient gene silencing method by introducing a tailor-designed uATG-containing element (ATGE) into the 5'UTR of genes in plants, occupying the original start site to act as a new pATG. Using base editing to introduce new uATGs failed to silence two of the tested three rice genes, indicating complex regulatory mechanisms. Precisely inserting an ATGE adjacent to pATG achieved significant target protein downregulation. Through extensive optimization, we demonstrated this strategy substantially and consistently downregulated target protein expression. By designing a bidirectional multifunctional ATGE4, we enabled tunable knockdown from 19% to 89% and observed expected phenotypes. Introducing ATGE into Waxy, which regulates starch synthesis, generated grains with lower amylose, revealing the value for crop breeding. Together, we have developed a programmable and robust method to knock down gene expression in plants, with potential for biological mechanism exploration and crop enhancement.

Dynamic regulation of alternative polyadenylation by PQBP1 during neurogenesis

Alternative polyadenylation (APA) is a critical post-transcriptional process that generates mRNA isoforms with distinct 3' untranslated regions (3' UTRs), thereby regulating mRNA localization, stability, and translational efficiency. Cell-type-specific APA extensively shapes the diversity of the cellular transcriptome, particularly during cell fate transition. Despite its recognized significance, the precise regulatory mechanisms governing cell-type-specific APA remain unclear. In this study, we uncover PQBP1 as an emerging APA regulator that actively maintains cell-specific APA profiles in neural progenitor cells (NPCs) and delicately manages the equilibrium between NPC proliferation and differentiation. Multi-omics analysis shows that PQBP1 directly interacts with the upstream UGUA elements, impeding the recruitment of the CFIm complex and influencing polyadenylation site selection within genes associated with the cell cycle. Our findings elucidate the molecular mechanism by which PQBP1 orchestrates dynamic APA changes during neurogenesis, providing valuable insights into the precise regulation of cell-type-specific APA and the underlying pathogenic mechanisms in neurodevelopmental disorders.

Hominini-specific regulation of the cell cycle by stop codon readthrough of FEM1B

FEM1B is a substrate-recognition component of the CRL2 E3 ubiquitin-protein ligase. This multi-protein complex targets specific proteins for ubiquitylation, which leads to their degradation. Here, we demonstrate the regulation of FEM1B expression by stop codon readthrough (SCR). In this process, translating ribosomes readthrough the stop codon of FEM1B to generate a C-terminally extended isoform that is highly unstable. A total of 81 nucleotides in the proximal 3'UTR of FEM1B constitute the necessary and sufficient cis-signal for SCR. Also, they encode the amino acid sequence responsible for the degradation of the SCR product. CRISPR-edited cells lacking this region, and therefore SCR of FEM1B, showed increased FEM1B expression. This in turn resulted in reduced expression of SLBP (a target of FEM1B-mediated degradation) and replication-dependent histones (target of SLBP for mRNA stability), causing cell cycle delay. Evolutionary analysis revealed that this phenomenon is specific to the genus Pan and Homo (Hominini). Overall, we show a relatively recently evolved SCR process that relieves the cell cycle from the negative regulation by FEM1B.

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

Quail GHRL and LEAP2 gene cloning, polymorphism detection, phylogenetic analysis, tissue expression profiling and its association analysis with feed intake

The GHRL, LEAP2, and GHSR system have recently been identified as important regulators of feed intake in mammals and chickens. However, the complete cloning of the quail GHRL (qGHRL) and quail LEAP2 (qLEAP2) genes, as well as their association with feed intake, remains unclear. This study cloned the entire qGHRL and qLEAP2 cDNA sequence in Chinese yellow quail (Coturnix japonica), including the 5' and 3' untranslated regions. Sanger sequencing analysis revealed no missense mutations in the coding region of qGHRL and qLEAP2. Subsequently, phylogenetic analysis and protein homology alignment were conducted on the qGHRL and qLEAP2 in major poultry species. The findings of this research indicated that the qGHRL and qLEAP2 sequences exhibit a high degree of similarity with those of chicken and turkey. Specifically, the N-terminal 6 amino acids of GHRL mature peptides and all the mature peptide sequence of LEAP2 exhibited consistent patterns across all species examined. The analysis of tissue gene expression profiles indicated that qGHRL was primarily expressed in the proventriculus and brain tissue, whereas qLEAP2 exhibited higher expression levels in the intestinal tissue, kidney, and liver tissue, differing slightly from previous studies conducted on chicken. It is necessary to investigate the significance of elevated expression of qGHRL in brain and qLEAP2 in kidney in the future. Further research has shown that the expression of qLEAP2 can quickly respond to changes in different energy states, whereas qGHRL does not exhibit the same capability. Overall, this study successfully cloned the complete cDNA sequences of qGHRL and qLEAP2, and conducted a comprehensive examination of their tissue expression profiles and gene expression levels in the main expressing organs across different energy states. Our current findings suggested that qLEAP2 is highly expressed in the liver, intestine, and kidney, and its expression level is regulated by feed intake.

Respiratory Syncytial Virus (RSV) optimizes the translational landscape during infection

Viral infection often triggers eukaryotic initiator factor 2α (eIF2α) phosphorylation, leading to global 5'-cap-dependent translation inhibition. RSV encodes messenger RNAs (mRNAs) mimicking 5'-cap structures of host mRNAs and thus inhibition of cap-dependent translation initiation would likely also reduce viral translation. We confirmed that RSV limits widespread translation initiation inhibition and unexpectedly found that the fraction of ribosomes within polysomes increases during infection, indicating higher ribosome loading on mRNAs during infection. We found that AU-rich host transcripts that are less efficiently translated under normal conditions become more efficient at recruiting ribosomes, similar to RSV transcripts. Viral transcripts are transcribed in cytoplasmic inclusion bodies, where the viral AU-rich binding protein M2-1 has been shown to bind viral transcripts and shuttle them into the cytoplasm. We further demonstrated that M2-1 is found on polysomes, and that M2-1 might deliver host AU-rich transcripts for translation.

How antisense transcripts can evolve to encode novel proteins

Protein coding features can emerge de novo in non coding transcripts, resulting in emergence of new protein coding genes. Studies across many species show that a large fraction of evolutionarily novel non-coding RNAs have an antisense overlap with protein coding genes. The open reading frames (ORFs) in these antisense RNAs could also overlap with existing ORFs. In this study, we investigate how the evolution an ORF could be constrained by its overlap with an existing ORF in three different reading frames. Using a combination of mathematical modeling and genome/transcriptome data analysis in two different model organisms, we show that antisense overlap can increase the likelihood of ORF emergence and reduce the likelihood of ORF loss, especially in one of the three reading frames. In addition to rationalising the repeatedly reported prevalence of de novo emerged genes in antisense transcripts, our work also provides a generic modeling and an analytical framework that can be used to understand evolution of antisense genes.

Untranslated yet indispensable-UTRs act as key regulators in the environmental control of gene expression

To survive and thrive in a dynamic environment, plants must continuously monitor their surroundings and adjust their development and physiology accordingly. Changes in gene expression underlie these developmental and physiological adjustments, and are traditionally attributed to widespread transcriptional reprogramming. Growing evidence, however, suggests that post-transcriptional mechanisms also play a vital role in tailoring gene expression to a plant's environment. Untranslated regions (UTRs) act as regulatory hubs for post-transcriptional control, harbouring cis-elements that affect an mRNA's processing, localization, translation, and stability, and thereby tune the abundance of the encoded protein. Here, we review recent advances made in understanding the critical function UTRs exert in the post-transcriptional control of gene expression in the context of a plant's abiotic environment. We summarize the molecular mechanisms at play, present examples of UTR-controlled signalling cascades, and discuss the potential that resides within UTRs to render plants more resilient to a changing climate.

Multilevel gene expression changes in lineages containing adaptive copy number variants

Copy-number variants (CNVs) are an important class of recurrent variants that mediate adaptive evolution. While CNVs can increase the relative fitness of the organism, they can also incur a cost. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition, and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- conditions. We find CNV-amplified genes correlate with higher RNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of expression by the abundance of the preceding step we were able to identify widespread divergence in the efficiency of each step in the gene in the efficiency of each step in gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA binding sites for SSD1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which adaptive evolution shapes gene expression.