Short RNA Sequences Research Articles

Knockdown of CD146 promotes endothelial-to-mesenchymal transition via Wnt/β-catenin pathway.

PurposeCardiac fibrosis is characterized by the excessive deposition of extracellular matrix (ECM) proteins and leads to the maladaptive changes in myocardium. Endothelial cells (ECs) undergoing mesenchymal transition contributes to the occurrence and development of cardiac fibrosis. CD146 is an adhesion molecule highly expressed in ECs. The present study was performed to explore the role of CD146 in modulating endothelial to mesenchymal transition (EndMT).MethodsC57BL/6 mice were subjected to subcutaneous implantation of osmotic minipump infused with angiotensin II (Ang Ⅱ). Adenovirus carrying CD146 short hairpin RNA (shRNA) or CD146 encoding sequence were infected into cultured human umbilical vein endothelial cells (HUVECs) followed by stimulation with Ang II or transforming growth factor-β1 (TGF-β1). Differentially expressed genes were revealed by RNA-sequencing (RNA-Seq) analysis. Gene expression was measured by quantitative real-time PCR, and protein expression and distribution were determined by Western blot and immunofluorescence staining, respectively.ResultsCD146 was predominantly expressed by ECs in normal mouse hearts. CD146 was upregulated in ECs but not fibroblasts and myocytes in hearts of Ang II-infused mice and in HUVECs stimulated with Ang Ⅱ. RNA-Seq analysis revealed the differentially expressed genes related to EndMT and Wnt/β-catenin signaling pathway. CD146 knockdown and overexpression facilitated and attenuated, respectively, EndMT induced by Ang II or TGF-β1. CD146 knockdown upregulated Wnt pathway-related genes including Wnt4, LEF1, HNF4A, FOXA1, SOX6, and CCND3, and increased the protein level and nuclear translocation of β-catenin.ConclusionsKnockdown of CD146 exerts promotional effects on EndMT via activating Wnt/β-catenin pathway and the upregulation of CD146 might play a protective role against EndMT and cardiac fibrosis.

Cellular and Molecular Profiling of Tumor Microenvironment and Early-Stage Lung Cancer.

Lung cancers are broadly divided into two categories: non-small-cell lung carcinoma (NSCLC), which accounts for 80–85% of all cancer cases, and small-cell lung carcinoma (SCLC), which covers the remaining 10–15%. Recent advances in cancer biology and genomics research have allowed an in-depth characterization of lung cancers that have revealed new therapy targets (EGFR, ALK, ROS, and KRAS mutations) and have the potential of revealing even more biomarkers for diagnostic, prognostic, and targeted therapies. A new source of biomarkers is represented by non-coding RNAs, especially microRNAs (miRNAs). MiRNAs are short non-coding RNA sequences that have essential regulatory roles in multiple cancers. Therefore, we aim to investigate the tumor microenvironment (TME) and miRNA tumor profile in a subset of 51 early-stage lung cancer samples (T1 and T2) to better understand early tumor and TME organization and molecular dysregulation. We analyzed the immunohistochemistry expression of CD4 and CD8 as markers of the main TME immune populations, E-cadherin to evaluate early-stage epithelial-to-mesenchymal transition (EMT), and p53, the main altered tumor suppressor gene in lung cancer. Starting from these 4 markers, we identified and validated 4 miRNAs that target TP53 and regulate EMT that can be further investigated as potential early-stage lung cancer biomarkers.

Detecting RNA dynamics in live mammalian cells with fluorescence lifetime‐based sensors

RNA function is closely linked with subcellular localization and interaction with other biological molecules such as other RNAs and proteins. Acute cellular perturbations and disease mutations may alter RNA dynamics and hence function. Therefore, visualizing RNAs in live cells is critically important for a mechanistic understanding. However, robust fluorescence labels for specific RNA molecules to achieve quantitative detection in the complex cellular environment for diverse RNA species are under‐developed. Typically, fluorescent RNA tags are engineered where a genetically encoded short RNA binds a small molecule that conveys fluorescence properties. These efforts aim to generate large fluorescence intensity changes of the fluorophore upon binding the RNA. Here, we use fluorescence lifetime as an orthogonal approach for a fluorescence RNA sensor. We used the Riboglow platform, an RNA tagging platform that relies on binding of a small fluorescent molecule to a short RNA sequence. First, we compared probe fluorescence signal in live cells in the presence and absence of the RNA tag. We found that fluorescence lifetime provides substantially improved contrast compared with fluorescence intensity. Second, we assessed sensitivity of monitoring subcellular RNA localization dynamics live via fluorescence lifetime. Third, we are exploring modularity of the system by varying the RNA‐ and probe component to enable multi‐parameter RNA detection. Together, we find that fluorescence lifetime offers advantages to track RNAs and quantify localization versus fluorescence intensity. We demonstrated the versatility of this system in diverse cellular examples. Our system is the first platform that exploits fluorescence lifetime for live RNA tracking in cells. We plan to expand this imaging modality further to quantitatively assess RNA dynamics.

Sirt1 attenuates astrocyte activation via modulating Dnajb1 and chaperone-mediated autophagy after closed head injury.

Our previous study indicates that Silent information regulator 1 (Sirt1) is involved in macroautophagy by upregulating light chain 3 (LC3) expression in astrocyte to exert a neuroprotective effect. Chaperon-mediated autophagy (CMA), another form of autophagy, is also upregulated after brain injury. However, little is known about the role of Sirt1 in regulation of the CMA. In the present study, an in vivo model of closed head injury (CHI) and an in vitro model of primary cortical astrocyte stimulated with interleukin-1β were employed to mimic the astrocyte activation induced by traumatic brain injury. Lentivirus carrying target complementary DNA (cDNA) or short hairpin RNA (shRNA) sequence was used to overexpress Sirt1 or knockdown DnaJ heat shock protein family member B1 (Dnajb1) (a molecular chaperone). We found that Sirt1 overexpression ameliorated neurological deficits, reduced tissue loss, and attenuated astrocyte activation after CHI, which was reversed by Dnajb1-shRNA administration. The upregulation of CMA activity induced by CHI in vivo and in vitro was inhibited after Dnajb1 knockdown. Sirt1 potently promoted CMA activity via upregulating Dnajb1 expression. Mechanically, Sirt1 could interact with Dnajb1 and modulate the deacetylation and ubiquitination of Dnajb1. These findings collectively suggest that Sirt1 plays a protective role against astrocyte activation, which may be associated with the regulation of the CMA activity via modulating the deacetylation and ubiquitination of Dnajb1 after CHI.

Clinical significance and role of CXCL16 in anti-neutrophil cytoplasmic autoantibody-associated vasculitis

Anti-neutrophil cytoplasmic autoantibody (ANCA)-associated vasculitis (AAV) is a group of potentially life-threatening autoimmune diseases. The kidney and lung are the most common and most severely affected organs. Previous studies have shown that the chemokine ligand CXCL16 and its receptor CXCR6 play an important role in kidney disease. However, whether CXCL16/CXCR6 is involved in the pathogenesis of AAV remains elusive. In this study, the levels of CXCL16 and its specific receptor CXCR6 were investigated. According to kidney outcome, patients were divided into two groups, specifically one with high CXCL16 levels and one with low CXCL16 levels, by cut-off values using receiver operating characteristic (ROC) curves. The clinical parameters and histological features were further compared between the two groups. The ability of CXCL16 to induce neutrophil chemotaxis was analysed using a Transwell migration assay in a coculture system of conditional immortalized human glomerular endothelial cells (ciGEnCs) and neutrophils. We observed that the levels of CXCL16 were significantly increased in the circulation, along with the expression in renal tissue of AAV patients compared to healthy controls (HCs). CXCR6 expression on neutrophils was significantly higher in patients with AAV than in HCs. There were positive correlations between the levels of CXCL16 and serum creatinine, IL-6, CRP, and TNF-α and negative correlations with eGFR. The serum levels of CXCL16 could act as a predictive biomarker of renal outcome in AAV. CXCL16 secretion was upregulated in ciGEnCs treated with AAV serum. CXCL16 released from ciGEnCs contributed to neutrophil migration. Furthermore, neutrophil migration was attenuated by silencing CXCL16 expression via transfection with short hairpin RNA (shRNA) sequences and lentivirus. Taken together, these data suggest that the inhibition of the CXCL16/CXCR6 axis may provide new therapeutic strategies targeting AAV.

The genome of the Paleogene relic tree Bretschneidera sinensis: insights into trade-offs in gene family evolution, demographic history, and adaptive SNPs.

Among relic species, genomic information may provide the key to inferring their long-term survival. Therefore, in this study, we investigated the genome of the Paleogene relic tree species, Bretschneidera sinensis, which is a rare endemic species within southeastern Asia. Specifically, we assembled a high-quality genome for B. sinensis using PacBio high-fidelity and high-throughput chromosome conformation capture reads and annotated it with long and short RNA sequencing reads. Using the genome, we then detected a trade-off between active and passive disease defences among the gene families. Gene families involved in salicylic acid and MAPK signalling pathways expanded as active defence mechanisms against disease, but families involved in terpene synthase activity as passive defences contracted. When inferring the long evolutionary history of B. sinensis, we detected population declines corresponding to historical climate change around the Eocene–Oligocene transition and to climatic fluctuations in the Quaternary. Additionally, based on this genome, we identified 388 single nucleotide polymorphisms (SNPs) that were likely under selection, and showed diverse functions in growth and stress responses. Among them, we further found 41 climate-associated SNPs. The genome of B. sinensis and the SNP dataset will be important resources for understanding extinction/diversification processes using comparative genomics in different lineages.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA.

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

De novo Transcriptome Assembly of Senna occidentalis Sheds Light on the Anthraquinone Biosynthesis Pathway.

Senna occidentalis is an annual leguminous herb that is rich in anthraquinones, which have various pharmacological activities. However, little is known about the genetics of S. occidentalis, particularly its anthraquinone biosynthesis pathway. To broaden our understanding of the key genes and regulatory mechanisms involved in the anthraquinone biosynthesis pathway, we used short RNA sequencing (RNA-Seq) and long-read isoform sequencing (Iso-Seq) to perform a spatial and temporal transcriptomic analysis of S. occidentalis. This generated 121,592 RNA-Seq unigenes and 38,440 Iso-Seq unigenes. Comprehensive functional annotation and classification of these datasets using public databases identified unigene sequences related to major secondary metabolite biosynthesis pathways and critical transcription factor families (bHLH, WRKY, MYB, and bZIP). A tissue-specific differential expression analysis of S. occidentalis and measurement of the amount of anthraquinones revealed that anthraquinone accumulation was related to the gene expression levels in the different tissues. In addition, the amounts and types of anthraquinones produced differ between S. occidentalis and S. tora. In conclusion, these results provide a broader understanding of the anthraquinone metabolic pathway in S. occidentalis.

MADS-box transcription factors determine the duration of temporary winter dormancy in closely related evergreen and deciduous Iris spp.

Winter dormancy (WD) is a crucial strategy for plants coping with potentially deadly environments. In recent decades, this process has been extensively studied in economically important perennial eudicots due to changing climate. However, in evergreen monocots with no chilling requirements, dormancy processes are so far a mystery. In this study, we compared the WD process in closely related evergreen (Iris japonica) and deciduous (I. tectorum) iris species across crucial developmental time points. Both iris species exhibit a 'temporary' WD process with distinct durations, and could easily resume growth under warm conditions. To decipher transcriptional changes, full-length sequencing for evergreen iris and short read RNA sequencing for deciduous iris were applied to generate respective reference transcriptomes. Combining results from a multipronged approach, SHORT VEGETATIVE PHASE and FRUITFULL (FUL) from MADS-box was associated with a dormancy- and a growth-related module, respectively. They were co-expressed with genes involved in phytohormone signaling, carbohydrate metabolism, and environmental adaptation. Also, gene expression patterns and physiological changes in the above pathways highlighted potential abscisic acid and jasmonic acid antagonism in coordinating growth and stress responses, whereas differences in carbohydrate metabolism and reactive oxygen species scavenging might lead to species-specific WD durations. Moreover, a detailed analysis of MIKCCMADS-box in irises revealed common features described in eudicots as well as possible new roles for monocots during temporary WD, such as FLOWERING LOCUS C and FUL. In essence, our results not only provide a portrait of temporary WD in perennial monocots but also offer new insights into the regulatory mechanism underlying WD in plants.

Endurance exercise training-responsive miR-19b-3p improves skeletal muscle glucose metabolism

Skeletal muscle is a highly adaptable tissue and remodels in response to exercise training. Using short RNA sequencing, we determine the miRNA profile of skeletal muscle from healthy male volunteers before and after a 14-day aerobic exercise training regime. Among the exercise training-responsive miRNAs identified, miR-19b-3p was selected for further validation. Overexpression of miR-19b-3p in human skeletal muscle cells increases insulin signaling, glucose uptake, and maximal oxygen consumption, recapitulating the adaptive response to aerobic exercise training. Overexpression of miR-19b-3p in mouse flexor digitorum brevis muscle enhances contraction-induced glucose uptake, indicating that miR-19b-3p exerts control on exercise training-induced adaptations in skeletal muscle. Potential targets of miR-19b-3p that are reduced after aerobic exercise training include KIF13A, MAPK6, RNF11, and VPS37A. Amongst these, RNF11 silencing potentiates glucose uptake in human skeletal muscle cells. Collectively, we identify miR-19b-3p as an aerobic exercise training-induced miRNA that regulates skeletal muscle glucose metabolism.

Tagging and Capturing of Lentiviral Vectors Using Short RNAs.

Lentiviral (LV) vectors have emerged as powerful tools for transgene delivery ex vivo but in vivo gene therapy applications involving LV vectors have faced a number of challenges, including the low efficiency of transgene delivery, a lack of tissue specificity, immunogenicity to both the product encoded by the transgene and the vector, and the inactivation of the vector by the human complement cascade. To mitigate these issues, several engineering approaches, involving the covalent modification of vector particles or the incorporation of specific protein domains into the vector’s envelope, have been tested. Short synthetic oligonucleotides, including aptamers bound to the surface of LV vectors, may provide a novel means with which to retarget LV vectors to specific cells and to shield these vectors from neutralization by sera. The purpose of this study was to develop strategies to tether nucleic acid sequences, including short RNA sequences, to LV vector particles in a specific and tight fashion. To bind short RNA sequences to LV vector particles, a bacteriophage lambda N protein-derived RNA binding domain (λN), fused to the measles virus hemagglutinin protein, was used. The λN protein bound RNA sequences bearing a boxB RNA hairpin. To test this approach, we used an RNA aptamer specific to the human epidermal growth factor receptor (EGFR), which was bound to LV vector particles via an RNA scaffold containing a boxB RNA motif. The results obtained confirmed that the EGFR-specific RNA aptamer bound to cells expressing EGFR and that the boxB containing the RNA scaffold was bound specifically to the λN RNA binding domain attached to the vector. These results show that LV vectors can be equipped with nucleic acid sequences to develop improved LV vectors for in vivo applications.

44 Current Challenges in miRNomics.

Mature microRNAs (miRNAs) are short RNA sequences about 18-24 nucleotide long, which provide the recognition key within RISC for the posttranscriptional regulation of target RNAs. Considering the canonical pathway, mature miRNAs are produced via a multistep process. Their transcription (pri-miRNAs) and first processing step via the microprocessor complex (pre-miRNAs) occur in the nucleus. Then they are exported into the cytosol, processed again by Dicer (dsRNA) and finally a single strand (mature miRNA) is incorporated into RISC (miRISC). The sequence of the incorporated miRNA provides the function of RNA target recognition via hybridization. Following binding of the target, the mRNA is either degraded or translation is inhibited, which ultimately leads to less protein production. Conversely, it has been shown that binding within the 5' UTR of the mRNA can lead to an increase in protein product. Regulation of homeostasis is very important for a cell; therefore, all steps in the miRNA-based regulation pathway, from transcription to the incorporation of the mature miRNA into RISC, are under tight control. While much research effort has been exerted in this area, the knowledgebase is not sufficient for accurately modelling miRNA regulation computationally. The computational prediction of miRNAs is, however, necessary because it is not feasible to investigate all possible pairs of a miRNA and its target, let alone miRNAs and their targets. We here point out open challenges important for computational modelling or for our general understanding of miRNA-based regulation and show how their investigation is beneficial. It is our hope that this collection of challenges will lead to their resolution in the near future.

Cow's milk may be delivering potentially harmful undetected cargoes to humans. Is it time to reconsider dairy recommendations?

Mammalian evolution has shaped milk into a species-specific vehicle for post-natal development, continuing what began within the mother's womb. Increased consumption of the mother's breast milk is associated with the most adequate metabolic programming and lowers the incidence of the diseases of civilization during adulthood. An abundance of short sequences of RNA, known as microRNA, exists in mammalian breast milk, enclosed within robust small extracellular vesicles known as exosomes. These microRNAs can epigenetically regulate over 60% of human genes. When cow's milk is consumed by humans, the bovine exosomes are transported through the gastrointestinal tract, detected intact in the blood stream, and taken up by target cells, where they alter protein expression. The aim of this review was to highlight the role of dairy exosomes and microRNA, and of the type of dairy product consumed, in human diseases. Given that microRNAs are involved in a vast array of physiological processes and associated with several diseases, perhaps caution should be practiced with regard to human consumption of dairy, particularly for individuals within developmentally critical time frames, such as pregnant and lactating mothers, and young children.

Differential RNA packaging into small extracellular vesicles by neurons and astrocytes

BackgroundSmall extracellular vesicles (sEVs) mediate intercellular communication by transferring RNA, proteins, and lipids to recipient cells. These cargo molecules are selectively loaded into sEVs and mirror the physiological state of the donor cells. Given that sEVs can cross the blood–brain barrier and their composition can change in neurological disorders, the molecular signatures of sEVs in circulation can be potential disease biomarkers. Characterizing the molecular composition of sEVs from different cell types is an important first step in determining which donor cells contribute to the circulating sEVs.MethodsCell culture supernatants from primary mouse cortical neurons and astrocytes were used to purify sEVs by differential ultracentrifugation and sEVs were characterized using nanoparticle tracking analysis, transmission electron microscopy and western blot. RNA sequencing was used to determine differential expression and loading patterns of miRNAs in sEVs released by primary neurons and astrocytes. Motif analysis was conducted on enriched miRNAs in sEVs and their respective donor cells.ResultsSequencing total cellular RNA, and miRNAs from sEVs isolated from culture media of postnatal mouse cortical neurons and astrocytes revealed a distinct profile between sEVs and their corresponding cells. Though the total number of detected miRNAs in astrocytes was greater than neurons, neurons expressed more sEV-associated miRNAs than astrocytes. Only 20.7% of astrocytic miRNAs were loaded into sEVs, while 41.0% of neuronal miRNAs were loaded into sEVs, suggesting differences in the cellular sorting mechanisms. We identified short RNA sequence motifs, or EXOmotifs, on the miRNAs that were differentially loaded or excluded from sEVs. A sequence motif GUAC was enriched in astrocytic sEVs. miRNAs preferably retained in neurons or astrocytes had a similar RNA motif CACACA, suggesting a cell-type-independent mechanism to maintain cellular miRNAs. mRNAs of five RNA-binding proteins associated with passive or active RNA sorting into sEVs were differentially expressed between neurons and astrocytes, one of which, major vault protein was higher in astrocytes than in neurons and detected in astrocytic sEVs.ConclusionsOur studies suggest differences in RNA sorting into sEVs. These differences in miRNA signatures can be used for determining the cellular sources of sEVs altered in neurological disorders.Emuizutjr-v5CYc2kJPkrbVideo abstract

The methyltransferase domain of the Respiratory Syncytial Virus L protein catalyzes cap N7 and 2'-O-methylation.

Respiratory syncytial virus (RSV) is a negative sense single-stranded RNA virus and one of the main causes of severe lower respiratory tract infections in infants and young children. RSV RNA replication/transcription and capping are ensured by the viral Large (L) protein. The L protein contains a polymerase domain associated with a polyribonucleotidyl transferase domain in its N-terminus, and a methyltransferase (MTase) domain followed by the C-terminal domain (CTD) enriched in basic amino acids at its C-terminus. The MTase-CTD of Mononegavirales forms a clamp to accommodate RNA that is subsequently methylated on the cap structure and depending on the virus, on internal positions. These enzymatic activities are essential for efficient viral mRNA translation into proteins, and to prevent the recognition of uncapped viral RNA by innate immunity sensors. In this work, we demonstrated that the MTase-CTD of RSV, as well as the full-length L protein in complex with phosphoprotein (P), catalyzes the N7- and 2’-O-methylation of the cap structure of a short RNA sequence that corresponds to the 5’ end of viral mRNA. Using different experimental systems, we showed that the RSV MTase-CTD methylates the cap structure with a preference for N7-methylation as first reaction. However, we did not observe cap-independent internal methylation, as recently evidenced for the Ebola virus MTase. We also found that at μM concentrations, sinefungin, a S-adenosylmethionine analogue, inhibits the MTase activity of the RSV L protein and of the MTase-CTD domain. Altogether, these results suggest that the RSV MTase domain specifically recognizes viral RNA decorated by a cap structure and catalyzes its methylation, which is required for translation and innate immune system subversion.

Classification of Precursor MicroRNAs from Different Species Based on K-mer Distance Features

MicroRNAs (miRNAs) are short RNA sequences that are actively involved in gene regulation. These regulators on the post-transcriptional level have been discovered in virtually all eukaryotic organisms. Additionally, miRNAs seem to exist in viruses and might also be produced in microbial pathogens. Initially, transcribed RNA is cleaved by Drosha, producing precursor miRNAs. We have previously shown that it is possible to distinguish between microRNA precursors of different clades by representing the sequences in a k-mer feature space. The k-mer representation considers the frequency of a k-mer in the given sequence. We further hypothesized that the relationship between k-mers (e.g., distance between k-mers) could be useful for classification. Three different distance-based features were created, tested, and compared. The three feature sets were entitled inter k-mer distance, k-mer location distance, and k-mer first–last distance. Here, we show that classification performance above 80% (depending on the evolutionary distance) is possible with a combination of distance-based and regular k-mer features. With these novel features, classification at closer evolutionary distances is better than using k-mers alone. Combining the features leads to accurate classification for larger evolutionary distances. For example, categorizing Homo sapiens versus Brassicaceae leads to an accuracy of 93%. When considering average accuracy, the novel distance-based features lead to an overall increase in effectiveness. On the contrary, secondary-structure-based features did not lead to any effective separation among clades in this study. With this line of research, we support the differentiation between true and false miRNAs detected from next-generation sequencing data, provide an additional viewpoint for confirming miRNAs when the species of origin is known, and open up a new strategy for analyzing miRNA evolution.

AN IMPROVED METHOD FOR MICRORNA-SEQUENCING ANALYSIS IDENTIFIES A SIGNATURE OF CIRCULATING MICRORNAS IN LATE-STAGE RADIOGRAPHIC KNEE OSTEOARTHRITIS

Purpose: While next generation sequencing enables comprehensive, unbiased, and absolute quantification of molecules in a biological sample, the way in which sequencing studies are designed, executed, and analyzed can dramatically impact the results. Among the promising applications of sequencing is biomarker discovery. Since there are currently no validated biomarkers for osteoarthritis (OA), sequencing represents a viable strategy for identifying new biomarkers in different disease stages. In particular, microRNAs (short noncoding RNA sequences) are excellent biomarker candidates because they are stable, easy to detect in circulation, and reflective of disease activity.

Characterization of Ocular Surface Microbial Profiles Revealed Discrepancies between Conjunctival and Corneal Microbiota.

The ocular microbiome composition has only been partially characterized. Here, we used RNA-sequencing (RNA-Seq) data to assess microbial diversity in human corneal tissue. Additionally, conjunctival swab samples were examined to characterize ocular surface microbiota. Short RNA-Seq reads, obtained from a previous transcriptome study of 50 corneal tissues, were mapped to the human reference genome GRCh38 to remove sequences of human origin. The unmapped reads were then used for taxonomic classification by comparing them with known bacterial, archaeal, and viral sequences from public databases. The components of microbial communities were identified and characterized using both conventional microbiology and polymerase chain reaction (PCR) techniques in 36 conjunctival swabs. The majority of ocular samples examined by conventional and molecular techniques showed very similar microbial taxonomic profiles, with most of the microorganisms being classified into Proteobacteria, Firmicutes, and Actinobacteria phyla. Only 50% of conjunctival samples exhibited bacterial growth. The PCR detection provided a broader overview of positive results for conjunctival materials. The RNA-Seq assessment revealed significant variability of the corneal microbial communities, including fastidious bacteria and viruses. The use of the combined techniques allowed for a comprehensive characterization of the eye microbiome’s elements, especially in aspects of microbiota diversity.

Comprehensive profiling of circular RNAs with nanopore sequencing and CIRI-long.

Reconstructing the sequence of circular RNAs (circRNAs) from short RNA sequencing reads has proved challenging given the similarity of circRNAs and their corresponding linear messenger RNAs. Previous sequencing methods were unable to achieve high-throughput detection of full-length circRNAs. Here we describe a protocol for enrichment and full-length sequencing of circRNA isoforms using nanopore technology. Circular reverse transcription and size selection achieves a 20-fold higher enrichment of circRNAs from total RNA compared to previous methods. We developed an algorithm, called circRNA identifier using long-read sequencing data (CIRI-long), to reconstruct the sequence of circRNAs. The workflow was validated with simulated data and by comparison to Illumina sequencing as well as quantitative real-time RT-PCR. We used CIRI-long to analyze adult mouse brain samples and systematically profile circRNAs, including mitochondria-derived and transcriptional read-through circRNAs. We identified a new type of intronic self-ligated circRNA that exhibits special splicing and expression patterns. Our method takes advantage of nanopore long reads and enables unbiased reconstruction of full-length circRNA sequences.

Affinity and Structural Analysis of the U1A RNA Recognition Motif with Engineered Methionines to Improve Experimental Phasing.

RNA plays a central role in all organisms and can fold into complex structures to orchestrate function. Visualization of such structures often requires crystallization, which can be a bottleneck in the structure-determination process. To promote crystallization, an RNA-recognition motif (RRM) of the U1A spliceosomal protein has been co-opted as a crystallization module. Specifically, the U1-snRNA hairpin II (hpII) single-stranded loop recognized by U1A can be transplanted into an RNA target to promote crystal contacts and to attain phase information via molecular replacement or anomalous diffraction methods using selenomethionine. Herein, we produced the F37M/F77M mutant of U1A to augment the phasing capability of this powerful crystallization module. Selenomethionine-substituted U1A(F37M/F77M) retains high affinity for hpII (KD of 59.7 ± 11.4 nM). The 2.20 Å resolution crystal structure reveals that the mutated sidechains make new S-π interactions in the hydrophobic core and are useful for single-wavelength anomalous diffraction. Crystals were also attained of U1A(F37M/F77M) in complex with a bacterial preQ1-II riboswitch. The F34M/F37M/F77M mutant was introduced similarly into a lab-evolved U1A variant (TBP6.9) that recognizes the internal bulged loop of HIV-1 TAR RNA. We envision that this short RNA sequence can be placed into non-essential duplex regions to promote crystallization and phasing of target RNAs. We show that selenomethionine-substituted TBP6.9(F34M/F37M/F77M) binds a TAR variant wherein the apical loop was replaced with a GNRA tetraloop (KD of 69.8 ± 2.9 nM), laying the groundwork for use of TBP6.9(F34M/F37M/F77M) as a crystallization module. These new tools are available to the research community.