Implementing FAIR data management within the German Network for Bioinformatics Infrastructure (de.NBI) exemplified by selected use cases.

This article describes some use case studies and self-assessments of FAIR status of de.NBI services to illustrate the challenges and requirements for the definition of the needs of adhering to the FAIR (findable, accessible, interoperable and reusable) data principles in a large distributed bioinformatics infrastructure. We address the challenge of heterogeneity of wet lab technologies, data, metadata, software, computational workflows and the levels of implementation and monitoring of FAIR principles within the different bioinformatics sub-disciplines joint in de.NBI. On the one hand, this broad service landscape and the excellent network of experts are a strong basis for the development of useful research data management plans. On the other hand, the large number of tools and techniques maintained by distributed teams renders FAIR compliance challenging.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

Genome sequence of the progenitor of the wheat D genome Aegilops tauschii.

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

Molecular landscape and clonal architecture of adult myelodysplastic/myeloproliferative neoplasms.

More than 90% of patients with myelodysplastic/myeloproliferative neoplasms (MDSs/MPNs) harbor somatic mutations in myeloid-related genes, but still, current diagnostic criteria do not include molecular data. We performed genome-wide sequencing techniques to characterize the mutational landscape of a large and clinically well-characterized cohort including 367 adults with MDS/MPN subtypes, including chronic myelomonocytic leukemia (CMML; n = 119), atypical chronic myeloid leukemia (aCML; n = 71), MDS/MPN with ring sideroblasts and thrombocytosis (MDS/MPN-RS-T; n = 71), and MDS/MPN unclassifiable (MDS/MPN-U; n = 106). A total of 30 genes were recurrently mutated in ≥3% of the cohort. Distribution of recurrently mutated genes and clonal architecture differed among MDS/MPN subtypes. Statistical analysis revealed significant correlations between recurrently mutated genes, as well as genotype-phenotype associations. We identified specific gene combinations that were associated with distinct MDS/MPN subtypes and that were mutually exclusive with most of the other MDSs/MPNs (eg, TET2-SRSF2 in CMML, ASXL1-SETBP1 in aCML, and SF3B1-JAK2 in MDS/MPN-RS-T). Patients with MDS/MPN-U were the most heterogeneous and displayed different molecular profiles that mimicked the ones observed in other MDS/MPN subtypes and that had an impact on the outcome of the patients. Specific gene mutations also had an impact on the outcome of the different MDS/MPN subtypes, which may be relevant for clinical decision-making. Overall, the results of this study help to elucidate the heterogeneity found in these neoplasms, which can be of use in the clinical setting of MDS/MPN.

Dark-matter matters: Discriminating subtle blood cancers using the darkest DNA.

The confluence of deep sequencing and powerful machine learning is providing an unprecedented peek at the darkest of the dark genomic matter, the non-coding genomic regions lacking any functional annotation. While deep sequencing uncovers rare tumor variants, the heterogeneity of the disease confounds the best of machine learning (ML) algorithms. Here we set out to answer if the dark-matter of the genome encompass signals that can distinguish the fine subtypes of disease that are otherwise genomically indistinguishable. We introduce a novel stochastic regularization, ReVeaL, that empowers ML to discriminate subtle cancer subtypes even from the same 'cell of origin'. Analogous to heritability, implicitly defined on whole genome, we use predictability (F1 score) definable on portions of the genome. In an effort to distinguish cancer subtypes using dark-matter DNA, we applied ReVeaL to a new WGS dataset from 727 patient samples with seven forms of hematological cancers and assessed the predictivity over several genomic regions including genic, non-dark, non-coding, non-genic, and dark. ReVeaL enabled improved discrimination of cancer subtypes for all segments of the genome. The non-genic, non-coding and dark-matter had the highest F1 scores, with dark-matter having the highest level of predictability. Based on ReVeaL's predictability of different genomic regions, dark-matter contains enough signal to significantly discriminate fine subtypes of disease. Hence, the agglomeration of rare variants, even in the hitherto unannotated and ill-understood regions of the genome, may play a substantial role in the disease etiology and deserve much more attention.

The pseudogenes of barley.

Pseudogenes have a reputation of being 'evolutionary relics' or 'junk DNA'. While they are well characterized in mammals, studies in more complex plant genomes have so far been hampered by the absence of reference genome sequences. Barley is one of the economically most important cereals and has a genome size of 5.1 Gb. With the first high-quality genome reference assembly available for a Triticeae crop, we conducted a whole-genome assessment of pseudogenes on the barley genome. We identified, characterized and classified 89 440 gene fragments and pseudogenes scattered along the chromosomes, with occasional hotspots and higher densities at the chromosome ends. Full-length pseudogenes (11 015) have preferentially retained their exon-intron structure. Retrotransposition of processed mRNAs only plays a marginal role in their creation. However, the distribution of retroposed pseudogenes reflects the Rabl configuration of barley chromosomes and thus hints at founding mechanisms. While parent genes related to the defense-response were found to be under-represented in cultivated barley, we detected several defense-related pseudogenes in wild barley accessions. The percentage of transcriptionally active pseudogenes is 7.2%, and these may potentially adopt new regulatory roles.The barley genome is rich in pseudogenes and small gene fragments mainly located towards chromosome tips or as tandemly repeated units. Our results indicate non-random duplication and pseudogenization preferences and improve our understanding of the dynamics of gene birth and death in large plant genomes and the mechanisms that lead to evolutionary innovations.

Combining RNA-seq data and homology-based gene prediction for plants, animals and fungi.

BACKGROUND: Genome annotation is of key importance in many research questions. The identification of protein-coding genes is often based on transcriptome sequencing data, ab-initio or homology-based prediction. Recently, it was demonstrated that intron position conservation improves homology-based gene prediction, and that experimental data improves ab-initio gene prediction. RESULTS: Here, we present an extension of the gene prediction program GeMoMa that utilizes amino acid sequence conservation, intron position conservation and optionally RNA-seq data for homology-based gene prediction. We show on published benchmark data for plants, animals and fungi that GeMoMa performs better than the gene prediction programs BRAKER1, MAKER2, and CodingQuarry, and purely RNA-seq-based pipelines for transcript identification. In addition, we demonstrate that using multiple reference organisms may help to further improve the performance of GeMoMa. Finally, we apply GeMoMa to four nematode species and to the recently published barley reference genome indicating that current annotations of protein-coding genes may be refined using GeMoMa predictions. CONCLUSIONS: GeMoMa might be of great utility for annotating newly sequenced genomes but also for finding homologs of a specific gene or gene family. GeMoMa has been published under GNU GPL3 and is freely available at http://www.jstacs.de/index.php/GeMoMa .

Towards a whole-genome sequence for rye (Secale cereale L.).

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Light and Plastid Signals Regulate Different Sets of Genes in the Albino Mutant Pap7-1.

Plants possessing dysfunctional plastids due to defects in pigment biosynthesis or translation are known to repress photosynthesis-associated nuclear genes via retrograde signals from the disturbed organelles toward the nucleus. These signals are thought to be essential for proper biogenesis and function of the plastid. Mutants lacking plastid-encoded RNA polymerase-associated proteins (PAPs) display a genetic arrest in eoplast-chloroplast transition leading to an albino phenotype in the light. Retrograde signaling in these mutants, therefore, could be expected to be similar as under conditions inducing plastid dysfunction. To answer this question, we performed plastome- and genomewide array analyses in the pap7-1 mutant of Arabidopsis (Arabidopsis thaliana). In parallel, we determined the potential overlap with light-regulated expression networks. To this end, we performed a comparative expression profiling approach using light- and dark-grown wild-type plants as relative control for the expression profiles obtained from light-grown pap7-1 mutants. Our data indicate a specific impact of retrograde signals on metabolism-related genes in pap7-1 mutants reflecting the starvation situation of the albino seedlings. In contrast, light regulation of PhANGs and other nuclear gene groups appears to be fully functional in this mutant, indicating that a block in chloroplast biogenesis per se does not repress expression of them as suggested by earlier studies. Only genes for light harvesting complex proteins displayed a significant repression indicating an exclusive retrograde impact on this gene family. Our results indicate that chloroplasts and arrested plastids each emit specific signals that control different target gene modules both in positive and negative manner.

Elevated dietary zinc oxide levels do not have a substantial effect on porcine reproductive and respiratory syndrome virus (PPRSV) vaccination and infection.

BACKGROUND: Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most important infectious agents for the swine industry worldwide. Zinc (Zn) salts, which are widely used as a dietary supplement in swine nutrition, have shown antiviral effects in vitro as well as in vivo. The purpose of this study was to determine the influence of dietary zinc oxide supplementation on vaccination and challenge infection with PRRSV. FINDINGS: The clinical course of PRRS and the success of vaccination with an experimental inactivated vaccine were compared between animals receiving a conventional diet (50 ppm Zn, control group) and diets supplemented with Zn oxide (ZnO) at final Zn concentrations of 150 or 2,500 ppm. Pigs receiving higher dietary Zn levels showed a tendency towards higher neutralizing antibody levels after infection, while dietary Zn levels did not substantially influence the number of antiviral IFN-gamma secreting cells (IFN-gamma-SC) or percentages of blood immune cell subsets after infection. Finally, feeding higher dietary Zn levels reduced neither clinical symptoms nor viral loads. CONCLUSIONS: Our results suggest that higher levels of dietary ZnO do not have the potential to stimulate or modulate systemic immune responses after vaccination and heterologous PRRSV infection to an extent that could improve the clinical and virological outcome.

High-dose dietary zinc oxide mitigates infection with transmissible gastroenteritis virus in piglets.

BACKGROUND: Zinc (Zn) supplementation has been shown to reduce the incidence of diarrhea and to protect animals from intestinal diseases, but the mechanisms of this protective effect against virus infection in vivo have not yet been elucidated. Transmissible gastroenteritis virus (TGEV) causes diarrhea in piglets with an age-dependent decrease of severity. RESULTS: We used 60 weaned piglets that were divided into three groups to evaluate the effect of different Zn levels added to a conventional diet (50 mg Zn/kg diet, Znlow, control group). The other groups received the diet supplemented with ZnO at final concentrations of 150 mg Zn/kg diet (Znmed), or 2,500 mg/kg diet (Znhigh). Oral challenge infection with TGEV was performed when the pigs had been fed for 1 week with the respective diet. Half of the piglets of each group were sacrificed at day 1 and 18 after challenge infection. Fecal consistency was improved and body weights increased in the Znhigh group when compared to the other groups, but no direct effect of Zn concentrations in the diet on fecal TGEV shedding and mucosal immune responses was detectable. However, in the Znhigh group, we found a prevention of villus atrophy and decreased caspase-3-mediated apoptosis of jejunal epithelium. Furthermore, pigs receiving high Zn diet showed a down-regulation of interferon (IFN)-α, oligoadenylate synthetase (OAS), Zn transporter SLC39A4 (ZIP4), but up-regulation of metallothionein-1 (MT1), as well as the Zn transporters SLC30A1 (ZnT1) and SLC30A5 (ZnT5). In addition, forskolin-induced chloride secretion and epithelial resistance were controlled at a physiological level in the Znhigh but not the other groups. Finally, in the Znhigh group, we documented an earlier and higher systemic TGEV-specific serum antibody response. CONCLUSIONS: These results suggest that high dietary Zn could provide enhanced protection in the intestinal tract and stimulate the systemic humoral immune response against TGEV infection.

Creating cloud platforms for supporting FAIR data management in biomedical research projects.

Biomedical research projects are becoming increasingly complex and require technological solutions that support all phases of the data lifecycle and application of the FAIR principles. At the Berlin Institute of Health (BIH), we have developed and established a flexible and cost-effective approach to building customized cloud platforms for supporting research projects. The approach is based on a microservice architecture and on the management of a portfolio of supported services. On this basis, we created and maintained cloud platforms for several international research projects. In this article, we present our approach and argue that building customized cloud platforms can offer multiple advantages over using multi-project platforms. Our approach is transferable to other research environments and can be easily adapted by other projects and other service providers.

Identification of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses.

Group A rotaviruses (RVAs) are an important cause of diarrhoeal illness in humans, as well as in mammalian and avian animal species. Previous sequence analyses indicated that avian RVAs are related only distantly to mammalian RVAs. Here, the complete genomes of RVA strain 03V0002E10 from turkey (Meleagris gallopavo) and RVA strain 10V0112H5 from pheasant (Phasianus colchicus) were analysed using a combination of 454 deep sequencing and Sanger sequencing technologies. An adenine-rich insertion similar to that found in the chicken RVA strain 02V0002G3, but considerably shorter, was found in the 3' NCR of the NSP1 gene of the pheasant strain. Most genome segments of both strains were related closely to those of avian RVAs. The novel genotype N10 was assigned to the NSP2 gene of the pheasant RVA, which is related most closely to genotype N6 found in avian RVAs. However, this virus contains a VP4 gene of the novel genotype P[37], which is related most closely to RVAs from pigs, dogs and humans. This strain either may represent an avian/mammalian rotavirus reassortant, or it carries an unusual avian rotavirus VP4 gene, thereby broadening the potential genetic and antigenic variability among RVAs.

Antiviral effects of a probiotic Enterococcus faecium strain against transmissible gastroenteritis coronavirus.

The enteropathogenic coronavirus transmissible gastroenteritis virus (TGEV) causes severe disease in young piglets. We have studied the protective effects of the probiotic Enterococcus faecium NCIMB 10415 (E. faecium), which is approved as a feed additive in the European Union, against TGEV infection. E. faecium was added to swine testicle (ST) cells before, concomitantly with, or after TGEV infection. Viability assays revealed that E. faecium led to a dose-dependent rescue of viability of TGEV-infected cells reaching nearly to complete protection. Virus yields of the E. faecium-treated cultures were reduced by up to three log10 units. Western blot analysis of purified TGEV revealed that the levels of all viral structural proteins were reduced after E. faecium treatment. Using transmission electron microscopy, we observed attachment of TGEV particles to the surface of E. faecium which might be a means to trap virus and to prevent infection. Increased production of nitric oxide in the cells treated with E. faecium and elevated expression of interleukin 6 and 8 pointed to stimulated cellular defense as a mechanism to fight TGEV infection.

Addressing Global Environmental Challenges to Mental Health Using Population Neuroscience: A Review.

Importance: Climate change, pollution, urbanization, socioeconomic inequality, and psychosocial effects of the COVID-19 pandemic have caused massive changes in environmental conditions that affect brain health during the life span, both on a population level as well as on the level of the individual. How these environmental factors influence the brain, behavior, and mental illness is not well known. Observations: A research strategy enabling population neuroscience to contribute to identify brain mechanisms underlying environment-related mental illness by leveraging innovative enrichment tools for data federation, geospatial observation, climate and pollution measures, digital health, and novel data integration techniques is described. This strategy can inform innovative treatments that target causal cognitive and molecular mechanisms of mental illness related to the environment. An example is presented of the environMENTAL Project that is leveraging federated cohort data of over 1.5 million European citizens and patients enriched with deep phenotyping data from large-scale behavioral neuroimaging cohorts to identify brain mechanisms related to environmental adversity underlying symptoms of depression, anxiety, stress, and substance misuse. Conclusions and Relevance: This research will lead to the development of objective biomarkers and evidence-based interventions that will significantly improve outcomes of environment-related mental illness.

Allelic variations in coding regions of the vitamin D receptor gene in dairy cows and potential susceptibility to periparturient hypocalcaemia.

Periparturient hypocalcaemia (milk fever) is a disorder of Ca metabolism in dairy cattle primarily affecting multiparous cows. The major reasons for the rapid decrease of blood Ca concentration after calving are the prompt increase of Ca secretion into the colostrum and the delayed activation of Ca regulation mechanisms including calcitriol, a metabolite of vitamin D. In man, vitamin D receptor (VDR) gene polymorphisms are reported to be associated with disturbances of Ca metabolism, whereas data confirming the same in dairy cows are still missing. Moreover, polymorphisms that only affect non-coding regions are sometimes difficult to ascribe to a specific disorder as pathways and unequivocal links remain elusive. Therefore, the idea of the present study was to investigate in a small group of dairy cows with documented clinical records whether polymorphisms in the coding regions of the VDR gene existed and whether these potentially found variations were correlated with the incidence of periparturient hypocalcaemia. For this purpose, blood DNA was isolated from 26 dairy cows in their 4th to 6th lactation, out of which 17 had experienced hypocalcaemia at least once, whereas 9 cows had never undergone periparturient hypocalcaemia in their lifetime. The 10 VDR exons and small parts of adjacent introns were sequenced and compared with the Bos taurus VDR sequence published on NCBI based on the DNA of one Hereford cow. In total, 8 sequence alterations were detected in the fragments, which were primarily heterozygous. However, only 4 of them were really located on exons thereby potentially causing changes of the encoded amino acid of the VDR protein, but were not correlated with the incidence of periparturient hypocalcaemia. Certainly, this lack of statistical correlation could be due to the small number of animals included; anyhow, it was not encouraging enough to initiate a larger study with hundreds of cows and document blood Ca levels post partum for at least four lactations.

Sustainable data analysis with Snakemake.

Data analysis often entails a multitude of heterogeneous steps, from the application of various command line tools to the usage of scripting languages like R or Python for the generation of plots and tables. It is widely recognized that data analyses should ideally be conducted in a reproducible way. Reproducibility enables technical validation and regeneration of results on the original or even new data. However, reproducibility alone is by no means sufficient to deliver an analysis that is of lasting impact (i.e., sustainable) for the field, or even just one research group. We postulate that it is equally important to ensure adaptability and transparency. The former describes the ability to modify the analysis to answer extended or slightly different research questions. The latter describes the ability to understand the analysis in order to judge whether it is not only technically, but methodologically valid. Here, we analyze the properties needed for a data analysis to become reproducible, adaptable, and transparent. We show how the popular workflow management system Snakemake can be used to guarantee this, and how it enables an ergonomic, combined, unified representation of all steps involved in data analysis, ranging from raw data processing, to quality control and fine-grained, interactive exploration and plotting of final results.

Complex landscape of alternative splicing in myeloid neoplasms.

Myeloid neoplasms are characterized by frequent mutations in at least seven components of the spliceosome that have distinct roles in the process of pre-mRNA splicing. Hotspot mutations in SF3B1, SRSF2, U2AF1 and loss of function mutations in ZRSR2 have revealed widely different aberrant splicing signatures with little overlap. However, previous studies lacked the power necessary to identify commonly mis-spliced transcripts in heterogeneous patient cohorts. By performing RNA-Seq on bone marrow samples from 1258 myeloid neoplasm patients and 63 healthy bone marrow donors, we identified transcripts frequently mis-spliced by mutated splicing factors (SF), rare SF mutations with common alternative splicing (AS) signatures, and SF-dependent neojunctions. We characterized 17,300 dysregulated AS events using a pipeline designed to predict the impact of mis-splicing on protein function. Meta-splicing analysis revealed a pattern of reduced levels of retained introns among disease samples that was exacerbated in patients with splicing factor mutations. These introns share characteristics with "detained introns," a class of introns that have been shown to promote differentiation by detaining pro-proliferative transcripts in the nucleus. In this study, we have functionally characterized 17,300 targets of mis-splicing by the SF mutations, identifying a common pathway by which AS may promote maintenance of a proliferative state.

Homozygous BCMA gene deletion in response to anti-BCMA CAR T cells in a patient with multiple myeloma.

B cell maturation antigen (BCMA) is a target for various immunotherapies and a biomarker for tumor load in multiple myeloma (MM). We report a case of irreversible BCMA loss in a patient with MM who was enrolled in the KarMMa trial ( NCT03361748 ) and progressed after anti-BCMA CAR T cell therapy. We identified selection of a clone with homozygous deletion of TNFRSF17 (BCMA) as the underlying mechanism of immune escape. Furthermore, we found heterozygous TNFRSF17 loss or monosomy 16 in 37 out of 168 patients with MM, including 28 out of 33 patients with hyperhaploid MM who had not been previously treated with BCMA-targeting therapies, suggesting that heterozygous TNFRSF17 deletion at baseline could theoretically be a risk factor for BCMA loss after immunotherapy.