Search

The extensive similarities between the genomes of human and model organisms are the foundation of much of modern biology, with model organism experimentation permitting valuable insights into biological function and the aetiology of human disease. In contrast, differences among genomes have received less attention. Yet these can be expected to govern the physiological and morphological distinctions apparent among species, especially if such differences are the result of evolutionary adaptation. A recent comparison of the draft sequences of mouse and human genomes has shed light on the selective forces that have predominated in their recent evolutionary histories. In particular, mouse-specific clusters of homologues associated with roles in reproduction, immunity and host defence appear to be under diversifying positive selective pressure, as indicated by high ratios of non-synonymous to synonymous substitution rates. These clusters are also frequently punctuated by homologous pseudogenes. They thus have experienced numerous gene death, as well as gene birth, events. These regions appear, therefore, to have borne the brunt of adaptive evolution that underlies physiological and behavioural innovation in mice. We predict that the availability of numerous animal genomes will give rise to a new field of genome zoology in which differences in animal physiology and ethology are illuminated by the study of genomic sequence variations.

Metals are important cofactors in the metabolic processes of cyanobacteria, including photosynthesis, cellular respiration, DNA replication, and the biosynthesis of primary and secondary metabolites. In adaptation to the marine environment, cyanobacteria use metallophores to acquire trace metals when necessary as well as to reduce potential toxicity from excessive metal concentrations. Leptochelins A-C were identified as structurally novel metallophores from three geographically dispersed cyanobacteria of the genus Leptothoe. Determination of the complex structures of these metabolites presented numerous challenges, but they were ultimately solved using integrated data from NMR, mass spectrometry and deductions from the biosynthetic gene cluster. The leptochelins are comprised of halogenated linear NRPS-PKS hybrid products with multiple heterocycles that have potential for hexadentate and tetradentate coordination with metal ions. The genomes of the three leptochelin producers were sequenced, and retrobiosynthetic analysis revealed one candidate biosynthetic gene cluster (BGC) consistent with the structure of leptochelin. The putative BGC is highly homologous in all three Leptothoe strains, and all possess genetic signatures associated with metallophores. Postcolumn infusion of metals using an LC-MS metabolomics workflow performed with leptochelins A and B revealed promiscuous binding of iron, copper, cobalt, and zinc, with greatest preference for copper. Iron depletion and copper toxicity experiments support the hypothesis that leptochelin metallophores may play key ecological roles in iron acquisition and in copper detoxification. In addition, the leptochelins possess significant cytotoxicity against several cancer cell lines.

The tropical marine cyanobacterium Moorena producens JHB is a prolific source of secondary metabolites with potential biomedical utility. Previous studies on this strain led to the discovery of several novel compounds such as hectochlorins and jamaicamides. However, bioinformatic analyses of its genome indicate the presence of numerous cryptic biosynthetic gene clusters that have yet to be characterized. To potentially stimulate the production of novel compounds from this strain, it was cocultured with Candida albicans. From this experiment, we observed the increased production of a new compound that we characterize here as hectoramide B. Bioinformatic analysis of the M. producens JHB genome enabled the identification of a putative biosynthetic gene cluster responsible for hectoramide B biosynthesis. This work demonstrates that coculture competition experiments can be a valuable method to facilitate the discovery of novel natural products from cyanobacteria.

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

Azaserine is a bacterial metabolite containing a biologically unusual and synthetically enabling α-diazoester functional group. Herein, we report the discovery of the azaserine (aza) biosynthetic gene cluster from Glycomyces harbinensis. Discovery of related gene clusters reveals previously unappreciated azaserine producers, and heterologous expression of the aza gene cluster confirms its role in azaserine assembly. Notably, this gene cluster encodes homologues of hydrazonoacetic acid (HYAA)-producing enzymes, implicating HYAA in α-diazoester biosynthesis. Isotope feeding and biochemical experiments support this hypothesis. These discoveries indicate that a 2-electron oxidation of a hydrazonoacetyl intermediate is required for α-diazoester formation, constituting a distinct logic for diazo biosynthesis. Uncovering this biological route for α-diazoester synthesis now enables the production of a highly versatile carbene precursor in cells, facilitating approaches for engineering complete carbene-mediated biosynthetic transformations in vivo.

With an ever-increasing amount of (meta)genomic data being deposited in sequence databases, (meta)genome mining for natural product biosynthetic pathways occupies a critical role in the discovery of novel pharmaceutical drugs, crop protection agents and biomaterials. The genes that encode these pathways are often organised into biosynthetic gene clusters (BGCs). In 2015, we defined the Minimum Information about a Biosynthetic Gene cluster (MIBiG): a standardised data format that describes the minimally required information to uniquely characterise a BGC. We simultaneously constructed an accompanying online database of BGCs, which has since been widely used by the community as a reference dataset for BGCs and was expanded to 2021 entries in 2019 (MIBiG 2.0). Here, we describe MIBiG 3.0, a database update comprising large-scale validation and re-annotation of existing entries and 661 new entries. Particular attention was paid to the annotation of compound structures and biological activities, as well as protein domain selectivities. Together, these new features keep the database up-to-date, and will provide new opportunities for the scientific community to use its freely available data, e.g. for the training of new machine learning models to predict sequence-structure-function relationships for diverse natural products. MIBiG 3.0 is accessible online at https://mibig.secondarymetabolites.org/.

Poly zinc finger proteins (ZFP) that contain a KRAB (Krüppel-associated box) domain represent the largest class of transcription factors in higher organisms, but their roles in development and pathogenesis are largely undefined. ZFP14 (also known as ZNF531) contains thirteen zinc fingers and is highly conserved across species. Notably, we found that ZFP14 is frequently down-regulated in a multitude of human cancers, which correlates with poor prognosis of patients. Since ZFP14 has never been characterized, we generated a Zfp14-deficient mouse model to investigate the role of ZFP14 in development and pathogenesis. We showed that the mice deficient in Zfp14 had a short lifespan and were prone to diffuse large B-cell lymphoma (DLBCL), hyperplasia in multiple organs, systemic chronic inflammation, liver steatosis, and pancreatitis. Additionally, several pro-inflammatory cytokines, including IL-1β, IL18, and TNFα, were highly expressed in inflamed Zfp14-/- mice spleens, livers, kidneys and lungs. To determine the underlying mechanism, RNA-seq analysis was performed and showed that the loss of ZFP14 led to increased expression of inflammatory and tumor-promoting genes. Out of the various tumor-promoting genes upregulated by ZFP14 loss, the HOXA gene cluster, which is known to promote lymphomagenesis and conserved between mouse and human, is consistently induced by loss of ZFP14. Moreover, we showed that the HOXA gene expression was inversely correlated with that of ZFP14 in human cancer patients and higher HOXA1 expression was correlated with poor patient prognosis. Together, we postulate that ZFP14 suppresses lymphomagenesis and abnormal inflammatory response by maintaining proper expression of the HOXA gene cluster.

Transcription factors (TFs) bind to the promoters of target genes to regulate gene expression in response to different stimuli. The functions and regulatory mechanisms of transcription factors (TFs) in Verticillium dahliae are, however, still largely unclear. This study showed that a C2H2-type zinc finger TF, VdCf2 (V. dahliae chorion transcription factor 2), plays key roles in V. dahliae growth, melanin production, and virulence. Transcriptome sequencing analysis showed that VdCf2 was involved in the regulation of expression of genes encoding secreted proteins, pathogen-host interaction (PHI) homologs, TFs, and G protein-coupled receptors (GPCRs). Furthermore, VdCf2 positively regulated the expression of VdPevD1 (VDAG_02735), a previously reported virulence factor. VdCf2 thus regulates the expression of several pathogenicity-related genes that also contribute to virulence in V. dahliae. VdCf2 also inhibited the transcription of the Vd276-280 gene cluster and interacted with two members encoding proteins (VDAG_07276 and VDAG_07278) in the gene cluster. IMPORTANCE Verticillium dahliae is an important soilborne phytopathogen which can ruinously attack numerous host plants and cause significant economic losses. Transcription factors (TFs) were reported to be involved in various biological processes, such as hyphal growth and virulence of pathogenic fungi. However, the functions and regulatory mechanisms of TFs in V. dahliae remain largely unclear. In this study, we identified a new transcription factor, VdCf2 (V. dahliae chorion transcription factor 2), based on previous transcriptome data, which participates in growth, melanin production, and virulence of V. dahliae. We provide evidence that VdCf2 regulates the expression of the pathogenicity-related gene VdPevD1 (VDAG_02735) and Vd276-280 gene cluster. VdCf2 also interacts with VDAG_07276 and VDAG_07278 in this gene cluster based on a yeast two-hybrid and bimolecular fluorescence complementation assay. These results revealed the regulatory mechanisms of a pivotal pathogenicity-related transcription factor, VdCf2 in V. dahliae.

Genome sequencing of filamentous fungi has demonstrated that most secondary metabolite biosynthetic gene clusters (BGCs) are silent under standard laboratory conditions. In this work, we have established an in vitro CRISPR-Cas9 system in Aspergillus wentii. To activate otherwise silent BGCs, we deleted the negative transcriptional regulator mcrA. Deletion of mcrA (mcrAΔ) resulted in differential production of 17 SMs in total when the strain was cultivated on potato dextrose media (PDA). Nine out of fifteen of these SMs were fully characterized, including emodin (1), physcion (2), sulochrin (3), physcion bianthrone (4), 14-O-demethylsulochrin (5), (trans/cis)-emodin bianthrone (6 and 7), and (trans/cis)-emodin physcion bianthrone (8 and 9). These compounds were all found to be produced by the same polyketide synthase (PKS) BGC. We then performed a secondary knockout targeting this PKS cluster in the mcrAΔ background. The metabolite profile of the dual-knockout strain revealed new metabolites that were not previously detected in the mcrAΔ parent strain. Two additional SMs were purified from the dual-knockout strain and were characterized as aspergillus acid B (16) and a structurally related but previously unidentified compound (17). For the first time, this work presents a facile genetic system capable of targeted gene editing in A. wentii. This work also illustrates the utility of performing a dual knockout to eliminate major metabolic products, enabling additional SM discovery.

Harmful cyanobacterial blooms (cyanoHABs) cause recurrent toxic events in global watersheds. Although public health agencies monitor the causal toxins of most cyanoHABs and scientists in the field continue developing precise detection and prediction tools, the potent anticholinesterase neurotoxin, guanitoxin, is not presently environmentally monitored. This is largely due to its incompatibility with widely employed analytical methods and instability in the environment, despite guanitoxin being among the most lethal cyanotoxins. Here, we describe the guanitoxin biosynthesis gene cluster and its rigorously characterized nine-step metabolic pathway from l-arginine in the cyanobacterium Sphaerospermopsis torques-reginae ITEP-024. Through environmental sequencing data sets, guanitoxin (gnt) biosynthetic genes are repeatedly detected and expressed in municipal freshwater bodies that have undergone past toxic events. Knowledge of the genetic basis of guanitoxin biosynthesis now allows for environmental, biosynthetic gene monitoring to establish the global scope of this neurotoxic organophosphate.

With the advent of genome sequencing and mining technologies, secondary metabolite biosynthetic gene clusters (BGCs) within bacterial genomes are becoming easier to predict. For subsequent BGC characterization, clustered regularly interspaced short palindromic repeats (CRISPR) has contributed to knocking out target genes and/or modulating their expression; however, CRISPR is limited to strains for which robust genetic tools are available. Here we present a strategy that combines CRISPR with chassis-independent recombinase-assisted genome engineering (CRAGE), which enables CRISPR systems in diverse bacteria. To demonstrate CRAGE-CRISPR, we select 10 polyketide/non-ribosomal peptide BGCs in Photorhabdus luminescens as models and create their deletion and activation mutants. Subsequent loss- and gain-of-function studies confirm 22 secondary metabolites associated with the BGCs, including a metabolite from a previously uncharacterized BGC. These results demonstrate that the CRAGE-CRISPR system is a simple yet powerful approach to rapidly perturb expression of defined BGCs and to profile genotype-phenotype relationships in bacteria.

Domoic acid (DA), the causative agent of amnesic shellfish poisoning, is produced by select organisms within two distantly related algal clades: planktonic diatoms and red macroalgae. The biosynthetic pathway to isodomoic acid A was recently solved in the harmful algal bloom-forming diatom Pseudonitzschia multiseries, establishing the genetic basis for the global production of this potent neurotoxin. Herein, we sequenced the 507-Mb genome of Chondria armata, the red macroalgal seaweed from which DA was first isolated in the 1950s, identifying several copies of the red algal DA (rad) biosynthetic gene cluster. The rad genes are organized similarly to the diatom DA biosynthesis cluster in terms of gene synteny, including a cytochrome P450 (CYP450) enzyme critical to DA production that is notably absent in red algae that produce the simpler kainoid neurochemical, kainic acid. The biochemical characterization of the N-prenyltransferase (RadA) and kainoid synthase (RadC) enzymes support a slightly altered DA biosynthetic model in C. armata via the congener isodomoic acid B, with RadC behaving more like the homologous diatom enzyme despite higher amino acid similarity to red algal kainic acid synthesis enzymes. A phylogenetic analysis of the rad genes suggests unique origins for the red macroalgal and diatom genes in their respective hosts, with native eukaryotic CYP450 neofunctionalization combining with the horizontal gene transfer of N-prenyltransferases and kainoid synthases to establish DA production within the algal lineages.

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Streptomycetes are highly metabolically gifted bacteria with the abilities to produce bioproducts that have profound economic and societal importance. These bioproducts are produced by metabolic pathways including those for the biosynthesis of secondary metabolites and catabolism of plant biomass constituents. Advancements in genome sequencing technologies have revealed a wealth of untapped metabolic potential from Streptomyces genomes. Here, we report the largest Streptomyces pangenome generated by using 205 complete genomes. Metabolic potentials of the pangenome and individual genomes were analyzed, revealing degrees of conservation of individual metabolic pathways and strains potentially suitable for metabolic engineering. Of them, Streptomyces bingchenggensis was identified as a potent degrader of plant biomass. Polyketide, non-ribosomal peptide, and gamma-butyrolactone biosynthetic enzymes are primarily strain specific while ectoine and some terpene biosynthetic pathways are highly conserved. A large number of transcription factors associated with secondary metabolism are strain-specific while those controlling basic biological processes are highly conserved. Although the majority of genes involved in morphological development are highly conserved, there are strain-specific varieties which may contribute to fine tuning the timing of cellular differentiation. Overall, these results provide insights into the metabolic potential, regulation and physiology of streptomycetes, which will facilitate further exploitation of these important bacteria.

The Antarctic marine ecosystem harbors a wealth of biological and chemical innovation that has risen in concert over millennia since the isolation of the continent and formation of the Antarctic circumpolar current. Scientific inquiry into the novelty of marine natural products produced by Antarctic benthic invertebrates led to the discovery of a bioactive macrolide, palmerolide A, that has specific activity against melanoma and holds considerable promise as an anticancer therapeutic. While this compound was isolated from the Antarctic ascidian Synoicum adareanum, its biosynthesis has since been hypothesized to be microbially mediated, given structural similarities to microbially produced hybrid nonribosomal peptide-polyketide macrolides. Here, we describe a metagenome-enabled investigation aimed at identifying the biosynthetic gene cluster (BGC) and palmerolide A-producing organism. A 74-kbp candidate BGC encoding the multimodular enzymatic machinery (hybrid type I-trans-AT polyketide synthase-nonribosomal peptide synthetase and tailoring functional domains) was identified and found to harbor key features predicted as necessary for palmerolide A biosynthesis. Surveys of ascidian microbiome samples targeting the candidate BGC revealed a high correlation between palmerolide gene targets and a single 16S rRNA gene variant (R = 0.83 to 0.99). Through repeated rounds of metagenome sequencing followed by binning contigs into metagenome-assembled genomes, we were able to retrieve a nearly complete genome (10 contigs) of the BGC-producing organism, a novel verrucomicrobium within the Opitutaceae family that we propose here as "Candidatus Synoicihabitans palmerolidicus." The refined genome assembly harbors five highly similar BGC copies, along with structural and functional features that shed light on the host-associated nature of this unique bacterium. IMPORTANCE Palmerolide A has potential as a chemotherapeutic agent to target melanoma. We interrogated the microbiome of the Antarctic ascidian, Synoicum adareanum, using a cultivation-independent high-throughput sequencing and bioinformatic strategy. The metagenome-encoded biosynthetic machinery predicted to produce palmerolide A was found to be associated with the genome of a member of the S. adareanum core microbiome. Phylogenomic analysis suggests the organism represents a new deeply branching genus, "Candidatus Synoicihabitans palmerolidicus," in the Opitutaceae family of the Verrucomicrobia phylum. The Ca. Synoicihabitans palmerolidicus 4.29-Mb genome encodes a repertoire of carbohydrate-utilizing and transport pathways, a chemotaxis system, flagellar biosynthetic capacity, and other regulatory elements enabling its ascidian-associated lifestyle. The palmerolide producer's genome also contains five distinct copies of the large palmerolide biosynthetic gene cluster that may provide structural complexity of palmerolide variants.

While specialized metabolites are thought to mediate ecological interactions, the evolutionary processes driving chemical diversification, particularly among closely related lineages, remain poorly understood. Here, we examine the evolutionary dynamics governing the distribution of natural product biosynthetic gene clusters (BGCs) among 118 strains representing all nine currently named species of the marine actinobacterial genus Salinispora. While much attention has been given to the role of horizontal gene transfer (HGT) in structuring BGC distributions, we find that vertical descent facilitates interspecies BGC diversification over evolutionary timescales. Moreover, we identified a distinct phylogenetic signal among Salinispora species at both the BGC and metabolite level, indicating that specialized metabolism represents a conserved phylogenetic trait. Using a combination of genomic analyses and liquid chromatography-high-resolution tandem mass spectrometry (LC-MS/MS) targeting nine experimentally characterized BGCs and their small molecule products, we identified gene gain/loss events, constrained interspecies recombination, and other evolutionary processes associated with vertical inheritance as major contributors to BGC diversification. These evolutionary dynamics had direct consequences for the compounds produced, as exemplified by species-level differences in salinosporamide production. Together, our results support the concept that specialized metabolites, and their cognate BGCs, can represent phylogenetically conserved functional traits with chemical diversification proceeding in species-specific patterns over evolutionary time frames. IMPORTANCE Microbial natural products are traditionally exploited for their pharmaceutical potential, yet our understanding of the evolutionary processes driving BGC evolution and compound diversification remain poorly developed. While HGT is recognized as an integral driver of BGC distributions, we find that the effects of vertical inheritance on BGC diversification had direct implications for species-level specialized metabolite production. As such, understanding the degree of genetic variation that corresponds to species delineations can enhance natural product discovery efforts. Resolving the evolutionary relationships between closely related strains and specialized metabolism can also facilitate our understanding of the ecological roles of small molecules in structuring the environmental distribution of microbes.

BackgroundGenetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human brain.MethodsRNA-seq data from 783 human brain samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples.ResultsIn the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A, a brain-wide synchronized 5N to 5A transition occurs between 24 post-conceptual weeks (2nd trimester) and 6 years of age. In mice, the equivalent 5N to 5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories.ConclusionsExon usage in SCN1A, SCN2A, SCN3A, and SCN8A changes dramatically during human brain development. These splice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

TMEM132D is a human gene identified with multiple risk alleles for panic disorders, anxiety and major depressive disorders. Defining a conserved family of transmembrane proteins, TMEM132D and its homologs are still of unknown molecular functions. By generating loss-of-function mutants of the sole TMEM132 ortholog in C. elegans, we identify abnormal morphologic phenotypes in the dopaminergic PDE neurons. Using a yeast two-hybrid screen, we find that NAP1 directly interacts with the cytoplasmic domain of human TMEM132D, and mutations in C. elegans tmem-132 that disrupt interaction with NAP1 cause similar morphologic defects in the PDE neurons. NAP1 is a component of the WAVE regulatory complex (WRC) that controls F-actin cytoskeletal dynamics. Decreasing activity of WRC rescues the PDE defects in tmem-132 mutants, whereas gain-of-function of TMEM132D in mammalian cells inhibits WRC, leading to decreased abundance of select WRC components, impaired actin nucleation and cell motility. We propose that metazoan TMEM132 family proteins play evolutionarily conserved roles in regulating NAP1 protein homologs to restrict inappropriate WRC activity, cytoskeletal and morphologic changes in the cell.

The continuous increase in global population prompts increased wheat production. Future wheat (Triticum aestivum L.) breeding will heavily rely on dissecting molecular and genetic bases of wheat yield and related traits which is possible through the discovery of quantitative trait loci (QTLs) in constructed populations, such as recombinant inbred lines (RILs). Here, we present an evaluation of 92 RILs in a bi-parental RIL mapping population (the International Triticeae Mapping Initiative Mapping Population [ITMI/MP]) using newly generated phenotypic data in 3-year experiments (2015), older phenotypic data (1997-2009), and newly created single nucleotide polymorphism (SNP) marker data based on 92 of the original RILs to search for novel and stable QTLs. Our analyses of more than 15 unique traits observed in multiple experiments included analyses of 46 traits in three environments in the USA, 69 traits in eight environments in Germany, 149 traits in 10 environments in Russia, and 28 traits in four environments in India (292 traits in 25 environments) with 7584 SNPs (292 × 7584 = 2 214 528 data points). A total of 874 QTLs were detected with limit of detection (LOD) scores of 2.01-3.0 and 432 QTLs were detected with LOD > 3.0. Moreover, 769 QTLs could be assigned to 183 clusters based on the common markers and relative proximity of related QTLs, indicating gene-rich regions throughout the A, B, and D genomes of common wheat. This upgraded genotype-phenotype information of ITMI/MP can assist breeders and geneticists who can make crosses with suitable RILs to improve or investigate traits of interest.

RNA, DNA, and protein molecules are highly organized within three-dimensional (3D) structures in the nucleus. Although RNA has been proposed to play a role in nuclear organization, exploring this has been challenging because existing methods cannot measure higher-order RNA and DNA contacts within 3D structures. To address this, we developed RNA & DNA SPRITE (RD-SPRITE) to comprehensively map the spatial organization of RNA and DNA. These maps reveal higher-order RNA-chromatin structures associated with three major classes of nuclear function: RNA processing, heterochromatin assembly, and gene regulation. These data demonstrate that hundreds of ncRNAs form high-concentration territories throughout the nucleus, that specific RNAs are required to recruit various regulators into these territories, and that these RNAs can shape long-range DNA contacts, heterochromatin assembly, and gene expression. These results demonstrate a mechanism where RNAs form high-concentration territories, bind to diffusible regulators, and guide them into compartments to regulate essential nuclear functions.

Fungi produce secondary metabolites that are not directly involved in their growth, but often contribute to their adaptation to extreme environmental stimuli and enable their survival. Conidial pigment or melanin is one of the secondary metabolites produced naturally by a polyketide synthesis (PKS) gene cluster in several filamentous fungi and is known to protect these fungi from extreme radiation conditions. Several pigmented or melanized fungi have been shown to grow under extreme radiation conditions at the Chernobyl nuclear accident site. Some of these fungi, including Paecilomyces variotii, were observed to grow towards the source of radiation. Therefore, in this study, we wanted to identify if the pigment produced by P. variotii, contributes to providing protection against radiation condition. We first identified the PKS gene responsible for synthesis of pigment in P. variotii and confirmed its role in providing protection against UV irradiation through CRISPR-Cas9 mediated gene deletion. This is the first report that describes the use of CRISPR methodology to create gene deletions in P. variotii. Further, we showed that the pigment produced by this fungus, was not inhibited by DHN-melanin pathway inhibitors, indicating that the fungus does not produce melanin. We then identified the pigment synthesized by the PKS gene of P. variotii, as a naptho-pyrone Ywa1, by heterologously expressing the gene in Aspergillus nidulans. The results obtained will further aid in understanding the mechanistic basis of radiation resistance.

Transposable elements (TEs) are mobile sequences that engender widespread mutations and thus are a major hazard that must be silenced. The most abundant active class of TEs in mammalian genomes is long interspersed element class 1 (LINE1). Here, we report that LINE1 transposition is suppressed in the male germline by transcription factors encoded by a rapidly evolving X-linked homeobox gene cluster. LINE1 transposition is repressed by many members of this RHOX transcription factor family, including those with different patterns of expression during spermatogenesis. One family member-RHOX10-suppresses LINE1 transposition during fetal development in vivo when the germline would otherwise be susceptible to LINE1 activation because of epigenetic reprogramming. We provide evidence that RHOX10 suppresses LINE transposition by inducing Piwil2, which encodes a key component in the Piwi-interacting RNA pathway that protects against TEs. The ability of RHOX transcription factors to suppress LINE1 is conserved in humans but is lost in RHOXF2 mutants from several infertile human patients, raising the possibility that loss of RHOXF2 causes human infertility by allowing uncontrolled LINE1 expression in the germline. Together, our results support a model in which the Rhox gene cluster is in an evolutionary arms race with TEs, resulting in expansion of the Rhox gene cluster to suppress TEs in different biological contexts.

While thousands of environmental metagenomes have been mined for the presence of novel biosynthetic gene clusters, such computational predictions do not provide evidence of their in vivo biosynthetic functionality. Using fluorescent in situ enzyme assay targeting carrier proteins common to polyketide (PKS) and nonribosomal peptide synthetases (NRPS), we applied fluorescence-activated cell sorting to tunicate microbiome to enrich for microbes with active secondary metabolic capabilities. Single-cell genomics uncovered the genetic basis for a wide biosynthetic diversity in the enzyme-active cells and revealed a member of marine Oceanospirillales harboring a novel NRPS gene cluster with high similarity to phylogenetically distant marine and terrestrial bacteria. Interestingly, this synthase belongs to a larger class of siderophore biosynthetic gene clusters commonly associated with pestilence and disease. This demonstrates activity-guided single-cell genomics as a tool to guide novel biosynthetic discovery.

Bacteria communicate by small-molecule chemicals that facilitate intra- and inter-species interactions. These extracellular signalling molecules mediate diverse processes including virulence, bioluminescence, biofilm formation, motility and specialized metabolism. The signalling molecules produced by members of the phylum Actinobacteria generally comprise γ-butyrolactones, γ-butenolides and furans. The best-known actinomycete γ-butyrolactone is A-factor, which triggers specialized metabolism and morphological differentiation in the genus Streptomyces . Salinipostins A–K are unique γ-butyrolactone molecules with rare phosphotriester moieties that were recently characterized from the marine actinomycete genus Salinispora . The production of these compounds has been linked to the nine-gene biosynthetic gene cluster (BGC) spt. Critical to salinipostin assembly is the γ-butyrolactone synthase encoded by spt9. Here, we report the surprising distribution of spt9 homologues across 12 bacterial phyla, the majority of which are not known to produce γ-butyrolactones. Further analyses uncovered a large group of spt-like gene clusters outside of the genus Salinispora , suggesting the production of new salinipostin-like diversity. These gene clusters show evidence of horizontal transfer and location-specific recombination among Salinispora strains. The results suggest that γ-butyrolactone production may be more widespread than previously recognized. The identification of new γ-butyrolactone BGCs is the first step towards understanding the regulatory roles of the encoded small molecules in Actinobacteria.

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates-namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose-methanol-choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases-we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq) creates new opportunities to dissect cell type-specific mechanisms of complex diseases. Since pancreatic islets are central to type 2 diabetes (T2D), we profiled 15,298 islet cells by using combinatorial barcoding snATAC-seq and identified 12 clusters, including multiple alpha, beta and delta cell states. We cataloged 228,873 accessible chromatin sites and identified transcription factors underlying lineage- and state-specific regulation. We observed state-specific enrichment of fasting glucose and T2D genome-wide association studies for beta cells and enrichment for other endocrine cell types. At T2D signals localized to islet-accessible chromatin, we prioritized variants with predicted regulatory function and co-accessibility with target genes. A causal T2D variant rs231361 at the KCNQ1 locus had predicted effects on a beta cell enhancer co-accessible with INS and genome editing in embryonic stem cell-derived beta cells affected INS levels. Together our findings demonstrate the power of single-cell epigenomics for interpreting complex disease genetics.

The genus Aspergillus encompasses human pathogens such as Aspergillus fumigatus and industrial powerhouses such as Aspergillus niger In both cases, Aspergillus biofilms have consequences for infection outcomes and yields of economically important products. However, the molecular components influencing filamentous fungal biofilm development, structure, and function remain ill defined. Macroscopic colony morphology is an indicator of underlying biofilm architecture and fungal physiology. A hypoxia-locked colony morphotype of A. fumigatus has abundant colony furrows that coincide with a reduction in vertically oriented hyphae within biofilms and increased low oxygen growth and virulence. Investigation of this morphotype has led to the identification of the causative gene, biofilm architecture factor A (bafA), a small cryptic open reading frame within a subtelomeric gene cluster. BafA is sufficient to induce the hypoxia-locked colony morphology and biofilm architecture in A. fumigatus Analysis across a large population of A. fumigatus isolates identified a larger family of baf genes, all of which have the capacity to modulate hyphal architecture, biofilm development, and hypoxic growth. Furthermore, introduction of A. fumigatusbafA into A. niger is sufficient to generate the hypoxia-locked colony morphology, biofilm architecture, and increased hypoxic growth. Together, these data indicate the potential broad impacts of this previously uncharacterized family of small genes to modulate biofilm architecture and function in clinical and industrial settings.IMPORTANCE The manipulation of microbial biofilms in industrial and clinical applications remains a difficult task. The problem is particularly acute with regard to filamentous fungal biofilms for which molecular mechanisms of biofilm formation, maintenance, and function are only just being elucidated. Here, we describe a family of small genes heterogeneously expressed across Aspergillus fumigatus strains that are capable of modifying colony biofilm morphology and microscopic hyphal architecture. Specifically, these genes are implicated in the formation of a hypoxia-locked colony morphotype that is associated with increased virulence of A. fumigatus Synthetic introduction of these gene family members, here referred to as biofilm architecture factors, in both A. fumigatus and A. niger additionally modulates low oxygen growth and surface adherence. Thus, these genes are candidates for genetic manipulation of biofilm development in aspergilli.

The DEFECTIVE EMBRYO AND MERISTEMS 1 (DEM1) gene encodes a protein of unknown biochemical function required for meristem formation and seedling development in tomato, but it was unclear whether DEM1's primary role was in cell division or alternatively, in defining the identity of meristematic cells. Genome sequence analysis indicates that flowering plants possess at least two DEM genes. Arabidopsis has two DEM genes, DEM1 and DEM2, which we show are expressed in developing embryos and meristems in a punctate pattern that is typical of genes involved in cell division. Homozygous dem1 dem2 double mutants were not recovered, and plants carrying a single functional DEM1 allele and no functional copies of DEM2, i.e. DEM1/dem1 dem2/dem2 plants, exhibit normal development through to the time of flowering but during male reproductive development, chromosomes fail to align on the metaphase plate at meiosis II and result in abnormal numbers of daughter cells following meiosis. Additionally, these plants show defects in both pollen and embryo sac development, and produce defective male and female gametes. In contrast, dem1/dem1 DEM2/dem2 plants showed normal levels of fertility, indicating that DEM2 plays a more important role than DEM1 in gamete viability. The increased importance of DEM2 in gamete viability correlated with higher mRNA levels of DEM2 compared to DEM1 in most tissues examined and particularly in the vegetative shoot apex, developing siliques, pollen and sperm. We also demonstrate that gamete viability depends not only on the number of functional DEM alleles inherited following meiosis, but also on the number of functional DEM alleles in the parent plant that undergoes meiosis. Furthermore, DEM1 interacts with RAS-RELATED NUCLEAR PROTEIN 1 (RAN1) in yeast two-hybrid and pull-down binding assays, and we show that fluorescent proteins fused to DEM1 and RAN1 co-localize transiently during male meiosis and pollen development. In eukaryotes, RAN is a highly conserved GTPase that plays key roles in cell cycle progression, spindle assembly during cell division, reformation of the nuclear envelope following cell division, and nucleocytoplasmic transport. Our results demonstrate that DEM proteins play an essential role in cell division in plants, most likely through an interaction with RAN1.

Non-Ribosomal Peptides (NRPs) represent a biomedically important class of natural products that include a multitude of antibiotics and other clinically used drugs. NRPs are not directly encoded in the genome but are instead produced by metabolic pathways encoded by biosynthetic gene clusters (BGCs). Since the existing genome mining tools predict many putative NRPs synthesized by a given BGC, it remains unclear which of these putative NRPs are correct and how to identify post-assembly modifications of amino acids in these NRPs in a blind mode, without knowing which modifications exist in the sample. To address this challenge, here we report NRPminer, a modification-tolerant tool for NRP discovery from large (meta)genomic and mass spectrometry datasets. We show that NRPminer is able to identify many NRPs from different environments, including four previously unreported NRP families from soil-associated microbes and NRPs from human microbiota. Furthermore, in this work we demonstrate the anti-parasitic activities and the structure of two of these NRP families using direct bioactivity screening and nuclear magnetic resonance spectrometry, illustrating the power of NRPminer for discovering bioactive NRPs.

Microbial natural products are important for the understanding of microbial interactions, chemical defense and communication, and have also served as an inspirational source for numerous pharmaceutical drugs. Tropical marine cyanobacteria have been highlighted as a great source of new natural products, however, few reports have appeared wherein a multi-omics approach has been used to study their natural products potential (i.e., reports are often focused on an individual natural product and its biosynthesis). This study focuses on describing the natural product genetic potential as well as the expressed natural product molecules in benthic tropical cyanobacteria. We collected from several sites around the world and sequenced the genomes of 24 tropical filamentous marine cyanobacteria. The informatics program antiSMASH was used to annotate the major classes of gene clusters. BiG-SCAPE phylum-wide analysis revealed the most promising strains for natural product discovery among these cyanobacteria. LCMS/MS-based metabolomics highlighted the most abundant molecules and molecular classes among 10 of these marine cyanobacterial samples. We observed that despite many genes encoding for peptidic natural products, peptides were not as abundant as lipids and lipopeptides in the chemical extracts. Our results highlight a number of highly interesting biosynthetic gene clusters for genome mining among these cyanobacterial samples.

Although most organisms synthesize methionine from homocysteine and methyl folates, some have "core" methionine synthases that lack folate-binding domains and use other methyl donors. In vitro, the characterized core synthases use methylcobalamin as a methyl donor, but in vivo, they probably rely on corrinoid (vitamin B12-binding) proteins. We identified four families of core methionine synthases that are distantly related to each other (under 30% pairwise amino acid identity). From the characterized enzymes, we identified the families MesA, which is found in methanogens, and MesB, which is found in anaerobic bacteria and archaea with the Wood-Ljungdahl pathway. A third uncharacterized family, MesC, is found in anaerobic archaea that have the Wood-Ljungdahl pathway and lack known forms of methionine synthase. We predict that most members of the MesB and MesC families accept methyl groups from the iron-sulfur corrinoid protein of that pathway. The fourth family, MesD, is found only in aerobic bacteria. Using transposon mutants and complementation, we show that MesD does not require 5-methyltetrahydrofolate or cobalamin. Instead, MesD requires an uncharacterized protein family (DUF1852) and oxygen for activity.

G-protein-coupled receptors (GPCRs) play an integral role in the neurobiology of psychiatric disorders. Almost all neurotransmitters involved in psychiatric disorders act through GPCRs, and GPCRs are the most common targets of therapeutic drugs currently used in the treatment of psychiatric disorders. However, the roles of GPCRs in the etiology and pathophysiology of psychiatric disorders are not fully understood. Using publically available datasets, we performed a comprehensive analysis of the transcriptomic signatures of G-protein-linked signaling across the major psychiatric disorders: autism spectrum disorder (ASD), schizophrenia (SCZ), bipolar disorder (BP), and major depressive disorder (MDD). We also used the BrainSpan transcriptomic dataset of the developing human brain to examine whether GPCRs that exhibit chronological age-associated expressions have a higher tendency to be dysregulated in psychiatric disorders than age-independent GPCRs. We found that most GPCR genes were differentially expressed in the four disorders and that the GPCR superfamily as a gene cluster was overrepresented in the four disorders. We also identified a greater amplitude of gene expression changes in GPCRs than other gene families in the four psychiatric disorders. Further, dysregulated GPCRs overlapped across the four psychiatric disorders, with SCZ exhibiting the highest overlap with the three other disorders. Finally, the results revealed a greater tendency of age-associated GPCRs to be dysregulated in ASD than random GPCRs. Our results substantiate the central role of GPCR signaling pathways in the etiology and pathophysiology of psychiatric disorders. Furthermore, our study suggests that common GPCRs' signaling may mediate distinct phenotypic presentations across psychiatric disorders. Consequently, targeting these GPCRs could serve as a common therapeutic strategy to treat specific clinical symptoms across psychiatric disorders.

Microbial biosynthetic gene clusters (BGCs) encoding secondary metabolites are thought to impact a plethora of biologically mediated environmental processes, yet their discovery and functional characterization in natural microbiomes remains challenging. Here we describe deep long-read sequencing and assembly of metagenomes from biological soil crusts, a group of soil communities that are rich in BGCs. Taking advantage of the unusually long assemblies produced by this approach, we recovered nearly 3,000 BGCs for analysis, including 712 full-length BGCs. Functional exploration through metatranscriptome analysis of a 3-day wetting experiment uncovered phylum-specific BGC expression upon activation from dormancy, elucidating distinct roles and complex phylogenetic and temporal dynamics in wetting processes. For example, a pronounced increase in BGC transcription occurs at night primarily in cyanobacteria, implicating BGCs in nutrient scavenging roles and niche competition. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional dynamics of BGCs in environmental processes and suggests a central role of secondary metabolites in maintaining phylogenetically conserved niches within biocrusts.

Filamentous marine cyanobacteria make a variety of bioactive molecules that are produced by polyketide synthases, nonribosomal peptide synthetases, and hybrid pathways that are encoded by large biosynthetic gene clusters. These cyanobacterial natural products represent potential drug leads; however, thorough pharmacological investigations have been impeded by the limited quantity of compound that is typically available from the native organisms. Additionally, investigations of the biosynthetic gene clusters and enzymatic pathways have been difficult due to the inability to conduct genetic manipulations in the native producers. Here we report a set of genetic tools for the heterologous expression of biosynthetic gene clusters in the cyanobacteria Synechococcus elongatus PCC 7942 and Anabaena (Nostoc) PCC 7120. To facilitate the transfer of gene clusters in both strains, we engineered a strain of Anabaena that contains S. elongatus homologous sequences for chromosomal recombination at a neutral site and devised a CRISPR-based strategy to efficiently obtain segregated double recombinant clones of Anabaena. These genetic tools were used to express the large 28.7 kb cryptomaldamide biosynthetic gene cluster from the marine cyanobacterium Moorena (Moorea) producens JHB in both model strains. S. elongatus did not produce cryptomaldamide; however, high-titer production of cryptomaldamide was obtained in Anabaena. The methods developed in this study will facilitate the heterologous expression of biosynthetic gene clusters isolated from marine cyanobacteria and complex metagenomic samples.

BackgroundCupriavidus strain STM 6070 was isolated from nickel-rich soil collected near Koniambo massif, New Caledonia, using the invasive legume trap host Mimosa pudica. STM 6070 is a heavy metal-tolerant strain that is highly effective at fixing nitrogen with M. pudica. Here we have provided an updated taxonomy for STM 6070 and described salient features of the annotated genome, focusing on heavy metal resistance (HMR) loci and heavy metal efflux (HME) systems.ResultsThe 6,771,773 bp high-quality-draft genome consists of 107 scaffolds containing 6118 protein-coding genes. ANI values show that STM 6070 is a new species of Cupriavidus. The STM 6070 symbiotic region was syntenic with that of the M. pudica-nodulating Cupriavidus taiwanensis LMG 19424T. In contrast to the nickel and zinc sensitivity of C. taiwanensis strains, STM 6070 grew at high Ni2+ and Zn2+ concentrations. The STM 6070 genome contains 55 genes, located in 12 clusters, that encode HMR structural proteins belonging to the RND, MFS, CHR, ARC3, CDF and P-ATPase protein superfamilies. These HMR molecular determinants are putatively involved in arsenic (ars), chromium (chr), cobalt-zinc-cadmium (czc), copper (cop, cup), nickel (nie and nre), and silver and/or copper (sil) resistance. Seven of these HMR clusters were common to symbiotic and non-symbiotic Cupriavidus species, while four clusters were specific to STM 6070, with three of these being associated with insertion sequences. Within the specific STM 6070 HMR clusters, three novel HME-RND systems (nieIC cep nieBA, czcC2B2A2, and hmxB zneAC zneR hmxS) were identified, which constitute new candidate genes for nickel and zinc resistance.ConclusionsSTM 6070 belongs to a new Cupriavidus species, for which we have proposed the name Cupriavidus neocaledonicus sp. nov.. STM6070 harbours a pSym with a high degree of gene conservation to the pSyms of M. pudica-nodulating C. taiwanensis strains, probably as a result of recent horizontal transfer. The presence of specific HMR clusters, associated with transposase genes, suggests that the selection pressure of the New Caledonian ultramafic soils has driven the specific adaptation of STM 6070 to heavy-metal-rich soils via horizontal gene transfer.

BACKGROUND:Whiteflies are a threat to cassava (Manihot esculenta), an important staple food in many tropical/subtropical regions. Understanding the molecular mechanisms regulating cassava's responses against this pest is crucial for developing control strategies. Pathogenesis-related (PR) protein families are an integral part of plant immunity. With the availability of whole genome sequences, the annotation and expression programs of the full complement of PR genes in an organism can now be achieved. An understanding of the responses of the entire complement of PR genes during biotic stress and to the defense hormones, salicylic acid (SA) and jasmonic acid (JA), is lacking. Here, we analyze the responses of cassava PR genes to whiteflies, SA, JA, and other biotic aggressors. RESULTS:The cassava genome possesses 14 of the 17 plant PR families, with a total of 447 PR genes. A cassava PR gene nomenclature is proposed. Phylogenetic relatedness of cassava PR proteins to each other and to homologs in poplar, rice and Arabidopsis identified cassava-specific PR gene family expansions. The temporal programs of PR gene expression in response to the whitefly (Aleurotrachelus socialis) in four whitefly-susceptible cassava genotypes showed that 167 of the 447 PR genes were regulated after whitefly infestation. While the timing of PR gene expression varied, over 37% of whitefly-regulated PR genes were downregulated in all four genotypes. Notably, whitefly-responsive PR genes were largely coordinately regulated by SA and JA. The analysis of cassava PR gene expression in response to five other biotic stresses revealed a strong positive correlation between whitefly and Xanthomonas axonopodis and Cassava Brown Streak Virus responses and negative correlations between whitefly and Cassava Mosaic Virus responses. Finally, certain associations between PR genes in cassava expansions and response to biotic stresses were observed among PR families. CONCLUSIONS:This study represents the first genome-wide characterization of PR genes in cassava. PR gene responses to six biotic stresses and to SA and JA are demonstrably different to other angiosperms. We propose that our approach could be applied in other species to fully understand PR gene regulation by pathogens, pests and the canonical defense hormones SA and JA.

Streptomyces are efficient producers of various bioactive compounds, which are mostly synthesized by their secondary metabolite biosynthetic gene clusters (smBGCs). The smBGCs are tightly controlled by complex regulatory systems at transcriptional and translational levels to effectively utilize precursors that are supplied by primary metabolism. Thus, dynamic changes in gene expression in response to cellular status at both the transcriptional and translational levels should be elucidated to directly reflect protein levels, rapid downstream responses, and cellular energy costs. In this study, RNA-Seq and ribosome profiling were performed for five industrially important Streptomyces species at different growth phases, for the deep sequencing of total mRNA, and only those mRNA fragments that are protected by translating ribosomes, respectively. Herein, 12.0 to 763.8 million raw reads were sufficiently obtained with high quality of more than 80% for the Phred score Q30 and high reproducibility. These data provide a comprehensive understanding of the transcriptional and translational landscape across the Streptomyces species and contribute to facilitating the rational engineering of secondary metabolite production.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Gene families underlie genetic innovation and phenotypic diversification. However, our understanding of the early genomic and functional evolution of tandemly arranged gene families remains incomplete as paralog sequence similarity hinders their accurate characterization. The Drosophila melanogaster-specific gene family Sdic is tandemly repeated and impacts sperm competition. We scrutinized Sdic in 20 geographically diverse populations using reference-quality genome assemblies, read-depth methodologies, and qPCR, finding that ∼90% of the individuals harbor 3-7 copies as well as evidence of population differentiation. In strains with reliable gene annotations, copy number variation (CNV) and differential transposable element insertions distinguish one structurally distinct version of the Sdic region per strain. All 31 annotated copies featured protein-coding potential and, based on the protein variant encoded, were categorized into 13 paratypes differing in their 3' ends, with 3-5 paratypes coexisting in any strain examined. Despite widespread gene conversion, the only copy present in all strains has functionally diverged at both coding and regulatory levels under positive selection. Contrary to artificial tandem duplications of the Sdic region that resulted in increased male expression, CNV in cosmopolitan strains did not correlate with expression levels, likely as a result of differential genome modifier composition. Duplicating the region did not enhance sperm competitiveness, suggesting a fitness cost at high expression levels or a plateau effect. Beyond facilitating a minimally optimal expression level, Sdic CNV acts as a catalyst of protein and regulatory diversity, showcasing a possible evolutionary path recently formed tandem multigene families can follow toward long-term consolidation in eukaryotic genomes.

The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.IMPORTANCE Differences in the relative frequencies of sexual and asexual reproduction have profound implications for the accumulation of deleterious mutations (Muller's ratchet), but little is known about how these differences impact the evolution of ecologically important phenotypes. Aspergillus flavus is the main producer of aflatoxin, a notoriously potent carcinogen that often contaminates food. We investigated if differences in the levels of production of aflatoxin by A. flavus could be explained by the accumulation of deleterious mutations due to a lack of recombination. Despite differences in the extent of recombination, variation in aflatoxin production is better explained by the demography and history of specific populations and may suggest important differences in the ecological roles of aflatoxin among populations. Furthermore, the association of aflatoxin production and populations provides a means of predicting the risk of aflatoxin contamination by determining the frequencies of isolates from low- and high-production populations.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Terpenoid metabolism plays vital roles in stress defense and the environmental adaptation of monocot crops. Here, we describe the identification of the terpene synthase (TPS) gene family of the panicoid food and bioenergy model crop foxtail millet (Setaria italica). The diploid S. italica genome contains 32 TPS genes, 17 of which were biochemically characterized in this study. Unlike other thus far investigated grasses, S. italica contains TPSs producing all three ent-, (+)- and syn-copalyl pyrophosphate stereoisomers that naturally occur as central building blocks in the biosynthesis of distinct monocot diterpenoids. Conversion of these intermediates by the promiscuous TPS SiTPS8 yielded different diterpenoid scaffolds. Additionally, a cytochrome P450 monooxygenase (CYP99A17), which genomically clustered with SiTPS8, catalyzes the C19 hydroxylation of SiTPS8 products to generate the corresponding diterpene alcohols. The presence of syntenic orthologs to about 19% of the S. italica TPSs in related grasses supports a common ancestry of selected pathway branches. Among the identified enzyme products, abietadien-19-ol, syn-pimara-7,15-dien-19-ol and germacrene-d-4-ol were detectable in planta, and gene expression analysis of the biosynthetic TPSs showed distinct and, albeit moderately, inducible expression patterns in response to biotic and abiotic stress. In vitro growth-inhibiting activity of abietadien-19-ol and syn-pimara-7,15-dien-19-ol against Fusarium verticillioides and Fusarium subglutinans may indicate pathogen defensive functions, whereas the low antifungal efficacy of tested sesquiterpenoids supports other bioactivities. Together, these findings expand the known chemical space of monocot terpenoid metabolism to enable further investigations of terpenoid-mediated stress resilience in these agriculturally important species.

Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface "self-associated molecular patterns" that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2-3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a "hotspot" in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human-Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal-human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution.

Despite the remarkable microbial diversity found within humans, our ability to link genes to phenotypes is based upon a handful of model microorganisms. We report a comparative genomics platform for Eggerthella lenta and other Coriobacteriia, a neglected taxon broadly relevant to human health and disease. We uncover extensive genetic and metabolic diversity and validate a tool for mapping phenotypes to genes and sequence variants. We also present a tool for the quantification of strains from metagenomic sequencing data, enabling the identification of genes that predict bacterial fitness. Competitive growth is reproducible under laboratory conditions and attributable to intrinsic growth rates and resource utilization. Unique signatures of in vivo competition in gnotobiotic mice include an adhesin enriched in poor colonizers. Together, these computational and experimental resources represent a strong foundation for the continued mechanistic dissection of the Coriobacteriia and a template that can be applied to study other genetically intractable taxa.

Bacteria isolated from soils are major sources of specialized metabolites, including antibiotics and other compounds with clinical value that likely shape interactions among microbial community members and impact biogeochemical cycles. Yet, isolated lineages represent a small fraction of all soil bacterial diversity. It remains unclear how the production of specialized metabolites varies across the phylogenetic diversity of bacterial species in soils and whether the genetic potential for production of these metabolites differs with soil depth and vegetation type within a geographic region. We sampled soils and saprolite from three sites in a northern California Critical Zone Observatory with various vegetation and bedrock characteristics and reconstructed 1,334 metagenome-assembled genomes containing diverse biosynthetic gene clusters (BGCs) for secondary metabolite production. We obtained genomes for prolific producers of secondary metabolites, including novel groups within the Actinobacteria, Chloroflexi, and candidate phylum "Candidatus Dormibacteraeota." Surprisingly, one genome of a candidate phyla radiation (CPR) bacterium coded for a ribosomally synthesized linear azole/azoline-containing peptide, a capacity we found in other publicly available CPR bacterial genomes. Overall, bacteria with higher biosynthetic potential were enriched in shallow soils and grassland soils, with patterns of abundance of BGC type varying by taxonomy.IMPORTANCE Microbes produce specialized compounds to compete or communicate with one another and their environment. Some of these compounds, such as antibiotics, are also useful in medicine and biotechnology. Historically, most antibiotics have come from soil bacteria which can be isolated and grown in the lab. Though the vast majority of soil bacteria cannot be isolated, we can extract their genetic information and search it for genes which produce these specialized compounds. These understudied soil bacteria offer a wealth of potential for the discovery of new and important microbial products. Here, we identified the ability to produce these specialized compounds in diverse and novel bacteria in a range of soil environments. This information will be useful to other researchers who wish to isolate certain products. Beyond their use to humans, understanding the distribution and function of microbial products is key to understanding microbial communities and their effects on biogeochemical cycles.

α/β hydrolases make up a large and diverse protein superfamily. In natural product biosynthesis, cis-acting thioesterase α/β hydrolases can terminate biosynthetic assembly lines and release products by hydrolyzing or cyclizing the biosynthetic intermediate. Thioesterases can also act in trans, removing aberrant intermediates and restarting stalled biosynthesis. Knockout of this "editing" function leads to reduced product titers. The borrelidin biosynthetic gene cluster from Streptomyces parvulus Tü4055 contains a hitherto uncharacterized stand-alone thioesterase, borB. In this work, we demonstrate that purified BorB cleaves acyl substrates with a preference for propionate, which supports the hypothesis that it is also an editing thioesterase. The crystal structure of BorB shows a wedgelike hydrophobic substrate binding crevice that limits substrate length. To investigate the structure-function relationship, we made chimeric BorB variants using loop regions from characterized homologues with different specificities. BorB chimeras slightly reduced activity, arguing that the modified region is a not major determinant of substrate preference. The structure-function relationships described here contribute to the process of elimination for understanding thioesterase specificity and, ultimately, engineering and applying trans-acting thioesterases in biosynthetic assembly lines.