Comparative Pangenomic Insights into the Distinct Evolution of Virulence Factors Among Grapevine Trunk Pathogens.

The permanent organs of grapevines (Vitis vinifera L.), like those of other woody perennials, are colonized by various unrelated pathogenic ascomycete fungi secreting cell wall-degrading enzymes and phytotoxic secondary metabolites that contribute to host damage and disease symptoms. Trunk pathogens differ in the symptoms they induce and the extent and speed of damage. Isolates of the same species often display a wide virulence range, even within the same vineyard. This study focuses on Eutypa lata, Neofusicoccum parvum, and Phaeoacremonium minimum, causal agents of Eutypa dieback, Botryosphaeria dieback, and Esca, respectively. We sequenced 50 isolates from viticulture regions worldwide and built nucleotide-level, reference-free pangenomes for each species. Through examination of genomic diversity and pangenome structure, we analyzed intraspecific conservation and variability of putative virulence factors, focusing on functions under positive selection and recent gene family dynamics of contraction and expansion. Our findings reveal contrasting distributions of putative virulence factors in the core, dispensable, and private genomes of each pangenome. For example, carbohydrate active enzymes (CAZymes) were prevalent in the core genomes of each pangenome, whereas biosynthetic gene clusters were prevalent in the dispensable genomes of E. lata and P. minimum. The dispensable fractions were also enriched in Gypsy transposable elements and virulence factors under positive selection (polyketide synthase genes in E. lata and P. minimum, glycosyltransferases in N. parvum). Our findings underscore the complexity of the genomic architecture in each species and provide insights into their adaptive strategies, enhancing our understanding of the underlying mechanisms of virulence. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative genomics of the Liberibacter genus reveals widespread diversity in genomic content and positive selection history.

'Candidatus Liberibacter' is a group of bacterial species that are obligate intracellular plant pathogens and cause Huanglongbing disease of citrus trees and Zebra Chip in potatoes. Here, we examined the extent of intra- and interspecific genetic diversity across the genus using comparative genomics. Our approach examined a wide set of Liberibacter genome sequences including five pathogenic species and one species not known to cause disease. By performing comparative genomics analyses, we sought to understand the evolutionary history of this genus and to identify genes or genome regions that may affect pathogenicity. With a set of 52 genomes, we performed comparative genomics, measured genome rearrangement, and completed statistical tests of positive selection. We explored markers of genetic diversity across the genus, such as average nucleotide identity across the whole genome. These analyses revealed the highest intraspecific diversity amongst the 'Ca. Liberibacter solanacearum' species, which also has the largest plant host range. We identified sets of core and accessory genes across the genus and within each species and measured the ratio of nonsynonymous to synonymous mutations (dN/dS) across genes. We identified ten genes with evidence of a history of positive selection in the Liberibacter genus, including genes in the Tad complex, which have been previously implicated as being highly divergent in the 'Ca. L. capsica' species based on high values of dN.

Sampling effect in predicting the evolutionary response of populations to climate change.

Genomic data and machine learning approaches have gained interest due to their potential to identify adaptive genetic variation across populations and to assess species vulnerability to climate change. By identifying gene-environment associations for putatively adaptive loci, these approaches project changes to adaptive genetic composition as a function of future climate change (genetic offsets), which are interpreted as measuring the future maladaptation of populations due to climate change. In principle, higher genetic offsets relate to increased population vulnerability and therefore can be used to set priorities for conservation and management. However, it is not clear how sensitive these metrics are to the intensity of population and individual sampling. Here, we use five genomic datasets with varying numbers of SNPs (NSNPs = 7006-1,398,773), sampled populations (Npop = 23-47) and individuals (Nind = 185-595) to evaluate the estimation sensitivity of genetic offsets to varying degrees of sampling intensity. We found that genetic offsets are sensitive to the number of populations being sampled, especially with less than 10 populations and when genetic structure is high. We also found that the number of individuals sampled per population had small effects on the estimation of genetic offsets, with more robust results when five or more individuals are sampled. Finally, uncertainty associated with the use of different future climate scenarios slightly increased estimation uncertainty in the genetic offsets. Our results suggest that sampling efforts should focus on increasing the number of populations, rather than the number of individuals per populations, and that multiple future climate scenarios should be evaluated to ascertain estimation sensitivity.

Multigenic resistance to Xylella fastidiosa in wild grapes (Vitis sps.) and its implications within a changing climate.

Xylella fastidiosa is a bacterium that infects crops like grapevines, coffee, almonds, citrus and olives. There is little understanding of the genes that contribute to plant resistance, the genomic architecture of resistance, and the potential role of climate in shaping resistance, in part because major crops like grapevines (Vitis vinifera) are not resistant to the bacterium. Here we study a wild grapevine species, V. arizonica, that segregates for resistance. Using genome-wide association, we identify candidate resistance genes. Resistance-associated kmers are shared with a sister species of V. arizonica but not with more distant species, suggesting that resistance evolved more than once. Finally, resistance is climate dependent, because individuals from low ( < 10 °C) temperature locations in the wettest quarter were typically susceptible to infection, likely reflecting a lack of pathogen pressure in colder climates. In fact, climate is as effective a predictor of resistance phenotypes as some genetic markers. We extend our climate observations to additional crops, predicting that increased pathogen pressure is more likely for grapevines and almonds than some other susceptible crops.

Insights into the domestication of avocado and potential genetic contributors to heterodichogamy.

The domestication history of the avocado (Persea americana) remains unclear. We created a reference genome from the Gwen varietal, which is closely related to the economically dominant Hass varietal. Our genome assembly had an N50 of 3.37 megabases, a BUSCO score of 91%, and was scaffolded with a genetic map, producing 12 pseudo-chromosomes with 49,450 genes. We used the Gwen genome as a reference to investigate population genomics, based on a sample of 34 resequenced accessions that represented the 3 botanical groups of P. americana. Our analyses were consistent with 3 separate domestication events; we estimated that the Mexican group diverged from the Lowland (formerly known as "West Indian") and Guatemalan groups >1 million years ago. We also identified putative targets of selective sweeps in domestication events; within the Guatemalan group, putative candidate genes were enriched for fruit development and ripening. We also investigated divergence between heterodichogamous flowering types, providing preliminary evidence for potential candidate genes involved in pollination and floral development.

Evaluating the persistence and utility of five wild Vitis species in the context of climate change.

Crop wild relatives (CWRs) have the capacity to contribute novel traits to agriculture. Given climate change, these contributions may be especially vital for the persistence of perennial crops, because perennials are often clonally propagated and consequently do not evolve rapidly. By studying the landscape genomics of samples from five Vitis CWRs (V. arizonica, V. mustangensis, V. riparia, V. berlandieri and V. girdiana) in the context of projected climate change, we addressed two goals. The first was to assess the relative potential of different CWR accessions to persist in the face of climate change. By integrating species distribution models with adaptive genetic variation, additional genetic features such as genomic load and a phenotype (resistance to Pierce's Disease), we predicted that accessions from one species (V. mustangensis) are particularly well-suited to persist in future climates. The second goal was to identify which CWR accessions may contribute to bioclimatic adaptation for grapevine (V. vinifera) cultivation. To do so, we evaluated whether CWR accessions have the allelic capacity to persist if moved to locations where grapevines are cultivated in the United States. We identified six candidates from V. mustangensis and hypothesized that they may prove useful for contributing alleles that can mitigate climate impacts on viticulture. By identifying candidate germplasm, this study takes a conceptual step toward assessing the genomic and bioclimatic characteristics of CWRs.

Using Genomes and Evolutionary Analyses to Screen for Host-Specificity and Positive Selection in the Plant Pathogen Xylella fastidiosa.

Xylella fastidiosa infects several economically important crops in the Americas, and it also recently emerged in Europe. Here, using a set of Xylella genomes reflective of the genus-wide diversity, we performed a pan-genome analysis based on both core and accessory genes for two purposes: (i) to test associations between genetic divergence and plant host species and (ii) to identify positively selected genes that are potentially involved in arms-race dynamics. For the former, tests yielded significant evidence for the specialization of X. fastidiosa to plant host species. This observation contributes to a growing literature suggesting that the phylogenetic history of X. fastidiosa lineages affects the host range. For the latter, our analyses uncovered evidence of positive selection across codons for 5.3% (67 of 1,257) of the core genes and 5.4% (201 of 3,691) of the accessory genes. These genes are candidates to encode interacting factors with plant and insect hosts. Most of these genes had unknown functions, but we did identify some tractable candidates, including nagZ_2, which encodes a beta-glucosidase that is important for Neisseria gonorrhoeae biofilm formation; cya, which modulates gene expression in pathogenic bacteria, and barA, a membrane associated histidine kinase that has roles in cell division, metabolism, and pili formation. IMPORTANCE Xylella fastidiosa causes devasting diseases to several critical crops. Because X. fastidiosa colonizes and infects many plant species, it is important to understand whether the genome of X. fastidiosa has genetic determinants that underlie specialization to specific host plants. We analyzed genome sequences of X. fastidiosa to investigate evolutionary relationships and to test for evidence of positive selection on specific genes. We found a significant signal between genome diversity and host plants, consistent with bacterial specialization to specific plant hosts. By screening for positive selection, we identified both core and accessory genes that may affect pathogenicity, including genes involved in biofilm formation.

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation.

BACKGROUND: Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks. RESULTS: By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors. CONCLUSIONS: Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Phylogenomics of Plant-Associated Botryosphaeriaceae Species.

The Botryosphaeriaceae is a fungal family that includes many destructive vascular pathogens of woody plants (e.g., Botryosphaeria dieback of grape, Panicle blight of pistachio). Species in the genera Botryosphaeria, Diplodia, Dothiorella, Lasiodiplodia, Neofusicoccum, and Neoscytalidium attack a range of horticultural crops, but they vary in virulence and their abilities to infect their hosts via different infection courts (flowers, green shoots, woody twigs). Isolates of seventeen species, originating from symptomatic apricot, grape, pistachio, and walnut were tested for pathogenicity on grapevine wood after 4 months of incubation in potted plants in the greenhouse. Results revealed significant variation in virulence in terms of the length of the internal wood lesions caused by these seventeen species. Phylogenomic comparisons of the seventeen species of wood-colonizing fungi revealed clade-specific expansion of gene families representing putative virulence factors involved in toxin production and mobilization, wood degradation, and nutrient uptake. Statistical analyses of the evolution of the size of gene families revealed expansions of secondary metabolism and transporter gene families in Lasiodiplodia and of secreted cell wall degrading enzymes (CAZymes) in Botryosphaeria and Neofusicoccum genomes. In contrast, Diplodia, Dothiorella, and Neoscytalidium generally showed a contraction in the number of members of these gene families. Overall, species with expansions of gene families, such as secreted CAZymes, secondary metabolism, and transporters, were the most virulent (i.e., were associated with the largest lesions), based on our pathogenicity tests and published reports. This study represents the first comparative phylogenomic investigation into the evolution of possible virulence factors from diverse, cosmopolitan members of the Botryosphaeriaceae.

Independent Whole-Genome Duplications Define the Architecture of the Genomes of the Devastating West African Cacao Black Pod Pathogen Phytophthora megakarya and Its Close Relative Phytophthora palmivora.

Phytophthora megakarya and P. palmivora are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While P. palmivora is a cosmopolitan pathogen, P. megakarya, which is more aggressive on cacao than P. palmivora, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule real-time sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. The P. megakarya genome is exceptionally large (222 Mbp) and nearly twice the size of P. palmivora (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). Previous reports pointed toward a whole-genome duplication (WGD) in P. palmivora In this study, we demonstrate that both species underwent independent and relatively recent WGD events. In P. megakarya we identified a unique combination of WGD and large-scale transposable element driven genome expansion, which places this genome in the upper range of Phytophthora genome sizes, as well as effector pools with 1,382 predicted RxLR effectors. Finally, this study provides evidence of adaptive evolution of effectors like RxLRs and Crinklers, and discusses the implications of effector expansion and diversification.

Profiling grapevine trunk pathogens in planta: a case for community-targeted DNA metabarcoding.

BACKGROUND: DNA metabarcoding, commonly used in exploratory microbial ecology studies, is a promising method for the simultaneous in planta-detection of multiple pathogens associated with disease complexes, such as the grapevine trunk diseases. Profiling of pathogen communities associated with grapevine trunk diseases is particularly challenging, due to the presence within an individual wood lesion of multiple co-infecting trunk pathogens and other wood-colonizing fungi, which span a broad range of taxa in the fungal kingdom. As such, we designed metabarcoding primers, using as template the ribosomal internal transcribed spacer of grapevine trunk-associated ascomycete fungi (GTAA) and compared them to two universal primer widely used in microbial ecology. RESULTS: We first performed in silico simulations and then tested the primers by high-throughput amplicon sequencing of (i) multiple combinations of mock communities, (ii) time-course experiments with controlled inoculations, and (iii) diseased field samples from vineyards under natural levels of infection. All analyses showed that GTAA had greater affinity and sensitivity, compared to those of the universal primers. Importantly, with GTAA, profiling of mock communities and comparisons with shotgun-sequencing metagenomics of field samples gave an accurate representation of genera of important trunk pathogens, namely Phaeomoniella, Phaeoacremonium, and Eutypa, the abundances of which were over- or under-estimated with universal primers. CONCLUSIONS: Overall, our findings not only demonstrate that DNA metabarcoding gives qualitatively and quantitatively accurate results when applied to grapevine trunk diseases, but also that primer customization and testing are crucial to ensure the validity of DNA metabarcoding results.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum.

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Lipopolysaccharide O-antigen delays plant innate immune recognition of Xylella fastidiosa.

Lipopolysaccharides (LPS) are among the known pathogen-associated molecular patterns (PAMPs). LPSs are potent elicitors of PAMP-triggered immunity (PTI), and bacteria have evolved intricate mechanisms to dampen PTI. Here we demonstrate that Xylella fastidiosa (Xf), a hemibiotrophic plant pathogenic bacterium, possesses a long chain O-antigen that enables it to delay initial plant recognition, thereby allowing it to effectively skirt initial elicitation of innate immunity and establish itself in the host. Lack of the O-antigen modifies plant perception of Xf and enables elicitation of hallmarks of PTI, such as ROS production specifically in the plant xylem tissue compartment, a tissue not traditionally considered a spatial location of PTI. To explore translational applications of our findings, we demonstrate that pre-treatment of plants with Xf LPS primes grapevine defenses to confer tolerance to Xf challenge.

Closed-reference metatranscriptomics enables in planta profiling of putative virulence activities in the grapevine trunk disease complex.

Grapevines, like other perennial crops, are affected by so-called 'trunk diseases', which damage the trunk and other woody tissues. Mature grapevines typically contract more than one trunk disease and often multiple grapevine trunk pathogens (GTPs) are recovered from infected tissues. The co-existence of different GTP species in complex and dynamic microbial communities complicates the study of the molecular mechanisms underlying disease development, especially under vineyard conditions. The objective of this study was to develop and optimize a community-level transcriptomics (i.e. metatranscriptomics) approach that could monitor simultaneously the virulence activities of multiple GTPs in planta. The availability of annotated genomes for the most relevant co-infecting GTPs in diseased grapevine wood provided the unprecedented opportunity to generate a multi-species reference for the mapping and quantification of DNA and RNA sequencing reads. We first evaluated popular sequence read mappers using permutations of multiple simulated datasets. Alignment parameters of the selected mapper were optimized to increase the specificity and sensitivity for its application to metagenomics and metatranscriptomics analyses. Initial testing on grapevine wood experimentally inoculated with individual GTPs confirmed the validity of the method. Using naturally infected field samples expressing a variety of trunk disease symptoms, we show that our approach provides quantitative assessments of species composition, as well as genome-wide transcriptional profiling of potential virulence factors, namely cell wall degradation, secondary metabolism and nutrient uptake for all co-infecting GTPs.

Condition-dependent co-regulation of genomic clusters of virulence factors in the grapevine trunk pathogen Neofusicoccum parvum.

The ascomycete Neofusicoccum parvum, one of the causal agents of Botryosphaeria dieback, is a destructive wood-infecting fungus and a serious threat to grape production worldwide. The capability to colonize woody tissue, combined with the secretion of phytotoxic compounds, is thought to underlie its pathogenicity and virulence. Here, we describe the repertoire of virulence factors and their transcriptional dynamics as the fungus feeds on different substrates and colonizes the woody stem. We assembled and annotated a highly contiguous genome using single-molecule real-time DNA sequencing. Transcriptome profiling by RNA sequencing determined the genome-wide patterns of expression of virulence factors both in vitro (potato dextrose agar or medium amended with grape wood as substrate) and in planta. Pairwise statistical testing of differential expression, followed by co-expression network analysis, revealed that physically clustered genes coding for putative virulence functions were induced depending on the substrate or stage of plant infection. Co-expressed gene clusters were significantly enriched not only in genes associated with secondary metabolism, but also in those associated with cell wall degradation, suggesting that dynamic co-regulation of transcriptional networks contributes to multiple aspects of N. parvum virulence. In most of the co-expressed clusters, all genes shared at least a common motif in their promoter region, indicative of co-regulation by the same transcription factor. Co-expression analysis also identified chromatin regulators with correlated expression with inducible clusters of virulence factors, suggesting a complex, multi-layered regulation of the virulence repertoire of N. parvum.

Phased diploid genome assembly with single-molecule real-time sequencing.

While genome assembly projects have been successful in many haploid and inbred species, the assembly of noninbred or rearranged heterozygous genomes remains a major challenge. To address this challenge, we introduce the open-source FALCON and FALCON-Unzip algorithms (https://github.com/PacificBiosciences/FALCON/) to assemble long-read sequencing data into highly accurate, contiguous, and correctly phased diploid genomes. We generate new reference sequences for heterozygous samples including an F1 hybrid of Arabidopsis thaliana, the widely cultivated Vitis vinifera cv. Cabernet Sauvignon, and the coral fungus Clavicorona pyxidata, samples that have challenged short-read assembly approaches. The FALCON-based assemblies are substantially more contiguous and complete than alternate short- or long-read approaches. The phased diploid assembly enabled the study of haplotype structure and heterozygosities between homologous chromosomes, including the identification of widespread heterozygous structural variation within coding sequences.

Developmental and Metabolic Plasticity of White-Skinned Grape Berries in Response to Botrytis cinerea during Noble Rot.

Noble rot results from exceptional infections of ripe grape (Vitis vinifera) berries by Botrytis cinerea. Unlike bunch rot, noble rot promotes favorable changes in grape berries and the accumulation of secondary metabolites that enhance wine grape composition. Noble rot-infected berries of cv Sémillon, a white-skinned variety, were collected over 3 years from a commercial vineyard at the same time that fruit were harvested for botrytized wine production. Using an integrated transcriptomics and metabolomics approach, we demonstrate that noble rot alters the metabolism of cv Sémillon berries by inducing biotic and abiotic stress responses as well as ripening processes. During noble rot, B. cinerea induced the expression of key regulators of ripening-associated pathways, some of which are distinctive to the normal ripening of red-skinned cultivars. Enhancement of phenylpropanoid metabolism, characterized by a restricted flux in white-skinned berries, was a common outcome of noble rot and red-skinned berry ripening. Transcript and metabolite analyses together with enzymatic assays determined that the biosynthesis of anthocyanins is a consistent hallmark of noble rot in cv Sémillon berries. The biosynthesis of terpenes and fatty acid aroma precursors also increased during noble rot. We finally characterized the impact of noble rot in botrytized wines. Altogether, the results of this work demonstrated that noble rot causes a major reprogramming of berry development and metabolism. This desirable interaction between a fruit and a fungus stimulates pathways otherwise inactive in white-skinned berries, leading to a greater accumulation of compounds involved in the unique flavor and aroma of botrytized wines.

Distinctive expansion of gene families associated with plant cell wall degradation, secondary metabolism, and nutrient uptake in the genomes of grapevine trunk pathogens.

BACKGROUND: Trunk diseases threaten the longevity and productivity of grapevines in all viticulture production systems. They are caused by distantly-related fungi that form chronic wood infections. Variation in wood-decay abilities and production of phytotoxic compounds are thought to contribute to their unique disease symptoms. We recently released the draft sequences of Eutypa lata, Neofusicoccum parvum and Togninia minima, causal agents of Eutypa dieback, Botryosphaeria dieback and Esca, respectively. In this work, we first expanded genomic resources to three important trunk pathogens, Diaporthe ampelina, Diplodia seriata, and Phaeomoniella chlamydospora, causal agents of Phomopsis dieback, Botryosphaeria dieback, and Esca, respectively. Then we integrated all currently-available information into a genome-wide comparative study to identify gene families potentially associated with host colonization and disease development. RESULTS: The integration of RNA-seq, comparative and ab initio approaches improved the protein-coding gene prediction in T. minima, whereas shotgun sequencing yielded nearly complete genome drafts of Dia. ampelina, Dip. seriata, and P. chlamydospora. The predicted proteomes of all sequenced trunk pathogens were annotated with a focus on functions likely associated with pathogenesis and virulence, namely (i) wood degradation, (ii) nutrient uptake, and (iii) toxin production. Specific patterns of gene family expansion were described using Computational Analysis of gene Family Evolution, which revealed lineage-specific evolution of distinct mechanisms of virulence, such as specific cell wall oxidative functions and secondary metabolic pathways in N. parvum, Dia. ampelina, and E. lata. Phylogenetically-informed principal component analysis revealed more similar repertoires of expanded functions among species that cause similar symptoms, which in some cases did not reflect phylogenetic relationships, thereby suggesting patterns of convergent evolution. CONCLUSIONS: This study describes the repertoires of putative virulence functions in the genomes of ubiquitous grapevine trunk pathogens. Gene families with significantly faster rates of gene gain can now provide a basis for further studies of in planta gene expression, diversity by genome re-sequencing, and targeted reverse genetic approaches. The functional validation of potential virulence factors will lead to a more comprehensive understanding of the mechanisms of pathogenesis and virulence, which ultimately will enable the development of accurate diagnostic tools and effective disease management.

Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator.

BACKGROUND: Powdery mildew, caused by the obligate biotrophic fungus Erysiphe necator, is an economically important disease of grapevines worldwide. Large quantities of fungicides are used for its control, accelerating the incidence of fungicide-resistance. Copy number variations (CNVs) are unbalanced changes in the structure of the genome that have been associated with complex traits. In addition to providing the first description of the large and highly repetitive genome of E. necator, this study describes the impact of genomic structural variation on fungicide resistance in Erysiphe necator. RESULTS: A shotgun approach was applied to sequence and assemble the genome of five E. necator isolates, and RNA-seq and comparative genomics were used to predict and annotate protein-coding genes. Our results show that the E. necator genome is exceptionally large and repetitive and suggest that transposable elements are responsible for genome expansion. Frequent structural variations were found between isolates and included copy number variation in EnCYP51, the target of the commonly used sterol demethylase inhibitor (DMI) fungicides. A panel of 89 additional E. necator isolates collected from diverse vineyard sites was screened for copy number variation in the EnCYP51 gene and for presence/absence of a point mutation (Y136F) known to result in higher fungicide tolerance. We show that an increase in EnCYP51 copy number is significantly more likely to be detected in isolates collected from fungicide-treated vineyards. Increased EnCYP51 copy numbers were detected with the Y136F allele, suggesting that an increase in copy number becomes advantageous only after the fungicide-tolerant allele is acquired. We also show that EnCYP51 copy number influences expression in a gene-dose dependent manner and correlates with fungal growth in the presence of a DMI fungicide. CONCLUSIONS: Taken together our results show that CNV can be adaptive in the development of resistance to fungicides by providing increasing quantitative protection in a gene-dosage dependent manner. The results of this work not only demonstrate the effectiveness of using genomics to dissect complex traits in organisms with very limited molecular information, but also may have broader implications for understanding genomic dynamics in response to strong selective pressure in other pathogens with similar genome architectures.

Genome-wide transcriptional profiling of Botrytis cinerea genes targeting plant cell walls during infections of different hosts.

Cell walls are barriers that impair colonization of host tissues, but also are important reservoirs of energy-rich sugars. Growing hyphae of necrotrophic fungal pathogens, such as Botrytis cinerea (Botrytis, henceforth), secrete enzymes that disassemble cell wall polysaccharides. In this work we describe the annotation of 275 putative secreted Carbohydrate-Active enZymes (CAZymes) identified in the Botrytis B05.10 genome. Using RNAseq we determined which Botrytis CAZymes were expressed during infections of lettuce leaves, ripe tomato fruit, and grape berries. On the three hosts, Botrytis expressed a common group of 229 potentially secreted CAZymes, including 28 pectin backbone-modifying enzymes, 21 hemicellulose-modifying proteins, 18 enzymes that might target pectin and hemicellulose side-branches, and 16 enzymes predicted to degrade cellulose. The diversity of the Botrytis CAZymes may be partly responsible for its wide host range. Thirty-six candidate CAZymes with secretion signals were found exclusively when Botrytis interacted with ripe tomato fruit and grape berries. Pectin polysaccharides are notably abundant in grape and tomato cell walls, but lettuce leaf walls have less pectin and are richer in hemicelluloses and cellulose. The results of this study not only suggest that Botrytis targets similar wall polysaccharide networks on fruit and leaves, but also that it may selectively attack host wall polysaccharide substrates depending on the host tissue.