Hybrid assemblies of microbiome Blastocystis protists reveal evolutionary diversification reflecting host ecology.

The most prevalent microbial eukaryote in the human gut is Blastocystis, an obligate commensal protist also common in many other vertebrates. Blastocystis is descended from free-living stramenopile ancestors; how it has adapted to thrive within humans and a wide range of hosts is unclear. Here, we cultivated six Blastocystis strains spanning the diversity of the genus and generated highly contiguous, annotated genomes with long-read DNA-seq, Hi-C, and RNA-seq. Comparative genomics between these strains and two closely related stramenopiles with different lifestyles, the lizard gut symbiont Proteromonas lacertae and the free-living marine flagellate Cafeteria burkhardae, reveal the evolutionary history of the Blastocystis genus. We find substantial gene content variability between Blastocystis strains. Blastocystis isolated from an herbivorous tortoise has many plant carbohydrate metabolizing enzymes, some horizontally acquired from bacteria, likely reflecting fermentation within the host gut. In contrast, human-isolated Blastocystis have gained many heat shock proteins, and we find numerous subtype-specific expansions of host-interfacing genes, including cell adhesion and cell surface glycan genes. In addition, we observe that human-isolated Blastocystis have substantial changes in gene structure, including shortened introns and intergenic regions, as well as genes lacking canonical termination codons. Finally, our data indicate that the common ancestor of Blastocystis lost nearly all ancestral genes for heterokont flagella morphology, including cilia proteins, microtubule motor proteins, and ion channel proteins. Together, these findings underscore the huge functional variability within the Blastocystis genus and provide candidate genes for the adaptations these lineages have undergone to thrive in the gut microbiomes of diverse vertebrates.

Evolutionary constraint and innovation across hundreds of placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Beyond the Biosynthetic Gene Cluster Paradigm: Genome-Wide Coexpression Networks Connect Clustered and Unclustered Transcription Factors to Secondary Metabolic Pathways.

Fungal secondary metabolites are widely used as therapeutics and are vital components of drug discovery programs. A major challenge hindering discovery of novel secondary metabolites is that the underlying pathways involved in their biosynthesis are transcriptionally silent under typical laboratory growth conditions, making it difficult to identify the transcriptional networks that they are embedded in. Furthermore, while the genes participating in secondary metabolic pathways are typically found in contiguous clusters on the genome, known as biosynthetic gene clusters (BGCs), this is not always the case, especially for global and pathway-specific regulators of pathways' activities. To address these challenges, we used 283 genome-wide gene expression data sets of the ascomycete cell factory Aspergillus niger generated during growth under 155 different conditions to construct two gene coexpression networks based on Spearman's correlation coefficients (SCCs) and on mutual rank-transformed Pearson's correlation coefficients (MR-PCCs). By mining these networks, we predicted six transcription factors, named MjkA to MjkF, to regulate secondary metabolism in A. niger. Overexpression of each transcription factor using the Tet-On cassette modulated the production of multiple secondary metabolites. We found that the SCC and MR-PCC approaches complemented each other, enabling the delineation of putative global (SCC) and pathway-specific (MR-PCC) transcription factors. These results highlight the potential of coexpression network approaches to identify and activate fungal secondary metabolic pathways and their products. More broadly, we argue that drug discovery programs in fungi should move beyond the BGC paradigm and focus on understanding the global regulatory networks in which secondary metabolic pathways are embedded. IMPORTANCE There is an urgent need for novel bioactive molecules in both agriculture and medicine. The genomes of fungi are thought to contain vast numbers of metabolic pathways involved in the biosynthesis of secondary metabolites with diverse bioactivities. Because these metabolites are biosynthesized only under specific conditions, the vast majority of the fungal pharmacopeia awaits discovery. To discover the genetic networks that regulate the activity of secondary metabolites, we examined the genome-wide profiles of gene activity of the cell factory Aspergillus niger across hundreds of conditions. By constructing global networks that link genes with similar activities across conditions, we identified six putative global and pathway-specific regulators of secondary metabolite biosynthesis. Our study shows that elucidating the behavior of the genetic networks of fungi under diverse conditions harbors enormous promise for understanding fungal secondary metabolism, which ultimately may lead to novel drug candidates.

Accurate and sensitive detection of microbial eukaryotes from whole metagenome shotgun sequencing.

BACKGROUND: Microbial eukaryotes are found alongside bacteria and archaea in natural microbial systems, including host-associated microbiomes. While microbial eukaryotes are critical to these communities, they are challenging to study with shotgun sequencing techniques and are therefore often excluded. RESULTS: Here, we present EukDetect, a bioinformatics method to identify eukaryotes in shotgun metagenomic sequencing data. Our approach uses a database of 521,824 universal marker genes from 241 conserved gene families, which we curated from 3713 fungal, protist, non-vertebrate metazoan, and non-streptophyte archaeplastida genomes and transcriptomes. EukDetect has a broad taxonomic coverage of microbial eukaryotes, performs well on low-abundance and closely related species, and is resilient against bacterial contamination in eukaryotic genomes. Using EukDetect, we describe the spatial distribution of eukaryotes along the human gastrointestinal tract, showing that fungi and protists are present in the lumen and mucosa throughout the large intestine. We discover that there is a succession of eukaryotes that colonize the human gut during the first years of life, mirroring patterns of developmental succession observed in gut bacteria. By comparing DNA and RNA sequencing of paired samples from human stool, we find that many eukaryotes continue active transcription after passage through the gut, though some do not, suggesting they are dormant or nonviable. We analyze metagenomic data from the Baltic Sea and find that eukaryotes differ across locations and salinity gradients. Finally, we observe eukaryotes in Arabidopsis leaf samples, many of which are not identifiable from public protein databases. CONCLUSIONS: EukDetect provides an automated and reliable way to characterize eukaryotes in shotgun sequencing datasets from diverse microbiomes. We demonstrate that it enables discoveries that would be missed or clouded by false positives with standard shotgun sequence analysis. EukDetect will greatly advance our understanding of how microbial eukaryotes contribute to microbiomes. Video abstract.

Pathogenic Allodiploid Hybrids of Aspergillus Fungi.

Interspecific hybridization substantially alters genotypes and phenotypes and can give rise to new lineages. Hybrid isolates that differ from their parental species in infection-relevant traits have been observed in several human-pathogenic yeasts and plant-pathogenic filamentous fungi but have yet to be found in human-pathogenic filamentous fungi. We discovered 6 clinical isolates from patients with aspergillosis originally identified as Aspergillus nidulans (section Nidulantes) that are actually allodiploid hybrids formed by the fusion of Aspergillus spinulosporus with an unknown close relative of Aspergillus quadrilineatus, both in section Nidulantes. Evolutionary genomic analyses revealed that these isolates belong to Aspergillus latus, an allodiploid hybrid species. Characterization of diverse infection-relevant traits further showed that A. latus hybrid isolates are genomically and phenotypically heterogeneous but also differ from A. nidulans, A. spinulosporus, and A. quadrilineatus. These results suggest that allodiploid hybridization contributes to the genomic and phenotypic diversity of filamentous fungal pathogens of humans.

A Robust Phylogenomic Time Tree for Biotechnologically and Medically Important Fungi in the Genera Aspergillus and Penicillium.

The filamentous fungal family Aspergillaceae contains >1,000 known species, mostly in the genera Aspergillus and Penicillium Several species are used in the food, biotechnology, and drug industries (e.g., Aspergillus oryzae and Penicillium camemberti), while others are dangerous human and plant pathogens (e.g., Aspergillus fumigatus and Penicillium digitatum). To infer a robust phylogeny and pinpoint poorly resolved branches and their likely underlying contributors, we used 81 genomes spanning the diversity of Aspergillus and Penicillium to construct a 1,668-gene data matrix. Phylogenies of the nucleotide and amino acid versions of this full data matrix as well as of several additional data matrices were generated using three different maximum likelihood schemes (i.e., gene-partitioned, unpartitioned, and coalescence) and using both site-homogenous and site-heterogeneous models (total of 64 species-level phylogenies). Examination of the topological agreement among these phylogenies and measures of internode certainty identified 11/78 (14.1%) bipartitions that were incongruent and pinpointed the likely underlying contributing factors, which included incomplete lineage sorting, hidden paralogy, hybridization or introgression, and reconstruction artifacts associated with poor taxon sampling. Relaxed molecular clock analyses suggest that Aspergillaceae likely originated in the lower Cretaceous and that the Aspergillus and Penicillium genera originated in the upper Cretaceous. Our results shed light on the ongoing debate on Aspergillus systematics and taxonomy and provide a robust evolutionary and temporal framework for comparative genomic analyses in Aspergillaceae More broadly, our approach provides a general template for phylogenomic identification of resolved and contentious branches in densely genome-sequenced lineages across the tree of life.IMPORTANCE Understanding the evolution of traits across technologically and medically significant fungi requires a robust phylogeny. Even though species in the Aspergillus and Penicillium genera (family Aspergillaceae, class Eurotiomycetes) are some of the most significant technologically and medically relevant fungi, we still lack a genome-scale phylogeny of the lineage or knowledge of the parts of the phylogeny that exhibit conflict among analyses. Here, we used a phylogenomic approach to infer evolutionary relationships among 81 genomes that span the diversity of Aspergillus and Penicillium species, to identify conflicts in the phylogeny, and to determine the likely underlying factors of the observed conflicts. Using a data matrix comprised of 1,668 genes, we found that while most branches of the phylogeny of the Aspergillaceae are robustly supported and recovered irrespective of method of analysis, a few exhibit various degrees of conflict among our analyses. Further examination of the observed conflict revealed that it largely stems from incomplete lineage sorting and hybridization or introgression. Our analyses provide a robust and comprehensive evolutionary genomic roadmap for this important lineage, which will facilitate the examination of the diverse technologically and medically relevant traits of these fungi in an evolutionary context.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards.

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

The birth, evolution and death of metabolic gene clusters in fungi.

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

An LaeA- and BrlA-Dependent Cellular Network Governs Tissue-Specific Secondary Metabolism in the Human Pathogen Aspergillus fumigatus.

Biosynthesis of many ecologically important secondary metabolites (SMs) in filamentous fungi is controlled by several global transcriptional regulators, like the chromatin modifier LaeA, and tied to both development and vegetative growth. In Aspergillus molds, asexual development is regulated by the BrlA > AbaA > WetA transcriptional cascade. To elucidate BrlA pathway involvement in SM regulation, we examined the transcriptional and metabolic profiles of ΔbrlA, ΔabaA, and ΔwetA mutant and wild-type strains of the human pathogen Aspergillus fumigatus. We find that BrlA, in addition to regulating production of developmental SMs, regulates vegetative SMs and the SrbA-regulated hypoxia stress response in a concordant fashion to LaeA. We further show that the transcriptional and metabolic equivalence of the ΔbrlA and ΔlaeA mutations is mediated by an LaeA requirement preventing heterochromatic marks in the brlA promoter. These results provide a framework for the cellular network regulating not only fungal SMs but diverse cellular processes linked to virulence of this pathogen. IMPORTANCE Filamentous fungi produce a spectacular variety of small molecules, commonly known as secondary or specialized metabolites (SMs), which are critical to their ecologies and lifestyles (e.g., penicillin, cyclosporine, and aflatoxin). Elucidation of the regulatory network that governs SM production is a major question of both fundamental and applied research relevance. To shed light on the relationship between regulation of development and regulation of secondary metabolism in filamentous fungi, we performed global transcriptomic and metabolomic analyses on mutant and wild-type strains of the human pathogen Aspergillus fumigatus under conditions previously shown to induce the production of both vegetative growth-specific and asexual development-specific SMs. We find that the gene brlA, previously known as a master regulator of asexual development, is also a master regulator of secondary metabolism and other cellular processes. We further show that brlA regulation of SM is mediated by laeA, one of the master regulators of SM, providing a framework for the cellular network regulating not only fungal SMs but diverse cellular processes linked to virulence of this pathogen.

Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species.

Filamentous fungi produce a diverse array of secondary metabolites (SMs) critical for defense, virulence, and communication. The metabolic pathways that produce SMs are found in contiguous gene clusters in fungal genomes, an atypical arrangement for metabolic pathways in other eukaryotes. Comparative studies of filamentous fungal species have shown that SM gene clusters are often either highly divergent or uniquely present in one or a handful of species, hampering efforts to determine the genetic basis and evolutionary drivers of SM gene cluster divergence. Here, we examined SM variation in 66 cosmopolitan strains of a single species, the opportunistic human pathogen Aspergillus fumigatus. Investigation of genome-wide within-species variation revealed 5 general types of variation in SM gene clusters: nonfunctional gene polymorphisms; gene gain and loss polymorphisms; whole cluster gain and loss polymorphisms; allelic polymorphisms, in which different alleles corresponded to distinct, nonhomologous clusters; and location polymorphisms, in which a cluster was found to differ in its genomic location across strains. These polymorphisms affect the function of representative A. fumigatus SM gene clusters, such as those involved in the production of gliotoxin, fumigaclavine, and helvolic acid as well as the function of clusters with undefined products. In addition to enabling the identification of polymorphisms, the detection of which requires extensive genome-wide synteny conservation (e.g., mobile gene clusters and nonhomologous cluster alleles), our approach also implicated multiple underlying genetic drivers, including point mutations, recombination, and genomic deletion and insertion events as well as horizontal gene transfer from distant fungi. Finally, most of the variants that we uncover within A. fumigatus have been previously hypothesized to contribute to SM gene cluster diversity across entire fungal classes and phyla. We suggest that the drivers of genetic diversity operating within a fungal species shown here are sufficient to explain SM cluster macroevolutionary patterns.

Nanoscale architecture of the Schizosaccharomyces pombe contractile ring.

The contractile ring is a complex molecular apparatus which physically divides many eukaryotic cells. Despite knowledge of its protein composition, the molecular architecture of the ring is not known. Here we have applied super-resolution microscopy and FRET to determine the nanoscale spatial organization of Schizosaccharomyces pombe contractile ring components relative to the plasma membrane. Similar to other membrane-tethered actin structures, we find proteins localize in specific layers relative to the membrane. The most membrane-proximal layer (0-80 nm) is composed of membrane-binding scaffolds, formin, and the tail of the essential myosin-II. An intermediate layer (80-160 nm) consists of a network of cytokinesis accessory proteins as well as multiple signaling components which influence cell division. Farthest from the membrane (160-350 nm) we find F-actin, the motor domains of myosins, and a major F-actin crosslinker. Circumferentially within the ring, multiple proteins proximal to the membrane form clusters of different sizes, while components farther from the membrane are uniformly distributed. This comprehensive organizational map provides a framework for understanding contractile ring function.

An alternative pathway to eusociality: Exploring the molecular and functional basis of fortress defense.

Some animals express a form of eusociality known as "fortress defense," in which defense rather than brood care is the primary social act. Aphids are small plant-feeding insects, but like termites, some species express division of labor and castes of aggressive juvenile "soldiers." What is the functional basis of fortress defense eusociality in aphids? Previous work showed that the acquisition of venoms might be a key innovation in aphid social evolution. We show that the lethality of aphid soldiers derives in part from the induction of exaggerated immune responses in insects they attack. Comparisons between closely related social and nonsocial species identified a number of secreted effector molecules that are candidates for immune modulation, including a convergently recruited protease described in unrelated aphid species with venom-like functions. These results suggest that aphids are capable of antagonizing conserved features of the insect immune response, and provide new insights into the mechanisms underlying the evolution of fortress defense eusociality in aphids.

Regulation of Secondary Metabolism by the Velvet Complex Is Temperature-Responsive in Aspergillus.

Sensing and responding to environmental cues is critical to the lifestyle of filamentous fungi. How environmental variation influences fungi to produce a wide diversity of ecologically important secondary metabolites (SMs) is not well understood. To address this question, we first examined changes in global gene expression of the opportunistic human pathogen, Aspergillus fumigatus, after exposure to different temperature conditions. We found that 11 of the 37 SM gene clusters in A. fumigatus were expressed at higher levels at 30° than at 37°. We next investigated the role of the light-responsive Velvet complex in environment-dependent gene expression by examining temperature-dependent transcription profiles in the absence of two key members of the Velvet protein complex, VeA and LaeA We found that the 11 temperature-regulated SM gene clusters required VeA at 37° and LaeA at both 30 and 37° for wild-type levels of expression. Interestingly, four SM gene clusters were regulated by VeA at 37° but not at 30°, and two additional ones were regulated by VeA at both temperatures but were substantially less so at 30°, indicating that the role of VeA and, more generally of the Velvet complex, in the regulation of certain SM gene clusters is temperature-dependent. Our findings support the hypothesis that fungal secondary metabolism is regulated by an intertwined network of transcriptional regulators responsive to multiple environmental factors.

Examining the evolution of the regulatory circuit controlling secondary metabolism and development in the fungal genus Aspergillus.

Filamentous fungi produce diverse secondary metabolites (SMs) essential to their ecology and adaptation. Although each SM is typically produced by only a handful of species, global SM production is governed by widely conserved transcriptional regulators in conjunction with other cellular processes, such as development. We examined the interplay between the taxonomic narrowness of SM distribution and the broad conservation of global regulation of SM and development in Aspergillus, a diverse fungal genus whose members produce well-known SMs such as penicillin and gliotoxin. Evolutionary analysis of the 2,124 genes comprising the 262 SM pathways in four Aspergillus species showed that most SM pathways were species-specific, that the number of SM gene orthologs was significantly lower than that of orthologs in primary metabolism, and that the few conserved SM orthologs typically belonged to non-homologous SM pathways. RNA sequencing of two master transcriptional regulators of SM and development, veA and mtfA, showed that the effects of deletion of each gene, especially veA, on SM pathway regulation were similar in A. fumigatus and A. nidulans, even though the underlying genes and pathways regulated in each species differed. In contrast, examination of the role of these two regulators in development, where 94% of the underlying genes are conserved in both species showed that whereas the role of veA is conserved, mtfA regulates development in the homothallic A. nidulans but not in the heterothallic A. fumigatus. Thus, the regulation of these highly conserved developmental genes is divergent, whereas-despite minimal conservation of target genes and pathways-the global regulation of SM production is largely conserved. We suggest that the evolution of the transcriptional regulation of secondary metabolism in Aspergillus represents a novel type of regulatory circuit rewiring and hypothesize that it has been largely driven by the dramatic turnover of the target genes involved in the process.

Variable suites of non-effector genes are co-regulated in the type III secretion virulence regulon across the Pseudomonas syringae phylogeny.

Pseudomonas syringae is a phylogenetically diverse species of Gram-negative bacterial plant pathogens responsible for crop diseases around the world. The HrpL sigma factor drives expression of the major P. syringae virulence regulon. HrpL controls expression of the genes encoding the structural and functional components of the type III secretion system (T3SS) and the type three secreted effector proteins (T3E) that are collectively essential for virulence. HrpL also regulates expression of an under-explored suite of non-type III effector genes (non-T3E), including toxin production systems and operons not previously associated with virulence. We implemented and refined genome-wide transcriptional analysis methods using cDNA-derived high-throughput sequencing (RNA-seq) data to characterize the HrpL regulon from six isolates of P. syringae spanning the diversity of the species. Our transcriptomes, mapped onto both complete and draft genomes, significantly extend earlier studies. We confirmed HrpL-regulation for a majority of previously defined T3E genes in these six strains. We identified two new T3E families from P. syringae pv. oryzae 1_6, a strain within the relatively underexplored phylogenetic Multi-Locus Sequence Typing (MLST) group IV. The HrpL regulons varied among strains in gene number and content across both their T3E and non-T3E gene suites. Strains within MLST group II consistently express the lowest number of HrpL-regulated genes. We identified events leading to recruitment into, and loss from, the HrpL regulon. These included gene gain and loss, and loss of HrpL regulation caused by group-specific cis element mutations in otherwise conserved genes. Novel non-T3E HrpL-regulated genes include an operon that we show is required for full virulence of P. syringae pv. phaseolicola 1448A on French bean. We highlight the power of integrating genomic, transcriptomic, and phylogenetic information to drive concise functional experimentation and to derive better insight into the evolution of virulence across an evolutionarily diverse pathogen species.

The fumagillin gene cluster, an example of hundreds of genes under veA control in Aspergillus fumigatus.

Aspergillus fumigatus is the causative agent of invasive aspergillosis, leading to infection-related mortality in immunocompromised patients. We previously showed that the conserved and unique-to-fungi veA gene affects different cell processes such as morphological development, gliotoxin biosynthesis and protease activity, suggesting a global regulatory effect on the genome of this medically relevant fungus. In this study, RNA sequencing analysis revealed that veA controls the expression of hundreds of genes in A. fumigatus, including those comprising more than a dozen known secondary metabolite gene clusters. Chemical analysis confirmed that veA controls the synthesis of other secondary metabolites in this organism in addition to gliotoxin. Among the secondary metabolite gene clusters regulated by veA is the elusive but recently identified gene cluster responsible for the biosynthesis of fumagillin, a meroterpenoid known for its anti-angiogenic activity by binding to human methionine aminopeptidase 2. The fumagillin gene cluster contains a veA-dependent regulatory gene, fumR (Afu8g00420), encoding a putative C6 type transcription factor. Deletion of fumR results in silencing of the gene cluster and elimination of fumagillin biosynthesis. We found expression of fumR to also be dependent on laeA, a gene encoding another component of the fungal velvet complex. The results in this study argue that veA is a global regulator of secondary metabolism in A. fumigatus, and that veA may be a conduit via which chemical development is coupled to morphological development and other cellular processes.