Post-transcriptional Regulation Research Articles

Expression profiles of circulating miRNAs in an endangered Piedmontese sheep breed during the estrus cycle

IntroductionThe preservation of locally endangered breeds is essential for maintaining ecosystem services that benefit both society and the environment. Reproductive fitness becomes a crucial consideration in this context. MicroRNAs (miRNAs) are small non-coding RNA molecules that play a key role in post-transcriptional regulation. Typically, they function within the tissues where they are produced. However, when they are released into extracellular fluid, they are referred to as circulating miRNAs (c-miRNAs). C-miRNAs may serve as potential biomarkers, whose profile changes under different physiological states. The purpose of this study is to establish a connection between distinctive variations in the expression of c-miRNAs and specific estrus cycle phases in Frabosana-Roaschina sheep, an endangered Piedmontese breed.MethodsTwo trials, each involving 20 ewes with different reproductive efficiencies (nulliparous in the first trial and pluriparous in the second trial), were sampled on alternate days after synchronization for blood, saliva, and feces. Ultrasound scans were performed during the induced estrus cycle. The animals’ behaviors were assessed through video recordings.ResultsIn the first trial, play behaviors were detected without sexual behaviors, whereas in the second trial, sexual behaviors were observed without play behaviors. Based on plasma trends of 17β-estradiol and progesterone and ultrasound images, two moments were identified for miRNAs analyses: the beginning of the follicular phase (day 2) and the beginning of the luteal phase (day 11). C-miRNAs of six representative animals from the second trial were sequenced. Analyses of the sequencing data have identified 12 c-miRNAs that were differentially expressed (DE) when comparing day 11 with day 2: five miRNAs were found to be upregulated, whereas seven miRNAs were downregulated. An enrichment analysis, based on predicted targets, using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases was performed. Many of these genes regulate reproductive pathways with the possible involvement of miRNAs. Finally, qRT-PCR was conducted to validate the DE miRNAs in all ewes. Differences in gene expression between the two sampling points and the two trials were observed, in line with existing literature.DiscussionInvestigating the role of these miRNAs in regulating estrus could improve the reproductive performance and welfare of Frabosana-Roaschina ewes.

Molecular mechanisms of acquired idiopathic generalized anhidrosis: differential gene expression and potential therapeutic targets

Eccrine sweat glands are essential for thermoregulation. Acquired idiopathic generalized anhidrosis (AIGA) is a rare disorder characterized by the absence of sweating, leading to severe complications like heatstroke and skin dryness. The underlying etiology remains unknown, complicating diagnosis and treatment. This study aimed to elucidate the genetic mechanisms of AIGA by analysing the gene expression in eccrine sweat glands using the GSE193125 dataset from the NCBI Gene Expression Omnibus (GEO). This study identified differentially expressed genes (DEGs) (genes that show statistically significant differences in expression levels in AIGA), with JUN and CCN1 significantly downregulated in anhidrotic lesions of AIGA patients. Receiver operating characteristic (ROC) curve demonstrated the accuracy of these genes in classification (JUN = 0.88 and CCN1 = 0.77). Gene and protein interaction networks which are defined as networks of physical interactions between genes and proteins, constructed using GeneMANIA and STRING, revealed substantial physical interactions among them. Gene set enrichment analysis highlighted the neuroinflammatory and estrogen receptor signaling pathways are associated with this condition, suggesting hormonal regulation’s role in the anhidrosis mechanism. Key miRNAs, including hsa-let-7b-5p, hsa-miR-10b-5p, and hsa-miR-124-3p, are associated with the DEGs, underscoring the importance of post-transcriptional and transcriptional regulation. Through the Drug Gene Interaction Database (DGIdb), ciprofibrate is identified as a potential therapeutic agent targeting the JUN gene. Molecular docking confirmed a strong binding affinity between ciprofibrate and JUN, suggesting ciprofibrate could modulate JUN activity and regulate AIGA symptoms. ADMET analysis, which refers to the assessment of absorption, distribution, metabolism, excretion, and toxicity of a drug or compound, is crucial for evaluating its pharmacokinetic and safety profile. The analysis indicated that ciprofibrate has favourable pharmacokinetic properties and a relatively safe toxicity profile. This comprehensive bioinformatics approach, integrating gene expression analysis, miRNA identification, and drug interaction, provides valuable insights into anhidrosis molecular pathology and highlights ciprofibrate’s potential as a targeted therapeutic strategy. Despite limitations, including a small sample size and reliance on computational predictions, this study paves the way for further research and potential therapeutic interventions for AIGA.

The miR-31b targets arylsulfatase B to regulate the ovarian development of Bactrocera dorsalis.

Reproduction is the basis of insect population growth and evolution, and encompasses ovarian development, reproductive behavior, and fecundity. Bactrocera dorsalis is a globally significant agricultural pest that is subject to quarantine, with mated females that can lay over 3000 eggs. The post-transcriptional regulation of ovarian development remains unclear. Here, miR-31b is shown to be involved in regulating Bactrocera dorsalis ovarian development. CRISPR/Cas9 was used to generate miR-31b loss-of-function mutations in Bactrocera dorsalis. The removal of miR-31b resulted in severely impaired ovarian development in adults, with phenotypes that included dramatically reduced egg production and hatching rates. The relationship between miR-31b and its target gene arylsulfatase B (ARSB) was subsequently identified using the methods of bioinformatics, transcriptomic sequencing, quantitative polymerase chain reaction (qPCR), RNA pull-down and dual-luciferase reporter assay. Finally, miR-31b was confirmed to bind the target gene arylsulfatase B to affect metabolism and thereby further hindered ovarian development of Bactrocera dorsalis. Overall, these results provide new insights into molecular mechanisms at the post-transcriptional level in regulating ovarian development and insect reproduction, consequently providing potential targets to control arthropod pests through the reproductive strategy. © 2024 Society of Chemical Industry.

Alternative polyadenylation shapes the molecular and clinical features of lung adenocarcinoma.

Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation that affects mRNA stability, localization and translation efficiency. Previous pan-cancer studies have revealed that APA is frequently disrupted in cancer and is associated with patient outcomes. Yet, little is known about cancer type-specific APA alterations. Here, we integrated RNA-sequencing data from a Korean cohort (GEO: GSE40419) and The Cancer Genome Atlas (TCGA) to comprehensively analyze APA alterations in lung adenocarcinomas (LUADs). Comparing expression levels of core genes involved in polyadenylation, we find that overall, the set of 28 of 31 genes are upregulated, with CSTF2 particularly upregulated. We observed broad and recurrent APA changes in LUAD growth-promoting genes. In addition, we find enrichment of APA events in genes associated with known LUAD pathways and an increased heterogeneity in polyadenylation (polyA) site usage of proliferation-associated genes. Upon further investigation, we report smoking-specific APA changes are also highly relevant to LUAD development. Overall, our in-depth analysis reveals APA as an important driver for the molecular and clinical features of lung adenocarcinoma.

P4PC: A Portal for Bioinformatics Resources of piRNAs and circRNAs

Background: PIWI-interacting RNAs (piRNAs) and circular RNAs (circRNAs) are two kinds of non-coding RNAs (ncRNAs) that play important roles in epigenetic regulation, transcriptional regulation, post-transcriptional regulation of many biological processes. Although there exist various resources, it is still challenging to select such resources for specific research projects on ncRNAs. Methods: In order to facilitate researchers in finding the appropriate bioinformatics sources for studying ncRNAs, we created a novel portal named P4PC that provides computational tools and data sources of piRNAs and circRNAs. Results: 249 computational tools, 126 databases and 420 papers are manually curated in P4PC. All entries in P4PC are classified in 5 groups and 26 subgroups. The list of resources is summarized in the first page of each group. Conclusion: According to their research proposes, users can quickly select proper resources for their research projects by viewing detail information and comments in P4PC. Database URL is http://www.ibiomedical.net/Portal4PC/ and https://43.138.46.5/Portal4PC/.

A Thorough Navigation of miRNA's Blueprint in Crafting Cardiovascular Fate.

Cardiovascular diseases contribute significantly to global morbidity and mortality. MicroRNAs are crucial in the development and progression of these diseases by regulating gene expression in various cells and tissues. Their roles in conditions like atherosclerosis, heart failure, myocardial infarction, and arrhythmias have been widely researched. The present study provides an overview of existing evidence regarding miRNAs' role in cardiovascular disease pathogenesis. Furthermore, the study examines current state-of-the-art technologies used in the study of miRNAs in cardiovascular disease. As a final point, we examine how miRNAs may serve as disease biomarkers, therapeutic targets, and prognostic indicators. In cardiology, microRNAs, small noncoding RNA molecules, are crucial to the posttranscriptional regulation of genes. Their role in regulating cardiac cell differentiation and maturation is critical during the development of the heart. They maintain the cardiac function of an adult heart by contributing to its electrical and contractile activity. By binding to messenger RNA molecules, they inhibit protein translation or degrade mRNA. Several cardiovascular diseases are associated with dysregulation of miRNAs, including arrhythmias, hypertension, atherosclerosis, and heart failure. miRNAs can be used as biomarkers to diagnose and predict diseases as well as therapeutic targets. A variety of state-of-the-art technologies have aided researchers in discovering, profiling, and analyzing miRNAs, including microarray analysis, next-generation sequencing, and others. Developing new diagnostics and therapeutic approaches is becoming more feasible as researchers refine their understanding of miRNA function. Ultimately, this will reduce the burden of cardiovascular disease around the world.

Heart failure with preserved ejection fraction in pigs causes shifts in posttranscriptional checkpoints.

Left ventricular pressure overload (LVPO) can lead to heart failure with a preserved ejection fraction (HFpEF) and LV chamber stiffness (LV Kc) is a hallmark. This project tested the hypothesis that the development of HFpEF due to an LVPO stimulus will alter posttranscriptional regulation, specifically microRNAs (miRs). LVPO was induced in pigs (n = 9) by sequential ascending aortic cuff and age- and weight-matched pigs (n = 6) served as controls. LV function was measured by echocardiography and LV Kc by speckle tracking. LV myocardial miRs were quantified using an 84-miR array. Treadmill testing and natriuretic peptide-A (NPPA) mRNA levels in controls and LVPO were performed (n = 10, n = 9, respectively). LV samples from LVPO and controls (n = 6, respectively) were subjected to RNA sequencing. LV mass and Kc increased by over 40% with LVPO (P < 0.05). A total of 30 miRs shifted with LVPO of which 11 miRs correlated to LV Kc (P < 0.05) that mapped to functional domains relevant to Kc such as fibrosis and calcium handling. LVPO resulted in reduced exercise tolerance (oxygen saturation, respiratory effort) and NPPA mRNA levels increased by fourfold (P < 0.05). RNA analysis identified several genes that mapped to specific miRs that were altered with LVPO. In conclusion, a specific set of miRs are changed in a large animal model consistent with the HFpEF phenotype, were related to LV stiffness properties, and several miRs mapped to molecular pathways that may hold relevance in terms of prognosis and therapeutic targets.NEW & NOTEWORTHY Heart failure with preserved ejection fraction (HFpEF) is an ever-growing cause for the HF burden. HFpEF is particularly difficult to treat as the mechanisms responsible for this specific form of HF are poorly understood. Using a relevant large animal model, this study uncovered a unique molecular signature with the development of HFpEF that regulates specific biological pathways relevant to the progression of this ever-growing cause of HF.

Unravelling alternative splicing patterns in susceptible and resistant Brassica napus lines in response to Xanthomonas campestris infection

BackgroundRapeseed (Brassica napus L.) is an important oil and industrial crop worldwide. Black rot caused by the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc) is an infectious vascular disease that leads to considerable yield losses in rapeseed. Resistance improvement through genetic breeding is an effective and sustainable approach to control black rot disease in B. napus. However, the molecular mechanisms underlying Brassica-Xcc interactions are not yet fully understood, especially regarding the impact of post-transcriptional gene regulation via alternative splicing (AS).ResultsIn this study, we compared the AS landscapes of a susceptible parental line and two mutagenized B. napus lines with contrasting levels of black rot resistance. Different types of AS events were identified in these B. napus lines at three time points upon Xcc infection, among which intron retention was the most common AS type. A total of 1,932 genes was found to show differential AS patterns between different B. napus lines. Multiple defense-related differential alternative splicing (DAS) hub candidates were pinpointed through an isoform-based co-expression network analysis, including genes involved in pathogen recognition, defense signalling, transcriptional regulation, and oxidation reduction.ConclusionThis study provides new insights into the potential effects of post-transcriptional regulation on immune responses in B. napus towards Xcc attack. These findings could be beneficial for the genetic improvement of B. napus to achieve durable black rot resistance in the future.

Comparative proteomic analysis provides insights into wood formation in immature xylem at different ages in Eucalyptus urophylla × Eucalyptus grandis

IntroductionWood formation is a crucial developmental stage in the life cycle of a woody plant; this process has substantial scientific research implications and practical applications. However, the mechanisms underlying woody plant development, especially the process of wood formation, remain poorly understood. As eucalyptus is one of the fastest growing tree species in the world, understanding the mechanism of wood formation in eucalyptus will greatly promote the development of molecular breeding technology for forest trees.ResultsIn this study, we investigated the proteomic profile of immature xylem at four different ages of Eucalyptus urophylla × Eucalyptus grandis (E. urograndis) using iTARQ technology. We identified 5236 proteins and 492 differentially abundant proteins (DAPs). The expression profiles of the DAPs corresponding to coding genes associated with wood formation were assessed using qRT-PCR. From the different expression profiles, it is inferred that the genes encoding kinesin, CDKD3, EXPA13, EXPA2, XTH27, EGases, UGT76E2, LAC, CCoAMT, CesA3, PAL, and CAD may undergo posttranscriptional regulation (PTR). Additionally, the genes encoding EIN2, ETR, MC4-like, and XCP may undergo posttranslational modifications (PTMs).ConclusionsWe investigated changes in wood formation-related proteins at the protein abundance level in the immature xylem of E. urograndis, thereby elucidating potential regulatory mechanisms of key proteins involved in eucalyptus wood formation. This study may provide theoretical guidance for further research on molecular breeding techniques and genetic improvement related to the cultivation of rapidly growing and high-quality trees.

Genetic Associations of Plasminogen Activator Inhibitor-1-Related miRNA Variants with Coronary Artery Disease

Coronary artery disease (CAD) is one of the most common types of cardiovascular disease and can lead to a heart attack as plaque gradually builds up inside the coronary arteries, blocking blood flow. Previous studies have shown that polymorphisms in the PAI-1 gene are associated with CAD; however, studies of the PAI-1 3′-untranslated region, containing a miRNA binding site, and the miRNAs that interact with it, are insufficient. To investigate the association between miRNA polymorphisms and CAD in the Korean population based on post-transcriptional regulation, we genotyped five polymorphisms in four miRNAs targeting the 3′-untranslated region of PAI-1 using real-time PCR and TaqMan assays. We found that the mutant genotype of miR-30c rs928508 A &gt; G was strongly associated with increased CAD susceptibility. In a genotype combination analysis, the combination of the homozygous mutant genotype (GG) of miR-30c rs928508 with the wild-type genotype (GG) of miR-143 rs41291957 resulted in increased risk for CAD. Also, in an allele combination analysis, the combination of the mutant allele (G) of miR-30c rs928508 and the wild-type allele (G) of miR-143 rs41291957 resulted in increased risk for CAD. Furthermore, metabolic syndrome and diabetes mellitus showed synergistic effects on CAD risk when combined with miR-30c rs928508. These results can be applied to identify CAD prognostic biomarkers among miRNA polymorphisms and various clinical factors.

Melatonin as an inducer of Th17 cell differentiation: new mechanisms

Cells of the immune system are sensitive to the action of melatonin, and the most obvious target of the hormone is the Th17 T helper subset: in addition to two high-affinity membrane receptors for melatonin, MT1 and MT2, these cells express a nuclear receptor, RORα. Among the secondary messengers involved in the transmission of signals from melatonin receptors, proteins of the sirtuin family are currently of particular interest. It is known that in non-tumor cells, sirtuin 1 (SIRT1) not only increases its expression in response to melatonin, but is also involved in the implementation of melatonin effects, as indicated by inhibitory assays using a specific SIRT1 blocker or corresponding siRNA/shRNA. This relates to the regulation of circadian oscillators, as well as the anti-inflammatory, antioxidant and anti-apoptotic effects of melatonin. Further mechanisms of SIRT1 activity in the cell include transcriptional and posttranscriptional regulation of gene expression through deacetylation of histones and non-histone proteins, and the apparent target of SIRT1 is the transcription factor RORα. It is this factor that mediates classical melatonin/SIRT1-dependent transcriptional control of key circadian regulators, the genes of which have ROR-binding sequences in their promoters. These data raise questions about the functions of SIRT1 in other cells expressing RORα, in particular in Th17 T helper cells, for which it is one of two key differentiation factors, along with RORγt. And if for RORα such a connection still remains hypothetical, for RORγt it has been convincingly demonstrated in studies both in vivo and in vitro: it has been shown that SIRT1 directly binds to RORγt and stimulates the development of Th17 cells, and blockade of sirtuin suppresses the differentiation of these cells in normal and prevents the development of Th17-associated pathology in mice. Summarizing these data, we can confidently predict the existence of a new mechanism for the regulation of the T helper Th17 population by melatonin, through the activation of sirtuin SIRT1, and this mechanism must be taken into account when interpreting data on the immunoregulatory activity of melatonin.

Identification of exosomal ceRNA networks as prognostic markers in clear cell renal cell carcinoma.

Aggressive clear cell renal cell carcinoma (ccRCC) has a bad prognosis. We seek new ccRCC biomarkers for diagnosis and treatment. We used exoRBase and The Cancer Genome Atlas Database to compare DEmRNAs, DEmiRNAs, DElncRNAs, and DEcircRNAs in ccRCC and normal renal tissues. CircRNAs and circRNAs targeting microRNAs (miRNAs) were anticipated and taken intersections, and several databases assessed the targeted link between common miRNAs and messenger RNAs (mRNAs). The Cancer Genome Atlas database was used to create a predictive mRNA signature that was validated in E-MTAB-1980. Finally, we examined competing endogenous RNA network miRNAs and long noncoding RNAs for ccRCC predictive biomarkers using overall survival analysis. We built the first competing endogenous RNA regulation network of circRNA-lncRNA-miRNA-mRNA and found that it substantially correlates with ccRCC prognosis. We unveiled ccRCC's posttranscriptional regulation mechanism in greater detail. Our findings identified novel biomarkers for ccRCC diagnosis, therapy, and prognosis.

TAF15 inhibits p53 nucleus translocation and promotes HCC cell 5-FU resistance via post-transcriptional regulation of UBE2N.

Chemotherapy resistance is an important factor responsible for the low 5-year survival rate of hepatocellular carcinoma (HCC) patients. Ubiquitin-conjugating enzyme E2N (UBE2N) is a cancer-associated ubiquitin-conjugating enzyme that is expressed in HCC tissues, and its high expression is associated with a poor prognosis. This study explored the role played by UBE2N in development of 5-fluorouracil (5-FU) resistance in HCC cells. Three HCC cell lines (HepG2 [p53 wild type], Huh7 [p53 point mutant type], Hep3B [p53 non-expression type]), and one normal liver cell line (MIHA) were used in our present study. The IC50 value of 5-FU was determined using a cell counting kit-8 (CCK-8) assay. Cell viability was assessed by colony formation assays. TUNEL assays and flow cytometry were used to analyze cell apoptosis. RNA pull-down and RNA immunoprecipitation (RIP) assays were performed to confirm the binding relationship between UBE2N mRNA and TAF15 protein. Our results showed that TAF15 and UBE2N were highly expressed in HCC cells. UBE2N inhibited the translocation of p53 protein into the cell nucleus to increase 5-FU resistance, as reflected by an increased IC50 value, an increase in cell viability, and a reduction in cell apoptosis. Overexpression of p53 reduced 5-FU resistance, but that effect could be reversed by UBE2N overexpression. TAF15 protein bound to and stabilized UBE2N mRNA, thereby inhibiting p53 translocation into the nucleus and promoting 5-FU resistance in HCC cells. Collectively, our present study identified a novel mechanism by which TAF15/UBE2N regulates p53 distribution to increase 5-FU resistance. Our results also suggest potential therapeutic strategies for treating HCC.

The YTHDF Proteins Shape the Brain Gene Signatures of Alzheimer's Disease.

The gene signatures of Alzheimer's Disease (AD) brains reflect an output of a complex interplay of genetic, epigenetic, epi-transcriptomic, and post-transcriptional regulations. To identify the most significant factor that shapes the AD brain signature, we developed a machine learning model (DEcode-tree) to integrate cellular and molecular factors explaining differential gene expression in AD. Our model indicates that YTHDF proteins, the canonical readers of N6-methyladenosine RNA modification (m6A), are the most influential predictors of the AD brain signature. We then show that protein modules containing YTHDFs are downregulated in human AD brains, and knocking out YTHDFs in iPSC-derived neural cells recapitulates the AD brain gene signature in vitro . Furthermore, eCLIP-seq analysis revealed that YTHDF proteins influence AD signatures through both m6A-dependent and independent pathways. These results indicate the central role of YTHDF proteins in shaping the gene signature of AD brains.

Evaluating the performance of multi-omics integration: a thyroid toxicity case study.

Multi-omics data integration has been repeatedly discussed as the way forward to more comprehensively cover the molecular responses of cells or organisms to chemical exposure in systems toxicology and regulatory risk assessment. In Canzler etal. (Arch Toxicol 94(2):371-388. https://doi.org/10.1007/s00204-020-02656-y ), we reviewed the state of the art in applying multi-omics approaches in toxicological research and chemical risk assessment. We developed best practices for the experimental design of multi-omics studies, omics data acquisition, and subsequent omics data integration. We found that multi-omics data sets for toxicological research questions were generally rare, with no datasets comprising more than two omics layers adhering to these best practices. Due to these limitations, we could not fully assess the benefits of different data integration approaches or quantitatively evaluate the contribution of various omics layers for toxicological research questions. Here, we report on a multi-omics study on thyroid toxicity that we conducted in compliance with these best practices. We induced direct and indirect thyroid toxicity through Propylthiouracil (PTU) and Phenytoin, respectively, in a 28-day plus 14-day recovery oral rat toxicity study. We collected clinical and histopathological data and six omics layers, including the long and short transcriptome, proteome, phosphoproteome, and metabolome from plasma, thyroid, and liver. We demonstrate that the multi-omics approach is superior to single-omics in detecting responses at the regulatory pathway level. We also show how combining omics data with clinical and histopathological parameters facilitates the interpretation of the data. Furthermore, we illustrate how multi-omics integration can hint at the involvement of non-coding RNAs in post-transcriptional regulation. Also, we show that multi-omics facilitates grouping, and we assess how much information individual and combinations of omics layers contribute to this approach.

Alternative splicing in the genome of HPV and its regulation.

Persistent infection with high-risk human papillomavirus (HR-HPV) is the main cause of cervical cancer. These chronic infections are characterized by high expression of the HPV E6 and E7 oncogenes and the absence of the L1 and L2 capsid proteins. The regulation of HPV gene expression plays a crucial role in both the viral life cycle and rare oncogenic events. Alternative splicing of HPV mRNA is a key mechanism in post-transcriptional regulation. Through alternative splicing, HPV mRNA is diversified into various splice isoforms with distinct coding potentials, encoding multiple proteins and influencing the expression of HPV genes. The spliced mRNAs derived from a donor splicing site within the E6 ORF and one of the different acceptor sites located in the early mRNA contain E6 truncated mRNAs, named E6*. E6* is one of the extensively studied splicing isoforms. However, the role of E6* proteins in cancer progression remains controversial. Here, we reviewed and compared the alternative splicing events occurring in the genomes of HR-HPV and LR-HPV. Recently, new HPV alternative splicing regulatory proteins have been continuously discovered, and we have updated the regulation of HPV alternative splicing. In addition, we summarized the functions of known splice isoforms from three aspects: anti-tumorigenic, tumorigenic, and other cancer-related functions, including not only E6*, but also E6^E7, E8^E2, and so on. Comprehending their contributions to cancer development enhances insights into the carcinogenic mechanisms of HPV and explores the potential utility of alternative splicing in the diagnosis and treatment of cervical cancer.

TGIRT-seq of Inflammatory Breast Cancer Tumor and Blood Samples Reveals Widespread Enhanced Transcription Impacting RNA Splicing and Intronic RNAs in Plasma.

Inflammatory breast cancer (IBC) is the most aggressive and lethal breast cancer subtype but lacks unequivocal genomic differences or robust biomarkers that differentiate it from non-IBC. Here, Thermostable Group II intron Reverse Transcriptase RNA-sequencing (TGIRT-seq) revealed myriad differences in tumor samples, Peripheral Blood Mononuclear Cells (PBMCs), and plasma that distinguished IBC from non-IBC patients and healthy donors across all tested receptor-based subtypes. These included numerous differentially expressed protein-coding gene and non-coding RNAs in all three sample types, a granulocytic immune response in IBC PBMCs, and over- expression of antisense RNAs, suggesting wide-spread enhanced transcription in both IBC tumors and PBMCs. By using TGIRT-seq to quantitate Intron-exon Depth Ratios (IDRs) and mapping reads to both genome and transcriptome reference sequences, we developed methods for parallel analysis of transcriptional and post-transcriptional gene regulation. This analysis identified numerous differentially and non-differentially expressed protein-coding genes in IBC tumors and PBMCs with high IDRs, the latter reflecting rate-limiting RNA splicing that negatively impacts mRNA production. Mirroring gene expression differences in tumors and PBMCs, over-represented protein-coding gene RNAs in IBC patient plasma were largely intronic RNAs, while those in non- IBC patients and healthy donor plasma were largely mRNA fragments. Potential IBC biomarkers in plasma included T-cell receptor pre-mRNAs and intronic, LINE-1, and antisense RNAs. Our findings provide new insights into IBC and set the stage for monitoring disease progression and response to treatment by liquid biopsy. The methods developed for parallel transcriptional and post- transcriptional gene regulation analysis have potentially broad RNA-seq and clinical applications.

RMVar 2.0: an updated database of functional variants in RNA modifications.

Evaluating the impact of genetic variants on RNA modifications (RMs) is crucial for identifying disease-associated variants and understanding the pathogenic mechanisms underlying human diseases. Previously, we developed a database called RMVar to catalog variants linked to RNA modifications in humans and mice. Here, we present an updated version RMVar 2.0 (http://rmvar.renlab.cn). In this updated version, we applied an enhanced analytical pipeline to the latest RNA modification datasets and genetic variant information to identify RM-associated variants. A notable advancement in RMVar 2.0 is our incorporation of allele-specific RNA modification analysis to identify RM-associated variants, a novel approach not utilized in RMVar 1.0 or other comparable databases. Furthermore, the database offers comprehensive annotations for various molecular events, including RNA-binding protein (RBP) interactions, RNA-RNA interactions, splicing events, and circular RNAs (circRNAs), which facilitate investigations into how RM-associated variants influence post-transcriptional regulation. Additionally, we provide disease-related information sourced from ClinVar and GWAS to help researchers explore the connections between RNA modifications and various diseases. We believe that RMVar 2.0 will significantly enhance our understanding of the functional implications of genetic variants affecting RNA modifications within the context of human disease research.

Pathological mechanisms and candidate therapeutic approaches in the hearing loss of mice carrying human MIR96 mutations.

Progressive hearing loss is a common problem in the human population with no effective therapeutics currently available. However, it has a strong genetic contribution, and investigating the genes and regulatory interactions underlying hearing loss offers the possibility of identifying therapeutic candidates. Mutations in regulatory genes are particularly useful for this, and an example is the microRNA miR-96, a post-transcriptional regulator which controls hair cell maturation. Mice and humans carrying mutations in miR-96 all exhibit hearing impairment, in homozygosis if not in heterozygosis, but different mutations result in different physiological, structural and transcriptional phenotypes. Here we present our characterisation of two lines of mice carrying different human mutations knocked-in to Mir96. We have carried out auditory brainstem response tests to examine their hearing with age and after noise exposure and have used confocal and scanning electron microscopy to examine the ultrastructure of the organ of Corti and hair cell synapses. Bulk RNA-seq was carried out on the organs of Corti of postnatal mice, followed by bioinformatic analyses to identify candidate targets. While mice homozygous for either mutation are profoundly deaf from 2weeks old, the heterozygous phenotypes differ markedly, with only one mutation resulting in hearing impairment in heterozygosis. Investigations of the structural phenotype showed that one mutation appears to lead to synaptic defects, while the other has a much more severe effect on the hair cell stereociliary bundles. Transcriptome analyses revealed a wide range of misregulated genes in both mutants which were notably dissimilar. We used the transcriptome analyses to investigate candidate therapeutics, and tested one, finding that it delayed the progression of hearing loss in heterozygous mice. Our work adds further support for the importance of the gain of novel targets in microRNA mutants and offers a proof of concept for the identification of pharmacological interventions to maintain hearing.

RNA-binding protein IGF2BP1 is required for spermatogenesis in an age-dependent manner.

Post-transcriptional regulation mediated by RNA binding proteins is crucial for male germline development. Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1), an RNA binding protein, is specifically expressed in human and mouse male gonads and is involved in manifold biological processes and tumorigenesis. However, the function of IGF2BP1 in mammalian spermatogenesis remains poorly understood. Herein, we generated an Igf2bp1 conditional knockout mouse model using Nanos3-Cre. Germ cell deficiency of Igf2bp1 in mice caused spermatogenic defects in an age-dependent manner, resulting in decreased numbers of undifferentiated spermatogonia and increased germ cell apoptosis. Immunoprecipitation-mass spectrometry analysis revealed that ELAV-like RNA binding protein 1, a well-recognized mRNA stabilizer, interacted with IGF2BP1. Single cell RNA-sequencing showed distinct mRNA profiles in spermatogonia from conditional knockout versus wide type mice. Further research showed that IGF2BP1 plays a vital role in the modulation of spermatogenesis by regulating Lin28a mRNA, which is essential for clonal expansion of undifferentiated spermatogonia. Thus, our results highlight the crucial effects of IGF2BP1 on spermatogonia for the long-term maintenance of spermatogenesis.