Cis-acting noncoding RNAs: friends and foes

Investigation of Some Long Noncoding RNAs (LncRNAs) in Pediatric Inflammatory Bowel Disease (IBD): An Iranian Study

Introduction: According to the importance of long noncoding RNAs (LncRNA) in the pathogenesis of inflammatory bowel disease (IBD) and also the lack of study for pediatric IBD in this regard, we investigated the expression of a selected panel of LncRNAs in Iranian pediatric cases of IBD compared to adult cases and healthy samples. Methods: In this gene expression study, blood samples were taken from the three groups of pediatric IBD cases, adult IBD cases, and pediatric healthy samples (for gene expression calibration). The investigated LncRNAs were UCA1, CCAT, IFNG-AS1, and CDKN2B. Real-time PCR was used and fold changes (FCs) were reported. Results: A total of 50 individuals were studied including 28 cases of pediatric IBD, 12 cases of controls, and 10 cases of adult IBD. UCA1 showed upregulation in adult IBD (FC = 10.56, p = 0.007). CCAT showed downregulations for pediatric IBD (FC = 0.01, p < 0.001) and adult IBD (FC = 0.10, p = 0.039). IFNG-AS1 showed downregulation in pediatric IBD (FC < 0.01, p < 0.001). CDKN2B showed upregulation in pediatric IBD (FC = 17.39, p < 0.001). The results were in contrast with the literature. Conclusion: It seems that these LncRNAs may have different roles in pediatric IBD. Further studies are needed on pediatric cases of IBD.

lncRNA CARINH regulates expression and function of innate immune transcription factor IRF1 in macrophages

The discovery of long non-coding RNAs (lncRNAs) has provided a new perspective on the centrality of RNA in gene regulation and genome organization. Here, we screened for lncRNAs with putative functions in the host response to single-stranded RNA respiratory viruses. We identify CARINH as a conserved cis-acting lncRNA up-regulated in three respiratory diseases to control the expression of its antisense gene IRF1, a key transcriptional regulator of the antiviral response. CARINH and IRF1 are coordinately increased in the circulation of patients infected with human metapneumovirus, influenza A virus, or SARS-CoV-2, and in macrophages in response to viral infection or TLR3 agonist treatment. Targeted depletion of CARINH or its mouse ortholog Carinh in macrophages reduces the expression of IRF1/Irf1 and their associated target gene networks, increasing susceptibility to viral infection. Accordingly, CRISPR-mediated deletion of Carinh in mice reduces antiviral immunity, increasing viral burden upon sublethal challenge with influenza A virus. Together, these findings identify a conserved role of lncRNA CARINH in coordinating interferon-stimulated genes and antiviral immune responses.

Identification of crucial LncRNAs associated with testicular development and LOC108635509 as a potential regulator in black goat spermatogenesis

The establishment and maintenance of spermatogenesis is a complex process involving a vast of regulatory pathways. There is growing evidence revealing that long noncoding RNAs (lncRNA) play important roles in regulating testicular development and spermatogenesis in a stage-specific way. However, our understanding of how lncRNA regulates testicular development and spermatogenesis in black goats is quite limited. In the current study, we screened the transcriptomes (lncRNA and mRNA) of testicular from Guangxi black goats before puberty (3 days old, D3; 30 days old, D30), puberty (90 days old, D90) and postpuberty (180 days old, D180), in order to identify the lncRNA interaction with mRNAs contributes to goat spermatogenesis. The RNA-sequencing (RNA-seq) analysis showed that there were 1211, 12,180, 834 differential lncRNAs and 1196, 8838,269 differential mRNAs at the ages of D30 vs. D3, D90 vs. D30, and D180 vs. D90. The lncRNAs showed the most significantly changes from D30 to D90, which indicated that D90 was a key node of lncRNAs participated in the regulation of testicular development and spermatogenesis in black goat. According to functional enrichment analysis of GO and KEGG, we found that differentially expressed lncRNAs (DE lncRNAs) and their target genes regulated spermatogenesis through signal pathways including MAPK, Ras, and PI3K-Akt. Using cis- and trans-acting, 39 DE lncRNAs-targeted genes were found to be enriched for male reproduction. Of these, LOC108635509, which specific expressed in testis and upregulated the expression levels at D90, was found participated in the regulation of testicular development through promoting the proliferation of Sertoli cells (SCs). Overall, this study provides new insight into the regulatory mechanisms that support spermatogenesis and testicular development in black goats.

LncRNA (BCO1-AS) regulate inflammatory responses in bacterial infection through caspase-1 in turbot (Scophthalmus maximus)

LncRNA plays key role in several biological processes, including transcriptional regulation, post transcriptional control and epigenetic regulation. However, research on the functional roles of lncRNAs in teleost species remains limited. Here, we discovered a lncRNA (BCO1-AS) with a critical role in antibacterial responses. Briefly, the full length of BCO1-AS was 2005 bp. Subsequently, BCO1-AS was distributed throughout the nucleus, where it may either trans- or cis-regulate the nearby genes. In addition, BCO1-AS was widely expressed in all the examined tissues with the highest expression level in intestine, while the lowest expression level was detected in muscle. Moreover, following Vibrio anguillarum challenge, BCO1-AS was significantly down-regulated in intestine, and up-regulated in gill and skin. In CHIRP experiment, BCO1-AS could effectively enrich RNA and might interact with several immune-related genes. Furthermore, we found that LPS could induce the expression of BCO1-AS. Finally, BCO1-AS could positively regulate caspase-1 at the mRNA and protein level. The BCO1-AS was speculated to inhibit the synthesis of inflammatory components. In summary, these results showed the roles of BCO1-AS in the regulation of inflammatory in turbot, which provided valuable information for further understanding the immune regulation network of lncRNA in teleost fish.

piRNA PROPER Suppresses DUSP1 Translation by Targeting N6-Methyladenosine-Mediated RNA Circularization to Promote Oncogenesis of Prostate Cancer

Genetic and epigenetic alterations occur in many physiological and pathological processes. The existing knowledge regarding the association of PIWI-interacting RNAs (piRNAs) and their genetic variants on risk and progression of prostate cancer (PCa) is limited. In this study, three genome-wide association study datasets are combined, including 85,707 PCa cases and 166,247 controls, to uncover genetic variants in piRNAs. Functional investigations involved manipulating piRNA expression in cellular and mouse models to study its oncogenetic role in PCa. A specific genetic variant, rs17201241 is identified, associated with increased expression of PROPER (piRNA overexpressed in prostate cancer) in tumors and are located within the gene, conferring an increased risk and malignant progression of PCa. Mechanistically, PROPER coupled with YTHDF2 to recognize N6-methyladenosine (m6A) and facilitated RNA-binding protein interactions between EIF2S3 at 5'-untranslated region (UTR) and YTHDF2/YBX3 at 3'-UTR to promote DUSP1 circularization. This m6A-dependent mRNA-looping pattern enhanced DUSP1 degradation and inhibited DUSP1 translation, ultimately reducing DUSP1 expression and promoting PCa metastasis via the p38 mitogen-activated protein kinase (MAPK) signaling pathway. Inhibition of PROPER expression using antagoPROPER effectively suppressed xenograft growth, suggesting its potential as a therapeutic target. Thus, targeting piRNA PROPER-mediated genetic and epigenetic fine control is a promising strategy for the concurrent prevention and treatment of PCa.

Cut from the same cloth: RNAs transcribed from regulatory elements

A certain degree of chromatin openness is necessary for the activity of transcription-regulating regions within the genome, facilitating accessibility to RNA polymerases and subsequent synthesis of regulatory element RNAs (regRNAs) from these regions. The rapidly increasing number of studies underscores the significance of regRNAs across diverse cellular processes and diseases, challenging the paradigm that these transcripts are non-functional transcriptional noise. This review explores the multifaceted roles of regRNAs in human cells, encompassing rather well-studied entities such as promoter RNAs and enhancer RNAs (eRNAs), while also providing insights into overshadowed silencer RNAs and insulator RNAs. Furthermore, we assess notable examples of shorter regRNAs, like miRNAs, snRNAs, and snoRNAs, playing important roles. Expanding our discourse, we deliberate on the potential usage of regRNAs as biomarkers and novel targets for cancer and other human diseases.

Generation of New β-Conglycinin-Deficient Soybean Lines by Editing the lincRNA lincCG1 Using the CRISPR/Cas9 System

Soybean β-conglycinin is a major allergen that adversely affects the nutritional properties of soybean. Soybean deficient in β-conglycinin is associated with low allergenicity and high nutritional value. Long intergenic noncoding RNAs (lincRNAs) regulate gene expression and are considered important regulators of essential biological processes. Despite increasing knowledge of the functions of lincRNAs, relatively little is known about the effects of lincRNAs on the accumulation of soybean β-conglycinin. The current study presents the identification of a lincRNA lincCG1 that was mapped to the intergenic noncoding region of the β-conglycinin α-subunit locus. The full-length lincCG1 sequence was cloned and found to regulate the expression of soybean seed storage protein (SSP) genes via both cis- and trans-acting regulatory mechanisms. Loss-of-function lincCG1 mutations generated using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system led to the deficiency of the allergenic α'-, α-, and β-subunits of soybean β-conglycinin as well as higher content of proteins, sulfur-containing amino acids, and free arginine. The dominant null allele LincCG1, and consequently, the β-conglycinin-deficient phenotype associated with the lincCG1-gene-edited line was stably inherited by the progenies in a Mendelian fashion. The dominant null allele LincCG1 may therefore be exploited for engineering/developing novel hypoallergenic soybean varieties. Furthermore, Cas9-free and β-conglycinin-deficient homozygous mutant lines were obtained in the T1 generation. This study is the first to employ the CRISPR/Cas9 technology for editing a lincRNA gene associated with the soybean allergenic protein β-conglycinin. Moreover, this study reveals that lincCG1 plays a crucial role in regulating the expression of the β-conglycinin subunit gene cluster, besides highlighting the efficiency of employing the CRISPR/Cas9 system for modulating lincRNAs, and thereby regulating soybean seed components.

SH-SY5Y human neuronal cells with mutations of the CDKN2B-AS1 gene are vulnerable under cultured conditions

Glaucoma is a common cause of blindness worldwide. Genetic effects are believed to contribute to the onset and progress of glaucoma, but the underlying pathological mechanisms are not fully understood. Here, we set out to introduce mutations into the CDKN2B-AS1 gene, which is known as being the closely associated with glaucoma, in a human neuronal cell line in vitro. We introduced gene mutations with CRISPR/Cas9 into exons and introns into the CDKN2B-AS1 gene. Both mutations strongly promoted neuronal cell death in normal culture conditions. RNA sequencing and pathway analysis revealed that the transcriptional factor Fos is a target molecule regulating CDKN2B-AS1 overexpression. We demonstrated that gene mutation of CDKN2B-AS1 is directly associated with neuronal cell vulnerability in vitro. Additionally, Fos, which is a downstream signaling molecule of CDKN2B-AS1, may be a potential source of new therapeutic targets for neuronal degeneration in diseases such as glaucoma.

m6A- and m5C- modified lncRNAs orchestrate the prognosis in cutaneous melanoma and m6A- modified LINC00893 regulates cutaneous melanoma cell metastasis

BACKGROUND

As the most important modifications on the RNA level, N6-methyladenosine (m6A-) and 5-methylcytosine (m5C-) modification could have a direct influence on the RNAs. Long non-coding RNAs (lncRNAs) could also be modified by methylcytosine modification. Compared with mRNAs, the function of lncRNAs could be more potent to some extent in biological processes like tumorigenesis. Until now, rare reports have been done associated with cutaneous melanoma. Herein, we wonder if the m6A- and m5C- modified lncRNAs could influence the immune landscape and prognosis in melanoma, and we also want to find some lncRNAs which could directly affect the malignant behaviors of melanoma.

METHODS

Systematically, we explored the expression pattern of m6A- and m5C- modified lncRNAs in melanoma from datasets including UCSC Xena and NCBI GEO, and the prognostic lncRNAs were selected. Then, according to the expression pattern of lncRNAs, melanoma samples from these datasets were divided into several subtypes. Prognostic model, nomogram survival model, drug sensitivity, GO, and KEGG pathway analysis were performed. Furthermore, among several selected lncRNAs, we identified one lncRNA named LINC00893 and investigated its expression pattern and its biological function in melanoma cell lines.

RESULTS

We identified 27 m6A- and m5C- related lncRNAs which were significantly associated with survival, and we made a subtype analysis of melanoma samples based on these 27 lncRNAs. Among the two subtypes, we found differences of immune cells infiltration between these two subtypes. Then, LASSO algorithm was used to screen the optimized lncRNAs combination including ZNF252P-AS1, MIAT, FAM13A-AS1, LINC-PINT, LINC00893, AGAP2-AS1, OIP5-AS1, and SEMA6A-AS1. We also found that there was a significant correlation between the different risk groups predicted based on RS model and the actual prognosis. The nomogram survival model based on independent survival prognostic factors was also constructed. Besides, sensitivity to chemotherapeutic agents, GO and KEGG analysis were performed. In different risk groups, a total of 14 drug molecules with different distributions were obtained, which included AZD6482, AZD7762, AZD8055, camptothecin, dasatinib, erlotinib, gefitinib, gemcitabine, GSK269962A, nilotinib, rapamycin, and sorafenib. A total of 55 significantly related biological processes and 17 KEGG signaling pathways were screened. At last, we noticed that LINC00893 had a relatively lower expression in melanoma tissue and cell lines compared with adjacent tissues and epidermal melanocyte, and down-regulation of LINC00893 could promote the malignant behavior of melanoma cells in A875 and MV3. In these two melanoma cell lines, down-regulation of m6A-related molecules like YTHDF3 and METTL3 could promote the expression of LINC00893.

CONCLUSION

We made an analysis of m6A- and m5C- related lncRNAs in melanoma samples and a prediction of these lncRNAs' role in prognosis, tumor microenvironment, immune infiltration, and clinicopathological features. We also found that LINC00893, which is potentially regulated by m6A modification, could serve as a tumor-suppressor in melanoma and play an inhibitory role in melanoma metastasis.

Novel insights on the positive correlation between sense and antisense pairs on gene expression

A considerable proportion of the eukaryotic genome undergoes transcription, leading to the generation of noncoding RNA molecules that lack protein-coding information and are not subjected to translation. These noncoding RNAs (ncRNAs) are well recognized to have essential roles in several biological processes. Long noncoding RNAs (lncRNAs) represent the most extensive category of ncRNAs found in the human genome. Much research has focused on investigating the roles of cis-acting lncRNAs in the regulation of specific target gene expression. In the majority of instances, the regulation of sense gene expression by its corresponding antisense pair occurs in a negative (discordant) manner, resulting in the suppression of the target genes. The notion that a negative correlation exists between sense and antisense pairings is, however, not universally valid. In fact, several recent studies have reported a positive relationship between corresponding cis antisense pairs within plants, budding yeast, and mammalian cancer cells. The positive (concordant) correlation between anti-sense and sense transcripts leads to an increase in the level of the sense transcript within the same genomic loci. In addition, mechanisms such as altering chromatin structure, the formation of R loops, and the recruitment of transcription factors can either enhance transcription or stabilize sense transcripts through their antisense pairs. The primary objective of this work is to provide a comprehensive understanding of both aspects of antisense regulation, specifically focusing on the positive correlation between sense and antisense transcripts in the context of eukaryotic gene expression, including its implications towards cancer progression. This article is categorized under: RNA Processing > 3' End Processing Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

KDM7A-DT induces genotoxic stress, tumorigenesis, and progression of p53 missense mutation-associated invasive breast cancer

Stress-induced promoter-associated and antisense lncRNAs (si-paancRNAs) originate from a reservoir of oxidative stress (OS)-specific promoters via RNAPII pausing-mediated divergent antisense transcription. Several studies have shown that the KDM7A divergent transcript gene (KDM7A-DT), which encodes a si-paancRNA, is overexpressed in some cancer types. However, the mechanisms of this overexpression and its corresponding roles in oncogenesis and cancer progression are poorly understood. We found that KDM7A-DT expression is correlated with highly aggressive cancer types and specific inherently determined subtypes (such as ductal invasive breast carcinoma (BRCA) basal subtype). Its regulation is determined by missense TP53 mutations in a subtype-specific context. KDM7A-DT transcribes several intermediate-sized ncRNAs and a full-length transcript, exhibiting distinct expression and localization patterns. Overexpression of KDM7A-DT upregulates TP53 protein expression and H2AX phosphorylation in nonmalignant fibroblasts, while in semi-transformed fibroblasts, OS superinduces KDM7A-DT expression in a TP53-dependent manner. KDM7A-DT knockdown and gene expression profiling in TP53-missense mutated luminal A BRCA variant, where it is abundantly expressed, indicate its significant role in cancer pathways. Endogenous over-expression of KDM7A-DT inhibits DNA damage response/repair (DDR/R) via the TP53BP1-mediated pathway, reducing apoptosis and promoting G2/M checkpoint arrest. Higher KDM7A-DT expression in BRCA is associated with KDM7A-DT locus gain/amplification, higher histologic grade, aneuploidy, hypoxia, immune modulation scores, and activation of the c-myc pathway. Higher KDM7A-DT expression is associated with relatively poor survival outcomes in patients with luminal A or Basal subtypes. In contrast, it is associated with favorable outcomes in patients with HER2+ER- or luminal B subtypes. KDM7A-DT levels are coregulated with critical transcripts and proteins aberrantly expressed in BRCA, including those involved in DNA repair via non-homologous end joining and epithelial-to-mesenchymal transition pathway. In summary, KDM7A-DT and its si-lncRNA exhibit several intrinsic biological and clinical characteristics that suggest important roles in invasive BRCA and its subtypes. KDM7A-DT-defined mRNA and protein subnetworks offer resources for identifying clinically relevant RNA-based signatures and prospective targets for therapeutic intervention.

Identification of Co-diagnostic Genes for Heart Failure and Hepatocellular Carcinoma Through WGCNA and Machine Learning Algorithms

This research delves into the intricate relationship between hepatocellular carcinoma (HCC) and heart failure (HF) by exploring shared genetic characteristics and molecular processes. Employing advanced methodologies such as differential analysis, weighted correlation network analysis (WGCNA), and algorithms like Random Forest (RF), Least Absolute Shrinkage Selection (LASSO), and XGBoost, we meticulously identified modular differential genes (DEGs) associated with both HF and HCC. Gene Set Variation Analysis (GSVA) and single sample gene set enrichment analysis (ssGSEA) were employed to unveil underlying biological mechanisms. The study revealed 88 core genes shared between HF and HCC, indicating a common mechanism. Enrichment analysis emphasized the roles of immune responses and inflammation in both diseases. Leveraging XGBoost, we crafted a robust multigene diagnostic model (including FCN3, MAP2K1, AP3M2, CDH19) with an area under the curve (AUC) > 0.9, showcasing exceptional predictive accuracy. GSVA and ssGSEA analyses unveiled the involvement of immune cells and metabolic pathways in the pathogenesis of HF and HCC. This research uncovers a pivotal interplay between HF and HCC, highlighting shared pathways and key genes, offering promising insights for future clinical treatments and experimental research endeavors.

Whole-transcriptome sequencing reveals a melanin-related ceRNA regulatory network in the breast muscle of Xichuan black-bone chicken

The economic losses incurred due to reduced muscle pigmentation highlight the crucial role of melanin-based coloration in the meat of black-bone chickens. Melanogenesis in the breast muscle of black-bone chickens is currently poorly understood in terms of molecular mechanisms. This study employed whole-transcriptome sequencing to analyze black and white breast muscle samples from black-bone chickens, leading to the identification of 367 differentially expressed (DE) mRNAs, 48 DElncRNAs, 104 DEcircRNAs, and 112 DEmiRNAs involved in melanin deposition. Based on these findings, a competitive endogenous RNA (ceRNA) network was developed to better understand the complex mechanisms of melanin deposition. Furthermore, our analysis revealed key DEmRNAs (TYR, DCT, EDNRB, MLPH and OCA2) regulated by DEmiRNAs (gga-miR-140-5p, gga-miR-1682, gga-miR-3529, gga-miR-499-3p, novel-m0012-3p, gga-miR-200b-5p, gga-miR-203a, gga-miR-6651-5p, gga-miR-7455-3p, gga-miR-31-5p, miR-140-x, miR-455-x, novel-m0065-3p, gga-miR-29b-1-5p, miR-455-y, novel-m0085-3p, and gga-miR-196-1-3p). These DEmiRNAs competitively interacted with DElncRNAs including MSTRG.2609.2, MSTRG.4185.1, LOC112530666, LOC112533366, LOC771030, LOC107054724, LOC121107411, LOC100859072, LOC101750037, LOC121108550, LOC121109224, LOC121110876, and LOC101749016, as well as DEcircRNAs, such as novel_circ_000158, novel_circ_000623, novel_001518, and novel_circ_003596. The findings from this study provide insight into the mechanisms that regulate lncRNA, circRNA, miRNA, and mRNA expression in chicken melanin deposition.

Deficiency of lncRNA MERRICAL abrogates macrophage chemotaxis and diabetes-associated atherosclerosis

Diabetes-associated atherosclerosis involves excessive immune cell recruitment and plaque formation. However, the mechanisms remain poorly understood. Transcriptomic analysis of the aortic intima in Ldlr-/- mice on a high-fat, high-sucrose-containing (HFSC) diet identifies a macrophage-enriched nuclear long noncoding RNA (lncRNA), MERRICAL (macrophage-enriched lncRNA regulates inflammation, chemotaxis, and atherosclerosis). MERRICAL expression increases by 249% in intimal lesions during progression. lncRNA-mRNA pair genomic mapping reveals that MERRICAL positively correlates with the chemokines Ccl3 and Ccl4. MERRICAL-deficient macrophages exhibit lower Ccl3 and Ccl4 expression, chemotaxis, and inflammatory responses. Mechanistically, MERRICAL guides the WDR5-MLL1 complex to activate CCL3 and CCL4 transcription via H3K4me3 modification. MERRICAL deficiency in HFSC diet-fed Ldlr-/- mice reduces lesion formation by 74% in the aortic sinus and 86% in the descending aorta by inhibiting leukocyte recruitment into the aortic wall and pro-inflammatory responses. These findings unveil a regulatory mechanism whereby a macrophage-enriched lncRNA potently inhibits chemotactic responses, alleviating lesion progression in diabetes.

V. S, Kaushik K, Sabharwal A, Lalwani MK, K. S, Singh N, Scaria V, Sivasubbu S. Forward genetic screen using a gene-breaking trap approach identifies a novel role of grin2bb-associated RNA transcript (grin2bbART) in zebrafish heart function

LncRNA-based control affects cardiac pathophysiologies like myocardial infarction, coronary artery disease, hypertrophy, and myotonic muscular dystrophy. This study used a gene-break transposon (GBT) to screen zebrafish (Danio rerio) for insertional mutagenesis. We identified three insertional mutants where the GBT captured a cardiac gene. One of the adult viable GBT mutants had bradycardia (heart arrhythmia) and enlarged cardiac chambers or hypertrophy; we named it "bigheart." Bigheart mutant insertion maps to grin2bb or N-methyl D-aspartate receptor (NMDAR2B) gene intron 2 in reverse orientation. Rapid amplification of adjacent cDNA ends analysis suggested a new insertion site transcript in the intron 2 of grin2bb. Analysis of the RNA sequencing of wild-type zebrafish heart chambers revealed a possible new transcript at the insertion site. As this putative lncRNA transcript satisfies the canonical signatures, we called this transcript grin2bb associated RNA transcript (grin2bbART). Using in situ hybridization, we confirmed localized grin2bbART expression in the heart, central nervous system, and muscles in the developing embryos and wild-type adult zebrafish atrium and bulbus arteriosus. The bigheart mutant had reduced Grin2bbART expression. We showed that bigheart gene trap insertion excision reversed cardiac-specific arrhythmia and atrial hypertrophy and restored grin2bbART expression. Morpholino-mediated antisense downregulation of grin2bbART in wild-type zebrafish embryos mimicked bigheart mutants; this suggests grin2bbART is linked to bigheart. Cardiovascular tissues use Grin2bb as a calcium-permeable ion channel. Calcium imaging experiments performed on bigheart mutants indicated calcium mishandling in the heart. The bigheart cardiac transcriptome showed differential expression of calcium homeostasis, cardiac remodeling, and contraction genes. Western blot analysis highlighted Camk2d1 and Hdac1 overexpression. We propose that altered calcium activity due to disruption of grin2bbART, a putative lncRNA in bigheart, altered the Camk2d-Hdac pathway, causing heart arrhythmia and hypertrophy in zebrafish.

Long noncoding RNA lncPostn links TGF-β and p53 signaling pathways to transcriptional regulation of cardiac fibrosis

Cardiac fibroblasts are essential for the homeostasis of the extracellular matrix, whose remodeling in many cardiovascular diseases leads to fibrosis. Long noncoding RNAs (lncRNAs) are associated with cardiac pathologies, but their functions in cardiac fibroblasts and contributions to cardiac fibrosis remain unclear. Here, we aimed to identify fibroblast-enriched lncRNAs essential in myocardial infarction (MI)-induced fibrosis and explore the molecular mechanisms responsible for their functions. Global lncRNA profiling was performed in post-MI mouse heart ventricles and transforming growth factor-β (TGF-β)-treated primary cardiac fibroblasts and confirmed in published data sets. We identified the cardiac fibroblast-enriched lncPostn, whose expression is stimulated in cardiac fibrosis induced by MI and the extracellular growth factor TGF-β. The promoter of lncPostn contains a functional TGF-β response element, and lncPostn knockdown suppresses TGF-β-stimulated cardiac fibroblast activation and improves cardiac functions post-MI. LncPostn stabilizes and recruits EP300 to the profibrotic periostin's promoter, representing a major mechanism for its transcriptional activation. Moreover, both MI and TGF-β enhance lncPostn expression while suppressing the cellular growth gatekeeper p53. TGF-β and p53 knockdown-induced profibrotic gene expression and fibrosis occur mainly through lncPostn and show additive effects. Finally, levels of serum lncPostn are significantly increased in patients' postacute MI and show a strong correlation with fibrosis markers, revealing a potential biomarker of cardiac fibrosis. Our findings identify the fibroblast-enriched lncPostn as a potent profibrotic factor, providing a transcriptional link between TGF-β and p53 signaling pathways to regulate fibrosis in cardiac fibroblasts.NEW & NOTEWORTHY Cardiac fibroblasts are essential for the homeostasis of the extracellular matrix, whose remodeling in many cardiovascular diseases leads to fibrosis. Long noncoding RNAs are functional and contribute to the biological processes of cardiovascular development and disorders. Our findings identify the fibroblast-enriched lncPostn as a potent profibrotic factor and demonstrate that serum lncPostn level may serve as a potential biomarker of human cardiac fibrosis postacute myocardial infarction.

Non-coding RNA in the gut of the blood-feeding parasitic worm, Haemonchus contortus

The intestine of Haemonchus contortus is an essential tissue that has been indicated to be a major target for the prevention of haemonchosis caused by this parasitic nematode of small ruminants. Biological peculiarities of the intestine warrant in-depth exploitation, which can be leveraged for future disease control efforts. Here, we determined the intestinal ncRNA (lncRNA, circRNA and miRNA) atlas using whole-transcriptome sequencing and bioinformatics approaches. In total, 4846 novel lncRNA, 982 circRNA, 96 miRNA (65 known and 31 novel) and 8821 mRNA were identified from the H. contortus intestine. The features of lncRNA, circRNA and miRNA were fully characterized. Comparison of miRNA from the intestines and extracellular vesicles supported the speculation that the miRNA from the latter were of intestinal origin in H. contortus. Further function analysis suggests that the cis-lncRNA targeted genes were involved in protein binding, intracellular anatomical structure, organelle and cellular process, whereas the circRNA parental genes were mainly enriched in molecular function categories, such as ribonucleotide binding, nucleotide binding, ATP binding and carbohydrate derivative binding. The miRNA target genes were related to the cellular process, cellular response to stimulus, cellular protein modification process and signal transduction. Moreover, competing endogenous RNA network analysis revealed that the majority of lncRNA, circRNA and mRNA only have one or two binding sites with specific miRNA. Lastly, randomly selected circRNA, lncRNA and miRNA were verified successfully using RT-PCR. Collectively, these data provide the most comprehensive compilation of intestinal transcripts and their functions, and it will be helpful to decipher the biological and molecular complexity of the intestine and lay the foundation for further functional research.

Histone lactylation-derived LINC01127 promotes the self-renewal of glioblastoma stem cells via the cis-regulating the MAP4K4 to activate JNK pathway

Gliomas are the most prevalent and aggressive brain tumors, exhibiting high proliferation, abnormal glycolysis, and poor prognosis. LncRNAs act as regulatory molecules and play crucial roles in the malignant behaviors of GBM cells, including cell self-renewal. However, the regulatory mechanisms involved are largely unknown. In this study, we performed bioinformatics analysis to explore NF-κB pathway-related lncRNAs. ECAR and qRT-PCR were used to measure the relationship between glycolytic activity and lncRNA expression. Assays such as RIP-PCR and ChIP-PCR were employed to reveal the regulatory mechanisms of the lncRNA. Neurosphere formation and limiting dilution assays were performed to evaluate the self-renewal capacity of GBM cells. In our study, we identified an NF-κB pathway-related lncRNA named LINC01127 in GBM, which was found to be associated with poor progression of GBM. Functionally, the NF-κB pathway promoted warburg effect, which, in turn, induced the lactylation of H3 histone and increased the expression of LINC01127. Consequently, this enhancement of LINC01127 expression led to the promotion of self-renewal in GBM cells. Furthermore, LINC01127 regulated MAP4K4 expression in cis by directly guiding POLR2A to the MAP4K4 promoter regions, thereby leading to JNK pathway activation, and ultimately modulating the self-renewal of GBM cells. Moreover, the activated JNK pathway promoted the phosphorylation of IκBα. Overall, targeting LINC01127-mediated axis impeded orthotopic tumor growth in GBM xenografts. Taken together these results revealed that warburg effect-induced histone lactylation drives NF-κB-related LINC01127 expression, thereby promoting the self-renewal of GBM cells through the MAP4K4/JNK/NF-κB axis, and providing substantial evidence that LINC01127 might provide a novel therapeutic strategy for GBM patients.

The role of ncRNA regulatory mechanisms in diseases-case on gestational diabetes

Non-coding RNAs (ncRNAs) are a class of RNA molecules that do not have the potential to encode proteins. Meanwhile, they can occupy a significant portion of the human genome and participate in gene expression regulation through various mechanisms. Gestational diabetes mellitus (GDM) is a pathologic condition of carbohydrate intolerance that begins or is first detected during pregnancy, making it one of the most common pregnancy complications. Although the exact pathogenesis of GDM remains unclear, several recent studies have shown that ncRNAs play a crucial regulatory role in GDM. Herein, we present a comprehensive review on the multiple mechanisms of ncRNAs in GDM along with their potential role as biomarkers. In addition, we investigate the contribution of deep learning-based models in discovering disease-specific ncRNA biomarkers and elucidate the underlying mechanisms of ncRNA. This might assist community-wide efforts to obtain insights into the regulatory mechanisms of ncRNAs in disease and guide a novel approach for early diagnosis and treatment of disease.

DNABERT-based explainable lncRNA identification in plant genome assemblies

Long non-coding ribonucleic acids (lncRNAs) have been shown to play an important role in plant gene regulation, involving both epigenetic and transcript regulation. LncRNAs are transcripts longer than 200 nucleotides that are not translated into functional proteins but can be translated into small peptides. Machine learning models have predominantly used transcriptome data with manually defined features to detect lncRNAs, however, they often underrepresent the abundance of lncRNAs and can be biased in their detection. Here we present a study using Natural Language Processing (NLP) models to identify plant lncRNAs from genomic sequences rather than transcriptomic data. The NLP models were trained to predict lncRNAs for seven model and crop species (Zea mays, Arabidopsis thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Glycine max and Oryza sativa) using publicly available genomic references. We demonstrated that lncRNAs can be accurately predicted from genomic sequences with the highest accuracy of 83.4% for Z. mays and the lowest accuracy of 57.9% for B. rapa, revealing that genome assembly quality might affect the accuracy of lncRNA identification. Furthermore, we demonstrated the potential of using NLP models for cross-species prediction with an average of 63.1% accuracy using target species not previously seen by the model. As more species are incorporated into the training datasets, we expect the accuracy to increase, becoming a more reliable tool for uncovering novel lncRNAs. Finally, we show that the models can be interpreted using explainable artificial intelligence to identify motifs important to lncRNA prediction and that these motifs frequently flanked the lncRNA sequence.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

Molecular regulatory mechanism of key LncRNAs in subclinical mastitic cows with folic acid supplementation

BACKGROUND

Folic acid is a water-soluble B vitamin (B9), which is closely related to the body's immune and other metabolic pathways. The folic acid synthesized by rumen microbes has been unable to meet the needs of high-yielding dairy cows. The incidence rate of subclinical mastitis in dairy herds worldwide ranged between 25%~65% with no obvious symptoms, but it significantly causes a decrease in lactation and milk quality. Therefore, this study aims at exploring the effects of folic acid supplementation on the expression profile of lncRNAs, exploring the molecular mechanism by which lncRNAs regulate immunity in subclinical mastitic dairy cows.

RESULTS

The analysis identified a total of 4384 lncRNA transcripts. Subsequently, differentially expressed lncRNAs in the comparison of two groups (SF vs. SC, HF vs. HC) were identified to be 84 and 55 respectively. Furthermore, the weighted gene co-expression network analysis (WGCNA) and the KEGG enrichment analysis result showed that folic acid supplementation affects inflammation and immune response-related pathways. The two groups have few pathways in common. One important lncRNA MSTRG.11108.1 and its target genes (ICAM1, CCL3, CCL4, etc.) were involved in immune-related pathways. Finally, through integrated analysis of lncRNAs with GWAS data and animal QTL database, we found that differential lncRNA and its target genes could be significantly enriched in SNPs and QTLs related to somatic cell count (SCC) and mastitis, such as MSTRG.11108.1 and its target gene ICAM1, CXCL3, GRO1.

CONCLUSIONS

For subclinical mastitic cows, folic acid supplementation can significantly affect the expression of immune-related pathway genes such as ICAM1 by regulating lncRNAs MSTRG.11108.1, thereby affecting related immune phenotypes. Our findings laid a ground foundation for theoretical and practical application for feeding folic acid supplementation in subclinical mastitic cows.

LncRNAs regulate the cyclic growth and development of hair follicles in Dorper sheep

Introduction

Hair follicles in Dorper sheep are characterized by seasonal cyclic growth and development, consequently resulting in hair shedding during spring. The cyclic growth and development of hair follicles are regulated by several influencing factors such as photoperiods, hormones, age of the animal, genes, long non-coding RNAs (lncRNAs), and signaling pathways.

Methods

In the present study, skin samples of five shedding sheep (S), used as experimental animals, and three non-shedding sheep (N), used as controls, were collected at three time points (September 27, 2019; January 3, 2020; and March 17, 2020) for RNA sequencing (RNA-seq) technology. Nine different groups (S1-vs-S2, S1-vs-S3, S2-vs-S3, N1- vs-N2, N1-vs-N3, N2-vs-N3, S1-vs-N1, S2-vs-N2, and S3-vs-N3) were compared using FDR < 0.05 and log 21 FC >as thresholds to assess the differences in the expression of lncRNAs.

Results and discussion

In total, 395 differentially expressed (DE) lncRNAs were screened. Cluster heatmap analysis identified two types of expression patterns, namely, high expression during the anagen phase (A pattern) and high expression during the telogen phase (T pattern). Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the target genes were largely enriched in the Estrogen signaling pathway, PI3K-Akt signaling pathway, Fc gamma R-mediated phagocytosis, and cell adhesion molecules (CAMs), which are associated with hair follicle cyclic growth and development-related pathways. In addition, 17 pairs of lncRNAs-target genes related to hair follicle cyclic growth and development were screened, and a regulatory network was constructed. Altogether, candidate lncRNAs and their regulated target genes were screened that contributed to sheep hair follicle cyclic growth and development. We believe these findings will provide useful insights into the underlying regulatory mechanisms.

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.

RESULTS

In this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.

CONCLUSION

In summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

Differential expression profiles and functional analysis of long non-coding RNAs in calcific aortic valve disease

AIM

To evaluate the expression profile of long non-coding RNAs (lncRNAs) in calcific aortic valve disease (CAVD) and explore their potential mechanism of action.

METHODS

The gene expression profiles (GSE153555, GSE148219, GSE199718) were downloaded from the Gene Expression Omnibus (GEO) database and FastQC was run for quality control checks. After filtering and classifying candidate lncRNAs by differentially expressed genes (DEGs) and weighted co-expression networks (WGCNA) in GSE153555, we predicted the potential cis- or trans-regulatory target genes of differentially expressed lncRNAs (DELs) by using FEELnc and established the competitive endogenous RNA (ceRNA) network by miRanda, more over functional enrichment was analyzed using the ClusterProfiler package in R Bioconductor. The hub cis- or trans-regulatory genes were verified in GSE148219 and GSE199718 respectively.

RESULTS

There were 340 up-regulated lncRNAs identified in AS group compared with the control group (|log₂Fold Change| ≥ 1.0 and P_adj ≤ 0.05), and 460 down-regulated lncRNAs. Based on target gene prediction and co-expression network construction, twelve Long non-coding RNAs (CDKN2B-AS1, AC244453.2, APCDD1L-DT, SLC12A5-AS1, TGFB3, AC243829.4, MIR4435-2HG, FAM225A, BHLHE40-AS1, LINC01614, AL356417.2, LINC01150) were identified as the hub cis- or trans-regulatory genes in the pathogenesis of CAVD which were validated in GSE148219 and GSE19971. Additionally, we found that MIR4435-2HG was the top hub trans-acting lncRNA which also plays a crucial role by ceRNA pattern.

CONCLUSION

LncRNAs may play an important role in CAVD and may provide a new perspective on the pathogenesis, diagnosis, and treatment of this disease. Further studies are required to illuminate the underlying mechanisms and provide potential therapeutic targets.

LncRNA109897-JrCCR4-JrTLP1b forms a positive feedback loop to regulate walnut resistance against anthracnose caused by Colletotrichum gloeosporioides

Walnut anthracnose induced by Colletotrichum gloeosporioides is a disastrous disease that severely restricts the development of the walnut industry in China. Long non-coding RNAs (lncRNAs) are involved in adaptive responses to disease, but their roles in the regulation of walnut anthracnose resistance response are not well defined. In this study, transcriptome analysis demonstrated that a C. gloeosporioides-induced lncRNA, lncRNA109897, located upstream from the target gene JrCCR4, upregulated the expression of JrCCR4. JrCCR4 interacted with JrTLP1b and promoted its transcriptional activity. In turn, JrTLP1b induced the transcription of lncRNA109897 to promote its expression. Meanwhile, transient expression in walnut leaves and stable transformation of Arabidopsis thaliana further proved that lncRNA, JrCCR4, and JrTLP1b improve the resistance of C. gloeosporioides. Collectively, these findings provide insights into the mechanism by which the lncRNA109897-JrCCR4-JrTLP1b transcriptional cascade regulates the resistance of walnut to anthracnose.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Modulation of long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression in the liver of Beagle dogs by Toxocara canis infection

BACKGROUND

Long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) play crucial roles in regulating various physiological and pathological processes. However, the role of lncRNAs and mRNAs in mediating the liver response during Toxocara canis infection remains incompletely understood.

METHODS

In the present study, the expression profile of lncRNAs and mRNAs was investigated in the liver of Beagle dogs infected by T. canis using high-throughput RNA sequencing.

RESULTS

Compared with the control groups, 876 differentially expressed (DE) lncRNAs and 288 DEmRNAs were identified at 12 h post-infection (hpi), 906 DElncRNAs and 261 DEmRNAs were identified at 24 hpi, and 876 DElncRNAs and 302 DEmRNAs were identified at 36 days post-infection (dpi). A total of 16 DEmRNAs (e.g. dpp4, crp and gnas) were commonly identified at the three infection stages. Enrichment and co-localization analyses identified several pathways involved in immune and inflammatory responses during T. canis infection. Some novel DElncRNAs, such as LNC_015756, LNC_011050 and LNC_011052, were also associated with immune and inflammatory responses. Also, LNC_005105 and LNC_005401 were associated with the secretion of anti-inflammatory cytokines, which may play a role in the healing of liver pathology at the late stage of infection.

CONCLUSIONS

Our data provided new insight into the regulatory roles of lncRNAs and mRNAs in the pathogenesis of T. canis and improved our understanding of the contribution of lncRNAs and mRNAs to the immune and inflammatory response of the liver during T. canis infection.

Identification of long noncoding RNAs involved in plumule-vernalization of Chinese cabbage

Vernalization is a phenomenon in which plants must undergo a period of continuous low temperatures to change from the vegetative growth stage to the reproductive growth stage. Chinese cabbage is a heading vegetable, and flowering time is an essential developmental trait. Premature vernalization leads to premature bolting, which causes a loss of product value and yield. While research into vernalization has provided a wealth of information, a complete understanding of the molecular mechanism for controlling vernalization requirements has not yet been elucidated. In this study, using high-throughput RNA sequencing, we analyzed the plumule-vernalization response of mRNA and long noncoding RNA in the bolting-resistant Chinese cabbage double haploid (DH) line 'Ju Hongxin' (JHX). A total of 3382 lncRNAs were identified, of which 1553 differentially expressed (DE) lncRNAs were characterized as plumule-vernalization responses. The ceRNA network revealed that 280 ceRNA pairs participated in the plumule-vernalization reaction of Chinese cabbage. Through identifying DE lncRNAs in Chinese cabbage and analyzing anti-, cis-, and trans-functional analysis, some candidate lncRNAs related to vernalization promoting flowering of Chinese cabbage and their regulated mRNA genes were found. Moreover, the expression of several critical lncRNAs and their targets was verified using qRT-PCR. Furthermore, we identified the candidate plumule-vernalization-related long noncoding RNAs that regulate BrFLCs in Chinese cabbage, which was interesting and different from previous studies and was a new discovery. Our findings expand the knowledge of lncRNAs in the vernalization of Chinese cabbage, and the identified lncRNAs provide rich resources for future comparative and functional studies.

Long Non-Coding RNA BNIP3 Inhibited the Proliferation of Bovine Intramuscular Preadipocytes via Cell Cycle

The intramuscular fat (or marbling fat) content is an essential economic trait of beef cattle and improves the flavor and palatability of meat. Several studies have highlighted the correlation between long non-coding RNAs (lncRNAs) and intramuscular fat development; however, the precise molecular mechanism remains unknown. Previously, through a high-throughput sequencing analysis, we found a lncRNA and named it a long non-coding RNA BNIP3 (lncBNIP3). The 5' RACE and 3' RACE explored 1945 bp total length of lncBNIP3, including 1621 bp of 5'RACE, and 464 bp of 3'RACE. The nucleoplasmic separation and FISH results explored the nuclear localization of lncBNIP3. Moreover, the tissue expression of lncBNIP3 was higher in the longissimus dorsi muscle, followed by intramuscular fat. Furthermore, down-regulation of lncBNIP3 increased the 5-Ethynyl-2'- deoxyuridine (EdU)-EdU-positive cells. The flow cytometry results showed that the number of cells in the S phase was significantly higher in preadipocytes transfected with si-lncBNIP3 than in the control group (si-NC). Similarly, CCK8 results showed that the number of cells after transfection of si-lncBNIP3 was significantly higher than in the control group. In addition, the mRNA expressions of proliferative marker genes CyclinB1 (CCNB1) and Proliferating Cell Nuclear Antigen (PCNA) in the si-lncBNIP3 group were significantly higher than in the control group. The Western Blot (WB) results also showed that the protein expression level of PCNA transfection of si-lncBNIP3 was significantly higher than in the control group. Similarly, the enrichment of lncBNIP3 significantly decreased the EdU-positive cells in the bovine preadipocytes. The results of flow cytometry and CCK8 assay also showed that overexpression of lncBNIP3 inhibited the proliferation of bovine preadipocytes. In addition, the overexpression of lncBNIP3 significantly inhibited the mRNA expressions of CCNB1 and PCNA. The WB results showed that the overexpression of lncBNIP3 significantly inhibited the expression of the CCNB1 protein level. To further explore the mechanism of lncBNIP3 on the proliferation of intramuscular preadipocytes, RNA-seq was performed after interference with si-lncBNIP3, and 660 differentially expressed genes (DEGs) were found, including 417 up-regulated DEGs and 243 down-regulated DEGs. The KEGG pathway analysis showed that the cell cycle was the most significant pathway for the functional enrichment of DEGs, followed by the DNA replication pathway. The RT-qPCR quantified the expression of twenty DEGs in the cell cycle. Therefore, we speculated that lncBNIP3 regulated intramuscular preadipocyte proliferation through the cell cycle and DNA replication pathways. To further confirm this hypothesis, the cell cycle inhibitor Ara-C was used to inhibit DNA replication of the S phase in intramuscular preadipocytes. Herein, Ara-C and si-lncBNIP3 were simultaneously added to the preadipocytes, and the CCK8, flow cytometry, and EdU assays were performed. The results showed that the si-lncBNIP3 could rescue the inhibitory effect of Ara-C in the bovine preadipocyte proliferation. In addition, lncBNIP3 could bind to the promoter of cell division control protein 6 (CDC6), and down-regulation of lncBNIP3 promoted the transcription activity and the expression of CDC6. Therefore, the inhibitory effect of lncBNIP3 on cell proliferation might be understood through the cell cycle pathway and CDC6 expression. This study provided a valuable lncRNA with functional roles in intramuscular fat accumulation and revealed new strategies for improving beef quality.

Long noncoding RNA GATA2-AS1 augments endothelial Hypoxia Inducible Factor 1-α induction and regulates hypoxic signaling

Vascular endothelial cells form the inner cellular lining of blood vessels and have myriad physiologic functions including angiogenesis and response to hypoxia. We recently identified a set of endothelial cell (EC)-enriched long noncoding RNAs (lncRNAs) in differentiated human primary cell types and described the role of the STEEL lncRNA in angiogenic patterning. We sought to further understand the role of EC-enriched lncRNAs in physiologic adaptation of the vascular endothelium. In this work, we describe an abundant, cytoplasmic, and EC-enriched lncRNA, GATA2-AS1, that is divergently transcribed from the EC-enriched developmental regulator, GATA2. While GATA2-AS1 is largely co-expressed with GATA2 in ECs, GATA2-AS1 and GATA2 appear to be complementary rather than synergistic as they have mostly distinct target genes. Common single nucleotide variants (SNVs) in GATA2-AS1 exons are associated with early onset coronary artery disease (CAD) and decreased expression of GATA2-AS1 in endothelial cell lines. In most cells, HIF1-α is central to the transcriptional response to hypoxia, while in ECs, both HIF1-α and HIF2-α are required to coordinate an acute and chronic response respectively. In this setting, GATA2-AS1 contributes to the "HIF switch" and augments HIF1-α induction in acute hypoxia to regulate HIF1-α/ HIF2-α balance. In hypoxia, GATA2-AS1 orchestrates HIF1-α-dependent induction of the glycolytic pathway, and HIF1-α-independent maintenance of mitochondrial biogenesis. Similarly, GATA2-AS1 coordinates both metabolism and "tip/stalk" cell signaling to regulate angiogenesis in hypoxic ECs. Furthermore, we find that GATA2-AS1 expression patterns are perturbed in atherosclerotic disease. Together, these results define a role for GATA2-AS1 in the EC-specific response to hypoxia.

Genome-wide analysis reveals the crucial role of lncRNAs in regulating the expression of genes controlling pollen development

The genomic location and stage-specific expression pattern of many long non-coding RNAs reveal their critical role in regulating protein-coding genes crucial in pollen developmental progression and male germ line specification. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with no apparent protein-coding potential. Multiple investigations have revealed high expression of lncRNAs in plant reproductive organs in a cell and tissue-specific manner. However, their potential role as essential regulators of molecular processes involved in sexual reproduction remains largely unexplored. We have used developing field mustard (Brassica rapa) pollen as a model system for investigating the potential role of lncRNAs in reproductive development. Reference-based transcriptome assembly performed to update the existing genome annotation identified novel expressed protein-coding genes and long non-coding RNAs (lncRNAs), including 4347 long intergenic non-coding RNAs (lincRNAs, 1058 expressed) and 2,045 lncRNAs overlapping protein-coding genes on the opposite strand (lncNATs, 780 expressed). The analysis of expression profiles reveals that lncRNAs are significant and stage-specific contributors to the gene expression profile of developing pollen. Gene co-expression networks accompanied by genome location analysis identified 38 cis-acting lincRNA, 31 cis-acting lncNAT, 7 trans-acting lincRNA and 14 trans-acting lncNAT to be substantially co-expressed with target protein-coding genes involved in biological processes regulating pollen development and male lineage specification. These findings provide a foundation for future research aiming at developing strategies to employ lncRNAs as regulatory tools for gene expression control during reproductive development.

Behind the scenes: How RNA orchestrates the epigenetic regulation of gene expression

Non-coding DNA accounts for approximately 98.5% of the human genome. Once labeled as "junk DNA", this portion of the genome has undergone a progressive re-evaluation and it is now clear that some of its transcriptional products, belonging to the non-coding RNAs (ncRNAs), are key players in cell regulatory networks. A growing body of evidence demonstrates the crucial impact of regulatory ncRNAs on mammalian gene expression. Here, we focus on the defined relationship between chromatin-interacting RNAs, particularly long non-coding RNA (lncRNA), enhancer RNA (eRNA), non-coding natural antisense transcript (ncNAT), and circular RNA (circRNA) and epigenome, a common ground where both protein and RNA species converge to regulate cellular functions. Through several examples, this review provides an overview of the variety of targets, interactors, and mechanisms involved in the RNA-mediated modulation of loci-specific epigenetic states, a fundamental evolutive strategy to orchestrate mammalian gene expression in a timely and reversible manner. We will discuss how RNA-mediated epigenetic regulation impacts development and tissue homeostasis and how its alteration contributes to the onset and progression of many different human diseases, particularly cancer.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Long non-coding RNA ZMIZ1-AS1 promotes osteosarcoma progression by stabilization of ZMIZ1

BACKGROUND

Osteosarcomas (OS) are frequent primary sarcomas of the bone in children and adolescents. The long non-coding RNAs (lncRNAs) can affect the progression of many cancers by their sense transcripts. The present study was designed to probe the role of ZMIZ1-AS1 and the downstream pathway in OS progression.

METHODS

Cell proliferation, invasion, and migration were detected by colony formation, transwell, and wound healing assays. The binding of SOX2 or MYC protein with ZMIZ1-AS1 promoter was explored by ChIP assay and dual-luciferase reporter assay. Interaction between PTBP1 protein and ZMIZ1-AS1 (or ZMIZ1 mRNA) was detected by RIP assay.

RESULTS

SOX2 and MYC are the downstream effectors of the Hippo pathway and transcriptionally activated ZMIZ1-AS1. Compared to the controls, OS tissues and cells contained higher ZMIZ1-AS1 expression. Silencing of ZMIZ1-AS1 repressed OS cell viability, proliferation, migration, and invasion. Our findings further showed that ZMIZ1-AS1 recruits RNA-binding protein PTBP1 to stabilize ZMIZ1 mRNA. PTBP1 or ZMIZ1 overexpression rescues the suppressive effects of silenced ZMIZ1-AS1 on OS cellular processes. Importantly, ZMIZ1-AS1 promotes OS growth in vivo by stabilization of ZMIZ1.

CONCLUSIONS

Long non-coding RNA ZMIZ1-AS1 promotes OS progression by stabilization of ZMIZ1. The Hippo pathway is inactivated in osteosarcoma. Transcriptional factors SOX2 and MYC downstream the Hippo pathway induce the upregulation of ZMIZ1-AS1 in osteosarcoma. ZMIZ1-AS1 recruits RNA binding protein PTBP1 that stabilizes ZMIZ1, the sense transcript of ZMIZ1-AS1. ZMIZ1-AS1 promotes osteosarcoma cell viability, proliferation, migration, and invasion by ZMIZ1 in a PTBP1 dependent manner.

mRNA, lncRNA, and circRNA expression profiles in a new aortic dissection murine model induced by hypoxia and Ang II

Background and aims

Aortic dissection (AD) is a cardiovascular emergency with degeneration of the aortic media. Mounting evidence indicates obstructive sleep apnea (OSA) as an independent risk factor for AD development with unknown mechanisms. This study aims to establish a stable murine model of OSA-related AD (OSA-AD) and uncover the potential changes in gene transcripts in OSA-AD.

Materials and methods

ApoE^–/– mice were exposed to the chronic intermittent hypoxia (CIH) system combined with Ang II administration to establish the OSA-AD model. Pathological staining was performed to exhibit the physiological structure of the mouse aorta. The SBC mouse ceRNA microarray was used to identify significantly differentially expressed (DE) mRNAs, DE long-non-coding RNAs (DElncRNAs), and DE circular RNAs (DEcircRNAs) in OSA-AD tissues. Subsequently, bioinformatics analysis, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genome (KEGG), and protein–protein interaction (PPI) analyses, were performed to evaluate the function of the significantly differentially expressed transcripts (DETs). The hub genes were confirmed using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

ApoE^–/– mice exposed to CIH and Ang II showed a high ratio of aortic accident (73.33%) and significant aortic diameter dilatation (1.96 ± 0.175 mm). A total of 1,742 mRNAs, 2,625 lncRNAs, and 537 circRNAs were identified as DETs (LogFC ≥ 1.5 or ≤ –1.5, P < 0.05). GO and KEGG analyses demonstrated that the differentially expressed mRNAs (DEmRNAs) were most enriched in cell proliferation, migration, apoptosis, inflammation, and hypoxia-related terms, which are closely related to aortic structural homeostasis. The PPI network contained 609 nodes and 934 connections, the hub genes were highlighted with the CytoHubba plugin and confirmed by qRT-PCR in AD tissues. KEGG pathway analysis revealed that the cis-regulated genes of DElncRNAs and circRNAs-host genes were enriched in aortic structural homeostasis-related pathways.

Conclusion

Our findings help establish a de novo OSA-AD animal model using ApoE^–/– mice. Many DEmRNAs, DElncRNAs, and DEcircRNAs were screened for the first time in OSA-AD tissues. Our findings provide useful bioinformatics data for understanding the molecular mechanism of OSA-AD and developing potential therapeutic strategies for OSA-AD.

Single-molecule real-time sequencing of the full-length transcriptome of Halophila beccarii

Ecologically, Halophila beccarii Asch. is considered as a colonizing or a pioneer seagrass species and a “tiny but mighty” seagrass species, since it may recover quickly from disturbance generally. The use of transcriptome technology can provide a better understanding of the physiological processes of seagrasses. To date, little is known about the genome and transcriptome information of H. beccarii. In this study, we used single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterize the transcriptome structure. A total of 11,773 of the 15,348 transcripts were successfully annotated in seven databases. In addition, 1573 long non-coding RNAs, 8402 simple sequence repeats and 2567 transcription factors were predicted in all the transcripts. A GO analysis showed that 5843 transcripts were divided into three categories, including biological process (BP), cellular component (CC) and molecular function (MF). In these three categories, metabolic process (1603 transcripts), protein-containing complex (515 transcripts) and binding (3233 transcripts) were the primary terms in BP, CC, and MF, respectively. The major types of transcription factors were involved in MYB-related and NF-YB families. To the best of our knowledge, this is the first report of the transcriptome of H. beccarii using SMRT sequencing technology.

Hypoxia-induced lncRNA STEAP3-AS1 activates Wnt/β-catenin signaling to promote colorectal cancer progression by preventing m6A-mediated degradation of STEAP3 mRNA

Background

Hypoxia, a typical hallmark of solid tumors, exhibits an essential role in the progression of colorectal cancer (CRC), in which the dysregulation of long non-coding RNAs (lncRNAs) is frequently observed. However, the underlying mechanisms are not clearly defined.

Methods

The TCGA database was analyzed to identify differential lncRNA expression involved in hypoxia-induced CRC progression. qRT-PCR was conducted to validate the upregulation of lncRNA STEAP3-AS1 in CRC cell lines and tumor-bearing mouse and zebrafish models under hypoxia. ChIP-qRT-PCR was used to detect the transcriptional activation of STEAP3-AS1 mediated by HIF-1α. RNA-seq, fluorescent in situ hybridization, RNA pulldown, RNA immunoprecipitation, co-immunoprecipitation, immunofluorescence and immunoblot experiments were used to ascertain the involved mechanisms. Functional assays were performed in both in vitro and in vivo models to investigate the regulatory role of STEAP3-AS1/STEAP3/Wnt/β-catenin axis in CRC proliferation and metastasis.

Results

Here, we identified a hypoxia-induced antisense lncRNA STEAP3-AS1 that was highly expressed in clinical CRC tissues and positively correlated with poor prognosis of CRC patients. Upregulation of lncRNA STEAP3-AS1, which was induced by HIF-1α-mediated transcriptional activation, facilitated the proliferation and metastasis of CRC cells both in vitro and in vivo. Mechanistically, STEAP3-AS1 interacted competitively with the YTH domain-containing family protein 2 (YTHDF2), a N⁶-methyladenosine (m⁶A) reader, leading to the disassociation of YTHDF2 with STEAP3 mRNA. This effect protected STEAP3 mRNA from m⁶A-mediated degradation, enabling the high expression of STEAP3 protein and subsequent production of cellular ferrous iron (Fe²⁺). Increased Fe²⁺ levels elevated Ser 9 phosphorylation of glycogen synthase kinase 3 beta (GSK3β) and inhibited its kinase activity, thus releasing β-catenin for nuclear translocation and subsequent activation of Wnt signaling to support CRC progression.

Conclusions

Taken together, our study highlights the mechanisms of lncRNA STEAP3-AS1 in facilitating CRC progression involving the STEAP3-AS1/STEAP3/Wnt/β-catenin axis, which may provide novel diagnostic biomarkers or therapeutic targets to benefit CRC treatment.

Graphical abstract

Hypoxia-induced HIF-1α transcriptionally upregulates the expression of lncRNA STEAP3-AS1, which interacts competitively with YTHDF2, thus upregulating mRNA stability of STEAP3 and consequent STEAP3 protein expression. The enhanced STEAP3 expression results in production of cellular ferrous iron (Fe²⁺), which induces the Ser 9 phosphorylation and inactivation of GSK3β, releasing β-catenin for nuclear translocation and contributing to subsequent activation of Wnt signaling to promote CRC progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-022-01638-1.

Novel genetic associations with five aesthetic facial traits: A genome-wide association study in the Chinese population

Background: The aesthetic facial traits are closely related to life quality and strongly influenced by genetic factors, but the genetic predispositions in the Chinese population remain poorly understood.

Methods: A genome-wide association studies (GWAS) and subsequent validations were performed in 26,806 Chinese on five facial traits: widow’s peak, unibrow, double eyelid, earlobe attachment, and freckles. Functional annotation was performed based on the expression quantitative trait loci (eQTL) variants, genome-wide polygenic scores (GPSs) were developed to represent the combined polygenic effects, and single nucleotide polymorphism (SNP) heritability was presented to evaluate the contributions of the variants.

Results: In total, 21 genetic associations were identified, of which ten were novel: GMDS-AS1 (rs4959669, p = 1.29 × 10⁻⁴⁹) and SPRED2 (rs13423753, p = 2.99 × 10⁻¹⁴) for widow’s peak, a previously unreported trait; FARSB (rs36015125, p = 1.96 × 10⁻²¹) for unibrow; KIF26B (rs7549180, p = 2.41 × 10⁻¹⁵), CASC2 (rs79852633, p = 4.78 × 10⁻¹¹⁾, RPGRIP1L (rs6499632, p = 9.15 × 10⁻¹¹), and PAX1 (rs147581439, p = 3.07 × 10⁻⁸) for double eyelid; ZFHX3 (rs74030209, p = 9.77 × 10⁻¹⁴) and LINC01107 (rs10211400, p = 6.25 × 10⁻¹⁰) for earlobe attachment; and SPATA33 (rs35415928, p = 1.08 × 10⁻⁸) for freckles. Functionally, seven identified SNPs tag the missense variants and six may function as eQTLs. The combined polygenic effect of the associations was represented by GPSs and contributions of the variants were evaluated using SNP heritability.

Conclusion: These identifications may facilitate a better understanding of the genetic basis of features in the Chinese population and hopefully inspire further genetic research on facial development.

Hypoxia-induced LncRNA DACT3-AS1 upregulates PKM2 to promote metastasis in hepatocellular carcinoma through the HDAC2/FOXA3 pathway

Growing evidence has revealed that hypoxia is involved in multiple stages of cancer development. However, there are limited reports on the effects of long noncoding RNAs (lncRNAs) on hepatocellular carcinoma (HCC) progression under hypoxia. The main purposes of this study were to analyze the effect of the novel lncRNA DACT3-AS1 on metastasis in HCC and to elucidate the related molecular mechanism. Bioinformatics tools were employed. RT–qPCR or western blot assays were conducted to detect RNA or protein expression. Clinical samples and in vivo assays were utilized to reveal the role of DACT3-AS1 in HCC. Other mechanism and functional analyses were specifically designed and performed as well. Based on the collected data, this study revealed that HIF-1α transcriptionally activates DACT3-AS1 expression under hypoxia. DACT3-AS1 was verified to promote metastasis in HCC. Mechanistically, DACT3-AS1 promotes the interaction between HDAC2 and FOXA3 to stimulate FOXA3 deacetylation, which consequently downregulates the FOXA3 protein. Furthermore, FOXA3 serves as a transcription factor that can bind to the PKM2 promoter region, thus hindering PKM2 expression. To summarize, this study uncovered that HIF-1α-induced DACT3-AS1 promotes metastasis in HCC and can upregulate PKM2 via the HDAC2/FOXA3 pathway in HCC cells.

Understanding the role of an RNA molecule involved in metastasis (spread) of liver cancer may suggest potential therapeutic targets. Hepatocarcinoma is a common primary liver cancer, and mortality remains high due to late diagnosis and the risk of metastasis. Scientists believe hypoxic (low oxygen) conditions in solid tumors may trigger metastasis by a mechanism involving long non-coding RNAs. Bin Li and co-workers at the Affiliated Hospital of Guilin Medical College, China, used patient tissue samples to examine the role of the long non-coding RNA molecule DACT3-AS1 in promoting hepatocarcinoma metastasis. Hypoxia triggers the overexpression of HIF-1α. This protein activated DACT3-AS1, which was then highly expressed in metastatic tissues. DACT3-AS1 interacted with a nearby gene and associated enzyme to promote cell migration and invasion, hinting at possible treatment options.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

SMRT sequencing of the full-length transcriptome of Gekko gecko

Tokay Gecko (Gekko gecko) is a rare and endangered medicinal animal in China. Its dry body has been used as an anti-asthmatic agent for two thousand years. To date, the genome and transcriptome of this species remain poorly understood. Here, we adopted single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterized the transcriptome structure. We identified 882,273 circular consensus (CCS) reads, including 746,317 full-length nonchimeric (FLNC) reads. The transcript cluster analysis revealed 212,964 consensus sequences, including 203,994 high-quality isoforms. In total, 111,372 of 117,888 transcripts were successfully annotated against eight databases (Nr, eggNOG, Swiss-Prot, GO, COG, KOG, Pfam and KEGG). Furthermore, 23,877 alternative splicing events, 169,128 simple sequence repeats (SSRs), 10,437 lncRNAs and 7,932 transcription factors were predicted across all transcripts. To our knowledge, this report is the first to document the G. gecko transcriptome using SMRT sequencing. The full-length transcript data might accelerate transcriptome research and lay the foundation for further research on G. gecko.

Comprehensive Analysis of Long Non-coding RNA Modulates Axillary Bud Development in Tobacco (Nicotiana tabacum L.)

Long non-coding RNAs (lncRNAs) regulate gene expression and are crucial for plant growth and development. However, the mechanisms underlying the effects of activated lncRNAs on axillary bud development remain largely unknown. By lncRNA transcriptomes of axillary buds in topped and untopped tobacco plants, we identified a total of 13,694 lncRNAs. LncRNA analysis indicated that the promoted growth of axillary bud by topping might be partially ascribed to the genes related to hormone signal transduction and glycometabolism, trans-regulated by differentially expressed lncRNAs, such as MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1. Metabolite profiling indicated that auxin, abscisic acid and gibberellin were decreased in axillary buds of topped tobacco lines, while cytokinin was increased, consistent with the expression levels of related lncRNAs. MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1 were shown to be influenced by hormones and sucrose treatments, and were associated with changes of axillary bud growth in the overexpression of NtCCD8 plants (with reduced axillary buds) and RNA interference of NtTB1 plants (with increased axillary buds). Moreover, MSTRG.28151.1 was identified as the antisense lncRNA of NtTB1. Silencing of MSTRG.28151.1 in tobacco significantly attenuated the expression of NtTB1 and resulted in larger axillary buds, suggesting the vital function of MSTRG.28151.1 axillary bud developmen by NtTB1. Our findings shed light on lncRNA-mRNA interactions and their functional roles in axillary bud growth, which would improve our understanding of lncRNAs as important regulators of axillary bud development and plant architecture.

Transcriptomic Analysis of Long Non-coding RNA-MicroRNA-mRNA Interactions in the Nucleus Accumbens Related to Morphine Addiction in Mice

Recent research suggest that some non-coding RNAs (ncRNAs) are important regulators of chromatin dynamics and gene expression in nervous system development and neurological diseases. Nevertheless, the molecular mechanisms of long non-coding RNAs (lncRNAs), acting as competing endogenous RNAs (ceRNAs), underlying morphine addiction are still unknown. In this research, RNA sequencing (RNA-seq) was used to examine the expression profiles of lncRNAs, miRNAs and mRNAs on the nucleus accumbens (NAc) tissues of mice trained with morphine or saline conditioned place preference (CPP), with differential expression of 31 lncRNAs, 393 miRNAs, and 371 mRNAs found. A ceRNA network was established for reciprocal interactions for 9 differentially expressed lncRNAs (DElncRNAs), 10 differentially expressed miRNAs (DEmiRNAs) and 12 differentially expressed mRNAs (DEmRNAs) based on predicted miRNAs shared by lncRNAs and mRNAs. KEGG pathway enrichment analyses were conducted to explore the potential functions of DEmRNAs interacting with lncRNAs in the ceRNA network. These DEmRNAs were enriched in synaptic plasticity-related pathways, including pyrimidine metabolism, ECM-receptor interaction, and focal adhesion. The correlation between the relative expression of lncRNAs, miRNAs and mRNAs was analyzed to further validate predicted ceRNA networks, and the Lnc15qD3-miR-139-3p-Lrp2 ceRNA regulatory interaction was determined. These results suggest that the comprehensive network represents a new insight into the lncRNA-mediated ceRNA regulatory mechanisms underlying morphine addiction and provide new potential diagnostic and prognostic biomarkers for morphine addiction.

USP30-AS1 contributes to mitochondrial quality control in glioblastoma cells

Glioblastoma is the most serious type of brain cancer with poor prognosis. Here, using the publicly available glioma database, we identified that USP30-AS1, an antisense lncRNA locating on the opposite strand of USP30 locus, is upregulated in human gliomas, particularly in high grade glioma. High level of USP30-AS1 is correlated with poor survival in both primary and recurrent glioma patients. USP30-AS1 regulates mitochondrial homeostasis and mitophagy in glioblastoma cells. Knockdown of USP30-AS1 decreases mitochondrial protein expression and mitochondrial mass, promotes mitochondrial uncoupler-induced mitophagy. However, USP30-AS1 does not regulate USP30 expression in a cis-regulatory manner. In summary, this study proposed that USP30-AS1 may serve as a valuable prognostic marker for gliomas. USP3-AS1 is a negative regulator of mitophagy and the regulatory effect is USP30-independent. USP30-AS1 mediated repression of mitophagy may contribute to the loss of mitochondrial homeostasis and tumor development in glioma.

Nitrogen use efficiency (NUE): elucidated mechanisms, mapped genes and gene networks in maize (Zea mays L.)

Nitrogen, the vital primary plant growth nutrient at deficit soil conditions, drastically affects the growth and yield of a crop. Over the years, excess use of inorganic nitrogenous fertilizers resulted in pollution, eutrophication and thereby demanding the reduction in the use of chemical fertilizers. Being a C₄ plant with fibrous root system and high NUE, maize can be deployed to be the best candidate for better N uptake and utilization in nitrogen deficient soils. The maize germplasm sources has enormous genetic variation for better nitrogen uptake contributing traits. Adoption of single cross maize hybrids as well as inherent property of high NUE has helped maize cultivars to achieve the highest growth rate among the cereals during last decade. Further, considering the high cost of nitrogenous fertilizers, adverse effects on soil health and environmental impact, maize improvement demands better utilization of existing genetic variation for NUE via introgression of novel allelic combinations in existing cultivars. Marker assisted breeding efforts need to be supplemented with introgression of genes/QTLs related to NUE in ruling varieties and thereby enhancing the overall productivity of maize in a sustainable manner. To achieve this, we need mapped genes and network of interacting genes and proteins to be elucidated. Identified genes may be used in screening ideal maize genotypes in terms of better physiological functionality exhibiting high NUE. Future genome editing may help in developing lines with increased productivity under low N conditions in an environment of optimum agronomic practices.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01113-z.

Epigenetic modifications in acute myeloid leukemia: The emerging role of circular RNAs (Review)

Canonical epigenetic modifications, which include histone modification, chromatin remodeling and DNA methylation, play key roles in numerous cellular processes. Epigenetics underlies how cells that posses DNA with similar sequences develop into different cell types with different functions in an organism. Earlier epigenetic research has primarily been focused at the chromatin level. However, the number of studies on epigenetic modifications of RNA, such as N1‑methyladenosine, 2'‑O‑ribosemethylation, inosine, 5‑methylcytidine, N6‑methyladenosine (m6A) and pseudouridine, has seen an increase. Circular RNAs (circRNAs), a type of RNA species that lacks a 5' cap or 3' poly(A) tail, are abundantly expressed in acute myeloid leukemia (AML) and may regulate disease progression. circRNAs possess various functions, including microRNA sponging, gene transcription regulation and RNA‑binding protein interaction. Furthermore, circRNAs are m6A methylated in other types of cancer, such as colorectal and hypopharyngeal squamous cell cancers. Therefore, the critical roles of circRNA epigenetic modifications, particularly m6A, and their possible involvement in AML are discussed in the present review. Epigenetic modification of circRNAs may become a diagnostic and therapeutic target for AML in the future.