Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

Integrating apaQTL and eQTL analysis identifies a potential causal variant associated with lung adenocarcinoma risk in the Chinese population

Alternative polyadenylation (APA) plays a crucial role in cancer biology. Here, we used data from the 3'aQTL-atlas, GTEx, and the China Nanjing Lung Cancer GWAS database to explore the association between apaQTL/eQTL-SNPs and the risk of lung adenocarcinoma (LUAD). The variant T allele of rs277646 in NIT2 is associated with an increased risk of LUAD (OR = 1.12, P = 0.015), lower PDUI values, and higher NIT2 expression. The 3'RACE experiment showed multiple poly (A) sites in NIT2, with the rs277646-T allele causing preferential use of the proximal poly (A) site, resulting in a shorter 3'UTR transcript. This leads to the loss of the hsa-miR-650 binding site, thereby affecting LUAD malignant phenotypes by regulating the expression level of NIT2. Our findings may provide new insights into understanding and exploring APA events in LUAD carcinogenesis.

Dynamic Alternative Polyadenylation during Litopenaeus Vannamei Metamorphosis Development

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

Utilizing non-coding RNA-mediated regulation of ATP binding cassette (ABC) transporters to overcome multidrug resistance to cancer chemotherapy

Multidrug resistance (MDR) is one of the primary factors that produces treatment failure in patients receiving cancer chemotherapy. MDR is a complex multifactorial phenomenon, characterized by a decrease or abrogation of the efficacy of a wide spectrum of anticancer drugs that are structurally and mechanistically distinct. The overexpression of the ATP-binding cassette (ABC) transporters, notably ABCG2 and ABCB1, are one of the primary mediators of MDR in cancer cells, which promotes the efflux of certain chemotherapeutic drugs from cancer cells, thereby decreasing or abolishing their therapeutic efficacy. A number of studies have suggested that non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), play a pivotal role in mediating the upregulation of ABC transporters in certain MDR cancer cells. This review will provide updated information about the induction of ABC transporters due to the aberrant regulation of ncRNAs in cancer cells. We will also discuss the measurement and biological profile of circulating ncRNAs in various body fluids as potential biomarkers for predicting the response of cancer patients to chemotherapy. Sequence variations, such as alternative polyadenylation of mRNA and single nucleotide polymorphism (SNPs) at miRNA target sites, which may indicate the interaction of miRNA-mediated gene regulation with genetic variations to modulate the MDR phenotype, will be reviewed. Finally, we will highlight novel strategies that could be used to modulate ncRNAs and circumvent ABC transporter-mediated MDR.

The miRNA-target interactions: An underestimated intricacy

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Hidden features of NAD-RNA epitranscriptome in Drosophila life cycle

Nicotinamide adenine dinucleotide (NAD), a nucleotide-containing metabolite, can be incorporated into the RNA 5'-terminus to result in NAD-capped RNA (NAD-RNA). Since NAD has been heightened as one of the most essential metabolites in cells, its linkage to RNA represents a critical but poorly studied modification at the epitranscriptomic level. Here, we design a highly sensitive method, DO-seq, to capture NAD-RNAs. Using Drosophila, we identify thousands of previously unexplored NAD-RNAs and their dynamics in the fly life cycle, from embryo to adult. We show the evidence that chromosomal clustering might be the structural basis by which co-expression can couple with NAD capping on physically and functionally linked genes. Furthermore, we note that NAD capping of cuticle genes inversely correlates with their gene expression. Combined, we propose NAD-RNA epitranscriptome as a hidden layer of regulation that underlies biological processes. DO-seq empowers the identification of NAD-capped RNAs, facilitating functional investigation into this modification.

Non-coding RNAs: Emerging biomarkers and therapeutic targets in ulcerative colitis

Ulcerative colitis (UC) is a persistent inflammatory condition affecting the colon's mucosal lining, leading to chronic bowel inflammation. Despite extensive research, the precise molecular mechanisms underlying UC pathogenesis remain elusive. NcRNAs form a category of functional RNA molecules devoid of protein-coding capacity. They have recently surfaced as pivotal modulators of gene expression and integral participants in various pathological processes, particularly those related to inflammatory disorders. The diverse classes of ncRNAs, encompassing miRNAs, circRNAs, and lncRNAs, have been implicated in UC. It highlights their involvement in key UC-related processes, such as immune cell activation, epithelial barrier integrity, and the production of pro-inflammatory mediators. ncRNAs have been identified as potential biomarkers for UC diagnosis and monitoring disease progression, offering promising avenues for personalized medicine. This approach may pave the way for novel, more specific treatments with reduced side effects, addressing the current limitations of conventional therapies. A comprehensive understanding of the interplay between ncRNAs and UC will advance our knowledge of the disease, potentially leading to more effective and personalized treatments for patients suffering from this debilitating condition. This review explores the pivotal role of ncRNAs in the context of UC, shedding light on their possible targets for diagnosis, prognosis, and therapeutic interventions.

Validation of an Integrated Genetic-Epigenetic Test for the Assessment of Coronary Heart Disease

BACKGROUND

Coronary heart disease (CHD) is the leading cause of death in the world. Unfortunately, many of the key diagnostic tools for CHD are insensitive, invasive, and costly; require significant specialized infrastructure investments; and do not provide information to guide postdiagnosis therapy. In prior work using data from the Framingham Heart Study, we provided in silico evidence that integrated genetic-epigenetic tools may provide a new avenue for assessing CHD.

METHODS AND RESULTS

In this communication, we use an improved machine learning approach and data from 2 additional cohorts, totaling 449 cases and 2067 controls, to develop a better model for ascertaining symptomatic CHD. Using the DNA from the 2 new cohorts, we translate and validate the in silico findings into an artificial intelligence-guided, clinically implementable method that uses input from 6 methylation-sensitive digital polymerase chain reaction and 10 genotyping assays. Using this method, the overall average area under the curve, sensitivity, and specificity in the 3 test cohorts is 82%, 79%, and 76%, respectively. Analysis of targeted cytosine-phospho-guanine loci shows that they map to key risk pathways involved in atherosclerosis that suggest specific therapeutic approaches.

CONCLUSIONS

We conclude that this scalable integrated genetic-epigenetic approach is useful for the diagnosis of symptomatic CHD, performs favorably as compared with many existing methods, and may provide personalized insight to CHD therapy. Furthermore, given the dynamic nature of DNA methylation and the ease of methylation-sensitive digital polymerase chain reaction methodologies, these findings may pave a pathway for precision epigenetic approaches for monitoring CHD treatment response.

Functional Consequences of Shifting Transcript Boundaries in Glucose Starvation

Glucose is a major source of carbon and essential for the survival of many organisms, ranging from yeast to human. A sudden 60-fold reduction of glucose in exponentially growing fission yeast induces transcriptome-wide changes in gene expression. This regulation is multilayered, and the boundaries of transcripts are known to vary, with functional consequences at the protein level. By combining direct RNA sequencing with 5'-CAGE and short-read sequencing, we accurately defined the 5'- and 3'-ends of transcripts that are both poly(A) tailed and 5'-capped in glucose starvation, followed by proteome analysis. Our results confirm previous experimentally validated loci with alternative isoforms and reveal several transcriptome-wide patterns. First, we show that sense-antisense gene pairs are more strongly anticorrelated when a time lag is taken into account. Second, we show that the glucose starvation response initially elicits a shortening of 3'-UTRs and poly(A) tails, followed by a shortening of the 5'-UTRs at later time points. These result in domain gains and losses in proteins involved in the stress response. Finally, the relatively poor overlap both between differentially expressed genes (DEGs), differential transcript usage events (DTUs), and differentially detected proteins (DDPs) highlight the need for further study on post-transcriptional regulation mechanisms in glucose starvation.

Non-coding RNAs in gynecologic cancer

The term "gynecologic cancer" pertains to neoplasms impacting the reproductive tissues and organs of women encompassing the endometrium, vagina, cervix, uterus, vulva, and ovaries. The progression of gynecologic cancer is linked to various molecular mechanisms. Historically, cancer research primarily focused on protein-coding genes. However, recent years have unveiled the involvement of non-coding RNAs (ncRNAs), including microRNAs, long non-coding RNAs (LncRNAs), and circular RNAs, in modulating cellular functions within gynecological cancer. Substantial evidence suggests that ncRNAs may wield a dual role in gynecological cancer, acting as either oncogenic or tumor-suppressive agents. Numerous clinical trials are presently investigating the roles of ncRNAs as biomarkers and therapeutic agents. These endeavors may introduce a fresh perspective on the diagnosis and treatment of gynecological cancer. In this overview, we highlight some of the ncRNAs associated with gynecological cancers.

Effect of mRNA-LNP components of two globally-marketed COVID-19 vaccines on efficacy and stability

During the COVID-19 pandemic, Pfizer-BioNTech and Moderna successfully developed nucleoside-modified mRNA lipid nanoparticle (LNP) vaccines. SARS-CoV-2 spike protein expressed by those vaccines are identical in amino acid sequence, but several key components are distinct. Here, we compared the effect of ionizable lipids, untranslated regions (UTRs), and nucleotide composition of the two vaccines, focusing on mRNA delivery, antibody generation, and long-term stability. We found that the ionizable lipid, SM-102, in Moderna's vaccine performs better than ALC-0315 in Pfizer-BioNTech's vaccine for intramuscular delivery of mRNA and antibody production in mice and long-term stability at 4 °C. Moreover, Pfizer-BioNTech's 5' UTR and Moderna's 3' UTR outperform their counterparts in their contribution to transgene expression in mice. We further found that varying N1-methylpseudouridine content at the wobble position of mRNA has little effect on vaccine efficacy. These findings may contribute to the further improvement of nucleoside-modified mRNA-LNP vaccines and therapeutics.

RNA-RNA competitive interactions: a molecular civil war ruling cell physiology and diseases

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

Amplifying gene expression with RNA-targeted therapeutics

Many diseases are caused by insufficient expression of mutated genes and would benefit from increased expression of the corresponding protein. However, in drug development, it has been historically easier to develop drugs with inhibitory or antagonistic effects. Protein replacement and gene therapy can achieve the goal of increased protein expression but have limitations. Recent discoveries of the extensive regulatory networks formed by non-coding RNAs offer alternative targets and strategies to amplify the production of a specific protein. In addition to RNA-targeting small molecules, new nucleic acid-based therapeutic modalities that allow highly specific modulation of RNA-based regulatory networks are being developed. Such approaches can directly target the stability of mRNAs or modulate non-coding RNA-mediated regulation of transcription and translation. This Review highlights emerging RNA-targeted therapeutics for gene activation, focusing on opportunities and challenges for translation to the clinic.

Non-coding regions of nuclear-DNA-encoded mitochondrial genes and intergenic sequences are targeted by autoantibodies in breast cancer

Autoantibodies against mitochondrial-derived antigens play a key role in chronic tissue inflammation in autoimmune disorders and cancers. Here, we identify autoreactive nuclear genomic DNA (nDNA)-encoded mitochondrial gene products (GAPDH, PKM2, GSTP1, SPATA5, MFF, TSPOAP1, PHB2, COA4, and HAGH) recognized by breast cancer (BC) patients’ sera as nonself, supporting a direct relationship of mitochondrial autoimmunity to breast carcinogenesis. Autoreactivity of multiple nDNA-encoded mitochondrial gene products was mapped to protein-coding regions, 3’ untranslated regions (UTRs), as well as introns. In addition, autoantibodies in BC sera targeted intergenic sequences that may be parts of long non-coding RNA (lncRNA) genes, including LINC02381 and other putative lncRNA neighbors of the protein-coding genes ERCC4, CXCL13, SOX3, PCDH1, EDDM3B, and GRB2. Increasing evidence indicates that lncRNAs play a key role in carcinogenesis. Consistent with this, our findings suggest that lncRNAs, as well as mRNAs of nDNA-encoded mitochondrial genes, mechanistically contribute to BC progression. This work supports a new paradigm of breast carcinogenesis based on a globally dysfunctional genome with altered function of multiple mitochondrial and non-mitochondrial oncogenic pathways caused by the effects of autoreactivity-induced dysregulation of multiple genes and their products. This autoimmunity-based model of carcinogenesis will open novel avenues for BC treatment.

Preimplantation embryo gene expression: 56 years of discovery, and counting

The biology of preimplantation embryo gene expression began 56 years ago with studies of the effects of protein synthesis inhibition and discovery of changes in embryo metabolism and related enzyme activities. The field accelerated rapidly with the emergence of embryo culture systems and progressively evolving methodologies that have allowed early questions to be re-addressed in new ways and in greater detail, leading to deeper understanding and progressively more targeted studies to discover ever more fine details. The advent of technologies for assisted reproduction, preimplantation genetic testing, stem cell manipulations, artificial gametes, and genetic manipulation, particularly in experimental animal models and livestock species, has further elevated the desire to understand preimplantation development in greater detail. The questions that drove enquiry from the earliest years of the field remain drivers of enquiry today. Our understanding of the crucial roles of oocyte-expressed RNA and proteins in early embryos, temporal patterns of embryonic gene expression, and mechanisms controlling embryonic gene expression has increased exponentially over the past five and a half decades as new analytical methods emerged. This review combines early and recent discoveries on gene regulation and expression in mature oocytes and preimplantation stage embryos to provide a comprehensive understanding of preimplantation embryo biology and to anticipate exciting future advances that will build upon and extend what has been discovered so far.

The Functional Meaning of 5'UTR in Protein-Coding Genes

As it is well known, messenger RNA has many regulatory regions along its sequence length. One of them is the 5' untranslated region (5'UTR), which itself contains many regulatory elements such as upstream ORFs (uORFs), internal ribosome entry sites (IRESs), microRNA binding sites, and structural components involved in the regulation of mRNA stability, pre-mRNA splicing, and translation initiation. Activation of the alternative, more upstream transcription start site leads to an extension of 5'UTR. One of the consequences of 5'UTRs extension may be head-to-head gene overlap. This review describes elements in 5'UTR of protein-coding transcripts and the functional significance of protein-coding genes 5' overlap with implications for transcription, translation, and disease.

3'UTR Diversity: Expanding Repertoire of RNA Alterations in Human mRNAs

Genomic information stored in the DNA is transcribed to the mRNA and translated to proteins. The 3' untranslated regions (3'UTRs) of the mRNA serve pivotal roles in posttranscriptional gene expression, regulating mRNA stability, translation, and localization. Similar to DNA mutations producing aberrant proteins, RNA alterations expand the transcriptome landscape and change the cellular proteome. Recent global analyses reveal that many genes express various forms of altered RNAs, including 3'UTR length variants. Alternative polyadenylation and alternative splicing are involved in diversifying 3'UTRs, which could act as a hidden layer of eukaryotic gene expression control. In this review, we summarize the functions and regulations of 3'UTRs and elaborate on the generation and functional consequences of 3'UTR diversity. Given that dynamic 3'UTR length control contributes to phenotypic complexity, dysregulated 3'UTR diversity might be relevant to disease development, including cancers. Thus, 3'UTR diversity in cancer could open exciting new research areas and provide avenues for novel cancer theragnostics.

A common whole-genome paleotetraploidization in Cucurbitales

Cucurbitales are an important order of flowering plants known for encompassing edible plants of economic and medicinal value and numerous ornamental plants of horticultural value. By reanalyzing the genomes of two representative families (Cucurbitaceae and Begoniaceae) in Cucurbitales, we found that the previously identified Cucurbitaceae common paleotetraploidization that occurred shortly after the core-eudicot-common hexaploidization event is shared by Cucurbitales, including Begoniaceae. We built a multigenome alignment framework for Cucurbitales by identifying orthologs and paralogs and systematically redating key evolutionary events in Cucurbitales. Notably, characterizing the gene retention levels and genomic fractionation patterns between subgenomes generated from different polyploidizations in Cucurbitales suggested the autopolyploid nature of the Begoniaceae common tetraploidization and the allopolyploid nature of the Cucurbitales common tetraploidization and the Cucurbita-specific tetraploidization. Moreover, we constructed the ancestral Cucurbitales karyotype comprising 17 proto-chromosomes, confirming that the most recent common ancestor of Cucurbitaceae contained 15 proto-chromosomes and rejecting the previous hypothesis for an ancestral Cucurbitaceae karyotype with 12 proto-chromosomes. In addition, we found that the polyploidization and tandem duplication events promoted the expansion of gene families involved in the cucurbitacin biosynthesis pathway; however, gene loss and chromosomal rearrangements likely limited the expansion of these gene families.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Isolation and characterization of avian leukosis virus subgroup J associated with hemangioma and myelocytoma in layer chickens in China

A strain of avian leukosis virus (ALV) belonging to a new envelope subgroup J (ALV-J) emerged in 1988 as a new subgroup of ALV and spread rapidly throughout the world. Due to the infection and spread of ALV-J, the global poultry industry experienced a significant loss. Although the disease had been prevented and controlled effectively by culling domestic chickens in the infected zone, a few field cases of ALV-J infection were reported in China in recent years. This study was conducted to characterize the genome and analyze the lesions and histopathology of the ALV-J strain named HB2020, which was isolated from layer chickens in Hubei Province, China. The full-length proviral genome sequence analysis of ALV-J HB2020 revealed that it was a recombinant strain of ev-1 and HPRS-103 in the gag gene in comparison to ALV-J prototype HPRS-103. In the 3′-untranslated region (3'UTR) of the nucleotide sequence, there were found 205-base pairs (bp) deletion, of which 175 were detected in the redundant transmembrane (rTM) region. Besides, the surface glycoprotein gene gp85 had five mutations in a conservative site, whereas the transmembrane protein gene gp37 was relatively conserved. The animal experiments conducted later on this strain have shown that HB2020 can cause various neoplastic lesions in chickens, including enlarged livers with hemangiomas and spleens with white nodules. Additionally, as the exposure time increased, the number of tumor cells that resembled myelocytes in the blood smears of infected chickens gradually increased. These results indicated that HB2020 on recombination with ALV subgroup E (ALV-E) and ALV-J could induce severe hemangiomas and myelocytomas. This inference might provide a molecular basis for further research about the pathogenicity of ALV and emphasize the need for control and prevention of avian leukosis.

Molecular Mechanisms Underlying the Elevated Expression of a Potentially Type 2 Diabetes Mellitus Associated SCD1 Variant

Disturbances in lipid metabolism related to excessive food intake and sedentary lifestyle are among major risk of various metabolic disorders. Stearoyl-CoA desaturase-1 (SCD1) has an essential role in these diseases, as it catalyzes the synthesis of unsaturated fatty acids, both supplying for fat storage and contributing to cellular defense against saturated fatty acid toxicity. Recent studies show that increased activity or over-expression of SCD1 is one of the contributing factors for type 2 diabetes mellitus (T2DM). We aimed to investigate the impact of the common missense rs2234970 (M224L) polymorphism on SCD1 function in transfected cells. We found a higher expression of the minor Leu224 variant, which can be attributed to a combination of mRNA and protein stabilization. The latter was further enhanced by various fatty acids. The increased level of Leu224 variant resulted in an elevated unsaturated: saturated fatty acid ratio, due to higher oleate and palmitoleate contents. Accumulation of Leu224 variant was found in a T2DM patient group, however, the difference was statistically not significant. In conclusion, the minor variant of rs2234970 polymorphism might contribute to the development of obesity-related metabolic disorders, including T2DM, through an increased intracellular level of SCD1.