Weak Effect of Gypsy Retrotransposon Bursts on Sonneratia alba Salt Stress Gene Expression

Transposable elements (TEs) are an important source of genetic diversity and can be co-opted for the regulation of host genes. However, to what extent the pervasive TE colonization of plant genomes has contributed to stress adaptation remains controversial. Plants inhabiting harsh environments in nature provide a unique opportunity to answer this question. We compared TE compositions and their evolutionary dynamics in the genomes of two mangrove species: the pioneer Sonneratia alba and its less salt-tolerant relative S. caseolaris. Age distribution, strength of purifying selection and the removal rate of LTR (long terminal repeat) retrotransposons were estimated. Phylogenetic analysis of LTR retrotransposons and their distribution in the genome of S. alba were surveyed. Small RNA sequencing and whole-genome bisulfite sequencing was conducted using leaves of S. alba. Expression pattern of LTR retrotransposons and their nearby genes were examined using RNA-seq data of S. alba under different salt treatments. S. alba possesses more TEs than S. caseolaris. Particularly, many more young Gypsy LTR retrotransposons have accumulated in S. alba than in S. caseolaris despite an increase in purifying selection against TE insertions. The top two most abundant Gypsy families in S. alba preferentially insert in gene-poor regions. They are under relaxed epigenetic repression, probably due to the presence of CHROMO domains in their 3′-ends. Although a considerable number of TEs in S. alba showed differential expression under salt stress, only four copies were significantly correlated with their nearby genes in expression levels. One such TE-gene pair involves Abscisic acid 8'-hydroxylase 3 functioning in abscisic acid catabolism. This study sheds light on the evolutionary dynamics and potential function of TEs in an extremophile. Our results suggest that the conclusion on co-option of TEs should be cautious even though activation of TEs by stress might be prevalent.

Research on lncRNA Related to Drought Resistance of Shanlan Upland Rice

Background Drought has become the major abiotic stress that causes losses in rice yields and consequently is one of the main environmental factors threatening food security. Long non-coding RNA (lncRNA) is known to play an important role in plant response to drought stress, while the mechanisms of competing endogenous RNA (ceRNA) in drought resistance in upland rice have been rarely reported. Results In our study, a total of 191 lncRNAs, 2115 mRNAs and 32 miRNAs (microRNAs) were found by strand-specific sequencing and small RNA sequencing to be differentially expressed in drought-stressed rice. Functional analysis of results indicate that they play important roles in hormone signal transduction, chlorophyll synthesis, protein synthesis and other pathways. Construction of a ceRNA network revealed that MSTRG.28732.3 may interact with miR171 in the chlorophyll biosynthesis pathway and affect the ability of plants to withstand drought stress by regulating Os02g0662700, Os02g0663100 and Os06g0105350. The accuracy of the regulatory network was verified by qRT-PCR. Conclusion Our results provide a theoretical basis for future studies on the potential function of lncRNA in plant drought resistance, and they provide new genetic resources for drought-resistant rice breeding.

RNA-sequencing of human post-mortem hypothalamus and nucleus accumbens identifies expression profiles associated with obesity.

Obesity, the accumulation of body fat to excess, may cause serious negative health effects, including increased risk of heart disease, type 2 diabetes, stroke and certain cancers. The biology of obesity is complex and not well understood, involving both environmental and genetic factors and affecting metabolic and endocrine mechanisms in tissues of the gut, adipose, and brain. Previous RNA sequencing studies have identified transcripts associated with obesity and body mass index in blood and fat, often using animal models, but RNA sequencing studies in human brain tissue related to obesity have not been previously undertaken. We conducted both large and small RNA sequencing of hypothalamus (207 samples) and nucleus accumbens (276 samples) from individuals defined as consistently obese (124 samples), consistently normal weight as controls (148 samples) or selected without respect to BMI and falling within neither case nor control definition (211 samples), based on longitudinal BMI measures. The samples were provided by three cohort studies with brain donation programs; the Framingham Heart Study (FHS), the Religious Orders Study (ROS) and the Rush Memory and Aging Project (MAP). For each brain region and large/small RNA sequencing set, differential expression of obesity, BMI, brain region and sex was performed. Analyses were done transcriptome-wide as well as with a priori defined sets of obesity or BMI-associated mRNAs and microRNAs (miRNAs). There are sixteen mRNAs and five microRNAs that are differentially expressed (adjusted p < 0.05) by obesity or BMI in these tissues, several of which were validated with qPCR data. The results include many that are BMI-associated, such as APOBR and CES1, as well as many associated with the immune system and some with addiction, such as the gene sets 'cytokine signaling in immune system' and 'opioid signaling'. In spite of the relatively large number of samples, our study was likely under-powered to detect other transcripts or miRNA with relevant but smaller effects.

Postnatal developmental trajectory of sex-biased gene expression in the mouse pituitary gland

The pituitary gland controls sexually dimorphic processes such as growth, pubertal onset, and lactation. However, the mechanisms underlying sex biases in pituitary gene regulation are not fully understood. To capture pituitary gene regulation dynamics during postnatal development, we ascertained gene and miRNA expression across five postnatal days that span the pubertal transition in mice. Using three prime untranslated region and small RNA sequencing, we observed over 900 instances of sex-biased gene expression, including 18 genes that were putative targets of 5 sex-biased miRNAs. In addition, by combining bulk RNA-seq with scRNA-seq pituitary data, we obtained evidence that cell-type proportion sex differences exist prior to puberty and contribute substantially to the observed sex-biased gene expression post-puberty. This work provides a resource for postnatal mouse pituitary gene regulation and highlights the importance of sex-biases in both cell-type composition and gene regulation when understanding the sexually dimorphic processes regulated by the pituitary gland.

Salivary miRNAs as non-invasive biomarkers of hepatocellular carcinoma: a pilot study

Background Improved detection of hepatocellular carcinoma (HCC) is needed, as current detection methods, such as alpha fetoprotein (AFP) and ultrasound, suffer from poor sensitivity. MicroRNAs (miRNAs) are small, non-coding RNAs that regulate many cellular functions and impact cancer development and progression. Notably, miRNAs are detectable in saliva and have shown potential as non-invasive biomarkers for a number of cancers including breast, oral, and lung cancers. Here, we present, to our knowledge, the first report of salivary miRNAs in HCC and compare these findings to patients with cirrhosis, a high-risk cohort for HCC. Methods We performed small RNA sequencing in 20 patients with HCC and 19 with cirrhosis. Eleven patients with HCC had chronic liver disease, and analyses were performed with these samples combined and stratified by the presence of chronic liver disease. P values were adjusted for multiple comparisons using a false discovery rate (FDR) approach and miRNA with FDR P < 0.05 were considered statistically significant. Differential expression of salivary miRNAs was compared to a previously published report of miRNAs in liver tissue of patients with HCC vs cirrhosis. Support vector machines and leave-one-out cross-validation were performed to determine if salivary miRNAs have predictive potential for detecting HCC. Results A total of 4,565 precursor and mature miRNAs were detected in saliva and 365 were significantly different between those with HCC compared to cirrhosis (FDR P < 0.05). Interestingly, 283 of these miRNAs were significantly downregulated in patients with HCC. Machine-learning identified a combination of 10 miRNAs and covariates that accurately classified patients with HCC (AUC = 0.87). In addition, we identified three miRNAs that were differentially expressed in HCC saliva samples and in a previously published study of miRNAs in HCC tissue compared to cirrhotic liver tissue. Conclusions This study demonstrates, for the first time, that miRNAs relevant to HCC are detectable in saliva, that salivary miRNA signatures show potential to be highly sensitive and specific non-invasive biomarkers of HCC, and that additional studies utilizing larger cohorts are needed.

CRISPR guides induce gene silencing in plants in the absence of Cas

Background RNA-targeting CRISPR-Cas can provide potential advantages over DNA editing, such as avoiding pleiotropic effects of genome editing, providing precise spatiotemporal regulation, and expanded function including antiviral immunity. Results Here, we report the use of CRISPR-Cas13 in plants to reduce both viral and endogenous RNA. Unexpectedly, we observe that crRNA designed to guide Cas13 could, in the absence of the Cas13 protein, cause substantial reduction in RNA levels as well. We demonstrate Cas13-independent guide-induced gene silencing (GIGS) in three plant species, including stable transgenic Arabidopsis. Small RNA sequencing during GIGS identifies the production of small RNA that extend beyond the crRNA expressed sequence in samples expressing multi-guide crRNA. Additionally, we demonstrate that mismatches in guide sequences at position 10 and 11 abolish GIGS. Finally, we show that GIGS is elicited by guides that lack the Cas13 direct repeat and can extend to Cas9 designed crRNA of at least 28 base pairs, indicating that GIGS can be elicited through a variety of guide designs and is not dependent on Cas13 crRNA sequences or design. Conclusions Collectively, our results suggest that GIGS utilizes endogenous RNAi machinery despite the fact that crRNA are unlike canonical triggers of RNAi such as miRNA, hairpins, or long double-stranded RNA. Given similar evidence of Cas13-independent silencing in an insect system, it is likely GIGS is active across many eukaryotes. Our results show that GIGS offers a novel and flexible approach to RNA reduction with potential benefits over existing technologies for crop improvement and functional genomics.

Analysis of Small RNA Sequencing Data in Plants

Over the past decades, next-generation sequencing (NGS) has been employed extensively for investigating the regulatory mechanisms of small RNAs. Several bioinformatics tools are available for aiding biologists to extract meaningful information from enormous amounts of data generated by NGS platforms. This chapter describes a detailed methodology for analyzing small RNA sequencing data using different open source tools. We elaborate on various steps involved in analysis, from processing the raw sequencing reads to identifying miRNAs, their targets, and differential expression studies.

The microRNA miR-132 is a key regulator of lymphatic vascular remodelling

Lymphangiogenesis (growth of new lymphatic vessels), and lymphatic remodelling more broadly, are important for disease progression in cancer, lymphedema and the pulmonary disease lymphangioleiomyomatosis. Multiple molecular pathways which signal for aspects of lymphangiogenesis are known but little is understood about their co-ordinate regulation in lymphatic endothelial cells (LECs). Small RNA molecules co-ordinately regulate complex biological processes, but knowledge about their involvement in lymphangiogenesis is limited. Here we used high-throughput small RNA sequencing of LECs to identify microRNAs (miRs) regulating lymphatic remodelling driven by the lymphangiogenic growth factors VEGF-C and VEGF-D. We identified miR-132 as up-regulated by both growth factors, and demonstrated that inhibiting miR-132 in LECs in vitro blocked cell proliferation and tube formation, key steps in lymphangiogenesis. We showed that miR-132 is expressed in human LECs in vivo in the lymphatics of human breast tumours expressing VEGF-D. Importantly, we demonstrated that inhibiting miR-132 in vivo blocked many aspects of lymphangiogenesis in mice. Finally, we identified mRNAs regulated by miR-132 in LECs, by sequencing after RNA-protein cross-linking and Argonaute immunoprecipitation, which demonstrated how miR-132 co-ordinately regulates signalling pathways in lymphangiogenesis. This study shows miR-132 is a critical regulator of lymphangiogenesis and a potential target for therapeutically manipulating lymphatic remodelling in disease.

Identification and Evaluation of Circulating Small Extracellular Vesicle Micrornas as Diagnostic Biomarkers for Patients with Indeterminate Pulmonary Nodules

Background: The identification of indeterminate pulmonary nodules (IPNs) following a low-dose computed tomography is a major challenge for early diagnosis of lung cancer. The inadequate assessment of IPNs’ malignancy risk results in a large number of unnecessary surgeries or an increased risk of cancer metastases. However, limited studies on non-invasive diagnosis of IPNs have been reported.Methods: In this study, we identified and evaluated the diagnostic value of circulating sEV miRNAs in patients with IPNs that had been newly detected using LDCT scanning and were scheduled for surgery. Out of 459 recruited patients, 109 eligible patients with IPNs were enrolled in the training cohort (n = 47) and the test cohort (n = 62). An external cohort (n=99) was used for validation. MiRNAs were extracted from plasma sEVs, and assessed using Small RNA sequencing. 490 lung adenocarcinoma samples and follow-up data were used to investigate the role of miRNAs in overall survival.Results: A circulating sEV miRNA (CirsEV-miR) model was constructed from five differentially expressed miRNAs (DEMs), showing 0.920 AUC in the training cohort (n = 47), and further identified in the test cohort (n = 62) and in an external validation cohort (n = 99). Among five DEMs of the CirsEV-miR model, miR-101-3p and miR-150-5p were significantly associated with better overall survival (p = 0.0001 and p = 0.0069). The CirsEV-miR scores were calculated, which significantly correlated with IPNs diameters (p < 0.05), and were able to discriminate between benign and malignant PNs (diameter ≤ 1 cm). The expression patterns of sEV miRNAs in the benign, adenocarcinoma in situ/minimally invasive adenocarcinoma, and invasive adenocarcinoma subgroups were found to gradually change with the increase in aggressiveness for the first time. Among all DEMs of the three subgroups, five miRNAs (miR-30c-5p, miR-30e-5p, miR-500a-3p, miR-125a-5p, and miR-99a-5p) were also significantly associated with overall survival of lung adenocarcinoma patients.Conclusions: Our results indicate that the CirsEV-miR model could help distinguish between benign and malignant PNs, providing insights into the feasibility of circulating sEV miRNAs in diagnostic biomarker development.Trial registration: Chinese Clinical Trials, ChiCTR1800019877. Registered 05 December 2018, https://www.chictr.org.cn/showproj.aspx?proj=31346.