Rapid Gene Loss Research Articles

A fungal plant pathogen overcomes mlo-mediated broad-spectrum disease resistance by rapid gene loss.

Hosts and pathogens typically engage in a coevolutionary arms race. This also applies to phytopathogenic powdery mildew fungi, which can rapidly overcome plant resistance and perform host jumps. Using experimental evolution, we show that the powdery mildew pathogen Blumeria hordei is capable of breaking the agriculturally important broad-spectrum resistance conditioned by barley loss-of-function mlo mutants. Partial mlo virulence of evolved B. hordei isolates is correlated with a distinctive pattern of adaptive mutations, including small-sized (c. 8-40 kb) deletions, of which one is linked to the de novo insertion of a transposable element. Occurrence of the mutations is associated with a transcriptional induction of effector protein-encoding genes that is absent in mlo-avirulent isolates on mlo mutant plants. The detected mutational spectrum comprises the same loci in at least two independently isolated mlo-virulent isolates, indicating convergent multigenic evolution. The mutational events emerged in part early (within the first five asexual generations) during experimental evolution, likely generating a founder population in which incipient mlo virulence was later stabilized by additional events. This work highlights the rapid dynamic genome evolution of an obligate biotrophic plant pathogen with a transposon-enriched genome.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution.

Since its initial discovery over 50years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

An angiosperm NLR Atlas reveals that NLR gene reduction is associated with ecological specialization and signal transduction component deletion

Nucleotide-binding leucine-rich-repeat (NLR) genes comprise the largest family of plant disease-resistance genes. Angiosperm NLR genes are phylogenetically divided into the TNL, CNL, and RNL subclasses. NLR copy numbers and subclass composition vary tremendously across angiosperm genomes. However, the evolutionary associations between genomic NLR content and ecological adaptation, or between NLR content and signal transduction components, are poorly characterized because of limited genome availability. In this study, we established an angiosperm NLR atlas (ANNA, https://biobigdata.nju.edu.cn/ANNA/) that includes NLR genes from over 300 angiosperm genomes. Using ANNA, we revealed that NLR copy numbers differ up to 66-fold among closely related species owing to rapid gene loss and gain. Interestingly, NLR contraction was associated with adaptations to aquatic, parasitic, and carnivorous lifestyles. The convergent NLR reduction in aquatic plants resembles the lack of NLR expansion during the long-term evolution of green algae before the colonization of land. A co-evolutionary pattern between NLR subclasses and plant immune pathway components was also identified, suggesting that immune pathway deficiencies may drive TNL loss. Finally, we identified a conserved TNL lineage that may function independently of the EDS1–SAG101–NRG1 module. Collectively, these findings provide new insights into the evolution of NLR genes in the context of ecological adaptation and genome content variation.

Evolution and comparative genomics of the most common Trichoderma species

BackgroundThe growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood.ResultsWe selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function.ConclusionsOur study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

Rapid genome reshaping by multiple-gene loss after whole-genome duplication in teleost fish suggested by mathematical modeling

Whole-genome duplication (WGD) is believed to be a significant source of major evolutionary innovation. Redundant genes resulting from WGD are thought to be lost or acquire new functions. However, the rates of gene loss and thus temporal process of genome reshaping after WGD remain unclear. The WGD shared by all teleost fish, one-half of all jawed vertebrates, was more recent than the two ancient WGDs that occurred before the origin of jawed vertebrates, and thus lends itself to analysis of gene loss and genome reshaping. Using a newly developed orthology identification pipeline, we inferred the post-teleost-specific WGD evolutionary histories of 6,892 protein-coding genes from nine phylogenetically representative teleost genomes on a time-calibrated tree. We found that rapid gene loss did occur in the first 60 My, with a loss of more than 70-80% of duplicated genes, and produced similar genomic gene arrangements within teleosts in that relatively short time. Mathematical modeling suggests that rapid gene loss occurred mainly by events involving simultaneous loss of multiple genes. We found that the subsequent 250 My were characterized by slow and steady loss of individual genes. Our pipeline also identified about 1,100 shared single-copy genes that are inferred to have become singletons before the divergence of clupeocephalan teleosts. Therefore, our comparative genome analysis suggests that rapid gene loss just after the WGD reshaped teleost genomes before the major divergence, and provides a useful set of marker genes for future phylogenetic analysis.

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini.

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

Hybrid origins and the earliest stages of diploidization in the highly successful recent polyploid Capsella bursa-pastoris

Whole-genome duplication (WGD) events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection, and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole-genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole-genome data provide strong support for recent hybrid origins of the tetraploid species within the past 100,000-300,000 y from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide due to the combined effects of demography, selfing, and genome redundancy from WGD. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that nonfunctionalization can receive a "head start" through a legacy of deleterious variants and differential expression originating in parental diploid populations.

AZFc Deletions and Spermatogenic Failure: A Population-Based Survey of 20,000 Y Chromosomes

Deletions involving the Y chromosome's AZFc region are the most common known genetic cause of severe spermatogenic failure (SSF). Six recurrent interstitial deletions affecting the region have been reported, but their population genetics are largely unexplored. We assessed the deletions' prevalence in 20,884 men in five populations and found four of the six deletions (presented here in descending order of prevalence): gr/gr, b2/b3, b1/b3, and b2/b4. One of every 27 men carried one of these four deletions. The 1.6 Mb gr/gr deletion, found in one of every 41 men, almost doubles the risk of SSF and accounts for ∼2% of SSF, although <2% of men with the deletion are affected. The 1.8 Mb b2/b3 deletion, found in one of every 90 men, does not appear to be a risk factor for SSF. The 1.6 Mb b1/b3 deletion, found in one of every 994 men, appears to increase the risk of SSF by a factor of 2.5, although <2% of men with the deletion are affected, and it accounts for only 0.15% of SSF. The 3.5 Mb b2/b4 deletion, found in one of every 2,320 men, increases the risk of SSF 145 times and accounts for ∼6% of SSF; the observed prevalence should approximate the rate at which the deletion arises anew in each generation. We conclude that a single rare variant of major effect (the b2/b4 deletion) and a single common variant of modest effect (the gr/gr deletion) are largely responsible for the AZFc region's contribution to SSF in the population.

Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes

The human X and Y chromosomes evolved from an ordinary pair of autosomes during the past 200–300 million years1–3. Due to genetic decay, the human MSY (male-specific region of Y chromosome) retains only three percent of the ancestral autosomes’ genes4,5. This evolutionary decay was driven by a series of five “stratification” events. Each event suppressed X-Y crossing over within a chromosome segment or “stratum”, incorporated that segment into the MSY, and subjected its genes to the erosive forces that attend the absence of crossing over2,6. The last of these events occurred 30 million years ago (mya), or 5 million years before the human and Old World monkey (OWM) lineages diverged. Although speculation abounds regarding ongoing decay and looming extinction of the human Y chromosome7–10, remarkably little is known about how many MSY genes were lost in the human lineage in the 25 million years that have followed its separation from the OWM lineage. To explore this question, we sequenced the MSY of the rhesus macaque, an OWM, and compared it to the human MSY. We discovered that, during the last 25 million years, MSY gene loss in the human lineage was limited to the youngest stratum (stratum 5), which comprises three percent of the human MSY. Within the older strata, which collectively comprise the bulk of the human MSY, gene loss evidently ceased more than 25 mya. Likewise, the rhesus MSY has not lost any older genes (from strata 1–4) during the past 25 million years, despite major structural differences from the human MSY. The rhesus MSY is simpler, with few amplified gene families or palindromes that might enable intrachromosomal recombination and repair. We present an empirical reconstruction of human MSY evolution in which each stratum transitioned from rapid, exponential loss of ancestral genes to strict conservation through purifying selection.

Gene expression in a paleopolyploid: a transcriptome resource for the ciliate Paramecium tetraurelia

BackgroundThe genome of Paramecium tetraurelia, a unicellular model that belongs to the ciliate phylum, has been shaped by at least 3 successive whole genome duplications (WGD). These dramatic events, which have also been documented in plants, animals and fungi, are resolved over evolutionary time by the loss of one duplicate for the majority of genes. Thanks to a low rate of large scale genome rearrangement in Paramecium, an unprecedented large number of gene duplicates of different ages have been identified, making this organism an outstanding model to investigate the evolutionary consequences of polyploidization. The most recent WGD, with 51% of pre-duplication genes still in 2 copies, provides a snapshot of a phase of rapid gene loss that is not accessible in more ancient polyploids such as yeast.ResultsWe designed a custom oligonucleotide microarray platform for P. tetraurelia genome-wide expression profiling and used the platform to measure gene expression during 1) the sexual cycle of autogamy, 2) growth of new cilia in response to deciliation and 3) biogenesis of secretory granules after massive exocytosis. Genes that are differentially expressed during these time course experiments have expression patterns consistent with a very low rate of subfunctionalization (partition of ancestral functions between duplicated genes) in particular since the most recent polyploidization event.ConclusionsA public transcriptome resource is now available for Paramecium tetraurelia. The resource has been integrated into the ParameciumDB model organism database, providing searchable access to the data. The microarray platform, freely available through NimbleGen Systems, provides a robust, cost-effective approach for genome-wide expression profiling in P. tetraurelia. The expression data support previous studies showing that at short evolutionary times after a whole genome duplication, gene dosage balance constraints and not functional change are the major determinants of gene retention.