tatDB: a database of Ago1-mediated targets of transfer RNA fragments.

tRNA-derived fragments (tRFs) are a class of emerging post-transcriptional regulators of gene expression likely binding to the transcripts of target genes. However, only a few tRFs targets have been experimentally validated, making it hard to extrapolate the functions or binding mechanisms of tRFs. The paucity of resources supporting the identification of the targets of tRFs creates a bottleneck in the fast-developing field. We have previously analyzed chimeric reads in crosslinked Argonaute1-RNA complexes to help infer the guide-target pairs and binding mechanisms of multiple tRFs based on experimental data in human HEK293 cells. To efficiently disseminate these results to the research community, we designed a web-based database tatDB (targets of tRFs DataBase) populated with close to 250 000 experimentally determined guide-target pairs with ∼23 000 tRF isoforms. tatDB has a user-friendly interface with flexible query options/filters allowing one to obtain comprehensive information on given tRFs (or targets). Modes of interactions are supported by secondary structures of potential guide-target hybrids and binding motifs, essential for understanding the targeting mechanisms of tRFs. Further, we illustrate the value of the database on an example of hypothesis-building for a tRFs potentially involved in the lifecycle of the SARS-CoV-2 virus. tatDB is freely accessible at https://grigoriev-lab.camden.rutgers.edu/tatdb.

Computational meta-analysis of ribosomal RNA fragments: potential targets and interaction mechanisms.

The most abundant cellular RNA species, ribosomal RNA (rRNA), appears to be a source of massive amounts of non-randomly generated fragments. We found rRNA fragments (rRFs) in immunoprecipitated Argonaute (Ago-IP) complexes in human and mouse cells and in small RNA sequencing datasets. In human Ago1-IP, guanine-rich rRFs were preferentially cut in single-stranded regions of mature rRNAs between pyrimidines and adenosine, and non-randomly paired with cellular transcripts in crosslinked chimeras. Numerous identical rRFs were found in the cytoplasm and nucleus in mouse Ago2-IP. We report specific interaction motifs enriched in rRF-target pairs. Locations of such motifs on rRFs were compatible with the Ago structural features and patterns of the Ago-RNA crosslinking in both species. Strikingly, many of these motifs may bind to double-stranded regions on target RNAs, suggesting a potential pathway for regulating translation by unwinding mRNAs. Occurring on either end of rRFs and matching intronic, untranslated or coding regions in targets, such interaction sites extend the concept of microRNA seed regions. Targeting both borders of certain short introns, rRFs may be involved in their biogenesis or function, facilitated by Ago. Frequently dismissed as noise, rRFs are poised to greatly enrich the known functional spectrum of small RNA regulation.

Structural variants in 3000 rice genomes.

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5' UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Impact of age-associated increase in 2'-O-methylation of miRNAs on aging and neurodegeneration in Drosophila.

MicroRNAs (miRNAs) are 20- to ∼24-nucleotide (nt) small RNAs that impact a variety of biological processes, from development to age-associated events. To study the role of miRNAs in aging, studies have profiled the levels of miRNAs with time. However, evidence suggests that miRNAs show heterogeneity in length and sequence in different biological contexts. Here, by examining the expression pattern of miRNAs by Northern blot analysis, we found that Drosophila miRNAs show distinct isoform pattern changes with age. Surprisingly, an increase of some miRNAs reflects increased 2'-O-methylation of select isoforms. Small RNA deep sequencing revealed a global increase of miRNAs loaded into Ago2, but not into Ago1, with age. Our data suggest increased loading of miRNAs into Ago2, but not Ago1, with age, indicating a mechanism for differential loading of miRNAs with age between Ago1 and Ago2. Mutations in Hen1 and Ago2, which lack 2'-O-methylation of miRNAs, result in accelerated neurodegeneration and shorter life span, suggesting a potential impact of the age-associated increase of 2'-O-methylation of small RNAs on age-associated processes. Our study highlights that miRNA 2'-O-methylation at the 3' end is modulated by differential partitioning of miRNAs between Ago1 and Ago2 with age and that this process, along with other functions of Ago2, might impact age-associated events in Drosophila.

Inferring targeting modes of Argonaute-loaded tRNA fragments.

Transfer RNA fragments (tRFs) are an emerging class of small RNA molecules derived from mature or precursor tRNAs. They are found across a wide range of organisms and tissues, in small RNA fraction or loaded to Argonaute in numbers comparable to microRNAs. Their functions and mechanisms of action are largely unknown, and results obtained on individual tRFs are often hard to generalize. Here we predicted binding mechanisms and specific target interaction sites of 26 human Argonaute-loaded tRFs of different types using large-scale meta-analyses of available experimental data. Strikingly, our findings matched all interaction sites detected in a recent experimental screen, confirming the validity of our computational approach. Such sites are primarily located on the 5' end of tRFs and often involve additional binding along the tRF length, similar to microRNAs. Indicative of multiple layers of regulation, diverse regulatory non-coding RNAs comprised a third of the tRF targets, with the rest being protein-coding transcripts. In the latter, coding sequence and 3' UTRs were the likely primary target regions, although we observed interactions of tRFs with 5' UTRs. Another novel phenomenon we report, a large number of putative interactions between tRFs and introns, is compatible with the roles of Argonaute in the nucleus. Further, observed tRF-intron binding modes suggest a mechanism of interaction of tRFs with Argonaute-dependent introns, and we predict here >20 candidate introns of this type. Taken together, these results present tRFs as regulatory molecules with a rich functional spectrum.

Prediction of recurrence and remission within 3 years in patients with Cushing disease after successful transnasal adenomectomy.

BACKGROUND: Some laboratory and clinical features are associated with a probability of recurrence after transnasal adenomectomy for Cushing disease (CD). However, there is no consensus on a set of predictors. Rules for prediction of recurrence were not proposed earlier. AIM: To develop prediction model of recurrence/remission after successful neurosurgical treatment for CD. METHODS: Retrospective single-site comparative study included 349 patients (52 men and 297 women) with a verified diagnosis of CD who underwent effective endoscopic transsphenoidal adenomectomy between 2007 and 2014. Clinical and laboratory parameters were evaluated. Laboratory tests were performed using immunochemiluminescent method. Time-to-event analysis and ROC-analysis were applied. Multivariate models were developed using logistic regression and artificial neural network (ANN). RESULTS: Postoperative cortisol and ACTH levels and their combinations cannot be used for prediction of recurrence. ANN for prediction of recurrence within 3 years after successful surgery was developed. Input variables are age, duration of the disease, MRI data on adenoma, morning postoperative levels of ACTH and cortisol, output variable is binary (recurrence/remission). Predictive value for remission is 93%, 95% CI [89%; 96%], and predictive value for recurrence is 85%, 95% CI [71%; 94%]. Web-calculator based on the model is developed and free for use. CONCLUSION: Effective method for prediction of recurrence and long-term remission within 3 years after successful endoscopic transsphenoidal adenomectomy is proposed.

Visualization of nucleotide substitutions in the (micro)transcriptome.

BACKGROUND: RNA-related applications of the next-generation sequencing (NGS) technologies require context-specific interpretations: e.g., sequence mismatches may indicate sites of RNA editing, or uneven read coverage often points to mature form of microRNA. Existing visualization tools traditionally show RNA molecules in two dimensions, with their base pairing and the resulting secondary structure. However, it is not straightforward to combine a linear NGS data display with the 2-D RNA depictions. RESULTS: We present a novel approach for interactive representation of nucleotide substitutions and modifications in the transcribed genome. With the focus on RNA secondary structure in the context of NGS data, it provides intuitive visualization of genomic environment, sequence reads, nucleotide polymorphisms and editing events integrated with the structural and functional elements of both coding and non-coding RNA molecules. Using our approach we present and discuss examples and general trends of polymorphisms and editing in the context of the secondary structure of microRNAs. As expected, most of the substitutions comprised A to G and C to T events, consistent with typical RNA editing patterns. However, we did not observe prevalence of editing in double-stranded regions of the microRNA stem-loop. We describe novel prominent editing event candidates, observed across several small RNA libraries of Drosophila melanogaster. CONCLUSIONS: In contrast to the existing general tools for NGS data visualization, the power of our approach is not only in the display of read alignments and their counts, but the integration of RNA secondary structure, sequencing depth, and rates/patterns of editing or other modifications. It provides a comprehensive picture, important for large-scale studies and detailed analyses, helping to gain insight into the intricate relationships between different events in RNA biogenesis.

Malignant Endo-suprasellar Glioma: A Case Report and a Brief Review of the Literature.

Introduction: Malignant sellar gliomas are very rare phenomena. To date, only few cases of sellar and suprasellar glioblastomas have been reported, most of which originate from the optic nerve or optic chiasm. Case Presentation: We present a 34-year-old woman with malignant endo-suprasellar glioma, originating from the pituitary stalk, which was initially classified as a macroprolactinoma. Conclusions: Although malignant sellar gliomas can mimic the clinical, endocrinological, and radiological features of pituitary macroadenomas, rapid progression without appropriate hormonal activity suggests their diagnosis. Considering the high malignant potential of sellar glioblastomas, it is important to discuss the specific features of these tumors and to investigate the possibility of differential diagnosis in the preoperative stage, which can be useful for early selection of the treatment plan.

A Pilot Study of Fluorescence-Guided Resection of Pituitary Adenomas with Chlorin e6 Photosensitizer.

Fluorescence diagnostics is one of the promising methods for intraoperative detection of brain tumor boundaries and helps in maximizing the extent of resection. This paper presents the results of a pilot study on the first use of the chlorin e6 photosensitizer and a two-channel video system for fluorescence-guided resection of pituitary adenomas. The study's clinical part involved two patients diagnosed with hormonally inactive pituitary macroadenomas and one patient with a hormonally active one. All neoplasms had different sizes and growth patterns. The data showed accumulation of chlorin e6 in tumor tissues in high concentrations: Patient 1: 2 mg/kg, Patient 2: 5 mg/kg, and Patient 3: 4 mg/kg. For Patient 1, the residual part of the tumor was not resected since it was intimately attached to the anterior genu of the internal carotid artery. For Patients 2 and 3, no regions of increased Ce6 accumulation were detected in the tumor foci after resection. Therefore, the use of the Ce6 and a two-channel video system helped to achieve a high degree of tumor resection in each case.

COVID-19 and neurodegeneration: The mitochondrial connection.

There is still a significant lack of knowledge regarding many aspects of the etiopathology and consequences of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in humans. For example, the variety of molecular mechanisms mediating this infection, and the long-term consequences of the disease remain poorly understood. It first seemed like the SARS-CoV-2 infection primarily caused a serious respiratory syndrome. However, over the last years, an increasing number of studies also pointed towards the damaging effects of this infection has on the central nervous system (CNS). In fact, evidence suggests a possible disruption of the blood-brain barrier and deleterious effects on the CNS, especially in patients who already suffer from other pathologies, such as neurodegenerative disorders. The molecular mechanisms behind these effects on the CNS could involve the dysregulation of mitochondrial physiology, a well-known early marker of neurodegeneration and a hallmark of aging. Moreover, mitochondria are involved in the activation of the inflammatory response, which has also been broadly described in the CNS in COVID-19. Here, we critically review the current bibliography regarding the presence of neurodegenerative symptoms in COVID-19 patients, with a special emphasis on the mitochondrial mechanisms of these disorders.

Large-Scale Computational Discovery of Binding Motifs in tRNA Fragments.

Accumulating evidence has suggested that tRNA-derived fragments (tRFs) could be loaded to Argonaute proteins and function as regulatory small RNAs. However, their mode of action remains largely unknown, and investigations of their binding mechanisms have been limited, revealing little more than microRNA-like seed regions in a handful of tRFs and a few targets. Here, we identified such regions of potential interaction on a larger scale, using in vivo formed hybrids of guides and targets in crosslinked chimeric reads in two orientations. We considered "forward pairs" (with guides located on the 5' ends and targets on the 3' ends of hybrids) and "reverse pairs" (opposite orientation) and compared them as independent sets of biological constructs. We observed intriguing differences between the two chimera orientations, including the paucity of tRNA halves and abundance of polyT-containing targets in forward pairs. We found a total of 197 quality-ranked motifs supported by ∼120,000 tRF-mRNA chimeras, with 103 interacting motifs common in forward and reverse pairs. By analyzing T→C conversions in human and mouse PAR-CLIP datasets, we detected Argonaute crosslinking sites in tRFs, conserved across species. We proposed a novel model connecting the formation of asymmetric pairs in two sets to the potential binding mechanisms of tRFs, involving the identified interaction motifs and crosslinking sites to Argonaute proteins. Our results suggest the way forward for further experimental elucidation of tRF-binding mechanisms.

Age-Related Argonaute Loading of Ribosomal RNA Fragments.

BACKGROUND: Accumulating evidence points to the functional roles of rRNA derived Fragments (rRFs), often considered degradation byproducts. Small RNAs, including miRNAs and tRNA-derived Fragments (tRFs), have been implicated in the aging process and we considered rRFs in this context. OBJECTIVE: We performed a computational analysis of Argonaute-loaded rRFs in Drosophila melanogaster to study rRF changes with age. We determined rRF abundance in Ago1 and Ago2 at 3 and 30 days to identify Ago1-guided and Ago2-guided fragments. We searched for putative seed sequences in rRFs based on frequent matches of sliding k-mer windows to the conserved regions of 12 Drosophila genomes. We predicted putative targets (containing matches to seeds identified in four rRFs) and studied their functional enrichments using Gene Ontology. RESULTS: We identified precise cleavage sites of distinct rRF isoforms from both nuclear and mitochondrial rRNAs. The most prominent rRF isoforms were enriched in Ago2 at 3 days and that loading strongly decreased with age. For less abundant rRFs, loading of Ago2-guided rRFs generally increased in Ago2, whereas Ago1-guided rRFs revealed diverse patterns. The distribution of seed matches in targets suggested that rRFs may bind to various conserved regions of many genes, possibly via miRNA-like seed-based mechanisms. CONCLUSION: Our observations suggest that rRFs may be functional molecules, with age-dependent Argonaute loading, comparable to that of miRNAs and tRFs. The putative rRF targets showed significant enrichment in developmental processes and biological regulation, similar to tRFs and consistent with a possible involvement of these newly identified small RNAs in the Drosophila aging.

[The rules for predicting remission in patients with Cushing disease after successful endoscopic transnasal adenomectomy].

BACKGROUND: The recurrence rate after successful transnasal adenomectomy in Cushing’s disease (CD) can reach 47%. We have previously shown that patients with ACTH levels less than 7 pg/ml recurred over 3 years 4.5 times less often than patients with higher levels of ACTH, patients with cortisol levels below 123 nmol/l — in 3.4 times less than at higher values of this hormone, however, these indicators are dissociated in 41% of cases, so it is not possible use them for prediction separately. AIM: To develop a method for managing patients after successful transnasal adenomectomy depending on prognosis. METHODS: A monocenter retrospective comparative study included 349 patients (52 men, 297 women) with a confirmed diagnosis of CD, who underwent effective endoscopic transsphenoidal adenomectomy in 2007−2014. Various combinations of postoperative morning levels of ACTH and cortisol were analyzed. RESULTS: Based on the developed forecasting methods and their best characteristics, the following rules were formulated. If postoperative morning ACTH is less than 7 pg/ml and/or postoperative morning cortisol is less than 123 nmol/l, then the patient will remain in remission for 1 year with probability of 99% (95% CI 97%–100%) and for 3 years with probability of 86% (95% CI 80%–91%). CONCLUSION: The rules for predicting remission for 1 and 3 years for patients after neurosurgical treatment for CD are proposed. These rules are based on combinations of ACTH and cortisol levels.

Circulating Plasma microRNA to Differentiate Cushing's Disease From Ectopic ACTH Syndrome.

Corticotropinomas and adrenocorticotropic hormone (ACTH)-secreting neuroendocrine tumors exhibit differential levels of some microRNAs (miRs) compared to normal tissue. Because miRs can be released from tissues into circulation, they offer promise as novel disease biomarkers. Objective: To evaluate whether miRs are differentially detected in plasma samples of patients with ACTH-dependent Cushing's syndrome (CS). Design: Case-control study. Methods: Morning fasting plasma samples were collected from 41 consecutive patients with confirmed ACTH-dependent CS and 11 healthy subjects and stored at -80°C. Twenty-one miRs previously reported to be differentially expressed in ACTH-secreting tumors vs. healthy tissue samples were quantified in plasma by qPCR. Results: Among enrolled subjects, 28 were confirmed to have Cushing's disease (CD), 13 had ectopic ACTH secretion (EAS) and 11 were healthy controls. We found statistically significant differences in the circulating levels of miR-16-5p [45.04 (95% CI 28.77-61.31) in CD vs. 5.26 (2.65-7.87) in EAS, P < 0.001; q = 0.001], miR-145-5p [0.097 (0.027-0.167) in CD vs. undetectable levels in EAS, P = 0.008; q = 0.087] and differences in miR-7g-5p [1.842 (1.283-2.400) in CD vs. 0.847 (0.187-1.507) in EAS, P = 0.02; q = 0.14]. The area under the receiver-operator (ROC) curve was 0.879 (95% CI 0.770-0.987), p < 0.001, when using miR-16-5p to distinguish between CD and EAS. Circulating levels of miR-16-5p in the healthy control group differed from that of both the CD and EAS groups. Conclusions: Plasma miR levels differ in patients with CD and EAS. In particular, miR-16-5p, miR-145-5p and miR-7g-5p are promising biomarkers for further research to differentiate ACTH-dependent CS.

ARIADNA: machine learning method for ancient DNA variant discovery.

Ancient DNA (aDNA) studies often rely on standard methods of mutation calling, optimized for high-quality contemporary DNA but not for excessive contamination, time- or environment-related damage of aDNA. In the absence of validated datasets and despite showing extreme sensitivity to aDNA quality, these methods have been used in many published studies, sometimes with additions of arbitrary filters or modifications, designed to overcome aDNA degradation and contamination problems. The general lack of best practices for aDNA mutation calling may lead to inaccurate results. To address these problems, we present ARIADNA (ARtificial Intelligence for Ancient DNA), a novel approach based on machine learning techniques, using specific aDNA characteristics as features to yield improved mutation calls. In our comparisons of variant callers across several ancient genomes, ARIADNA consistently detected higher-quality genome variants with fast runtimes, while reducing the false positive rate compared with other approaches.

Effects of active acromegaly on bone mRNA and microRNA expression patterns.

OBJECTIVE: To evaluate the response of bone to chronic long-term growth hormone (GH) and insulin-like growth factor-1 (IGF1) excess by measuring the expression of selected mRNA and microRNA (miR) in bone tissue samples of patients with active acromegaly. DESIGN: Case-control study. METHODS: Bone tissue samples were obtained during transsphenoidal adenomectomy from the sphenoid bone (sella turcica) from 14 patients with clinically and biochemically confirmed acromegaly and 10 patients with clinically non-functioning pituitary adenoma (NFPA) matched by sex and age. Expression of genes involved in the regulation of bone remodeling was studied using quantitative polymerase chain reaction (qPCR). RESULTS: Of the genes involved in osteoblast and osteoclast activity, only alkaline phosphatase (ALP) mRNA was 50% downregulated in patients with acromegaly. GH excess caused increased expression of the Wnt signaling antagonists (DKK1) and agonists (WNT10B) and changes in the levels of miR involved in mesenchymal stem cell commitment to chondrocytes (miR-199a-5p) or adipocytes (miR-27-5p, miR-125b-5p, miR-34a-5p, miR-188-3p) P < 0.05; q < 0.1. Relevant compensatory mechanisms were found through the changes in miR involved in osteoblastogenesis (miR-210-5p, miR-135a-5p, miR-211, miR-23a-3p, miR-204-5p), but the expression of TWIST1 was 50% downregulated and RUNX2 was unchanged. CONCLUSIONS: Acromegaly had minimal effects on tested mRNAs specific to osteoblast or osteoclast function except for downregulated ALP expression. The expressions of miR known to be involved in mesenchymal stem cell commitment and downregulated TWIST1 expression suggest acromegaly has a negative effect on osteoblastogenesis.

Phylogenetic Heatmaps Highlight Composition Biases in Sequenced Reads.

Due to advancements in sequencing technology, sequence data production is no longer a constraint in the field of microbiology and has made it possible to study uncultured microbes or whole environments using metagenomics. However, these new technologies introduce different biases in metagenomic sequencing, affecting the nucleotide distribution of resulting sequence reads. Here, we illustrate such biases using two methods. One is based on phylogenetic heatmaps (PGHMs), a novel approach for compact visualization of sequence composition differences between two groups of sequences containing the same phylogenetic groups. This method is well suited for finding noise and biases when comparing metagenomics samples. We apply PGHMs to detect noise and bias in the data produced with different DNA extraction protocols, different sequencing platforms and different experimental frameworks. In parallel, we use principal component analysis displaying different clustering of sequences from each sample to support our findings and illustrate the utility of PGHMs. We considered contributions of the read length and GC-content variation and observed that in most cases biases were generally due to the GC-content of the reads.

Lightning-fast genome variant detection with GROM.

Current human whole genome sequencing projects produce massive amounts of data, often creating significant computational challenges. Different approaches have been developed for each type of genome variant and method of its detection, necessitating users to run multiple algorithms to find variants. We present Genome Rearrangement OmniMapper (GROM), a novel comprehensive variant detection algorithm accepting aligned read files as input and finding SNVs, indels, structural variants (SVs), and copy number variants (CNVs). We show that GROM outperforms state-of-the-art methods on 7 validated benchmarks using 2 whole genome sequencing (WGS) data sets. Additionally, GROM boasts lightning-fast run times, analyzing a 50× WGS human data set (NA12878) on commonly available computer hardware in 11 minutes, more than an order of magnitude (up to 72 times) faster than tools detecting a similar range of variants. Addressing the needs of big data analysis, GROM combines in 1 algorithm SNV, indel, SV, and CNV detection, providing superior speed, sensitivity, and precision. GROM is also able to detect CNVs, SNVs, and indels in non-paired-read WGS libraries, as well as SNVs and indels in whole exome or RNA sequencing data sets.

Evolutionary adaptation revealed by comparative genome analysis of woolly mammoths and elephants.

Comparative genomics studies typically limit their focus to single nucleotide variants (SNVs) and that was the case for previous comparisons of woolly mammoth genomes. We extended the analysis to systematically identify not only SNVs but also larger structural variants (SVs) and indels and found multiple mammoth-specific deletions and duplications affecting exons or even complete genes. The most prominent SV found was an amplification of RNase L (with different copy numbers in different mammoth genomes, up to 9-fold), involved in antiviral defense and inflammasome function. This amplification was accompanied by mutations affecting several domains of the protein including the active site and produced different sets of RNase L paralogs in four mammoth genomes likely contributing to adaptations to environmental threats. In addition to immunity and defense, we found many other unique genetic changes in woolly mammoths that suggest adaptations to life in harsh Arctic conditions, including variants involving lipid metabolism, circadian rhythms, and skeletal and body features. Together, these variants paint a complex picture of evolution of the mammoth species and may be relevant in the studies of their population history and extinction.