LINE-1 mRNA 3' end dynamics shape its biology and retrotransposition potential.

LINE-1 (L1) retrotransposons are mobile genetic elements that create new genomic insertions by a copy-paste mechanism involving L1 RNA/RNP intermediates. L1 encodes two ORFs, of which L1-ORF2p nicks genomic DNA and reverse transcribes L1 mRNA using the nicked DNA as a primer which base-pairs with poly(A) tail of L1 mRNA. To better understand the importance of non-templated L1 3' ends' dynamics and the interplay between L1 3' and 5' ends, we investigated the effects of genomic knock-outs and temporal knock-downs of XRN1, DCP2, and other factors. We hypothesized that in the absence of XRN1, the major 5'→3' exoribonuclease, there would be more L1 mRNA and retrotransposition. Conversely, we observed that loss of XRN1 decreased L1 retrotransposition. This occurred despite slight stabilization of L1 mRNA, but with decreased L1 RNP formation. Similarly, loss of DCP2, the catalytic subunit of the decapping complex, lowered retrotransposition despite increased steady-state levels of L1 proteins. In both XRN1 and DCP2 depletions we observed shortening of L1 3' poly(A) tails and their increased uridylation by TUT4/7. We explain the observed reduction of L1 retrotransposition by the changed qualities of non-templated L1 mRNA 3' ends demonstrating the important role of L1 3' end dynamics in L1 biology.

Detecting life signatures with RNA sequence similarity measures.

The RNA World is currently the most plausible hypothesis for explaining the origins of life on Earth. The supporting body of evidence is growing and it comes from multiple areas, including astrobiology, chemistry, biology, mathematics, and, in particular, from computer simulations. Such methods frequently assume the existence of a hypothetical species on Earth, around three billion years ago, with a base sequence probably dissimilar from any in known genomes. However, it is often hard to verify whether or not a hypothetical sequence has the characteristics of biological sequences, and is thus likely to be functional. The primary objective of the presented research was to verify the possibility of building a computational 'life probe' for determining whether a given genetic sequence is biological, and assessing the sensitivity of such probes to the signatures of life present in known biological sequences. We have proposed decision algorithms based on the normalized compression distance (NCD) and Levenshtein distance (LD). We have validated the proposed method in the context of the RNA World hypothesis using short genetic sequences shorter than the error threshold value (i.e., 100 nucleotides). We have demonstrated that both measures can be successfully used to construct life probes that are significantly better than a random decision procedure, while varying from each other when it comes to detailed characteristics. We also observed that fragments of sequences related to replication have better discriminatory power than sequences having other molecular functions. In a broader context, this shows that the signatures of life in short RNA samples can be effectively detected using relatively simple means.

New in silico approach to assessing RNA secondary structures with non-canonical base pairs.

BACKGROUND: The function of RNA is strongly dependent on its structure, so an appropriate recognition of this structure, on every level of organization, is of great importance. One particular concern is the assessment of base-base interactions, described as the secondary structure, the knowledge of which greatly facilitates an interpretation of RNA function and allows for structure analysis on the tertiary level. The RNA secondary structure can be predicted from a sequence using in silico methods often adjusted with experimental data, or assessed from 3D structure atom coordinates. Computational approaches typically consider only canonical, Watson-Crick and wobble base pairs. Handling of non-canonical interactions, important for a full description of RNA structure, is still very difficult. RESULTS: We introduce our novel approach to assessing an extended RNA secondary structure, which characterizes both canonical and non-canonical base pairs, along with their type classification. It is based on predicting the RNA 3D structure from a user-provided sequence or a secondary structure that only describes canonical base pairs, and then deriving the extended secondary structure from atom coordinates. In our example implementation, this was achieved by integrating the functionality of two fully automated, high fidelity methods in a computational pipeline: RNAComposer for the 3D RNA structure prediction and RNApdbee for base-pair annotation. CONCLUSIONS: The presented methodology ties together existing applications for RNA 3D structure prediction and base-pair annotation. The example performance, applying RNAComposer and RNApdbee, reveals better accuracy in non-canonical base pair assessment than the compared methods that directly predict RNA secondary structure.

Sorting signal targeting mRNA into hepatic extracellular vesicles.

Intercellular communication mediated by extracellular vesicles has proved to play an important role in normal and pathological scenarios. However not too much information about the sorting mechanisms involved in loading the vesicles is available. Recently, our group has characterized the mRNA content of vesicles released by hepatic cellular systems, showing that a set of transcripts was particularly enriched in the vesicles in comparison with their intracellular abundance. In the current work, based on in silico bioinformatics tools, we have mapped a novel sequence of 12 nucleotides C[TA]G[GC][AGT]G[CT]C[AT]GG[GA], which is significantly enriched in the set of mRNAs that accumulate in extracellular vesicles. By including a 3'-UTR containing this sequence in a luciferase mRNA reporter, we have shown that in a hepatic cellular system this reporter mRNA was incorporated into extracellular vesicles. This study identifies a sorting signal in mRNAs that is involved in their enrichment in EVs, within a hepatic non-tumoral cellular model.

Gut Mycobiota Dysbiosis Is Associated with Melanoma and Response to Anti-PD-1 Therapy.

Recent research indicates that gut microbiota may be vital in the advancement of melanoma. In this study, we found that melanoma patients exhibited a distinct gut mycobiota structure compared with healthy participants. Candida albicans, Candida dubliniensis, and Neurospora crassa were more abundant in samples from patients with melanoma, whereas Saccharomyces cerevisiae and Debaryomyces hansenii were less abundant. During anti-PD-1 treatment, the relative amount of Malassezia restricta and C. albicans increased. A higher level of Saccharomyces paradoxus was associated with a positive response to anti-PD-1 treatment, whereas a higher level of Tetrapisispora blattae was associated with a lack of clinical benefits. High levels of M. restricta and C. albicans, elevated serum lactate dehydrogenase, and being overweight were linked to increased risk of melanoma progression and poorer response to anti-PD-1 treatment. Thus, this study has revealed melanoma-associated mycobiome dysbiosis, characterized by altered fungal composition and fungi species associated with a higher risk of melanoma progression, identifying a role for the gut mycobiome in melanoma progression.

The emerging role of the gut mycobiome in liver diseases.

In recent years, it has become clear that gut microbiota plays a major role in the human body, both in health and disease. Because of that, the gut microbiome and its impact on human well-being are getting wider and wider attention. Studies focused on the liver are not an exception. However, the majority of the analyses are concentrated on the bacterial part of the gut microbiota, while the fungi living in the human intestines are often omitted or underappreciated. This review is focused on the gut mycobiome as an important factor that should be taken into consideration regarding liver homeostasis and its perturbations. We have collected the findings in this field and we discuss their importance. We aim to emphasize the fungal compositional changes related to liver diseases and, by that, provide novel insights into the directions of liver research and gut microbiota as a therapeutic target for liver diseases.

Host Factors Associated with Gut Mycobiome Structure.

Recent studies revealed a significant role of the gut fungal community in human health. Here, we investigated the content and variation of gut mycobiota among subjects from the European population. We explored the interplay between gut fungi and various host-related sociodemographic, lifestyle, health, and dietary factors. The study included 923 participants. Fecal DNA samples were analyzed by whole-metagenome high-throughput sequencing. Subsequently, fungi taxonomic profiles were determined and accompanied by computational and statistical analyses of the association with 53 host-related factors. Fungal communities were characterized by a high prevalence of Saccharomyces, Candida, and Sporisorium. Ten factors were found to correlate significantly with the overall mycobiota variation. Most were diet related, including the consumption of chips, meat, sodas, sweetening, processed food, and alcohol, followed by age and marital status. Differences in α- and/or ß-diversity were also reported for other factors such as body mass index (BMI), job type, autoimmunological diseases, and probiotics. Differential abundance analysis revealed fungal species that exhibited different patterns of changes under specific conditions. The human gut mycobiota is dominated by yeast, including Saccharomyces, Malassezia, and Candida. Although intervolunteer variability was high, several fungal species persisted across most samples, which may be evidence that a core gut mycobiota exists. Moreover, we showed that host-related factors such as diet, age, and marital status influence the variability of gut mycobiota. To our knowledge, this is the first large and comprehensive study of the European cohort in terms of gut mycobiota associations with such an extensive and differentiated host-related set of factors. IMPORTANCE The human gut is inhabited by many organisms, including bacteria and fungi, that may affect human health. However, research on human gut mycobiome is still rare. Moreover, the large European-based cohort study is missing. Here, we analyzed the first large European cohort in terms of gut mycobiota associations with a differentiated host-related set of factors. Our results showed that chips, meat, sodas, sweetening, processed food, beer, alcohol consumption, age, and marital status were associated with the variability of gut mycobiota. Moreover, our analysis revealed changes in abundances at the fungal species level for many investigated factors. Our results can suggest potentially valuable paths for further, narrowly focused research on gut mycobiome and its impact on human health. In the coming era of gut microbiome-based precision medicine, further research into the relationship between different mycobial structures and host-related factors may result in new preventive approaches or therapeutic procedures.

RNAmap2D - calculation, visualization and analysis of contact and distance maps for RNA and protein-RNA complex structures.

BACKGROUND: The structures of biological macromolecules provide a framework for studying their biological functions. Three-dimensional structures of proteins, nucleic acids, or their complexes, are difficult to visualize in detail on flat surfaces, and algorithms for their spatial superposition and comparison are computationally costly. Molecular structures, however, can be represented as 2D maps of interactions between the individual residues, which are easier to visualize and compare, and which can be reconverted to 3D structures with reasonable precision. There are many visualization tools for maps of protein structures, but few for nucleic acids. RESULTS: We developed RNAmap2D, a platform-independent software tool for calculation, visualization and analysis of contact and distance maps for nucleic acid molecules and their complexes with proteins or ligands. The program addresses the problem of paucity of bioinformatics tools dedicated to analyzing RNA 2D maps, given the growing number of experimentally solved RNA structures in the Protein Data Bank (PDB) repository, as well as the growing number of tools for RNA 2D and 3D structure prediction. RNAmap2D allows for calculation and analysis of contacts and distances between various classes of atoms in nucleic acid, protein, and small ligand molecules. It also discriminates between different types of base pairing and stacking. CONCLUSIONS: RNAmap2D is an easy to use method to visualize, analyze and compare structures of nucleic acid molecules and their complexes with other molecules, such as proteins or ligands and metal ions. Its special features make it a very useful tool for analysis of tertiary structures of RNAs. RNAmap2D for Windows/Linux/MacOSX is freely available for academic users at http://iimcb.genesilico.pl/rnamap2d.html.

The standardisation of the approach to metagenomic human gut analysis: from sample collection to microbiome profiling.

In recent years, the number of metagenomic studies increased significantly. Wide range of factors, including the tremendous community complexity and variability, is contributing to the challenge in reliable microbiome community profiling. Many approaches have been proposed to overcome these problems making hardly possible to compare results of different studies. The significant differences between procedures used in metagenomic research are reflected in a variation of the obtained results. This calls for the need for standardisation of the procedure, to reduce the confounding factors originating from DNA isolation, sequencing and bioinformatics analyses in order to ensure that the differences in microbiome composition are of a true biological origin. Although the best practices for metagenomics studies have been the topic of several publications and the main aim of the International Human Microbiome Standard (IHMS) project, standardisation of the procedure for generating and analysing metagenomic data is still far from being achieved. To highlight the difficulties in the standardisation of metagenomics methods, we thoroughly examined each step of the analysis of the human gut microbiome. We tested the DNA isolation procedure, preparation of NGS libraries for next-generation sequencing, and bioinformatics analysis, aimed at identifying microbial taxa. We showed that the homogenisation time is the leading factor impacting sample diversity, with the recommendation for a shorter homogenisation time (10 min). Ten minutes of homogenisation allows for better reflection of the bacteria gram-positive/gram-negative ratio, and the obtained results are the least heterogenous in terms of beta-diversity of samples microbial composition. Besides increasing the homogenisation time, we observed further potential impact of the library preparation kit on the gut microbiome profiling. Moreover, our analysis revealed that the choice of the library preparation kit influences the reproducibility of the results, which is an important factor that has to be taken into account in every experiment. In this study, a tagmentation-based kit allowed for obtaining the most reproducible results. We also considered the choice of the computational tool for determining the composition of intestinal microbiota, with Kraken2/Bracken pipeline outperforming MetaPhlAn2 in our in silico experiments. The design of an experiment and a detailed establishment of an experimental protocol may have a serious impact on determining the taxonomic profile of the intestinal microbiome community. Results of our experiment can be helpful for a wide range of studies that aim to better understand the role of the gut microbiome, as well as for clinical purposes.

Profile of Basal Cell Carcinoma Mutations and Copy Number Alterations - Focus on Gene-Associated Noncoding Variants.

Basal cell carcinoma (BCC) of the skin is the most common cancer in humans, characterized by the highest mutation rate among cancers, and is mostly driven by mutations in genes involved in the hedgehog pathway. To date, almost all BCC genetic studies have focused exclusively on protein-coding sequences; therefore, the impact of noncoding variants on the BCC genome is unrecognized. In this study, with the use of whole-exome sequencing of 27 tumor/normal pairs of BCC samples, we performed an analysis of somatic mutations in both protein-coding sequences and gene-associated noncoding regions, including 5'UTRs, 3'UTRs, and exon-adjacent intron sequences. Separately, in each region, we performed hotspot identification, mutation enrichment analysis, and cancer driver identification with OncodriveFML. Additionally, we performed a whole-genome copy number alteration analysis with GISTIC2. Of the >80,000 identified mutations, ~50% were localized in noncoding regions. The results of the analysis generally corroborated the previous findings regarding genes mutated in coding sequences, including PTCH1, TP53, and MYCN, but more importantly showed that mutations were also clustered in specific noncoding regions, including hotspots. Some of the genes specifically mutated in noncoding regions were identified as highly potent cancer drivers, of which BAD had a mutation hotspot in the 3'UTR, DHODH had a mutation hotspot in the Kozak sequence in the 5'UTR, and CHCHD2 frequently showed mutations in the 5'UTR. All of these genes are functionally implicated in cancer-related processes (e.g., apoptosis, mitochondrial metabolism, and de novo pyrimidine synthesis) or the pathogenesis of UV radiation-induced cancers. We also found that the identified BAD and CHCHD2 mutations frequently occur in melanoma but not in other cancers via The Cancer Genome Atlas analysis. Finally, we identified a frequent deletion of chr9q, encompassing PTCH1, and unreported frequent copy number gain of chr9p, encompassing the genes encoding the immune checkpoint ligands PD-L1 and PD-L2. In conclusion, this study is the first systematic analysis of coding and noncoding mutations in BCC and provides a strong basis for further analyses of the variants in BCC and cancer in general.

Prebiotic Soup Components Trapped in Montmorillonite Nanoclay Form New Molecules: Car-Parrinello Ab Initio Simulations.

The catalytic effects of complex minerals or meteorites are often mentioned as important factors for the origins of life. To assess the possible role of nanoconfinement within a catalyst consisting of montmorillonite (MMT) and the impact of local electric field on the formation efficiency of the simple hypothetical precursors of nucleic acid bases or amino acids, we performed ab initio Car-Parrinello molecular dynamics simulations. We prepared four condensed-phase systems corresponding to previously suggested prototypes of a primordial soup. We monitored possible chemical reactions occurring within gas-like bulk and MMT-confined four simulation boxes on a 20-ps time scale at 1 atm and 300 K, 400 K, and 600 K. Elevated temperatures did not affect the reactivity of the elementary components of the gas-like boxes considerably; however, the presence of the MMT nanoclay substantially increased the formation probability of new molecules. Approximately 20 different new compounds were found in boxes containing carbon monoxide or formaldehyde molecules. This observation and an analysis of the atom-atom radial distribution functions indicated that the presence of Ca2+ ions at the surface of the internal MMT cavities may be an important factor in the initial steps of the formation of complex molecules at the early stages of the Earth's history.

Simulating the origins of life: The dual role of RNA replicases as an obstacle to evolution.

Despite years of study, it is still not clear how life emerged from inanimate matter and evolved into the complex forms that we observe today. One of the most recognized hypotheses for the origins of life, the RNA World hypothesis, assumes that life was sparked by prebiotic replicating RNA chains. In this paper, we address the problems caused by the interplay between hypothetical prebiotic RNA replicases and RNA parasitic species. We consider the coexistence of parasite RNAs and RNA replicases as well as the impact of parasites on the further evolution of replicases. For these purposes, we used multi-agent modeling techniques that allow for realistic assumptions regarding the movement and spatial interactions of modeled species. The general model used in this study is based on work by Takeuchi and Hogeweg. Our results confirm that the coexistence of parasite RNAs and replicases is possible in a spatially extended system, even if we take into consideration more realistic assumptions than Takeuchi and Hogeweg. However, we also showed that the presence of trade-off that takes into the account an RNA folding process could still pose a serious obstacle to the evolution of replication. We conclude that this might be a cause for one of the greatest transitions in life that took place early in evolution-the separation of the function between DNA templates and protein enzymes, with a central role for RNA species.