DRESDEN-concept Genome Center

Abstract High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species 1–4 . To address this issue, the international Genome 10K (G10K) consortium 5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Abstract Fast and reliable detection of patients with severe and heterogeneous illnesses is a major goal of precision medicine 1,2 . Patients with leukaemia can be identified using machine learning on the basis of their blood transcriptomes 3 . However, there is an increasing divide between what is technically possible and what is allowed, because of privacy legislation 4,5 . Here, to facilitate the integration of any medical data from any data owner worldwide without violating privacy laws, we introduce Swarm Learning—a decentralized machine-learning approach that unites edge computing, blockchain-based peer-to-peer networking and coordination while maintaining confidentiality without the need for a central coordinator, thereby going beyond federated learning. To illustrate the feasibility of using Swarm Learning to develop disease classifiers using distributed data, we chose four use cases of heterogeneous diseases (COVID-19, tuberculosis, leukaemia and lung pathologies). With more than 16,400 blood transcriptomes derived from 127 clinical studies with non-uniform distributions of cases and controls and substantial study biases, as well as more than 95,000 chest X-ray images, we show that Swarm Learning classifiers outperform those developed at individual sites. In addition, Swarm Learning completely fulfils local confidentiality regulations by design. We believe that this approach will notably accelerate the introduction of precision medicine.

Abstract High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are only available for a few non-microbial species 1–4 . To address this issue, the international Genome 10K (G10K) consortium 5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling the most accurate and complete reference genomes to date. Here we summarize these developments, introduce a set of quality standards, and present lessons learned from sequencing and assembling 16 species representing major vertebrate lineages (mammals, birds, reptiles, amphibians, teleost fishes and cartilaginous fishes). We confirm that long-read sequencing technologies are essential for maximizing genome quality and that unresolved complex repeats and haplotype heterozygosity are major sources of error in assemblies. Our new assemblies identify and correct substantial errors in some of the best historical reference genomes. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an effort to generate high-quality, complete reference genomes for all ~70,000 extant vertebrate species and help enable a new era of discovery across the life sciences.

Biallelic mutations of the CEBPA gene (CEBPAbi) define a distinct entity associated with favorable prognosis; however, the role of monoallelic mutations (CEBPAsm) is poorly understood. We retrospectively analyzed 4708 adults with acute myeloid leukemia (AML) who had been recruited into the Study Alliance Leukemia trials, to investigate the prognostic impact of CEBPAsm. CEBPA mutations were identified in 240 patients (5.1%): 131 CEBPAbi and 109 CEBPAsm (60 affecting the N-terminal transactivation domains [CEBPAsmTAD] and 49 the C-terminal DNA-binding or basic leucine zipper region [CEBPAsmbZIP]). Interestingly, patients carrying CEBPAbi or CEBPAsmbZIP shared several clinical factors: they were significantly younger (median, 46 and 50 years, respectively) and had higher white blood cell (WBC) counts at diagnosis (median, 23.7 × 109/L and 35.7 × 109/L) than patients with CEBPAsmTAD (median age, 63 years, median WBC 13.1 × 109/L; P < .001). Co-mutations were similar in both groups: GATA2 mutations (35.1% CEBPAbi; 36.7% CEBPAsmbZIP vs 6.7% CEBPAsmTAD; P < .001) or NPM1 mutations (3.1% CEBPAbi; 8.2% CEBPAsmbZIP vs 38.3% CEBPAsmTAD; P < .001). CEBPAbi and CEBPAsmbZIP, but not CEBPAsmTAD were associated with significantly improved overall (OS; median 103 and 63 vs 13 months) and event-free survival (EFS; median, 20.7 and 17.1 months vs 5.7 months), in univariate and multivariable analyses. Additional analyses revealed that the clinical and molecular features as well as the favorable survival were confined to patients with in-frame mutations in bZIP (CEBPAbZIP-inf). When patients were classified according to CEBPAbZIP-inf and CEBPAother (including CEBPAsmTAD and non-CEBPAbZIP-inf), only patients bearing CEBPAbZIP-inf showed superior complete remission rates and the longest median OS and EFS, arguing for a previously undefined prognostic role of this type of mutation.

Macrophages populate every organ during homeostasis and disease, displaying features of tissue imprinting and heterogeneous activation. The disconnected picture of macrophage biology that has emerged from these observations is a barrier for integration across models or with in vitro macrophage activation paradigms. We set out to contextualize macrophage heterogeneity across mouse tissues and inflammatory conditions, specifically aiming to define a common framework of macrophage activation. We built a predictive model with which we mapped the activation of macrophages across 12 tissues and 25 biological conditions, finding a notable commonality and finite number of transcriptional profiles, in particular among infiltrating macrophages, which we modeled as defined stages along four conserved activation paths. These activation paths include a "phagocytic" regulatory path, an "inflammatory" cytokine-producing path, an "oxidative stress" antimicrobial path, or a "remodeling" extracellular matrix deposition path. We verified this model with adoptive cell transfer experiments and identified transient RELMɑ expression as a feature of monocyte-derived macrophage tissue engraftment. We propose that this integrative approach of macrophage classification allows the establishment of a common predictive framework of monocyte-derived macrophage activation in inflammation and homeostasis.

The neural stem cell (NSC) reservoir can be harnessed for stem cell-based regenerative therapies. Zebrafish remarkably regenerate their brain by inducing NSC plasticity in a Amyloid-β-42 (Aβ42)-induced experimental Alzheimer's disease (AD) model. Interleukin-4 (IL-4) is also critical for AD-induced NSC proliferation. However, the mechanisms of this response have remained unknown. Using single-cell transcriptomics in the adult zebrafish brain, we identify distinct subtypes of NSCs and neurons and differentially regulated pathways and their gene ontologies and investigate how cell-cell communication is altered through ligand-receptor pairs in AD conditions. Our results propose the existence of heterogeneous and spatially organized stem cell populations that react distinctly to amyloid toxicity. This resource article provides an extensive database for the molecular basis of NSC plasticity in the AD model of the adult zebrafish brain. Further analyses of stem cell heterogeneity and neuro-regenerative ability at single-cell resolution could yield drug targets for mobilizing NSCs for endogenous neuro-regeneration in humans.

Cellular redox states regulate the balance between stem cell maintenance and activation. Increased levels of intracellular reactive oxygen species (ROS) are linked to proliferation and lineage specification. In contrast to this general principle, we here show that in the hippocampus of adult mice, quiescent neural precursor cells (NPCs) maintain the highest ROS levels (hiROS). Classifying NPCs on the basis of cellular ROS content identified distinct functional states. Shifts in ROS content primed cells for a subsequent state transition, with lower ROS content marking proliferative activity and differentiation. Physical activity, a physiological activator of adult hippocampal neurogenesis, recruited hiROS NPCs into proliferation via a transient Nox2-dependent ROS surge. In the absence of Nox2, baseline neurogenesis was unaffected, but the activity-induced increase in proliferation disappeared. These results provide a metabolic classification of NPC functional states and describe a mechanism linking the modulation of cellular ROS by behavioral cues to the activation of adult NPCs.

Objective The intestinal epithelium is a rapidly renewing tissue which plays central roles in nutrient uptake, barrier function and the prevention of intestinal inflammation. Control of epithelial differentiation is essential to these processes and is dependent on cell type-specific activity of transcription factors which bind to accessible chromatin. Here, we studied the role of SET Domain Bifurcated Histone Lysine Methyltransferase 1, also known as ESET (SETDB1), a histone H3K9 methyltransferase, in intestinal epithelial homeostasis and IBD. Design We investigated mice with constitutive and inducible intestinal epithelial deletion of Setdb1 , studied the expression of SETDB1 in patients with IBD and mouse models of IBD, and investigated the abundance of SETDB1 variants in healthy individuals and patients with IBD. Results Deletion of intestinal epithelial Setdb1 in mice was associated with defects in intestinal epithelial differentiation, barrier disruption, inflammation and mortality. Mechanistic studies showed that loss of SETDB1 leads to de-silencing of endogenous retroviruses, DNA damage and intestinal epithelial cell death. Predicted loss-of-function variants in human SETDB1 were considerably less frequently observed than expected, consistent with a critical role of SETDB1 in human biology. While the vast majority of patients with IBD showed unimpaired mucosal SETDB1 expression, comparison of IBD and non-IBD exomes revealed over-representation of individual rare missense variants in SETDB1 in IBD, some of which are predicted to be associated with loss of function and may contribute to the pathogenesis of intestinal inflammation. Conclusion SETDB1 plays an essential role in intestinal epithelial homeostasis. Future work is required to investigate whether rare variants in SETDB1 contribute to the pathogenesis of IBD.

ABSTRACT Zebrafish display widespread and pronounced adult neurogenesis, which is fundamental for their regeneration capability after central nervous system injury. However, the cellular identity and the biological properties of adult newborn neurons are elusive for most brain areas. Here, we have used short-term lineage tracing of radial glia progeny to prospectively isolate newborn neurons from the her4.1+ radial glia lineage in the homeostatic adult forebrain. Transcriptome analysis of radial glia, newborn neurons and mature neurons using single cell sequencing identified distinct transcriptional profiles, including novel markers for each population. Specifically, we detected two separate newborn neuron types, which showed diversity of cell fate commitment and location. Further analyses showed that these cell types are homologous to neurogenic cells in the mammalian brain, identified neurogenic commitment in proliferating radial glia and indicated that glutamatergic projection neurons are generated in the adult zebrafish telencephalon. Thus, we prospectively isolated adult newborn neurons from the adult zebrafish forebrain, identified markers for newborn and mature neurons in the adult brain, and revealed intrinsic heterogeneity among adult newborn neurons and their homology with mammalian adult neurogenic cell types.

Abstract Hematopoietic stem cells (HSCs) produce highly diverse cell lineages. Here, we chart native lineage pathways emanating from HSCs and define their physiological regulation by computationally integrating experimental approaches for fate mapping, mitotic tracking, and single-cell RNA sequencing. We find that lineages begin to split when cells leave the tip HSC population, marked by high Sca-1 and CD201 expression. Downstream, HSCs either retain high Sca-1 expression and the ability to generate lymphocytes, or irreversibly reduce Sca-1 level and enter into erythro-myelopoiesis or thrombopoiesis. Thrombopoiesis is the sum of two pathways that make comparable contributions in steady state, a long route via multipotent progenitors and CD48 hi megakaryocyte progenitors (MkPs), and a short route from HSCs to developmentally distinct CD48 −/lo MkPs. Enhanced thrombopoietin signaling differentially accelerates the short pathway, enabling a rapid response to increasing demand. In sum, we provide a blueprint for mapping physiological differentiation fluxes from HSCs and decipher two functionally distinct pathways of native thrombopoiesis.

Pro-inflammatory autoantigen-specific CD4 + T helper (auto-Th) cells are central orchestrators of autoimmune diseases (AIDs). We aimed to characterize these cells in human AIDs with defined autoantigens by combining human leukocyte antigen (HLA)-tetramer-based and activation-based multidimensional ex vivo analyses. In aquaporin4-antibody-positive neuromyelitis optica spectrum disorder (AQP4-NMOSD) patients, auto-Th cells expressed CD154, but proliferative capacity and pro-inflammatory cytokines were strongly reduced. Instead, exhaustion-associated co-inhibitory receptors were expressed together with FOXP3, the canonical regulatory T cell (Treg) transcription factor. Auto-Th cells responded in vitro to checkpoint inhibition and provided potent B cell help. Cells with the same exhaustion-like (ThEx) phenotype were identified in soluble liver antigen (SLA)-antibody-autoimmune hepatitis and BP180-antibody-positive bullous pemphigoid, AIDs of the liver and skin, respectively. While originally described in cancer and chronic infection, our data point to T cell exhaustion as a common mechanism of adaptation to chronic (self-)stimulation across AID types and link exhausted CD4 + T cells to humoral autoimmune responses, with implications for therapeutic targeting.

The decline of brain function during aging is associated with epigenetic changes, including DNA methylation. Lifestyle interventions can improve brain function during aging, but their influence on age-related epigenetic changes is unknown. Using genome-wide DNA methylation sequencing, we here show that experiencing a stimulus-rich environment counteracts age-related DNA methylation changes in the hippocampal dentate gyrus of mice. Specifically, environmental enrichment prevented the aging-induced CpG hypomethylation at target sites of the methyl-CpG-binding protein Mecp2, which is critical to neuronal function. The genes at which environmental enrichment counteracted aging effects have described roles in neuronal plasticity, neuronal cell communication and adult hippocampal neurogenesis and are dysregulated with age-related cognitive decline in the human brain. Our results highlight the stimulating effects of environmental enrichment on hippocampal plasticity at the level of DNA methylation and give molecular insights into the specific aspects of brain aging that can be counteracted by lifestyle interventions.

Inferring historical population sizes is key to identifying drivers of ecological and evolutionary change and crucial to predicting the future of species on our rapidly changing planet. The pairwise sequentially Markovian coalescent (PSMC) method provided a revolutionary framework to reconstruct species' demographic histories over millions of years based on the genome sequence of a single individual. Here, we detected and solved a common artifact in PSMC and related methods: recent population peaks followed by population collapses. Combining real and simulated genomes, we show that these peaks do not represent true population dynamics. Instead, ill-set default parameters cause false peaks in our own and published data, which can be avoided by adjusting parameter settings. Furthermore, we show that certain changes in population structure can cause similar patterns. Newer methods, like Beta-PSMC, perform better but do not always avoid this artifact. Our results suggest testing multiple parameters that split the first time window before interpreting recent population peaks followed by collapses and call for the development of robust methods.

Zoonoses are infectious diseases transmitted from animals to humans. Bats have been suggested to harbour more zoonotic viruses than any other mammalian order1. Infections in bats are largely asymptomatic2,3, indicating limited tissue-damaging inflammation and immunopathology. To investigate the genomic basis of disease resistance, the Bat1K project generated reference-quality genomes of ten bat species, including potential viral reservoirs. Here we describe a systematic analysis covering 115 mammalian genomes that revealed that signatures of selection in immune genes are more prevalent in bats than in other mammalian orders. We found an excess of immune gene adaptations in the ancestral chiropteran branch and in many descending bat lineages, highlighting viral entry and detection factors, and regulators of antiviral and inflammatory responses. ISG15, which is an antiviral gene contributing to hyperinflammation during COVID-19 (refs. 4,5), exhibits key residue changes in rhinolophid and hipposiderid bats. Cellular infection experiments show species-specific antiviral differences and an essential role of protein conjugation in antiviral function of bat ISG15, separate from its role in secretion and inflammation in humans. Furthermore, in contrast to humans, ISG15 in most rhinolophid and hipposiderid bats has strong anti-SARS-CoV-2 activity. Our work reveals molecular mechanisms that contribute to viral tolerance and disease resistance in bats. A systematic analysis of 115 mammalian genomes, including 10 new bat genomes, reveals prevalent positive selection in immune genes in bats and shows key adaptations in the antiviral gene ISG15 that aid disease resistance in bats, including to coronaviruses.

Conditional gene inactivation is a powerful tool to determine gene function when constitutive mutations result in detrimental effects. The most commonly used technique to achieve conditional gene inactivation employs the Cre/loxP system and its ability to delete DNA sequences flanked by two loxP sites. However, targeting a gene with two loxP sites is time and labor consuming. Here, we show Cre-Controlled CRISPR (3C) mutagenesis to circumvent these issues. 3C relies on gRNA and Cre-dependent Cas9-GFP expression from the same transgene. Exogenous or transgenic supply of Cre results in Cas9-GFP expression and subsequent mutagenesis of the gene of interest. The recombined cells become fluorescently visible enabling their isolation and subjection to various omics techniques. Hence, 3C mutagenesis provides a valuable alternative to the production of loxP-flanked alleles. It might even enable the conditional inactivation of multiple genes simultaneously and should be applicable to other model organisms amenable to single integration transgenesis.

The influence of aging on intestinal stem cells and their niche can explain underlying causes for perturbation in their function observed during aging. Molecular mechanisms for such a decrease in the functionality of intestinal stem cells during aging remain largely undetermined. Using transcriptome-wide approaches, our study demonstrates that aging intestinal stem cells strongly upregulate antigen presenting pathway genes and over-express secretory lineage marker genes resulting in lineage skewed differentiation into the secretory lineage and strong upregulation of MHC class II antigens in the aged intestinal epithelium. Mechanistically, we identified an increase in proinflammatory cells in the lamina propria as the main source of elevated interferon gamma (IFNγ) in the aged intestine, that leads to the induction of Stat1 activity in intestinal stem cells thus priming the aberrant differentiation and elevated antigen presentation in epithelial cells. Of note, systemic inhibition of IFNγ-signaling completely reverses these aging phenotypes and reinstalls regenerative capacity of the aged intestinal epithelium.

Loss of limbs evolved many times in squamate reptiles. Here we investigated the genomic basis of convergent limb loss in reptiles. We sequenced the genomes of a closely related pair of limbless-limbed gymnophthalmid lizards and performed a comparative genomic analysis including five snakes and the limbless glass lizard. Our analysis of these three independent limbless lineages revealed that signatures of shared sequence or transcription factor binding site divergence in individual limb regulatory elements are generally rare. Instead, shared divergence occurs more often at the level of signaling pathways, involving different regulatory elements associated with the same limb genes (such as Hand2 or Hox) and/or patterning mechanisms (such as Shh signaling). Interestingly, although snakes are known to have mutations in the Shh ZRS limb enhancer, this enhancer lacks relevant mutations in limbless lizards. Thus, different mechanisms could contribute to limb loss, and there are likely multiple evolutionary paths to limblessness in reptiles.

A genomic database of all Earth's eukaryotic species could contribute to many scientific discoveries; however, only a tiny fraction of species have genomic information available. In 2018, scientists across the world united under the Earth BioGenome Project (EBP), aiming to produce a database of high-quality reference genomes containing all ~1.5 million recognized eukaryotic species. As the European node of the EBP, the European Reference Genome Atlas (ERGA) sought to implement a new decentralised, equitable and inclusive model for producing reference genomes. For this, ERGA launched a Pilot Project establishing the first distributed reference genome production infrastructure and testing it on 98 eukaryotic species from 33 European countries. Here we outline the infrastructure and explore its effectiveness for scaling high-quality reference genome production, whilst considering equity and inclusion. The outcomes and lessons learned provide a solid foundation for ERGA while offering key learnings to other transnational, national genomic resource projects and the EBP.

The immaturity of human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) is a major limitation for their use in drug screening to identify pro-arrhythmogenic or cardiotoxic molecules. Here, we demonstrate an approach that combines lipid-enriched maturation medium with a high concentration of calcium, nanopatterning of culture surfaces and electrostimulation to generate iPSC-CMs with advanced electrophysiological, structural and metabolic phenotypes. Systematic testing reveals that electrostimulation is the key driver of enhanced mitochondrial development and metabolic maturation and improved electrophysiological properties of iPSC-CMs. Increased calcium concentration strongly promotes electrophysiological maturation, while nanopatterning primarily facilitates sarcomere organisation with minor effect on electrophysiological properties. Transcriptome analysis reveals that activation of HMCES and TFAM targets contributes to mitochondrial development, whereas downregulation of MAPK/PI3K and SRF targets is associated with iPSC-CM polyploidy. These findings provide mechanistic insights into iPSC-CM maturation, paving the way for pharmacological responses that more closely resemble those of adult CMs.

CD8+ T cells are important effectors of adaptive immunity against pathogens, tumors and self antigens. Here, we asked how human cognate antigen-responsive CD8+ T cells and their receptors could be identified in unselected single-cell gene expression data. Single-cell RNA sequencing and qPCR of dye-labelled antigen-specific cells identified large gene sets that were congruently up- or downregulated in virus-responsive CD8+ T cells under different antigen presentation conditions. Combined expression of TNFRSF9, XCL1, XCL2, and CRTAM was the most distinct marker of virus-responsive cells on a single-cell level. Using transcriptomic data, we developed a machine learning-based classifier that provides sensitive and specific detection of virus-responsive CD8+ T cells from unselected populations. Gene response profiles of CD8+ T cells specific for the autoantigen islet-specific glucose-6-phosphatase catalytic subunit-related protein differed markedly from virus-specific cells. These findings provide single-cell gene expression parameters for comprehensive identification of rare antigen-responsive cells and T cell receptors.