Comparison of gene repertoires and patterns of evolutionary rates in eight aphid species that differ by reproductive mode

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Consequences of Asexuality in Natural Populations: Insights from Stick Insects

Recombination is a fundamental process with significant impacts on genome evolution. Predicted consequences of the loss of recombination include a reduced effectiveness of selection, changes in the amount of neutral polymorphisms segregating in populations, and an arrest of GC-biased gene conversion. Although these consequences are empirically well documented for nonrecombining genome portions, it remains largely unknown if they extend to the whole genome scale in asexual organisms. We identify the consequences of asexuality using de novo transcriptomes of five independently derived, obligately asexual lineages of stick insects, and their sexual sister-species. We find strong evidence for higher rates of deleterious mutation accumulation, lower levels of segregating polymorphisms and arrested GC-biased gene conversion in asexuals as compared with sexuals. Taken together, our study conclusively shows that predicted consequences of genome evolution under asexuality can indeed be found in natural populations.

No signal of deleterious mutation accumulation in conserved gene sequences of extant asexual hexapods

Loss of sex and recombination is generally assumed to impede the effectiveness of purifying selection and to result in the accumulation of slightly deleterious mutations. Empirical evidence for this has come from several studies investigating mutational load in a small number of individual genes. However, recent whole transcriptome based studies have yielded inconsistent results, hence questioning the validity of the assumption of mutational meltdown in asexual populations. Here, we study the effectiveness of purifying selection in eight asexual hexapod lineages and their sexual relatives, as present in the 1 K Insect Transcriptome Evolution (1KITE) project, covering eight hexapod groups. We analyse the accumulation of slightly deleterious nonsynonymous and synonymous point mutations in 99 single copy orthologue protein-coding loci shared among the investigated taxa. While accumulation rates of nonsynonymous mutations differed between genes and hexapod groups, we found no effect of reproductive mode on the effectiveness of purifying selection acting at nonsynonymous and synonymous sites. Although the setup of this study does not fully rule out nondetection of subtle effects, our data does not support the established consensus of asexual lineages undergoing ‘mutational meltdown’.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Effective purifying selection in ancient asexual oribatid mites

Sex is beneficial in the long term because it can prevent mutational meltdown through increased effectiveness of selection. This idea is supported by empirical evidence of deleterious mutation accumulation in species with a recent transition to asexuality. Here, we study the effectiveness of purifying selection in oribatid mites which have lost sex millions of years ago and diversified into different families and species while reproducing asexually. We compare the accumulation of deleterious nonsynonymous and synonymous mutations between three asexual and three sexual lineages using transcriptome data. Contrasting studies of young asexual lineages, we find evidence for strong purifying selection that is more effective in asexual as compared to sexual oribatid mite lineages. Our results suggest that large populations likely sustain effective purifying selection and facilitate the escape of mutational meltdown in the absence of sex. Thus, sex per se is not a prerequisite for the long-term persistence of animal lineages.

Asexual reproduction is thought to be an evolutionary dead end in eukaryotes because deleterious mutations will not be purged effectively. Here, Brandt and colleagues show that anciently asexual oribatid mites in fact have reduced accumulation of deleterious mutations compared to their sexual relatives.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Comparative transcriptomic analyses revealed divergences of two agriculturally important aphid species

BACKGROUND

Grain aphid (Sitobion avenae F) and pea aphid (Acyrthosiphon pisum) are two agriculturally important pest species, which cause significant yield losses to crop plants each year by inflicting damage both through the direct effects of feeding and by vectoring debilitating plant viruses. Although a close phylogenetic relationship between grain aphid and pea aphid was proposed, the biological variations between these two aphid species are obvious. While the host ranges of grain aphid is restricted to cereal crops and in particular wheat, that of pea aphid is wider, mainly colonizing leguminous plant species. Until now, the genetic factors underlying the divergence between grain aphid and pea aphid still remain unclear due to the limited genomic data of grain aphid available in public databases.

RESULTS

Based on a set of transcriptome data of grain aphid generated by using Roche 454 GS-FLX pyrosequencing, comparative analysis between this set of transcriptome data of grain aphid and mRNA sequences of pea aphid available in the public databases was performed. Compared with mRNA sequences of pea aphid, 4,857 unigenes were found to be specifically presented in the transcriptome of grain aphid under the rearing conditions described in this study. Furthermore, 3,368 orthologous pairs which could be calculated with both nonsynonymous (Ka) and synonymous (Ks) substitutions were used to infer their sequence divergences. The average differences in the coding, 5' and 3' untranslated regions of these orthologs were 10.53%, 21.29% and 18.96%, respectively. Moreover, of 340 orthologs which were identified to have evolved in response to positive selection based on the rates of Ka and Ks substitutions, 186 were predicted to be involved in secondary metabolism and xenobiotic metabolisms which might contribute to the divergence of these two aphid species.

CONCLUSIONS

The comprehensive transcriptome divergent sequence analysis between grain aphid and pea aphid provides an invaluable resource for the investigation of genes involved in host plant adaptation and evolution. Moreover, the demonstration of divergent transcriptome sequences between grain aphid and pea aphid pave the way for the investigation of the molecular mechanisms underpinning the biological variations of these two agriculturally important aphid species.

Mating systems and selection efficacy: a test using chloroplastic sequence data in Angiosperms

Selfing is assumed to reduce selection efficacy, especially purifying selection. This can be tested using molecular data, for example by comparing the Dn/Ds ratio between selfing and outcrossing lineages. So far, little evidence of relaxed selection against weakly deleterious mutations (as inferred by a higher Dn/Ds ratio) in selfers as compared to outcrossers has been found, contrary to the pattern often observed between asexual and sexual lineages. However, few groups have been studied to date. To further test this hypothesis, we compiled and analysed chloroplastic sequence data sets in several plant groups. We found a general trend towards relaxed selection in selfers in our data sets but with weak statistical support. Simulations suggested that the results were compatible with weak-to-moderate Dn/Ds ratio differences in selfing lineages. Simple theoretical predictions also showed that the ability to detect relaxed selection in selfers could strongly depend on the distribution of the effects of deleterious mutations on fitness. Our results are compatible with a recent origin of selfing lineages whereby deleterious mutations potentially have a strong impact on population extinction or with a more ancient origin but without a marked effect of deleterious mutations on the extinction dynamics.

Aphid protein effectors promote aphid colonization in a plant species-specific manner

Microbial pathogens and pests produce effectors to modulate host processes. Aphids are phloem-feeding insects, which introduce effectors via saliva into plant cells. However, it is not known if aphid effectors have adapted to modulate processes in specific plant species. Myzus persicae is a polyphagous insect that colonizes Arabidopsis thaliana and Nicotiana benthamiana, while the pea aphid Acyrthosiphon pisum specializes on colonizing plant species of the family Fabaceae. We found that M. persicae reproduction increased on transgenic Arabidopsis, producing the M. persicae effectors C002, PIntO1 (Mp1), and PIntO2 (Mp2), whereas reproduction of M. persicae did not increase on Arabidopsis producing the A. pisum orthologs of these three proteins. Plant-mediated RNA interference experiments showed that c002- and PIntO2-silenced M. persicae produce less progeny on Arabidopsis and N. benthamiana than nonsilenced aphids. Orthologs of c002, PIntO1, and PIntO2 were identified in multiple aphid species with dissimilar plant host ranges. We revealed high nonsynonymous versus synonymous nucleotide substitution rates within the effector orthologs, indicating that the effectors are fast evolving. Application of maximum likelihood methods identified specific sites with high probabilities of being under positive selection in PIntO1, whereas those of C002 and PIntO2 may be located in alignment gaps. In support of the latter, a M. persicae c002 mutant without the NDNQGEE repeat region, which overlaps with an alignment gap in C002, does not promote M. persicae colonization on Arabidopsis. Taken together, these results provide evidence that aphid effectors are under positive selection to promote aphid colonization on specific plant species.

Characterization of transcriptomes from sexual and asexual lineages of a New Zealand snail (Potamopyrgus antipodarum)

Understanding the evolution and maintenance of sexual reproduction is one of the central challenges of evolutionary biology, yet we know very little about how sex influences molecular evolution. The New Zealand freshwater snail Potamopyrgus antipodarum is ideally suited to address this knowledge gap because obligately sexual individuals often coexist with multiple independently derived obligately asexual lineages. This unusual situation allows direct comparisons both between sexual and asexual P. antipodarum and across populations that differ in the relative frequency of sexual individuals. As such, P. antipodarum has received a great deal of attention as a model system for the maintenance of sex in nature and is also used as a model for environmental toxicology and biological invasions. Molecular genetic resources for P. antipodarum will thus be useful to investigators in a variety of biological fields. We used 454 sequencing of cDNA libraries to generate transcriptomes from two sexual and two asexual P. antipodarum lineages. A de novo assembly of 116.7 Mb of sequence reads produced 41 396 contigs, and sequence similarity-based Gene Ontology annotations were obtained for 3740 contigs. We detected 408 315 SNP loci and 7315 microsatellite loci, which together represent the first genome-scale resource available for P. antipodarum. Raw 454 read sequences, contig sequences, annotation data and polymorphism data are publicly available in a searchable online database and for download at http://www.biology.uiowa.edu/neiman/transcriptome.php.