A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici.

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Quantitative disease resistance to the bacterial pathogen Xanthomonas campestris involves an Arabidopsis immune receptor pair and a gene of unknown function.

Although quantitative disease resistance (QDR) is a durable and broad-spectrum form of resistance in plants, the identification of the genes underlying QDR is still in its infancy. RKS1 (Resistance related KinaSe1) has been reported recently to confer QDR in Arabidopsis thaliana to most but not all races of the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc). We therefore explored the genetic bases of QDR in A. thaliana to diverse races of X. campestris (Xc). A nested genome-wide association mapping approach was used to finely map the genomic regions associated with QDR to Xcc12824 (race 2) and XccCFBP6943 (race 6). To identify the gene(s) implicated in QDR, insertional mutants (T-DNA) were selected for the candidate genes and phenotyped in response to Xc. We identified two major QTLs that confer resistance specifically to Xcc12824 and XccCFBP6943. Although QDR to Xcc12824 is conferred by At5g22540 encoding for a protein of unknown function, QDR to XccCFBP6943 involves the well-known immune receptor pair RRS1/RPS4. In addition to RKS1, this study reveals that three genes are involved in resistance to Xc with strikingly different ranges of specificity, suggesting that QDR to Xc involves a complex network integrating multiple response pathways triggered by distinct pathogen molecular determinants.

Natural genetic variation of Xanthomonas campestris pv. campestris pathogenicity on arabidopsis revealed by association and reverse genetics.

ABSTRACT The pathogenic bacterium Xanthomonas campestris pv. campestris, the causal agent of black rot of Brassicaceae, manipulates the physiology and the innate immunity of its hosts. Association genetic and reverse-genetic analyses of a world panel of 45 X. campestris pv. campestris strains were used to gain understanding of the genetic basis of the bacterium's pathogenicity to Arabidopsis thaliana. We found that the compositions of the minimal predicted type III secretome varied extensively, with 18 to 28 proteins per strain. There were clear differences in aggressiveness of those X. campestris pv. campestris strains on two Arabidopsis natural accessions. We identified 3 effector genes (xopAC, xopJ5, and xopAL2) and 67 amplified fragment length polymorphism (AFLP) markers that were associated with variations in disease symptoms. The nature and distribution of the AFLP markers remain to be determined, but we observed a low linkage disequilibrium level between predicted effectors and other significant markers, suggesting that additional genetic factors make a meaningful contribution to pathogenicity. Mutagenesis of type III effectors in X. campestris pv. campestris confirmed that xopAC functions as both a virulence and an avirulence gene in Arabidopsis and that xopAM functions as a second avirulence gene on plants of the Col-0 ecotype. However, we did not detect the effect of any other effector in the X. campestris pv. campestris 8004 strain, likely due to other genetic background effects. These results highlight the complex genetic basis of pathogenicity at the pathovar level and encourage us to challenge the agronomical relevance of some virulence determinants identified solely in model strains. IMPORTANCE The identification and understanding of the genetic determinants of bacterial virulence are essential to be able to design efficient protection strategies for infected plants. The recent availability of genomic resources for a limited number of pathogen isolates and host genotypes has strongly biased our research toward genotype-specific approaches. Indeed, these do not consider the natural variation in both pathogens and hosts, so their applied relevance should be challenged. In our study, we exploited the genetic diversity of Xanthomonas campestris pv. campestris, the causal agent of black rot on Brassicaceae (e.g., cabbage), to mine for pathogenicity determinants. This work evidenced the contribution of known and unknown loci to pathogenicity relevant at the pathovar level and identified these virulence determinants as prime targets for breeding resistance to X. campestris pv. campestris in Brassicaceae.

Does segregating variation in sexual or microhabitat preferences lead to non-random mating within a population of Drosophila melanogaster?

Variation in female choice for mates has implications for the maintenance of genetic variation and the evolution of male traits. Yet, estimates of population-level variation in male mating success owing to female genotype are rare. Here, we used a panel of recombinant inbred lines to estimate the strength of selection at many genetic loci in a single generation and attempt to assess differences between females with respect to the males they mated with. We performed selection assays in a complex environment to allow differences in habitat or social group preference to be expressed. We detected directional selection at loci across the genome, but are unable to provide support for differential male success because of variation in female genotype.

Quantitative trait loci affecting phenotypic plasticity and the allometric relationship of ovariole number and thorax length in Drosophila melanogaster.

Environmental factors during juvenile growth such as temperature and nutrition have major effects on adult morphology and life-history traits. In Drosophila melanogaster, ovary size, measured as ovariole number, and body size, measured as thorax length, are developmentally plastic traits with respect to larval nutrition. Herein we investigated the genetic basis for plasticity of ovariole number and body size, as well the genetic basis for their allometric relationship using recombinant inbred lines (RILs) derived from a natural population in Winters, California. We reared 196 RILs in four yeast concentrations and measured ovariole number and body size. The genetic correlation between ovariole number and thorax length was positive, but the strength of this correlation decreased with increasing yeast concentration. Genetic variation and genotype-by-environment (G x E) interactions were observed for both traits. We identified quantitative trait loci (QTL), epistatic, QTL-by-environment, and epistatic-by-environment interactions for both traits and their scaling relationships. The results are discussed in the context of multivariate trait evolution.

Cis and trans regulatory effects contribute to natural variation in transcriptome of Drosophila melanogaster.

The dissection of intraspecific variation in transcriptome is a central theme of many recent quantitative genomic analyses. Transcript level variation has been attributed to factors at the gene itself (cis) and elsewhere in the genome (trans). Previous analyses of Drosophila intraspecific transcriptome variation pointed toward a larger contribution of trans factors. However, data from other genera, and from interspecific comparisons within Drosophila, are more consistent with a major role for cis factors. We investigated the relative amount of cis and trans variation in Drosophila melanogaster, using whole-genome expression from an oligonucleotide microarray in the 2 extensively studied genotypes Ore and 2b3, and 6 recombinant inbred (RI) lines derived from these parents. We examined 2 types of models to decompose cis and trans contributions to genetic variation in transcript level: 1) an infinitesimal model assuming that the transcription variation is highly polygenic and due to many small effects and 2) contrast models assuming that a few large effects contribute to the transcriptional variation. We explicitly fitted cis-by-trans interactions and extended our analyses to consider regulation of alternatively spliced transcripts. We estimated that approximately 10% of the transcriptome was differentially regulated among the lines. We were able to identify cis and trans effects that contribute to this differential regulation for 1,340 genes. Our analyses revealed numerous cis effects (90%) but much fewer trans effects, perhaps due to reduced power of detection for trans effects. In addition, we identified 15 genes that have alternative splice variants differentially regulated in cis.

How repeatable are associations between polymorphisms in achaete-scute and bristle number variation in Drosophila?

Currently, the relevance of common genetic variants--particularly those significantly associated with phenotypic variation in laboratory studies--to standing phenotypic variation in the wild is poorly understood. To address this, we quantified the relationship between achaete-scute complex (ASC) polymorphisms and Drosophila bristle number phenotypes in several new population samples. MC22 is a biallelic, nonrepetitive-length polymorphism 97 bp downstream of the scute transcript. It has been previously shown to be associated with sternopleural bristle number variation in both sexes in a set of isogenic lines. We replicated this association in a large cohort of wild-caught Drosophila melanogaster. We also detected a significant association at MC22 in an outbred population maintained under laboratory conditions for approximately 25 years, but the phenotypic effects in this sample were opposite from the direction estimated in the initial study. Finally, no significant associations were detected in a second large wild-caught cohort or in a set of 134 nearly isogenic lines. Our ability to repeat the initial association in wild samples suggests that it was not spurious. Nevertheless, inconsistent results from the other three panels suggest that the relationship between polymorphic genetic markers and loci contributing to continuous variation is not a simple one.

Sex-specific expression of alternative transcripts in Drosophila.

BACKGROUND: Many genes produce multiple transcripts due to alternative splicing or utilization of alternative transcription initiation/termination sites. This 'transcriptome expansion' is thought to increase phenotypic complexity by allowing a single locus to produce several functionally distinct proteins. However, sex, genetic and developmental variation in the representation of alternative transcripts has never been examined systematically. Here, we describe a genome-wide analysis of sex-specific expression of alternative transcripts in Drosophila melanogaster. RESULTS: We compared transcript profiles in males and females from eight Drosophila lines (OregonR and 2b, and 6 RIL) using a newly designed 60-mer oligonucleotide microarray that allows us to distinguish a large proportion of alternative transcripts. The new microarray incorporates 7,207 oligonucleotides, satisfying stringent binding and specificity criteria that target both the common and the unique regions of 2,768 multi-transcript genes, as well as 12,912 oligonucleotides that target genes with a single known transcript. We estimate that up to 22% of genes that produce multiple transcripts show a sex-specific bias in the representation of alternative transcripts. Sexual dimorphism in overall transcript abundance was evident for 53% of genes. The X chromosome contains a significantly higher proportion of genes with female-biased transcription than the autosomes. However, genes on the X chromosome are no more likely to have a sexual bias in alternative transcript representation than autosomal genes. CONCLUSION: Widespread sex-specific expression of alternative transcripts in Drosophila suggests that a new level of sexual dimorphism at the molecular level exists.

A low-cost open-source SNP genotyping platform for association mapping applications.

Association mapping aimed at identifying DNA polymorphisms that contribute to variation in complex traits entails genotyping a large number of single-nucleotide polymorphisms (SNPs) in a very large panel of individuals. Few technologies, however, provide inexpensive high-throughput genotyping. Here, we present an efficient approach developed specifically for genotyping large fixed panels of diploid individuals. The cost-effective, open-source nature of our methodology may make it particularly attractive to those working in nonmodel systems.

Genetics of resistance to transgenic Bacillus thuringiensis poplars in Chrysomela tremulae (Coleoptera: Chrysomelidae).

The area under genetically engineered plants producing Bacillus thuringiensis (Bt) toxins is steadily increasing. This increase has magnified the risk of alleles conferring resistance to these toxins being selected in natural populations of target insect pests. The speed at which this selection is likely to occur depends on the genetic characteristics of Bt resistance. We selected a strain of the beetle Chrysomela tremulae Fabricius on a transgenic Bt poplar clone Populus tremula L. x Populus tremuloides Michx producing high levels of B. thuringiensis Cry3Aa toxin. This strain was derived from an isofemale line that generated some F2 offspring that actively fed on this Bt poplar clone. The resistance ratio of the strain was >6400. Susceptibility had decreased to such an extent that the mortality of beetles of the strain fed Bt poplar leaves was similar to that of beetles fed nontransgenic poplar leaves. Genetic crosses between susceptible, resistant, and F1 hybrids showed that resistance to the Cry3Aa toxin was almost completely recessive (D(LC) = 0.07) and conferred by a single autosomal gene. The concentration of Cry3Aa produced in the transgenic Bt poplar used in this study was 6.34 times higher than the LC99 of the F1 hybrids, accounting for the complete recessivity (D(ML) = 0) of survival on Bt poplar leaves. Overall, the genetic characteristics of the resistance of C. tremulae to the Cry3Aa toxin are consistent with the assumptions underlying the high-dose refuge strategy, which aims to decrease the selection of Bt resistance alleles in natural target pest populations.

No evidence for an association between common nonsynonymous polymorphisms in delta and bristle number variation in natural and laboratory populations of Drosophila melanogaster.

We test the hypothesis that naturally occurring nonsynonymous variants in the Delta ligand of the Notch signaling pathway contribute to standing variation in sternopleural and/or abdominal bristle number in Drosophila melanogaster, for both a large cohort of wild-caught flies and previously described laboratory lines. We sequenced the transcribed region of Delta for 16 naturally occurring chromosomes and 65 SNPs, including 7 nonsynonymous SNPs (nsSNPs), were observed. Identified nsSNPs and 6 additional common SNPs, all located in exon 6 and the 3' UTR, were genotyped in 2060 wild-caught flies using an OLA-based methodology and genotyped in 38 additional natural chromosomes via DNA sequencing. None of the genotyped nsSNPs were significantly associated with natural variation in bristle number as assessed by a permutation test. A 95% upper bound on the additive genetic variance attributable to each genotyped SNP in the large natural cohort is <2% of the total phenotypic variation. Results suggest that two previously detected genotype/phenotype associations between bristle number and variants in the introns of Delta cannot be explained by linkage disequilibrium between these variants and nearby nonsynonymous variants. Unidentified regulatory variants more parsimoniously explain previous observations.

Initial frequency of alleles conferring resistance to Bacillus thuringiensis poplar in a field population of Chrysomela tremulae.

Globally, the estimated total area planted with transgenic plants producing Bacillus thuringiensis (Bt) toxins was 12 million hectares in 2001. The risk of target pests becoming resistant to these toxins has led to the implementation of resistance-management strategies. The efficiency and sustainability of these strategies, including the high-dose plus refuge strategy currently recommended for North American maize, depend on the initial frequency of resistance alleles. In this study, we estimated the initial frequencies of alleles conferring resistance to transgenic Bt poplars producing Cry3A in a natural population of the poplar pest Chrysomela tremulae (Coleoptera: Chrysomelidae). We used the F(2) screen method developed for detecting resistance alleles in natural pest populations. At least three parents of the 270 lines tested were heterozygous for a major Bt resistance allele. We estimated mean resistance-allele frequency for the period 1999-2001 at 0.0037 (95% confidence interval = 0.00045-0.0080) with a detection probability of 90%. These results demonstrate that (i) the F(2) screen method can be used to detect major alleles conferring resistance to Bt-producing plants in insects and (ii) the initial frequency of alleles conferring resistance to Bt toxin can be close to the highest theoretical values that are expected prior to the use of Bt plants if considering fitness costs and typical mutation rates.