Reconstituting the genome of a young allopolyploid crop,Brassica napus,with its related species

Brassica napus (AⁿAⁿCⁿCⁿ) is an important worldwide oilseed crop, but it is a young allotetraploid with a short evolutionary history and limited genetic diversity. To significantly broaden its genetic diversity and create a novel heterotic population for sustainable rapeseed breeding, this study reconstituted the genome of B. napus by replacing it with the subgenomes from 122 accessions of Brassica rapa (A^rA^r) and 74 accessions of Brassica carinata (B^cB^cC^cC^c) and developing a novel gene pool of B. napus through five rounds of extensive recurrent selection. When compared with traditional B. napus using SSR markers and high‐throughput SNP/Indel markers through genotyping by sequencing, the newly developed gene pool and its homozygous progenies exhibited a large genetic distance, rich allelic diversity, new alleles and exotic allelic introgression across all 19 AC chromosomes. In addition to the abundant genomic variation detected in the AC genome, we also detected considerable introgression from the eight chromosomes of the B genome. Extensive trait variation and some genetic improvements were present from the early recurrent selection to later generations. This novel gene pool produced equally rich phenotypic variation and should be valuable for rapeseed genetic improvement. By reconstituting the genome of B. napus by introducing subgenomic variation within and between the related species using intense selection and recombination, the whole genome could be substantially reorganized. These results serve as an example of the manipulation of the genome of a young allopolyploid and provide insights into its rapid genome evolution affected by interspecific and intraspecific crosses.     

Mapping of homoeologous chromosome exchanges influencing quantitative trait variation in Brassica napus

Genomic rearrangements arising during polyploidization are an important source of genetic and phenotypic variation in the recent allopolyploid crop Brassica napus. Exchanges among homoeologous chromosomes, due to interhomoeologue pairing, and deletions without compensating homoeologous duplications are observed in both natural B. napus and synthetic B. napus. Rearrangements of large or small chromosome segments induce gene copy number variation (CNV) and can potentially cause phenotypic changes. Unfortunately, complex genome restructuring is difficult to deal with in linkage mapping studies. Here, we demonstrate how high‐density genetic mapping with codominant, physically anchored SNP markers can detect segmental homoeologous exchanges (HE) as well as deletions and accurately link these to QTL. We validated rearrangements detected in genetic mapping data by whole‐genome resequencing of parental lines along with cytogenetic analysis using fluorescence in situ hybridization with bacterial artificial chromosome probes (BAC‐FISH) coupled with PCR using primers specific to the rearranged region. Using a well‐known QTL region influencing seed quality traits as an example, we confirmed that HE underlies the trait variation in a DH population involving a synthetic B. napus trait donor, and succeeded in narrowing the QTL to a small defined interval that enables delineation of key candidate genes.     

SNP markers‐based map construction and genome‐wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag‐Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non‐SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous‐Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag‐Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high‐resolution mapping of loci in B. napus.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Bottlenecks in erucic acid accumulation in genetically engineered ultrahigh erucic acid Crambe abyssinica

Erucic acid is a valuable industrial fatty acid with many applications. The main producers of this acid are today high erucic rapeseed (Brassica napus) and mustard (Brassica juncea), which have 45%–50% of erucic acid in their seed oils. Crambe abyssinica is an alternative promising producer of this acid as it has 55%–60% of erucic acid in its oil. Through genetic modification (GM) of three genes, we have previously increased the level of erucic acid to 71% (68 mol%) in Crambe seed oil. In this study, we further investigated different aspects of oil biosynthesis in the developing GM Crambe seeds in comparison with wild‐type (Wt) Crambe, rapeseed and safflower (Carthamus tinctorius). We show that Crambe seeds have very low phosphatidylcholine‐diacylglycerol interconversion, suggesting it to be the main reason why erucic acid is limited in the membrane lipids during oil biosynthesis. We further show that GM Crambe seeds have slower seed development than Wt, accompanied by slower oil accumulation during the first 20 days after flowering (DAF). Despite low accumulation of erucic acid during early stages of GM seed development, nearly 86 mol% of all fatty acids accumulated between 27 and 50 DAF was erucic acid, when 40% of the total oil is laid down. Likely bottlenecks in the accumulation of erucic acid during early stages of GM Crambe seed development are discussed.     

Quantitative trait loci for growth and body size in the nine‐spined stickleback Pungitius pungitius L.

Body size is an ecologically important trait shown to be genetically variable both within and among different animal populations as revealed by quantitative genetic studies. However, few studies have looked into underlying genetic architecture of body size variability in the wild using genetic mapping methods. With the aid of quantitative trait loci (QTL) analyses based on 226 microsatellite markers, we mapped body size and growth rate traits in the nine‐spined stickleback (Pungitius pungitius) using an F₂‐intercross (n = 283 offspring) between size‐divergent populations. In total, 17 QTL locations were detected. The proportion of phenotypic variation explained by individual body size‐related QTL ranged from 3% to 12% and those related to growth parameters and increments from 3% to 10%. Several of the detected QTL affected either early or late growth. These results provide a solid starting point for more in depth investigations of structure and function of genomic regions involved in determination of body size in this popular model of ecological and evolutionary research.     

Development of iFOX‐hunting as a functional genomic tool and demonstration of its use to identify early senescence‐related genes in the polyploid Brassica napus

Functional genomic studies of many polyploid crops, including rapeseed (Brassica napus), are constrained by limited tool sets. Here we report development of a gain‐of‐function platform, termed ‘iFOX (inducible Full‐length cDNA OvereXpressor gene)‐Hunting’, for inducible expression of B. napus seed cDNAs in Arabidopsis. A Gateway‐compatible plant gene expression vector containing a methoxyfenozide‐inducible constitutive promoter for transgene expression was developed. This vector was used for cloning of random cDNAs from developing B. napus seeds and subsequent Agrobacterium‐mediated transformation of Arabidopsis. The inducible promoter of this vector enabled identification of genes upon induction that are otherwise lethal when constitutively overexpressed and to control developmental timing of transgene expression. Evaluation of a subset of the resulting ~6000 Arabidopsis transformants revealed a high percentage of lines with full‐length B. napus transgene insertions. Upon induction, numerous iFOX lines with visible phenotypes were identified, including one that displayed early leaf senescence. Phenotypic analysis of this line (rsl‐1327) after methoxyfenozide induction indicated high degree of leaf chlorosis. The integrated B. napuscDNA was identified as a homolog of an Arabidopsis acyl‐CoA binding protein (ACBP) gene designated BnACBP1‐like. The early senescence phenotype conferred by BnACBP1‐like was confirmed by constitutive expression of this gene in Arabidopsis and B. napus. Use of the inducible promoter in the iFOX line coupled with RNA‐Seq analyses allowed mechanistic clues and a working model for the phenotype associated with BnACBP1‐like expression. Our results demonstrate the utility of iFOX‐Hunting as a tool for gene discovery and functional characterization of Brassica napus genome.     

Induction of telomere‐mediated chromosomal truncation and behavior of truncated chromosomes in Brassica napus

Engineered minichromosomes could be stably inherited and serve as a platform for simultaneously transferring and stably expressing multiple genes. Chromosomal truncation mediated by repeats of telomeric sequences is a promising approach for the generation of minichromosomes. In the present work, direct repetitive sequences of Arabidopsis telomere were used to study telomere‐mediated truncation of chromosomes in Brassica napus. Transgenes containing alien Arabidopsis telomere were successfully obtained, and Southern blotting and fluorescence in situ hybridization (FISH) results show that the transgenes resulted in successful chromosomal truncation in B. napus. In addition, truncated chromosomes were inherited at rates lower than that predicted by Mendelian rules. To determine the potential manipulations and applications of the engineered chromosomes, such as the stacking of multiple transgenes and the Cre/lox and FRT/FLP recombination systems, both amenable to genetic manipulations through site‐specific recombination in somatic cells, were tested for their ability to undergo recombination in B. napus. These results demonstrate that alien Arabidopsis telomere is able to mediate chromosomal truncation in B. napus. This technology would be feasible for chromosomal engineering and for studies on chromosome structure and function in B. napus.     

The Brassica napus receptor‐like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus

Leucine‐rich repeat receptor‐like proteins (LRR‐RLPs) are highly adaptable parts of the signalling apparatus for extracellular detection of plant pathogens. Resistance to blackleg disease of Brassica spp. caused by Leptosphaeria maculans is largely governed by host race‐specific R‐genes, including the LRR‐RLP gene LepR3. The blackleg resistance gene Rlm2 was previously mapped to the same genetic interval as LepR3. In this study, the LepR3 locus of the Rlm2 Brassica napus line ‘Glacier DH24287’ was cloned, and B. napus transformants were analysed for recovery of the Rlm2 phenotype. Multiple B. napus, B. rapa and B. juncea lines were assessed for sequence variation at the locus. Rlm2 was found to be an allelic variant of the LepR3 LRR‐RLP locus, conveying race‐specific resistance to L. maculans isolates harbouring AvrLm2. Several defence‐related LRR‐RLPs have previously been shown to associate with the RLK SOBIR1 to facilitate defence signalling. Bimolecular fluorescence complementation (BiFC) and co‐immunoprecipitation of RLM2‐SOBIR1 studies revealed that RLM2 interacts with SOBIR1 of Arabidopsis thaliana when co‐expressed in Nicotiana benthamiana. The interaction of RLM2 with AtSOBIR1 is suggestive of a conserved defence signalling pathway between B. napus and its close relative A. thaliana.     

The impact of reducing fatty acid desaturation on the composition and thermal stability of rapeseed oil

Oilseed rape (Brassica napus) is the third largest source of vegetable oil globally. In addition to food uses, there are industrial applications that exploit the ability of the species to accumulate the very‐long‐chain fatty acid (VLCFA) erucic acid in its seed oil, controlled by orthologues of FATTY ACID ELONGASE 1 (Bna.FAE1.A8 and Bna.FAE1.C3). The proportion of polyunsaturated fatty acids (PUFAs) in rapeseed oil is predicted to affect its thermal stability and is controlled by orthologues of FATTY ACID DESATURASE 2, particularly Bna.FAD2.C5. Our aim was to develop rapeseed lines combining high erucic and low PUFA characters and to assess the impact on thermal stability of the oil they produce. The new type of rapeseed oil (high erucic low polyunsaturate; HELP) contained a substantially greater proportion of erucic acid (54%) compared with high erucic rapeseed oil (46%). Although the total VLCFA content was greater in oil from HELP lines (64%) than from high erucic rapeseed (57%), analysis of triacylglycerol composition showed negligible incorporation of VLCFAs into the sn‐2 position. Rancimat analysis showed that the thermal stability of rapeseed oil was improved greatly as a consequence of reduction of PUFA content, from 3.8 and 4.2 h in conventional low erucic and high erucic rapeseed oils, respectively, to 11.3 and 16.4 h in high oleic low PUFA (HOLP) and HELP oils, respectively. Our results demonstrate that engineering of the lipid biosynthetic pathway of rapeseed, using traditional approaches, enables the production of renewable industrial oils with novel composition and properties.     

Association of gene-linked SSR markers to seed glucosinolate content in oilseed rape (Brassica napus ssp. napus)

Breeding of oilseed rape (Brassica napus ssp. napus) has evoked a strong bottleneck selection towards double-low (00) seed quality with zero erucic acid and low seed glucosinolate content. The resulting reduction of genetic variability in elite 00-quality oilseed rape is particularly relevant with regard to the development of genetically diverse heterotic pools for hybrid breeding. In contrast, B. napus genotypes containing high levels of erucic acid and seed glucosinolates (++ quality) represent a comparatively genetically divergent source of germplasm. Seed glucosinolate content is a complex quantitative trait, however, meaning that the introgression of novel germplasm from this gene pool requires recurrent backcrossing to avoid linkage drag for high glucosinolate content. Molecular markers for key low-glucosinolate alleles could potentially improve the selection process. The aim of this study was to identify potentially gene-linked markers for important seed glucosinolate loci via structure-based allele-trait association studies in genetically diverse B. napus genotypes. The analyses included a set of new simple-sequence repeat (SSR) markers whose orthologs in Arabidopsis thaliana are physically closely linked to promising candidate genes for glucosinolate biosynthesis. We found evidence that four genes involved in the biosynthesis of indole, aliphatic and aromatic glucosinolates might be associated with known quantitative trait loci for total seed glucosinolate content in B. napus. Markers linked to homoeologous loci of these genes in the paleopolyploid B. napus genome were found to be associated with a significant effect on the seed glucosinolate content. This example shows the potential of Arabidopsis-Brassica comparative genome analysis for synteny-based identification of gene-linked SSR markers that can potentially be used in marker-assisted selection for an important trait in oilseed rape. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Distribution of mating‐type alleles and genetic variability in field populations of Leptosphaeria maculans in western Canada

Leptosphaeria maculans is the most important fungal pathogen of canola (Brassica napus, oilseed rape) that causes the devastating stem canker in canola fields of western Canada. The population genetic structure of L. maculans, represented by nine subpopulations from a 6‐year period and three different provinces in western Canada, was determined using ten minisatellite markers. Isolates collected at different locations in six consecutive years had an even distribution of MAT1‐1 and MAT1‐2 across the nine subpopulations. All subpopulations of L. maculans exhibited a moderate gene diversity (H = 0.356–0.585). The majority of the genetic variation occurred within subpopulations. Approximately 8% and 4% of the variations were distributed between sampling year and location, respectively. Genetic distance (F_ST) results, using analysis of molecular variation (AMOVA), indicated that subpopulation pairing within isolates by year ranged from F_ST = 0.010 to 0.109, and the location subpopulation ranged from F_ST = 0.038 to 0.085. Bayesian clustering analyses of multiloci inferred two distinct clusters in all the subpopulations examined. This study indicates a relatively high degree of gene exchange between the different L. maculans isolates. Our results suggest that this can occur in the wide growing areas of canola fields in western Canada. This gene exchange produced different gene allele frequencies and divergence between populations.     

Time‐specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage‐linkage disequilibrium (linkage‐LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full‐sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole‐genome re‐sequencing of the 435‐member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10^?7) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time‐specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time‐specific effects showed species‐specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.     

The identification and mapping of candidate genes and QTL involved in the fatty acid desaturation pathway in <Emphasis Type="Italic">Brassica napus</Emphasis>

We constructed a linkage map for the population QDH, which was derived from a cross between an oilseed rape cultivar and a resynthesised Brassica napus. The linkage map included ten markers linked to loci orthologous to those encoding fatty acid biosynthesis genes in Arabidopsis thaliana. The QDH population contains a high level of allelic variation, particularly in the C genome. We conducted quantitative trait locus (QTL) analyses, using field data obtained over 3 years, for the fatty acid composition of seed oil. The population segregates for the two major loci controlling erucic acid content, on linkage groups A8 and C3, which quantitatively affect the content of other fatty acids and is a problem generally encountered when crossing “wild” germplasm with cultivated “double low” oilseed rape cultivars. We assessed three methods for QTL analysis, interval mapping, multiple QTL mapping and single marker regression analysis of the subset of lines with low erucic acid. We found the third of these methods to be most appropriate for our main purpose, which was the study of the genetic control of the desaturation of 18-carbon fatty acids. This method enabled us to decouple the effect of the segregation of the erucic acid-controlling loci and identify 34 QTL for fatty acid content of seed oil, 14 in the A genome and 20 in the C genome. The QTL indicate the presence of 13 loci with novel alleles inherited from the progenitors of the resynthesised B. napus that might be useful for modulating the content or extent of desaturation of polyunsaturated fatty acids, only one of which coincides with the anticipated position of a candidate gene, an orthologue of FAD2.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Dissection of the genetic architecture of three seed‐quality traits and consequences for breeding in Brassica napus

Genome‐wide association studies (GWASs) combining high‐throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed‐quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well‐characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed‐quality traits, which would be useful for genetic improvement of B. napus.     

Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor‐bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R‐genes harbour more SNPs than variable genes. More nucleotide binding site‐leucine‐rich repeat (NBS‐LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics‐based improvement and improved disease resistance.     

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2‐year common‐garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non‐native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.     

Strong specificity in the interaction between parasitoids and symbiont‐protected hosts

Coevolution between hosts and parasites may promote the maintenance of genetic variation in both antagonists by negative frequency‐dependence if the host–parasite interaction is genotype‐specific. Here we tested for specificity in the interaction between parasitoids (Lysiphlebus fabarum) and aphid hosts (Aphis fabae) that are protected by a heritable defensive endosymbiont, the γ‐proteobacterium Hamiltonella defensa. Previous studies reported a lack of genotype specificity between unprotected aphids and parasitoids, but suggested that symbiont‐conferred resistance might exhibit a higher degree of specificity. Indeed, in addition to ample variation in host resistance as well as parasitoid infectivity, we found a strong aphid clone‐by‐parasitoid line interaction on the rates of successful parasitism. This genotype specificity appears to be mediated by H. defensa, highlighting the important role that endosymbionts can play in host–parasite coevolution.