Accelerated evolution of sex chromosomes in aphids, an x0 system

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Evolutionary consequences of loss of sexual reproduction on male-related traits in parthenogenetic lineages of the pea aphid

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Gene Editing and Genetic Control of Hemipteran Pests: Progress, Challenges and Perspectives

The origin of the order Hemiptera can be traced to the late Permian Period more than 230 MYA, well before the origin of flowering plants 100 MY later in during the Cretaceous period. Hemipteran species consume their liquid diets using a sucking proboscis; for phytophagous hemipterans their mouthparts (stylets) are elegant structures that enable voracious feeding from plant xylem or phloem. This adaptation has resulted in some hemipteran species becoming globally significant pests of agriculture resulting in significant annual crop losses. Due to the reliance on chemical insecticides for the control of insect pests in agricultural settings, many hemipteran pests have evolved resistance to insecticides resulting in an urgent need to develop new, species-specific and environmentally friendly methods of pest control. The rapid advances in CRISPR/Cas9 technologies in model insects such as Drosophila melanogaster, Tribolium castaneum, Bombyx mori, and Aedes aegypti has spurred a new round of innovative genetic control strategies in the Diptera and Lepidoptera and an increased interest in assessing genetic control technologies for the Hemiptera. Genetic control approaches in the Hemiptera have, to date, been largely overlooked due to the problems of introducing genetic material into the germline of these insects. The high frequency of CRISPR-mediated mutagenesis in model insect species suggest that, if the delivery problem for Hemiptera could be solved, then gene editing in the Hemiptera might be quickly achieved. Significant advances in CRISPR/Cas9 editing have been realized in nine species of Hemiptera over the past 4 years. Here we review progress in the Hemiptera and discuss the challenges and opportunities for extending contemporary genetic control strategies into species in this agriculturally important insect orderr.

Detection of sex in adults and larvae of Leptinotarsa decemlineata on principle of copy number variation

The identification of sex in larvae of insects is usually challenging or even impossible, while in adults the sexual dimorphism is usually evident. Here, we used copy number analysis to develop a method of sex detection in Colorado potato beetle (Leptinotarsa decemlineata), which has an X0 sex determination system. The X linked gene LdVssc and autosomal gene LdUBE3B were identified as appropriate target and reference loci, respectively. The copy numbers (CNV) of LdVssc in males and females were estimated using standard droplet digital PCR (ddPCR) and real-time PCR (qPCR). With both methods, CNVs were bimodally distributed (BA_ddPCR = 0.709 and BA_qPCR = 0.683) with 100% ability to distinguish females from males. The use of qPCR-based sex detection in a broad collection of 448 random CPB adults showed a perfect association (Phi = 1.0, p < 0.05) with the true sexes of adults, with mean CNV in females of 2.032 (SD = 0.227) and 0.989 in males (SD = 0.147). In the collection of 50 random 4th instar larvae, 27 females and 23 males were identified, consistent with the expected 1:1 sex ratio (p = 0.689). The method is suitable for sexing in all stages of ontogenesis. The optimal cost-effective application of the method in large populations requires the DNA extraction using CTAB, the qPCR assay in one biological replicate and three technical replicates of each marker, and the use of one randomly chosen male per run to calibrate calculation of CNV.

Comparative analysis of chromosome numbers and sex chromosome systems in Paraneoptera (Insecta)

This article is part (the 4th article) of the themed issue (a monograph) "Aberrant cytogenetic and reproductive patterns in the evolution of Paraneoptera". The purpose of this article is to consider chromosome structure and evolution, chromosome numbers and sex chromosome systems, which all together constitute the chromosomal basis of reproduction and are essential for reproductive success. We are based on our own observations and literature data available for all major lineages of Paraneoptera including Zoraptera (angel insects), Copeognatha (=Psocoptera; bark lice), Parasita (=Phthiraptera s. str; true lice), Thysanoptera (thrips), Homoptera (scale insects, aphids, jumping plant-lice, whiteflies, and true hoppers), Heteroptera (true bugs), and Coleorrhyncha (moss bugs). Terminology, nomenclature, classification, and the study methods are given in the first paper of the issue (Gavrilov-Zimin et al. 2021).

No evidence for genetic sex determination in Daphnia magna

Mechanisms of sex determination (SD) differ widely across the tree of life. In genotypic sex determination (GSD), genetic elements determine whether individuals are male or female, while in environmental sex determination (ESD), external cues control the sex of the offspring. In cyclical parthenogens, females produce mostly asexual daughters, but environmental stimuli such as crowding, temperature or photoperiod may cause them to produce sons. In aphids, sons are induced by ESD, even though GSD is present, with females carrying two X chromosomes and males only one (X0 SD system). By contrast, although ESD exists in Daphnia, the two sexes were suggested to be genetically identical, based on a 1972 study on Daphnia magna (2n=20) that used three allozyme markers. This study cannot, however, rule out an X0 system, as all three markers may be located on autosomes. Motivated by the life cycle similarities of Daphnia and aphids, and the absence of karyotype information for Daphnia males, we tested for GSD (homomorphic sex chromosomes and X0) systems in D. magna using a whole-genome approach by comparing males and females of three genotypes. Our results confirm the absence of haploid chromosomes or haploid genomic regions in D. magna males as well as the absence of sex-linked genomic regions and sex-specific single-nucleotide polymorphisms. Within the limitations of the three studied populations here and the methods used, we suggest that our results make the possibility of genetic differences among sexes in the widely used Daphnia model system very unlikely.

Male-Biased microRNA Discovery in the Pea Aphid

Epigenetic mechanisms modulate gene expression levels during development, shaping how a single genome produces a diversity of phenotypes. Here, we begin to explore the epigenetic regulation of sexual dimorphism in pea aphids (Acyrthosiphon pisum) by focusing on microRNAs. Previous analyses of microRNAs in aphids have focused solely on females, so we performed deep sequencing of a sample containing early-stage males. We used this sample, plus samples from Genbank, to find 207 novel pea aphid microRNA coding loci. We localized microRNA loci to a chromosome-level assembly of the pea aphid genome and found that those on the X chromosome have lower overall expression compared to those on autosomes. We then identified a set of 19 putative male-biased microRNAs and found them enriched on the X chromosome. Finally, we performed protein-coding RNA-Seq of first instar female and male pea aphids to identify genes with lower expression in males. 10 of these genes were predicted targets of the 19 male-biased microRNAs. Our study provides the most complete set of microRNAs in the pea aphid to date and serves as foundational work for future studies on the epigenetic control of sexual dimorphism.

The Aphid X Chromosome Is a Dangerous Place for Functionally Important Genes: Diverse Evolution of Hemipteran Genomes Based on Chromosome-Level Assemblies

Different evolutionary forces shape gene content and sequence evolution on autosomes versus sex chromosomes. Location on a sex chromosome can favor male-beneficial or female-beneficial mutations depending on the sex determination system and selective pressure on different sexual morphs. An X0 sex determination can lead to autosomal enrichment of male-biased genes, as observed in some hemipteran insect species. Aphids share X0 sex determination; however, models predict the opposite pattern, due to their unusual life cycles, which alternate between all-female asexual generations and a single sexual generation. Predictions include enrichment of female-biased genes on autosomes and of male-biased genes on the X, in contrast to expectations for obligately sexual species. Robust tests of these models require chromosome-level genome assemblies for aphids and related hemipterans with X0 sex determination and obligate sexual reproduction. In this study, we built the first chromosome-level assembly of a psyllid, an aphid relative with X0 sex determination and obligate sexuality, and compared it with recently resolved chromosome-level assemblies of aphid genomes. Aphid and psyllid X chromosomes differ strikingly. In aphids, female-biased genes are strongly enriched on autosomes and male-biased genes are enriched on the X. In psyllids, male-biased genes are enriched on autosomes. Furthermore, functionally important gene categories of aphids are enriched on autosomes. Aphid X-linked genes and male-biased genes are under relaxed purifying selection, but gene content and order on the X is highly conserved, possibly reflecting constraints imposed by unique chromosomal mechanisms associated with the unusual aphid life cycle.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Analysis of the extent of synteny and conservation in the gene order in aphids: A first glimpse from the Aphis glycines genome

In the last decade several insect genomes have been sequenced, but for most the chromosomal mapping of the identified scaffolds/annotated genes is not available. The lack of this information makes it difficult to analyse various genetic aspects, including the presence of genome rearrangements and the extent of synteny within and across species. We mapped five multigenic DNA families (major and minor rDNAs, histone gene cluster, esterases and carotenoid desaturases) and seven scaffolds corresponding to 9 Mb of the soybean aphid, Aphis glycines, genome and identified loci spanning the four soybean aphid chromosomes. A comparative analysis of the localization of the annotated A. glycines genes with respect to the peach potato aphid, Myzus persicae, and the fly, Drosophila melanogaster, evidenced a lower degree of synteny between the two aphid species than in the aphid-fly comparison. Only 1.4 genes per syntenic block were observed in aphids in contrast to 2.3 genes per block in flies. This higher chromosomal rearrangement rate in aphids could be explained considering that they possess holocentric chromosomes that can favour the stabilization and inheritance of chromosomal rearrangements. Lastly, our experiments did not detect the presence of chimeric assemblies in the newly available A. glycines biotype 1 genome, differently from what reported in assembled genome of other aphid species, suggesting that chromosomal mapping can be used to ascertain the quality of assembled genomes.

Evidence for Faster X Chromosome Evolution in Spiders

In species with chromosomal sex determination, X chromosomes are predicted to evolve faster than autosomes because of positive selection on recessive alleles or weak purifying selection. We investigated X chromosome evolution in Stegodyphus spiders that differ in mating system, sex ratio, and population dynamics. We assigned scaffolds to X chromosomes and autosomes using a novel method based on flow cytometry of sperm cells and reduced representation sequencing. We estimated coding substitution patterns (dN/dS) in a subsocial outcrossing species (S. africanus) and its social inbreeding and female-biased sister species (S. mimosarum), and found evidence for faster-X evolution in both species. X chromosome-to-autosome diversity (piX/piA) ratios were estimated in multiple populations. The average piX/piA estimates of S. africanus (0.57 [95% CI: 0.55-0.60]) was lower than the neutral expectation of 0.75, consistent with more hitchhiking events on X-linked loci and/or a lower X chromosome mutation rate, and we provide evidence in support of both. The social species S. mimosarum has a significantly higher piX/piA ratio (0.72 [95% CI: 0.65-0.79]) in agreement with its female-biased sex ratio. Stegodyphus mimosarum also have different piX/piA estimates among populations, which we interpret as evidence for recurrent founder events. Simulations show that recurrent founder events are expected to decrease the piX/piA estimates in S. mimosarum, thus underestimating the true effect of female-biased sex ratios. Finally, we found lower synonymous divergence on X chromosomes in both species, and the male-to-female substitution ratio to be higher than 1, indicating a higher mutation rate in males.

Disentangling the Causes for Faster-X Evolution in Aphids

The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima’s D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution.

Coalescence and Linkage Disequilibrium in Facultatively Sexual Diploids

Under neutrality, linkage disequilibrium results from physically linked sites having nonindependent coalescent histories. In obligately sexual organisms, meiotic recombination is the dominant force separating linked variants from one another, and thus in determining the decay of linkage disequilibrium with physical distance. In facultatively sexual diploid organisms that principally reproduce clonally, mechanisms of mitotic exchange are expected to become relatively more important in shaping linkage disequilibrium. Here we outline mathematical and computational models of a facultative-sex coalescent process that includes meiotic and mitotic recombination, via both crossovers and gene conversion, to determine how linkage disequilibrium is affected with facultative sex. We demonstrate that the degree to which linkage disequilibrium is broken down by meiotic recombination simply scales with the probability of sex if it is sufficiently high (much greater than [Formula: see text] for population size N). However, with very rare sex (occurring with frequency on the order of [Formula: see text]), mitotic gene conversion plays a particularly important and complicated role because it both breaks down associations between sites and removes within-individual diversity. Strong population structure under rare sex leads to lower average linkage disequilibrium values than in panmictic populations, due to the influence of low-frequency polymorphisms created by allelic sequence divergence acting in individual subpopulations. These analyses provide information on how to interpret observed linkage disequilibrium patterns in facultative sexuals and to determine what genomic forces are likely to shape them.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

Dosage compensation and sex-specific epigenetic landscape of the X chromosome in the pea aphid

Background

Heterogametic species display a differential number of sex chromosomes resulting in imbalanced transcription levels for these chromosomes between males and females. To correct this disequilibrium, dosage compensation mechanisms involving gene expression and chromatin accessibility regulations have emerged throughout evolution. In insects, these mechanisms have been extensively characterized only in Drosophila but not in insects of agronomical importance. Aphids are indeed major pests of a wide range of crops. Their remarkable ability to switch from asexual to sexual reproduction during their life cycle largely explains the economic losses they can cause. As heterogametic insects, male aphids are X0, while females (asexual and sexual) are XX.

Results

Here, we analyzed transcriptomic and open chromatin data obtained from whole male and female individuals to evaluate the putative existence of a dosage compensation mechanism involving differential chromatin accessibility of the pea aphid’s X chromosome. Transcriptomic analyses first showed X/AA and XX/AA expression ratios for expressed genes close to 1 in males and females, respectively, suggesting dosage compensation in the pea aphid. Analyses of open chromatin data obtained by Formaldehyde-Assisted Isolation of Regulatory Elements (FAIRE-seq) revealed a X chromosome chromatin accessibility globally and significantly higher in males than in females, while autosomes’ chromatin accessibility is similar between sexes. Moreover, chromatin environment of X-linked genes displaying similar expression levels in males and females—and thus likely to be compensated—is significantly more accessible in males.

Conclusions

Our results suggest the existence of an underlying epigenetic mechanism enhancing the X chromosome chromatin accessibility in males to allow X-linked gene dose correction between sexes in the pea aphid, similar to Drosophila. Our study gives new evidence into the comprehension of dosage compensation in link with chromatin biology in insects and newly in a major crop pest, taking benefits from both transcriptomic and open chromatin data.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0137-1) contains supplementary material, which is available to authorized users.

Chromosome Painting in Triatomine Insects Reveals Shared Sequences Between X Chromosomes and Autosomes

In order to provide a broad picture on the origin and evolution of holocentric X chromosomes in heteropteran species, we prepared a sex chromosome painting probe by microdissection of the X₁ and X₂ chromosomes from a kissing bug Mepraia spinolai (Hemiptera: Reduviidae: Triatominae). Fluorescence in situ hybridization on four species of the Triatomini having different amounts of autosomal heterochromatin and sex chromosome systems show that the Xs probe hybridizes on the euchromatin, located both on autosomes and X chromosomes. The heterochromatic Y chromosome and autosomal heterochromatic regions always appear free of hybridization signals. The hybridization results of the Xs probe on Rhodnius prolixus (Rhodniini) is completely different to that observed in Triatomini species. The hybridization signals are small and scattered on all euchromatin, without specific regions including the X chromosome. These results are in accordance with previous data obtained by genomic in situ hybridization and fluorescent banding, suggesting a clear differentiation in the repeat sequence composition of both sex chromosomes between Triatomini and Rhodniini tribes. These results also support that each sex chromosome in Triatomini has evolved independently from different autosomal pairs of a common ancestor, as described in other insect orders.

The microRNA toolkit of insects

Is there a correlation between miRNA diversity and levels of organismic complexity? Exhibiting extraordinary levels of morphological and developmental complexity, insects are the most diverse animal class on earth. Their evolutionary success was in particular shaped by the innovation of holometabolan metamorphosis in endopterygotes. Previously, miRNA evolution had been linked to morphological complexity, but astonishing variation in the currently available miRNA complements of insects made this link unclear. To address this issue, we sequenced the miRNA complement of the hemimetabolan Blattella germanica and reannotated that of two other hemimetabolan species, Locusta migratoria and Acyrthosiphon pisum, and of four holometabolan species, Apis mellifera, Tribolium castaneum, Bombyx mori and Drosophila melanogaster. Our analyses show that the variation of insect miRNAs is an artefact mainly resulting from poor sampling and inaccurate miRNA annotation, and that insects share a conserved microRNA toolkit of 65 families exhibiting very low variation. For example, the evolutionary shift toward a complete metamorphosis was accompanied only by the acquisition of three and the loss of one miRNA families.

Contrasting Levels of Molecular Evolution on the Mouse X Chromosome

The mammalian X chromosome has unusual evolutionary dynamics compared to autosomes. Faster-X evolution of spermatogenic protein-coding genes is known to be most pronounced for genes expressed late in spermatogenesis, but it is unclear if these patterns extend to other forms of molecular divergence. We tested for faster-X evolution in mice spanning three different forms of molecular evolution-divergence in protein sequence, gene expression, and DNA methylation-across different developmental stages of spermatogenesis. We used FACS to isolate individual cell populations and then generated cell-specific transcriptome profiles across different stages of spermatogenesis in two subspecies of house mice (Mus musculus), thereby overcoming a fundamental limitation of previous studies on whole tissues. We found faster-X protein evolution at all stages of spermatogenesis and faster-late protein evolution for both X-linked and autosomal genes. In contrast, there was less expression divergence late in spermatogenesis (slower late) on the X chromosome and for autosomal genes expressed primarily in testis (testis-biased). We argue that slower-late expression divergence reflects strong regulatory constraints imposed during this critical stage of sperm development and that these constraints are particularly acute on the tightly regulated sex chromosomes. We also found slower-X DNA methylation divergence based on genome-wide bisulfite sequencing of sperm from two species of mice (M. musculus and M. spretus), although it is unclear whether slower-X DNA methylation reflects development constraints in sperm or other X-linked phenomena. Our study clarifies key differences in patterns of regulatory and protein evolution across spermatogenesis that are likely to have important consequences for mammalian sex chromosome evolution, male fertility, and speciation.

The X Chromosome of Hemipteran Insects: Conservation, Dosage Compensation and Sex-Biased Expression

Insects of the order Hemiptera (true bugs) use a wide range of mechanisms of sex determination, including genetic sex determination, paternal genome elimination, and haplodiploidy. Genetic sex determination, the prevalent mode, is generally controlled by a pair of XY sex chromosomes or by an XX/X0 system, but different configurations that include additional sex chromosomes are also present. Although this diversity of sex determining systems has been extensively studied at the cytogenetic level, only the X chromosome of the model pea aphid Acyrthosiphon pisum has been analyzed at the genomic level, and little is known about X chromosome biology in the rest of the order.

In this study, we take advantage of published DNA- and RNA-seq data from three additional Hemiptera species to perform a comparative analysis of the gene content and expression of the X chromosome throughout this clade. We find that, despite showing evidence of dosage compensation, the X chromosomes of these species show female-biased expression, and a deficit of male-biased genes, in direct contrast to the pea aphid X. We further detect an excess of shared gene content between these very distant species, suggesting that despite the diversity of sex determining systems, the same chromosomal element is used as the X throughout a large portion of the order.

Parallel faster-X evolution of gene expression and protein sequences in Drosophila: beyond differences in expression properties and protein interactions

Population genetics models predict that the X (or Z) chromosome will evolve at faster rates than the autosomes in XY (or ZW) systems. Studies of molecular evolution using large datasets in multiple species have provided evidence supporting this faster-X effect in protein-coding sequences and, more recently, in transcriptomes. However, X-linked and autosomal genes differ significantly in important properties besides hemizygosity in males, including gene expression levels, tissue specificity in gene expression, and the number of interactions in which they are involved (i.e., protein-protein or DNA-protein interactions). Most important, these properties are known to correlate with rates of evolution, which raises the question of whether differences between the X chromosome and autosomes in gene properties, rather than hemizygosity, are sufficient to explain faster-X evolution. Here I investigate this possibility using whole genome sequences and transcriptomes of Drosophila yakuba and D. santomea and show that this is not the case. Additional factors are needed to account for faster-X evolution of both gene expression and protein-coding sequences beyond differences in gene properties, likely a higher incidence of positive selection in combination with the accumulation of weakly deleterious mutations.

Genetic control of contagious asexuality in the pea aphid

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of such shifts toward asexuality remain largely unknown. We addressed this issue in an aphid species where both sexual and obligate asexual lineages coexist in natural populations. These sexual and asexual lineages may occasionally interbreed because some asexual lineages maintain a residual production of males potentially able to mate with the females produced by sexual lineages. Hence, this species is an ideal model to study the genetic basis of the loss of sexual reproduction with quantitative genetic and population genomic approaches. Our analysis of the co-segregation of ∼ 300 molecular markers and reproductive phenotype in experimental crosses pinpointed an X-linked region controlling obligate asexuality, this state of character being recessive. A population genetic analysis (>400-marker genome scan) on wild sexual and asexual genotypes from geographically distant populations under divergent selection for reproductive strategies detected a strong signature of divergent selection in the genomic region identified by the experimental crosses. These population genetic data confirm the implication of the candidate region in the control of reproductive mode in wild populations originating from 700 km apart. Patterns of genetic differentiation along chromosomes suggest bidirectional gene flow between populations with distinct reproductive modes, supporting contagious asexuality as a prevailing route to permanent parthenogenesis in pea aphids. This genetic system provides new insights into the mechanisms of coexistence of sexual and asexual aphid lineages.

Genomic regions repeatedly involved in divergence among plant-specialized pea aphid biotypes

Understanding the genetic bases of biological diversification is a long-standing goal in evolutionary biology. Here, we investigate whether replicated cases of adaptive divergence involve the same genomic regions in the pea aphid, Acyrthosiphon pisum, a large complex of genetically differentiated biotypes, each specialized on different species of legumes. A previous study identified genomic regions putatively involved in host-plant adaptation and/or reproductive isolation by performing a hierarchical genome scan in three biotypes. This led to the identification of 11 F(ST) outliers among 390 polymorphic microsatellite markers. In this study, the outlier status of these 11 loci was assessed in eight biotypes specialized on other host plants. Four of the 11 previously identified outliers showed greater genetic differentiation among these additional biotypes than expected under the null hypothesis of neutral evolution (α < 0.01). Whether these hotspots of genomic divergence result from adaptive events, intrinsic barriers or reduced recombination is discussed.

The faster-X effect: integrating theory and data

Population genetics theory predicts that X (or Z) chromosomes could play disproportionate roles in speciation and evolutionary divergence, and recent genome-wide analyses have identified situations in which X or Z-linked divergence exceeds that on the autosomes (the so-called 'faster-X effect'). Here, we summarize the current state of both the theory and data surrounding the study of faster-X evolution. Our survey indicates that the faster-X effect is pervasive across a taxonomically diverse array of evolutionary lineages. These patterns could be informative of the dominance or recessivity of beneficial mutations and the nature of genetic variation acted upon by natural selection. We also identify several aspects of disagreement between these empirical results and the population genetic models used to interpret them. However, there are clearly delineated aspects of the problem for which additional modeling and collection of genomic data will address these discrepancies and provide novel insights into the population genetics of adaptation.

Masculinization of the x chromosome in the pea aphid

Evolutionary theory predicts that sexually antagonistic mutations accumulate differentially on the X chromosome and autosomes in species with an XY sex-determination system, with effects (masculinization or feminization of the X) depending on the dominance of mutations. Organisms with alternative modes of inheritance of sex chromosomes offer interesting opportunities for studying sexual conflicts and their resolution, because expectations for the preferred genomic location of sexually antagonistic alleles may differ from standard systems. Aphids display an XX/X0 system and combine an unusual inheritance of the X chromosome with the alternation of sexual and asexual reproduction. In this study, we first investigated theoretically the accumulation of sexually antagonistic mutations on the aphid X chromosome. Our results show that i) the X is always more favourable to the spread of male-beneficial alleles than autosomes, and should thus be enriched in sexually antagonistic alleles beneficial for males, ii) sexually antagonistic mutations beneficial for asexual females accumulate preferentially on autosomes, iii) in contrast to predictions for standard systems, these qualitative results are not affected by the dominance of mutations. Under the assumption that sex-biased gene expression evolves to solve conflicts raised by the spread of sexually antagonistic alleles, one expects that male-biased genes should be enriched on the X while asexual female-biased genes should be enriched on autosomes. Using gene expression data (RNA-Seq) in males, sexual females and asexual females of the pea aphid, we confirm these theoretical predictions. Although other mechanisms than the resolution of sexual antagonism may lead to sex-biased gene expression, we argue that they could hardly explain the observed difference between X and autosomes. On top of reporting a strong masculinization of the aphid X chromosome, our study highlights the relevance of organisms displaying an alternative mode of sex chromosome inheritance to understanding the forces shaping chromosome evolution.

Widespread selection across coding and noncoding DNA in the pea aphid genome

Genome-wide patterns of diversity and selection are critical measures for understanding how evolution has shaped the genome. Yet, these population genomic estimates are available for only a limited number of model organisms. Here we focus on the population genomics of the pea aphid (Acyrthosiphon pisum). The pea aphid is an emerging model system that exhibits a range of intriguing biological traits not present in classic model systems. We performed low-coverage genome resequencing of 21 clonal pea aphid lines collected from alfalfa host plants in North America to characterize genome-wide patterns of diversity and selection. We observed an excess of low-frequency polymorphisms throughout coding and noncoding DNA, which we suggest is the result of a founding event and subsequent population expansion in North America. Most gene regions showed lower levels of Tajima's D than synonymous sites, suggesting that the majority of the genome is not evolving neutrally but rather exhibits significant constraint. Furthermore, we used the pea aphid's unique manner of X-chromosome inheritance to assign genomic scaffolds to either autosomes or the X chromosome. Comparing autosomal vs. X-linked sequence variation, we discovered that autosomal genes show an excess of low frequency variants indicating that purifying selection acts more efficiently on the X chromosome. Overall, our results provide a critical first step in characterizing the genetic diversity and evolutionary pressures on an aphid genome.

An excess of gene expression divergence on the X chromosome in Drosophila embryos: implications for the faster-X hypothesis

The X chromosome is present as a single copy in the heterogametic sex, and this hemizygosity is expected to drive unusual patterns of evolution on the X relative to the autosomes. For example, the hemizgosity of the X may lead to a lower chromosomal effective population size compared to the autosomes, suggesting that the X might be more strongly affected by genetic drift. However, the X may also experience stronger positive selection than the autosomes, because recessive beneficial mutations will be more visible to selection on the X where they will spend less time being masked by the dominant, less beneficial allele--a proposal known as the faster-X hypothesis. Thus, empirical studies demonstrating increased genetic divergence on the X chromosome could be indicative of either adaptive or non-adaptive evolution. We measured gene expression in Drosophila species and in D. melanogaster inbred strains for both embryos and adults. In the embryos we found that expression divergence is on average more than 20% higher for genes on the X chromosome relative to the autosomes; but in contrast, in the inbred strains, gene expression variation is significantly lower on the X chromosome. Furthermore, expression divergence of genes on Muller's D element is significantly greater along the branch leading to the obscura sub-group, in which this element segregates as a neo-X chromosome. In the adults, divergence is greatest on the X chromosome for males, but not for females, yet in both sexes inbred strains harbour the lowest level of gene expression variation on the X chromosome. We consider different explanations for our results and conclude that they are most consistent within the framework of the faster-X hypothesis.

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome-wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host-adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.