What do spring migrants reveal about sex and host selection in the melon aphid?

Chromosome-level genome assemblies of two cotton-melon aphid Aphis gossypii biotypes unveil mechanisms of host adaption

The cotton-melon aphid Aphis gossypii is a sap-sucking insect that is considered a serious global pest. The species is distributed over a large geographical range and uses a wide variety of hosts, with some populations being specialized to attack different plant species. Here, we provide de novo chromosome-level genome assemblies of a cotton specialist population (Hap1) and a cucurbit specialist population (Hap3). We achieved this by using a combination of third-generation sequencing platforms, namely Illumina and Hi-C sequencing technologies. We were able to anchor a total of 334.89 Mb (scaffold N50 of 89.13 Mb) and 359.95 Mb (scaffold N50 of 68.88 Mb) to four chromosomes for Hap1 and Hap3, respectively. Moreover, our results showed that the X-chromosome of Hap3 (113.01 Mb) was significantly longer than that of Hap1 (100.26 Mb), with a high level of sequence conservation between the aphid species. We also report variation in the number of protein-coding genes and repeat sequences between Hap1 and Hap3. In particular, olfactory and gustatory receptor genes underwent a high level of gene duplication and expansion events in A. gossypii, including between Hap1 and Hap3. Moreover, we identified two glutathione S-transferase genes which underwent single gene duplications in Hap3, and tandem duplication and inversion events affecting the cytochrome P450 monooxygenase between Hap1 and Hap3, all of which include the CYP3 family. Our results illustrate the variance in the genomic composition of two specialized A. gossypii populations and provide a helpful resource for the study of aphid population evolution, host adaption and insecticide resistance.

Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

Aphis gossypii/Aphis frangulae collected worldwide: Microsatellite markers data and genetic cluster assignment

Aphis gossypii is a cosmopolitan aphid species able to colonize hundreds of plant species from various families [1]. It causes serious damage to a wide range of crops and it is considered a major pest of cucurbits and cotton [2]. It reproduces clonally, by obligate parthenogenesis, on secondary hosts present throughout the year in the intertropical area. At higher latitude, some lineages clonally overwinter but part of the population may have a sexual reproduction in autumn on primary host such as Hibiscus syriacus, to generate cold resistant overwintering eggs [3]. It is highly challenging to distinguish A. gossypii from its sister species Aphis frangulae as both are colonizing solanaceous plants as secondary hosts but the primary host of A. frangulae is Frangula alnus[4]. This paper describes a worldwide collection of both species from December 1989 to September 2019. Aphids were collected individually on plants (19 families) or in traps. The location, the morph type and the botanical family of the host plant were registered. DNA was extracted from each aphid and amplified at 8 microsatellite loci [5]. Amplicons were analysed with ABI technology and their size was defined with Genemapper software. We named each unique combination of alleles, called a multilocus genotype (MLG), and then each individual was given its MLG. The matrix of alleles of all MLGs was run for a Bayesian analysis to describe the genetic structure of the diversity collected and then each MLG had a probability to belong to a genetic group [6,7]. Probability of assignation to each genetic group revealed by the analysis was reported to each individual according to its MLG. This dataset can be used to analyze host plant specificities in A. gossypii, genetic diversity in A. gossypii and relative incidence of variants in diverse geographical regions, admixture between two sister species (Aphis gossypii and Aphis frangulae).

Genetic diversity of melon aphids Aphis gossypii associated with landscape features

Despite increasing evidence that landscape features strongly influence the abundance and dispersal of insect populations, landscape composition has seldom been explicitly linked to genetic structure. We conducted a genetic study of the melon aphid, Aphis gossypii, in two counties of Beijing, China during spring migration using samples from watermelon. We performed aphid genetic analysis using restriction site associated DNA sequencing (2b-RAD) and investigated the relationship between land cover and the genetic diversity. The percentage area of land cover (cropland, vegetable, orchard, grassland, woodland) was quantified in each particular scale (ranging from 0.5 km to 3 km) and was used as a predictor variable in our generalized linear models. We found a moderate level of genetic differentiation among nine sampled populations. Geographic distance and genetic distance were not significantly associated, indicating that geographic location was not a barrier to migration. These nine populations could be clustered depending on their level of genetic diversity (high and low). The genetic diversity (Shannon's information index) was positively correlated with grassland at the spatial scales of 1 and 2 km and negatively with orchard and vegetable at 0.5 and 1 km. Genetic diversity was best predicted by the grassland + orchard + vegetable model at a spatial scale of 1 km. Based on the method of relative weights, orchard land had the greatest relative importance, followed by grassland and vegetable land, in that order. This study contributes to our understanding of the genetic variation of aphids in agricultural landscapes.

Vat, an Amazing Gene Conferring Resistance to Aphids and Viruses They Carry: From Molecular Structure to Field Effects

We review half a century of research on Cucumis melo resistance to Aphis gossypii from molecular to field levels. The Vat gene is unique in conferring resistance to both A. gossypii and the viruses it transmits. This double phenotype is aphid clone-dependent and has been observed in 25 melon accessions, mostly from Asia. It is controlled by a cluster of genes including CC-NLR, which has been characterized in detail. Copy-number polymorphisms (for the whole gene and for a domain that stands out in the LLR region) and single-nucleotide polymorphisms have been identified in the Vat cluster. The role of these polymorphisms in plant/aphid interactions remains unclear. The Vat gene structure suggests a functioning with separate recognition and response phases. During the recognition phase, the VAT protein is thought to interact (likely indirectly) with an aphid effector introduced during cell puncture by the aphid. A few hours later, several miRNAs are upregulated in Vat plants. Peroxidase activity increases, and callose and lignin are deposited in the walls of the cells adjacent to the stylet path, disturbing aphid behavior. In aphids feeding on Vat plants, Piwi-interacting RNA-like sequences are abundant and the levels of other miRNAs are modified. At the plant level, resistance to aphids is quantitative (aphids escape the plant and display low rates of reproduction). Resistance to viruses is qualitative and local. Durability of NLR genes is highly variable. A. gossypii clones are adapted to Vat resistance, either by introducing a new effector that interferes with the deployment of plant defenses, or by adapting to the defenses it triggered. Viruses transmitted in a non-persistent manner cannot adapt to Vat resistance. At population level, Vat reduces aphid density and genetic diversity. The durability of Vat resistance to A. gossypii populations depends strongly on the agro-ecosystem, including, in particular, the presence of other cucurbit crops serving as alternative hosts for adapted clones in fall and winter. At the crop level, Vat resistance decreases the intensity of virus epidemics when A. gossypii is the main aphid vector in the crop environment.

Insight into the durability of plant resistance to aphids from a demo-genetic study of Aphis gossypii in melon crops

Resistance breakdown has been observed following the deployment of plant cultivars resistant to pests. Assessing the durability of a resistance requires long-term experiments at least at a regional scale. We collected such data for melon resistance conferred by the Vat gene cluster to melon aphids. We examined landscape-level populations of Aphis gossypii collected in 2004-2015, from melon-producing regions with and without the deployment of Vat resistance and with different climates. We conducted demo-genetic analyses of the aphid populations on Vat and non-Vat plants during the cropping seasons. The Vat resistance decreased the density of aphid populations in all areas and changed the genetic structure and composition of these populations. Two bottlenecks were identified in the dynamics of adapted clones, due to the low levels of production of dispersal morphs and winter extinction. Our results suggest that (i) Vat resistance will not be durable in the Lesser Antilles, where no bottleneck affected the dynamics of adapted clones, (ii) Vat resistance will be durable in south-west France, where both bottlenecks affected the dynamics of adapted clones and (iii) Vat resistance will be less durable in south-east France, where only one of the two bottlenecks was observed.

The Distribution and Host Shifts of Cotton-Melon Aphids in Northern China

Aphis gossypii Glover (Hemiptera: Aphididae) is a serious pest of cotton in northern China. A microsatellite analysis was used to characterize the genetic structure of A. gossypii populations from different geographic, host plant, and seasonal populations in 2014. Among 906 individuals, 507 multilocus genotypes were identified, with genotypic richness values of 0.07–1.00 for the populations. We observed moderate levels of genetic differentiation among geographic populations (F_ST = 0.103; 95% confidence interval: 0.065–0.145) and host plant populations (F_ST = 0.237; 95% confidence interval: 0.187–0.296). A Mantel test of isolation by distance revealed no significant correlations between Slatkin’s linearized F_ST and the natural logarithm of geographic distance. A Bayesian analysis of population genetic structures identified three clusters. An analysis of molecular variance revealed significant differences among the three clusters (F = 0.26596, P < 0.0001), among seasons (F = 0.04244, P = 0.00381), and among host populations (F = 0.12975, P = 0.0029). Thus, the A. gossypii populations in northern China exhibit considerable genotypic diversity. Additionally, our findings indicated that the 31 analyzed populations could be classified as one of three host biotypes (i.e., cotton, cucumber, and pomegranate biotypes). There were also clear seasonal effects on population genetic structure diversity among aphids collected from Anyang.

Identification of Aphis gossypii Glover (Hemiptera: Aphididae) Biotypes from Different Host Plants in North China

Background

The cotton-melon aphid, Aphis gossypii Glover (Hemiptera: Aphididae), is a polyphagous species with a worldwide distribution and a variety of biotypes. North China is a traditional agricultural area with abundant winter and summer hosts of A. gossypii. While the life cycles of A. gossypii on different plants have been well studied, those of the biotypes of North China are still unclear.

Results

Host transfer experiments showed that A. gossypii from North China has two host-specialized biotypes: cotton and cucumber. Based on complete mitochondrial sequences, we identified a molecular marker with five single-nucleotide polymorphisms to distinguish the biotypes. Using this marker, a large-scale study of biotypes on primary winter and summer hosts was conducted. All A. gossypii collected from three primary hosts—hibiscus, pomegranate, and Chinese prickly ash—were cotton biotypes, with more cotton-melon aphids found on hibiscus than the other two species. In May, alate cotton and cucumber biotypes coexisted on cotton and cucumber seedlings, but each preferred its natal host. Both biotypes existed on zucchini, although the cucumber biotype was more numerous. Aphids on muskmelon were all cucumber biotypes, whereas most aphids on kidney bean were cotton biotypes. Aphids on seedlings of potato and cowpea belong to other species. In August, aphids on cotton and cucumber were the respective biotypes, with zucchini still hosting both biotypes as before. Thus, the biotypes had different fitnesses on different host plants.

Conclusions

Two host-specialized biotypes (cotton and cucumber) are present in North China. Hibiscus, pomegranate, and Chinese prickly ash can serve as winter hosts for the cotton biotype but not the cucumber biotype in North China. The fitnesses of the two host-specialized biotypes differ on various summer hosts. When alate aphids migrate to summer hosts, they cannot accurately land on the corresponding plant.

The Vat locus encodes for a CC-NBS-LRR protein that confers resistance to Aphis gossypii infestation and A. gossypii-mediated virus resistance

Aphis gossypii is a polyphagous sucking aphid and a vector for many viruses. In Cucumis melo, a dominant locus, Vat, confers a high level of resistance to Aphis gossypii infestation and to viruses transmitted by this vector. To investigate the mechanism underlying this double resistance, we first genetically dissected the Vat locus. We delimited the double resistance to a single gene that encodes for a coiled-coil-nucleotide-binding-site-leucine-rich repeat (CC-NBS-LRR) protein type. To validate the genetic data, transgenic lines expressing the Vat gene were generated and assessed for the double resistance. In this analysis, Vat-transgenic plants were resistant to A. gossypii infestation as well as A. gossypii-mediated virus transmission. When the plants were infected mechanically, virus infection occurred on both transgenic and non-transgenic control plants. These results confirmed that the cloned CC-NBS-LRR gene mediates both resistance to aphid infestation and virus infection using A. gossypii as a vector. This resistance also invokes a separate recognition and response phases in which the recognition phase involves the interaction of an elicitor molecule from the aphid and Vat from the plant. The response phase is not specific and blocks both aphid infestation and virus infection. Sequence analysis of Vat alleles suggests a major role of an unusual conserved LRR repeat in the recognition of A. gossypii.

Low genetic diversity in Melanaphis sacchari aphid populations at the worldwide scale

Numerous studies have examined the genetic diversity and genetic structure of invading species, with contrasting results concerning the relative roles of genetic diversity and phenotypic plasticity in the success of introduced populations. Increasing evidence shows that asexual lineages of aphids are able to occupy a wide geographical and ecological range of habitats despite low genetic diversity. The anholocyclic aphid Melanaphis sacchari is a pest of sugarcane and sorghum which originated in the old world, was introduced into the Americas, and is now distributed worldwide. Our purpose was to assess the genetic diversity and structuring of populations of this species according to host and locality. We used 10 microsatellite markers to genotype 1333 individuals (57 samples, 42 localities, 15 countries) collected mainly on sugarcane or sorghum. Five multilocus lineages (MLL) were defined, grouping multilocus genotypes (MLG) differing by only a few mutations or scoring errors. Analysis of a 658 bp sequence of mitochondrial COI gene on 96 individuals revealed five haplotypes, with a mean divergence of only 0.19 %. The distribution of MLL appeared to be strongly influenced by geography but not by host plant. Each of the five MLL grouped individuals from (A) Africa, (B) Australia, (C) South America, the Caribbean and the Indian Ocean including East Africa, (D) USA, and (E) China. The MLL A and C, with a wide geographic distribution, matched the definition of superclone. Among aphids, M. sacchari has one of the lowest known rates of genetic diversity for such a wide geographical distribution.

Regular or covert sex defines two lineages and worldwide superclones within the leaf-curl plum aphid (Brachycaudus helichrysi, Kaltenbach)

Asexual reproduction occurs widely in plants and animals, particularly in insects. Aphid species usually reproduce by cyclic parthenogenesis, but many species include obligate asexual lineages. We recently showed that the leaf-curl plum aphid, Brachycaudus helichrysi, actually encompasses two lineages, B. helichrysi H1 and H2. Ecological data suggest that these lineages have different life cycles. We conducted a large population genetics study, based on 14 microsatellite loci, to infer their respective life cycles and investigate their population structure and geographical distribution. Brachycaudus helichrysi H1 displayed the genetic signature of cyclical parthenogenesis, using plum trees as primary hosts for sexual reproduction, as classically described for B. helichrysi. This global survey showed that the Central Asian population of H1 was clearly differentiated from American-European populations. By contrast, B. helichrysi H2 displayed the typical signature of obligate asexual reproduction. H2 encompassed at least eight highly successful genotypes or superclones. This lack of ability to undergo sexual reproduction was confirmed for one of the superclones by sex induction experiments. We found only one B. helichrysi H2 population that underwent sexual reproduction, which was collected from peach trees, in Northern India. Our results confirm that H1 and H2 have different life cycles. Brachycaudus helichrysi H1 is clearly heteroecious using plum trees as primary hosts, while B. helichrysi H2 encompasses several anholocyclic lineages, and some heteroecious populations that until now have only been found associated with peach trees as primary hosts. We discuss implications of these findings for the pest status of B. helichrysi lineages.