TeloBase: a community-curated database of telomere sequences across the tree of life.

Discoveries over the recent decade have demonstrated the unexpected diversity of telomere DNA motifs in nature. However, currently available resources, 'Telomerase database' and 'Plant rDNA database', contain just fragments of all relevant literature published over decades of telomere research as they have a different primary focus and limited updates. To fill this gap, we gathered data about telomere DNA sequences from a thorough literature screen as well as by analysing publicly available NGS data, and we created TeloBase (http://cfb.ceitec.muni.cz/telobase/) as a comprehensive database of information about telomere motif diversity. TeloBase is supplemented by internal taxonomy utilizing popular on-line taxonomic resources that enables in-house data filtration and graphical visualisation of telomere DNA evolutionary dynamics in the form of heat tree plots. TeloBase avoids overreliance on administrators for future data updates by having a simple form and community-curation system for application and approval, respectively, of new telomere sequences by users, which should ensure timeliness of the database and topicality. To demonstrate TeloBase utility, we examined telomere motif diversity in species from the fungal genus Aspergillus, and discovered (TTTATTAGGG)n sequence as a putative telomere motif in the plant family Chrysobalanaceae. This was bioinformatically confirmed by analysing template regions of identified telomerase RNAs.

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction.

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

HITS-CLIP analysis of human ALKBH8 reveals interactions with fully processed substrate tRNAs and with specific noncoding RNAs.

Transfer RNAs acquire a large plethora of chemical modifications. Among those, modifications of the anticodon loop play important roles in translational fidelity and tRNA stability. Four human wobble U-containing tRNAs obtain 5-methoxycarbonylmethyluridine (mcm5U34) or 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U34), which play a role in decoding. This mark involves a cascade of enzymatic activities. The last step is mediated by alkylation repair homolog 8 (ALKBH8). In this study, we performed a transcriptome-wide analysis of the repertoire of ALKBH8 RNA targets. Using a combination of HITS-CLIP and RIP-seq analyses, we uncover ALKBH8-bound RNAs. We show that ALKBH8 targets fully processed and CCA modified tRNAs. Our analyses uncovered the previously known set of wobble U-containing tRNAs. In addition, both our approaches revealed ALKBH8 binding to several other types of noncoding RNAs, in particular C/D box snoRNAs.

Hormone-regulated expansins: Expression, localization, and cell wall biomechanics in Arabidopsis root growth.

Expansins facilitate cell expansion by mediating pH-dependent cell wall (CW) loosening. However, the role of expansins in controlling CW biomechanical properties in specific tissues and organs remains elusive. We monitored hormonal responsiveness and spatial specificity of expression and localization of expansins predicted to be the direct targets of cytokinin signaling in Arabidopsis (Arabidopsis thaliana). We found EXPANSIN1 (EXPA1) homogenously distributed throughout the CW of columella/lateral root cap, while EXPA10 and EXPA14 localized predominantly at 3-cell boundaries in the epidermis/cortex in various root zones. EXPA15 revealed cell-type-specific combination of homogenous vs. 3-cell boundaries localization. By comparing Brillouin frequency shift and AFM-measured Young's modulus, we demonstrated Brillouin light scattering (BLS) as a tool suitable for non-invasive in vivo quantitative assessment of CW viscoelasticity. Using both BLS and AFM, we showed that EXPA1 overexpression upregulated CW stiffness in the root transition zone (TZ). The dexamethasone-controlled EXPA1 overexpression induced fast changes in the transcription of numerous CW-associated genes, including several EXPAs and XYLOGLUCAN:XYLOGLUCOSYL TRANSFERASEs (XTHs), and associated with rapid pectin methylesterification determined by in situ Fourier-transform infrared spectroscopy in the root TZ. The EXPA1-induced CW remodeling is associated with the shortening of the root apical meristem, leading to root growth arrest. Based on our results, we propose that expansins control root growth by a delicate orchestration of CW biomechanical properties, possibly regulating both CW loosening and CW remodeling.

Transcriptome analysis of thermomorphogenesis in ovules and during early seed development in Brassica napus.

BACKGROUND: Plant sexual reproduction is highly sensitive to elevated ambient temperatures, impacting seed development and production. We previously phenotyped this effect on three rapeseed cultivars (DH12075, Topas DH4079, and Westar). This work describes the transcriptional response associated with the phenotypic changes induced by heat stress during early seed development in Brassica napus. RESULTS: We compared the differential transcriptional response in unfertilized ovules and seeds bearing embryos at 8-cell and globular developmental stages of the three cultivars exposed to high temperatures. We identified that all tissues and cultivars shared a common transcriptional response with the upregulation of genes linked to heat stress, protein folding and binding to heat shock proteins, and the downregulation of cell metabolism. The comparative analysis identified an enrichment for a response to reactive oxygen species (ROS) in the heat-tolerant cultivar Topas, correlating with the phenotypic changes. The highest heat-induced transcriptional response in Topas seeds was detected for genes encoding various peroxidases, temperature-induced lipocalin (TIL1), or protein SAG21/LEA5. On the contrary, the transcriptional response in the two heat-sensitive cultivars, DH12075 and Westar, was characterized by heat-induced cellular damages with the upregulation of genes involved in the photosynthesis and plant hormone signaling pathways. Particularly, the TIFY/JAZ genes involved in jasmonate signaling were induced by stress, specifically in ovules of heat-sensitive cultivars. Using a weighted gene co-expression network analysis (WGCNA), we identified key modules and hub genes involved in the heat stress response in studied tissues of either heat-tolerant or sensitive cultivars. CONCLUSIONS: Our transcriptional analysis complements a previous phenotyping analysis by characterizing the growth response to elevated temperatures during early seed development and reveals the molecular mechanisms underlying the phenotypic response. The results demonstrated that response to ROS, seed photosynthesis, and hormonal regulation might be the critical factors for stress tolerance in oilseed rape.

WALTER: an easy way to online evaluate telomere lengths from terminal restriction fragment analysis.

BACKGROUND: Telomeres, nucleoprotein structures comprising short tandem repeats and delimiting the ends of linear eukaryotic chromosomes, play an important role in the maintenance of genome stability. Therefore, the determination of the length of telomeres is of high importance for many studies. Over the last years, new methods for the analysis of the length of telomeres have been developed, including those based on PCR or analysis of NGS data. Despite that, terminal restriction fragment (TRF) method remains the gold standard to this day. However, this method lacks universally accepted and precise tool capable to analyse and statistically evaluate TRF results. RESULTS: To standardize the processing of TRF results, we have developed WALTER, an online toolset allowing rapid, reproducible, and user-friendly analysis including statistical evaluation of the data. Given its web-based nature, it provides an easily accessible way to analyse TRF data without any need to install additional software. CONCLUSIONS: WALTER represents a major upgrade from currently available tools for the image processing of TRF scans. This toolset enables a rapid, highly reproducible, and user-friendly evaluation of almost any TRF scan including in-house statistical evaluation of the data. WALTER platform together with user manual describing the evaluation of TRF scans in detail and presenting tips and troubleshooting, as well as test data to demo the software are available at https://www.ceitec.eu/chromatin-molecular-complexes-jiri-fajkus/rg51/tab?tabId=125#WALTER and the source code at https://github.com/mlyc93/WALTER .

GLASS: assisted and standardized assessment of gene variations from Sanger sequence trace data.

MOTIVATION: Sanger sequencing is still being employed for sequence variant detection by many laboratories, especially in a clinical setting. However, chromatogram interpretation often requires manual inspection and in some cases, considerable expertise. RESULTS: We present GLASS, a web-based Sanger sequence trace viewer, editor, aligner and variant caller, built to assist with the assessment of variations in 'curated' or user-provided genes. Critically, it produces a standardized variant output as recommended by the Human Genome Variation Society. AVAILABILITY AND IMPLEMENTATION: GLASS is freely available at http://bat.infspire.org/genomepd/glass/ with source code at https://github.com/infspiredBAT/GLASS. CONTACT: nikos.darzentas@gmail.com or malcikova.jitka@fnbrno.cz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ARResT/Interrogate: an interactive immunoprofiler for IG/TR NGS data.

Motivation: The study of immunoglobulins and T cell receptors using next-generation sequencing has finally allowed exploring immune repertoires and responses in their immense variability and complexity. Unsurprisingly, their analysis and interpretation is a highly convoluted task. Results: We thus implemented ARResT/Interrogate, a web-based, interactive application. It can organize and filter large amounts of immunogenetic data by numerous criteria, calculate several relevant statistics, and present results in the form of multiple interconnected visualizations. Availability and Implementation: ARResT/Interrogate is implemented primarily in R, and is freely available at http://bat.infspire.org/arrest/interrogate/ Contact: nikos.darzentas@gmail.com Supplementary Information: Supplementary data are available at Bioinformatics online.

Trematode Diplostomum pseudospathaceum inducing differential immune gene expression in sexual and gynogenetic gibel carp (Carassius gibelio): parasites facilitating the coexistence of two reproductive forms of the invasive species.

Introduction: Parasite-mediated selection is considered one of the potential mechanisms contributing to the coexistence of asexual-sexual complexes. Gibel carp (Carassius gibelio), an invasive fish species in Europe, often forms populations composed of gynogenetic and sexual specimens. Methods: The experimental infection was induced in gynogenetic and sexual gibel carp using eye-fluke Diplostomum pseudospathaceum (Trematoda), and the transcriptome profile of the spleen as a major immune organ in fish was analyzed to reveal the differentially expressed immunity-associated genes related to D. pseudospathaceum infection differing between gynogenetic and sexual gibel carp. Results: High parasite infection was found in gynogenetic fish when compared to genetically diverse sexuals. Although metacercariae of D. pseudospathaceum are situated in an immune-privileged organ, our results show that eye trematodes may induce a host immune response. We found differential gene expression induced by eye-fluke infection, with various impacts on gynogenetic and sexual hosts, documenting for the majority of DEGs upregulation in sexuals, and downregulation in asexuals. Differences in gene regulation between gynogenetic and sexual gibel carp were evidenced in many immunity-associated genes. GO analyses revealed the importance of genes assigned to the GO terms: immune function, the Notch signaling pathway, MAP kinase tyrosine/threonine/phosphatase activity, and chemokine receptor activity. KEGG analyses revealed the importance of the genes involved in 12 immunity-associated pathways - specifically, FoxO signaling, adipocytokine signaling, TGF-beta signaling, apoptosis, Notch signaling, C-type lectin receptor signaling, efferocytosis, intestinal immune network for IgA production, insulin signaling, virion - human immunodeficiency virus, Toll-like receptor signaling, and phosphatidylinositol signaling system. Discussion: Our study indicates the limited potential of asexual fish to cope with higher parasite infection (likely a loss of capacity to induce an effective immune response) and highlights the important role of molecular mechanisms associated with immunity for the coexistence of gynogenetic and sexual gibel carp, potentially contributing to its invasiveness.

Transcriptome Profile Analyses of Head Kidney in Roach (Rutilus rutilus), Common Bream (Abramis brama) and Their Hybrids: Does Infection by Monogenean Parasites in Freshwater Fish Reveal Differences in Fish Vigour among Parental Species and Their Hybrids?

Hybrid generations usually face either a heterosis advantage or a breakdown, that can be expressed by the level of parasite infection in hybrid hosts. Hybrids are less infected by parasites than parental species (especially F1 generations) or more infected than parental species (especially post-F1 generations). We performed the experiment with blood-feeding gill parasite Paradiplozoon homoion (Monogenea) infecting leuciscid species, Abramis brama and Rutilus rutilus, their F1 generation and two backcross generations. Backcross generations tended to be more parasitized than parental lines and the F1 generation. The number of differentially expressed genes (DEGs) was lower in F1 hybrids and higher in backcross hybrids when compared to each of the parental lines. The main groups of DEGs were shared among lines; however, A. brama and R. rutilus differed in some of the top gene ontology (GO) terms. DEG analyses revealed the role of heme binding and erythrocyte differentiation after infection by blood-feeding P. homoion. Two backcross generations shared some of the top GO terms, representing mostly downregulated genes associated with P. homoion infection. KEGG analysis revealed the importance of disease-associated pathways; the majority of them were shared by two backcross generations. Our study revealed the most pronounced DEGs associated with blood-feeding monogeneans in backcross hybrids, potentially (but not exclusively) explainable by hybrid breakdown. The lower DEGs reported in F1 hybrids being less parasitized than backcross hybrids is in line with the hybrid advantage.

Computational Modelling of Metabolic Burden and Substrate Toxicity in Escherichia coli Carrying a Synthetic Metabolic Pathway.

In our previous work, we designed and implemented a synthetic metabolic pathway for 1,2,3-trichloropropane (TCP) biodegradation in Escherichia coli. Significant effects of metabolic burden and toxicity exacerbation were observed on single cell and population levels. Deeper understanding of mechanisms underlying these effects is extremely important for metabolic engineering of efficient microbial cell factories for biotechnological processes. In this paper, we present a novel mathematical model of the pathway. The model addresses for the first time the combined effects of toxicity exacerbation and metabolic burden in the context of bacterial population growth. The model is calibrated with respect to the real data obtained with our original synthetically modified E. coli strain. Using the model, we explore the dynamics of the population growth along with the outcome of the TCP biodegradation pathway considering the toxicity exacerbation and metabolic burden. On the methodological side, we introduce a unique computational workflow utilising algorithmic methods of computer science for the particular modelling problem.