Royal Netherlands Academy of Arts and Sciences

Fungi play major roles in ecosystem processes, but the determinants of fungal diversity and biogeographic patterns remain poorly understood. Using DNA metabarcoding data from hundreds of globally distributed soil samples, we demonstrate that fungal richness is decoupled from plant diversity. The plant-to-fungus richness ratio declines exponentially toward the poles. Climatic factors, followed by edaphic and spatial variables, constitute the best predictors of fungal richness and community composition at the global scale. Fungi show similar latitudinal diversity gradients to other organisms, with several notable exceptions. These findings advance our understanding of global fungal diversity patterns and permit integration of fungi into a general macroecological framework.

Microbial communities play a pivotal role in the functioning of plants by influencing their physiology and development. While many members of the rhizosphere microbiome are beneficial to plant growth, also plant pathogenic microorganisms colonize the rhizosphere striving to break through the protective microbial shield and to overcome the innate plant defense mechanisms in order to cause disease. A third group of microorganisms that can be found in the rhizosphere are the true and opportunistic human pathogenic bacteria, which can be carried on or in plant tissue and may cause disease when introduced into debilitated humans. Although the importance of the rhizosphere microbiome for plant growth has been widely recognized, for the vast majority of rhizosphere microorganisms no knowledge exists. To enhance plant growth and health, it is essential to know which microorganism is present in the rhizosphere microbiome and what they are doing. Here, we review the main functions of rhizosphere microorganisms and how they impact on health and disease. We discuss the mechanisms involved in the multitrophic interactions and chemical dialogues that occur in the rhizosphere. Finally, we highlight several strategies to redirect or reshape the rhizosphere microbiome in favor of microorganisms that are beneficial to plant growth and health.

UNLABELLED: Sambamba is a high-performance robust tool and library for working with SAM, BAM and CRAM sequence alignment files; the most common file formats for aligned next generation sequencing data. Sambamba is a faster alternative to samtools that exploits multi-core processing and dramatically reduces processing time. Sambamba is being adopted at sequencing centers, not only because of its speed, but also because of additional functionality, including coverage analysis and powerful filtering capability. AVAILABILITY AND IMPLEMENTATION: Sambamba is free and open source software, available under a GPLv2 license. Sambamba can be downloaded and installed from http://www.open-bio.org/wiki/Sambamba.Sambamba v0.5.0 was released with doi:10.5281/zenodo.13200.

Pancreatic stellate cells (PSCs) differentiate into cancer-associated fibroblasts (CAFs) that produce desmoplastic stroma, thereby modulating disease progression and therapeutic response in pancreatic ductal adenocarcinoma (PDA). However, it is unknown whether CAFs uniformly carry out these tasks or if subtypes of CAFs with distinct phenotypes in PDA exist. We identified a CAF subpopulation with elevated expression of α-smooth muscle actin (αSMA) located immediately adjacent to neoplastic cells in mouse and human PDA tissue. We recapitulated this finding in co-cultures of murine PSCs and PDA organoids, and demonstrated that organoid-activated CAFs produced desmoplastic stroma. The co-cultures showed cooperative interactions and revealed another distinct subpopulation of CAFs, located more distantly from neoplastic cells, which lacked elevated αSMA expression and instead secreted IL6 and additional inflammatory mediators. These findings were corroborated in mouse and human PDA tissue, providing direct evidence for CAF heterogeneity in PDA tumor biology with implications for disease etiology and therapeutic development.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly ‘housekeeping’, whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research. A study from the FANTOM consortium using single-molecule cDNA sequencing of transcription start sites and their usage in human and mouse primary cells, cell lines and tissues reveals insights into the specificity and diversity of transcription patterns across different mammalian cell types. FANTOM5 (standing for functional annotation of the mammalian genome 5) is the fifth major stage of a major international collaboration that aims to dissect the transcriptional regulatory networks that define every human cell type. Two Articles in this issue of Nature present some of the project's latest results. The first paper uses the FANTOM5 panel of tissue and primary cell samples to define an atlas of active, in vivo bidirectionally transcribed enhancers across the human body. These authors show that bidirectional capped RNAs are a signature feature of active enhancers and identify more than 40,000 enhancer candidates from over 800 human cell and tissue samples. The enhancer atlas is used to compare regulatory programs between different cell types and identify disease-associated regulatory SNPs, and will be a resource for studies on cell-type-specific enhancers. In the second paper, single-molecule sequencing is used to map human and mouse transcription start sites and their usage in a panel of distinct human and mouse primary cells, cell lines and tissues to produce the most comprehensive mammalian gene expression atlas to date. The data provide a plethora of insights into open reading frames and promoters across different cell types in addition to valuable annotation of mammalian cell-type-specific transcriptomes.

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

The mammalian intestine is covered by a single layer of epithelial cells that is renewed every 4-5 days. This high cell turnover makes it a very attractive and comprehensive adult organ system for the study of cell proliferation and differentiation. The intestine is composed of proliferative crypts, which contain intestinal stem cells, and villi, which contain differentiated specialized cell types. Through the recent identification of Lgr5, an intestinal stem cell marker, it is now possible to visualize stem cells and study their behavior and differentiation in a much broader context. In this review we describe the identification of intestinal stem cells. We also discuss genetic studies that have helped to elucidate those signals important for progenitor cells to differentiate into one of the specialized intestinal epithelial cell types. These studies describe a genetic hierarchy responsible for cell fate commitment in normal gut physiology. Where relevant we also mention aberrant deregulation of these molecular pathways that results in colon cancer.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can cause coronavirus disease 2019 (COVID-19), an influenza-like disease that is primarily thought to infect the lungs with transmission through the respiratory route. However, clinical evidence suggests that the intestine may present another viral target organ. Indeed, the SARS-CoV-2 receptor angiotensin-converting enzyme 2 (ACE2) is highly expressed on differentiated enterocytes. In human small intestinal organoids (hSIOs), enterocytes were readily infected by SARS-CoV and SARS-CoV-2, as demonstrated by confocal and electron microscopy. Enterocytes produced infectious viral particles, whereas messenger RNA expression analysis of hSIOs revealed induction of a generic viral response program. Therefore, the intestinal epithelium supports SARS-CoV-2 replication, and hSIOs serve as an experimental model for coronavirus infection and biology.

To understand organ function, it is important to have an inventory of its cell types and of their corresponding marker genes. This is a particularly challenging task for human tissues like the pancreas, because reliable markers are limited. Hence, transcriptome-wide studies are typically done on pooled islets of Langerhans, obscuring contributions from rare cell types and of potential subpopulations. To overcome this challenge, we developed an automated platform that uses FACS, robotics, and the CEL-Seq2 protocol to obtain the transcriptomes of thousands of single pancreatic cells from deceased organ donors, allowing in silico purification of all main pancreatic cell types. We identify cell type-specific transcription factors and a subpopulation of REG3A-positive acinar cells. We also show that CD24 and TM4SF4 expression can be used to sort live alpha and beta cells with high purity. This resource will be useful for developing a deeper understanding of pancreatic biology and pathophysiology of diabetes mellitus.

Stem cells fuel tissue development, renewal, and regeneration, and these activities are controlled by the local stem cell microenvironment, the "niche." Wnt signals emanating from the niche can act as self-renewal factors for stem cells in multiple mammalian tissues. Wnt proteins are lipid-modified, which constrains them to act as short-range cellular signals. The locality of Wnt signaling dictates that stem cells exiting the Wnt signaling domain differentiate, spatially delimiting the niche in certain tissues. In some instances, stem cells may act as or generate their own niche, enabling the self-organization of patterned tissues. In this Review, we discuss the various ways by which Wnt operates in stem cell control and, in doing so, identify an integral program for tissue renewal and regeneration.

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

The phylum Acidobacteria is one of the most widespread and abundant on the planet, yet remarkably our knowledge of the role of these diverse organisms in the functioning of terrestrial ecosystems remains surprisingly rudimentary. This blatant knowledge gap stems to a large degree from the difficulties associated with the cultivation of these bacteria by classical means. Given the phylogenetic breadth of the Acidobacteria, which is similar to the metabolically diverse Proteobacteria, it is clear that detailed and functional descriptions of acidobacterial assemblages are necessary. Fortunately, recent advances are providing a glimpse into the ecology of members of the phylum Acidobacteria. These include novel cultivation and enrichment strategies, genomic characterization and analyses of metagenomic DNA from environmental samples. Here, we couple the data from these complementary approaches for a better understanding of their role in the environment, thereby providing some initial insights into the ecology of this important phylum. All cultured acidobacterial type species are heterotrophic, and members of subdivisions 1, 3, and 4 appear to be more versatile in carbohydrate utilization. Genomic and metagenomic data predict a number of ecologically relevant capabilities for some acidobacteria, including the ability to: use of nitrite as N source, respond to soil macro-, micro nutrients and soil acidity, express multiple active transporters, degrade gellan gum and produce exopolysaccharide (EPS). Although these predicted properties allude to a competitive life style in soil, only very few of these prediction shave been confirmed via physiological studies. The increased availability of genomic and physiological information, coupled to distribution data in field surveys and experiments, should direct future progress in unraveling the ecology of this important but still enigmatic phylum.

Abstract Pancreatic cancer is the most lethal common solid malignancy. Systemic therapies are often ineffective, and predictive biomarkers to guide treatment are urgently needed. We generated a pancreatic cancer patient–derived organoid (PDO) library that recapitulates the mutational spectrum and transcriptional subtypes of primary pancreatic cancer. New driver oncogenes were nominated and transcriptomic analyses revealed unique clusters. PDOs exhibited heterogeneous responses to standard-of-care chemotherapeutics and investigational agents. In a case study manner, we found that PDO therapeutic profiles paralleled patient outcomes and that PDOs enabled longitudinal assessment of chemosensitivity and evaluation of synchronous metastases. We derived organoid-based gene expression signatures of chemosensitivity that predicted improved responses for many patients to chemotherapy in both the adjuvant and advanced disease settings. Finally, we nominated alternative treatment strategies for chemorefractory PDOs using targeted agent therapeutic profiling. We propose that combined molecular and therapeutic profiling of PDOs may predict clinical response and enable prospective therapeutic selection. Significance: New approaches to prioritize treatment strategies are urgently needed to improve survival and quality of life for patients with pancreatic cancer. Combined genomic, transcriptomic, and therapeutic profiling of PDOs can identify molecular and functional subtypes of pancreatic cancer, predict therapeutic responses, and facilitate precision medicine for patients with pancreatic cancer. Cancer Discov; 8(9); 1112–29. ©2018 AACR. See related commentary by Collisson, p. 1062. This article is highlighted in the In This Issue feature, p. 1047

Richard O'Connell and colleagues report the genomes and transcriptomes of two Colletotrichum plant fungal pathogens. C. higginsianum infects Arabidopsis thaliana, and C. graminicola infects maize (Zea mays); comparative genomics in both species lead to molecular insights into the transition from biotrophic to necrotrophic life stages. Colletotrichum species are fungal pathogens that devastate crop plants worldwide. Host infection involves the differentiation of specialized cell types that are associated with penetration, growth inside living host cells (biotrophy) and tissue destruction (necrotrophy). We report here genome and transcriptome analyses of Colletotrichum higginsianum infecting Arabidopsis thaliana and Colletotrichum graminicola infecting maize. Comparative genomics showed that both fungi have large sets of pathogenicity-related genes, but families of genes encoding secreted effectors, pectin-degrading enzymes, secondary metabolism enzymes, transporters and peptidases are expanded in C. higginsianum. Genome-wide expression profiling revealed that these genes are transcribed in successive waves that are linked to pathogenic transitions: effectors and secondary metabolism enzymes are induced before penetration and during biotrophy, whereas most hydrolases and transporters are upregulated later, at the switch to necrotrophy. Our findings show that preinvasion perception of plant-derived signals substantially reprograms fungal gene expression and indicate previously unknown functions for particular fungal cell types.

Cognitive functions rely on the coordinated activity of neurons in many brain regions, but the interactions between cortical areas are not yet well understood. Here we investigated whether low-frequency (α) and high-frequency (γ) oscillations characterize different directions of information flow in monkey visual cortex. We recorded from all layers of the primary visual cortex (V1) and found that γ-waves are initiated in input layer 4 and propagate to the deep and superficial layers of cortex, whereas α-waves propagate in the opposite direction. Simultaneous recordings from V1 and downstream area V4 confirmed that γ- and α-waves propagate in the feedforward and feedback direction, respectively. Microstimulation in V1 elicited γ-oscillations in V4, whereas microstimulation in V4 elicited α-oscillations in V1, thus providing causal evidence for the opposite propagation of these rhythms. Furthermore, blocking NMDA receptors, thought to be involved in feedback processing, suppressed α while boosting γ. These results provide new insights into the relation between brain rhythms and cognition.

Organoids are microscopic self-organizing, three-dimensional structures that are grown from stem cells in vitro. They recapitulate many structural and functional aspects of their in vivo counterpart organs. This versatile technology has led to the development of many novel human cancer models. It is now possible to create indefinitely expanding organoids starting from tumor tissue of individuals suffering from a range of carcinomas. Alternatively, CRISPR-based gene modification allows the engineering of organoid models of cancer through the introduction of any combination of cancer gene alterations to normal organoids. When combined with immune cells and fibroblasts, tumor organoids become models for the cancer microenvironment enabling immune-oncology applications. Emerging evidence indicates that organoids can be used to accurately predict drug responses in a personalized treatment setting. Here, we review the current state and future prospects of the rapidly evolving tumor organoid field.

In order to identify patterns in bacterial community composition in freshwater habitats, we analyzed the available database of 16S rDNA sequences from freshwater plankton, including 24 new sequences from Parker River (Massachusetts, USA), 42 from Lake Soyang (South Korea) and 148 from Lake IJssel (The Netherlands). At this point, combined diversity studies using random cloning have deposited 689 bacterial and 75 plastid 16S rDNA sequences from the water column of rivers and lakes in North America, Europe and Asia. Systematic comparisons with the global database showed that the majority of the bacterial sequences were most closely related to other freshwater clones or isolates, while relatively few were closest to sequences recovered from soils or marine habitats. This habitat-specific clustering suggests that the clustered 16S rDNA sequences represent species or groups of species that are indigenous to freshwater. We have discerned 34 phylogenetic clusters of closely related sequences that are either restricted to freshwater or dominated by freshwater sequences. Of these clusters, 23 contained no cultivated organisms. These putative freshwater clusters were found among the alpha-, beta-and gamma-Proteobacteria, the Cytophaga-Flavobacterium-Bacteroides group, the Cyanobacteria, the Actinobacteria, the Verrucomicrobia, the green non-sulfur bacteria and candidate division OP10. This study shows that rivers and lakes have a specific planktonic bacterial community distinct from bacteria in neighboring environments such as soil and sediments. It also points out that these planktonic bacteria are distributed in diverse freshwater ecosystems around the world.

Phenotypic variation is ubiquitous in biology and is often traceable to underlying genetic and environmental variation. However, even genetically identical cells in identical environments display variable phenotypes. Stochastic gene expression, or gene expression "noise," has been suggested as a major source of this variability, and its physiological consequences have been topics of intense research for the last decade. Several recent studies have measured variability in protein and messenger RNA levels, and they have discovered strong connections between noise and gene regulation mechanisms. When integrated with discrete stochastic models, measurements of cell-to-cell variability provide a sensitive "fingerprint" with which to explore fundamental questions of gene regulation. In this review, we highlight several studies that used gene expression variability to develop a quantitative understanding of the mechanisms and dynamics of gene regulation.

Cytotoxic T lymphocyte (CTL) peptide epitopes can be used for immunization of mice against lethal virus infection. To study whether this approach can be successful against virus-induced tumors we generated a B6 (H-2b) tumorigenic cell line transformed by human papillomavirus (HPV). This virus is detected in over 90% of all human cervical cancers. To identify vaccine candidates, we generated a set of 240 overlapping peptides derived from the HPV type 16 (HPV16) oncogenes E6 and E7. These peptides were tested for their ability to bind H-2Kb and H-2Db MHC class I molecules. Binding peptides were compared with the presently known peptide-binding motifs for H-2Kb and H-2Db and the predictive value of these motifs is shortly discussed. The high-affinity H-2Db-binding peptide and putative CTL epitope E7 49-57 (RAHYNIVTF) was used in vaccination studies against HPV 16-transformed tumor cells. Immunization with peptide E7 49-57 rendered mice insensitive to a subsequent challenge with HPV 16-transformed tumor cells in vivo, and induced a CTL response which lysed the tumor cells in vitro.

There is a clear and unmet clinical need for biomarkers to predict responsiveness to chemotherapy for cancer. We developed an in vitro test based on patient-derived tumor organoids (PDOs) from metastatic lesions to identify nonresponders to standard-of-care chemotherapy in colorectal cancer (CRC). In a prospective clinical study, we show the feasibility of generating and testing PDOs for evaluation of sensitivity to chemotherapy. Our PDO test predicted response of the biopsied lesion in more than 80% of patients treated with irinotecan-based therapies without misclassifying patients who would have benefited from treatment. This correlation was specific to irinotecan-based chemotherapy, however, and the PDOs failed to predict outcome for treatment with 5-fluorouracil plus oxaliplatin. Our data suggest that PDOs could be used to prevent cancer patients from undergoing ineffective irinotecan-based chemotherapy.