Ryan D. Bickel

Developmental Origin and Evolution of Bacteriocytes in the Aphid–Buchnera Symbiosis

Symbiotic relationships between bacteria and insect hosts are common. Although the bacterial endosymbionts have been subjected to intense investigation, little is known of the host cells in which they reside, the bacteriocytes. We have studied the development and evolution of aphid bacteriocytes, the host cells that contain the endosymbiotic bacteria Buchnera aphidicola. We show that bacteriocytes of Acyrthosiphon pisum express several gene products (or their paralogues): Distal-less, Ultrabithorax/Abdominal-A, and Engrailed. Using these markers, we find that a subpopulation of the bacteriocytes is specified prior to the transmission of maternal bacteria to the embryo. In addition, we discovered that a second population of cells is recruited to the bacteriocyte fate later in development. We experimentally demonstrate that bacteriocyte induction and proliferation occur independently of B. aphidicola. Major features of bacteriocyte development, including the two-step recruitment of bacteriocytes, have been conserved in aphids for 80–150 million years. Furthermore, we have investigated two cases of evolutionary loss of bacterial symbionts: in one case, where novel extracellular, eukaryotic symbionts replaced the bacteria, the bacteriocyte is maintained; in another case, where symbionts are absent, the bacteriocytes are initiated but not maintained. The bacteriocyte represents an evolutionarily novel cell fate, which is developmentally determined independently of the bacteria. Three of five transcription factors we examined show novel expression patterns in bacteriocytes, suggesting that bacteriocytes may have evolved to express many additional transcription factors. The evolutionary transition to a symbiosis in which bacteria and an aphid cell form a functional unit, similar to the origin of plastids, has apparently involved extensive molecular adaptations on the part of the host cell.

Comprehensive survey of developmental genes in the pea aphid, Acyrthosiphon pisum: frequent lineage‐specific duplications and losses of developmental genes

Aphids exhibit unique attributes, such as polyphenisms and specialized cells to house endosymbionts, that make them an interesting system for studies at the interface of ecology, evolution and development. Here we present a comprehensive characterization of the developmental genes in the pea aphid, Acyrthosiphon pisum, and compare our results to other sequenced insects. We investigated genes involved in fundamental developmental processes such as establishment of the body plan and organogenesis, focusing on transcription factors and components of signalling pathways. We found that most developmental genes were well conserved in the pea aphid, although many lineage‐specific gene duplications and gene losses have occurred in several gene families. In particular, genetic components of transforming growth factor beta (TGFβ) Wnt, JAK/STAT (Janus kinase/signal transducer and activator of transcription) and EGF (Epidermal Growth Factor) pathways appear to have been significantly modified in the pea aphid.

Composite Effects of Polymorphisms near Multiple Regulatory Elements Create a Major-Effect QTL

Many agriculturally, evolutionarily, and medically important characters vary in a quantitative fashion. Unfortunately, the genes and sequence variants accounting for this variation remain largely unknown due to a variety of biological and technical challenges. Drosophila melanogaster contains high levels of sequence variation and low linkage disequilibrium, allowing us to dissect the effects of many causative variants within a single locus. Here, we take advantage of these features to identify and characterize the sequence polymorphisms that comprise major effect QTL alleles segregating at the bric-a-brac locus. We show that natural bric-a-brac alleles with large effects on cuticular pigmentation reflect a cumulative impact of polymorphisms that affect three functional regions: a promoter, a tissue-specific enhancer, and a Polycomb response element. Analysis of allele-specific expression at the bric-a-brac locus confirms that these polymorphisms modulate transcription at the cis-regulatory level. Our results establish that a single QTL can act through a confluence of multiple molecular mechanisms and that sequence variation in regions flanking experimentally validated functional elements can have significant quantitative effects on transcriptional activity and phenotype. These findings have important design and conceptual implications for basic and medical genomics.

Allele-specific expression assays using Solexa

BackgroundAllele-specific expression (ASE) assays can be used to identify cis, trans, and cis-by-trans regulatory variation. Understanding the source of expression variation has important implications for disease susceptibility, phenotypic diversity, and adaptation. While ASE is commonly measured via relative fluorescence at a SNP, next generation sequencing provides an opportunity to measure ASE in an accurate and high-throughput manner using read counts.ResultsWe introduce a Solexa-based method to perform large numbers of ASE assays using only a single lane of a Solexa flowcell. In brief, transcripts of interest, which contain a known SNP, are PCR enriched and barcoded to enable multiplexing. Then high-throughput sequencing is used to estimate allele-specific expression using sequencing counts. To validate this method, we measured the allelic bias in a dilution series and found high correlations between measured and expected values (r>0.9, p < 0.001). We applied this method to a set of 5 genes in a Drosophila simulans parental mix, F1 and introgression and found that for these genes the majority of expression divergence can be explained by cis-regulatory variation.ConclusionWe present a new method with the capacity to measure ASE for large numbers of assays using as little as one lane of a Solexa flowcell. This will be a valuable technique for molecular and population genetic studies, as well as for verification of genome-wide data sets.

Accelerated Evolution of Morph-Biased Genes in Pea Aphids

Phenotypic plasticity, the production of alternative phenotypes (or morphs) from the same genotype due to environmental factors, results in some genes being expressed in a morph-biased manner. Theoretically, these morph-biased genes experience relaxed selection, the consequence of which is the buildup of slightly deleterious mutations at these genes. Over time, this is expected to result in increased protein divergence at these genes between species and a signature of relaxed purifying selection within species. Here we test these theoretical expectations using morph-biased genes in the pea aphid, a species that produces multiple morphs via polyphenism. We find that morph-biased genes exhibit faster rates of evolution (in terms of dN/dS) relative to unbiased genes and that divergence generally increases with increasing morph bias. Further, genes with expression biased toward rarer morphs (sexual females and males) show faster rates of evolution than genes expressed in the more common morph (asexual females), demonstrating that the amount of time a gene spends being expressed in a morph is associated with its rate of evolution. And finally, we show that genes expressed in the rarer morphs experience decreased purifying selection relative to unbiased genes, suggesting that it is a relaxation of purifying selection that contributes to their faster rates of evolution. Our results provide an important empirical look at the impact of phenotypic plasticity on gene evolution.

Abiotic and biotic stressors causing equivalent mortality induce highly variable transcriptional responses in the soybean aphid.

Environmental stress affects basic organismal functioning and can cause physiological, developmental, and reproductive impairment. However, in many nonmodel organisms, the core molecular stress response remains poorly characterized and the extent to which stress-induced transcriptional changes differ across qualitatively different stress types is largely unexplored. The current study examines the molecular stress response of the soybean aphid (Aphis glycines) using RNA sequencing and compares transcriptional responses to multiple stressors (heat, starvation, and plant defenses) at a standardized stress level (27% adult mortality). Stress-induced transcriptional changes showed remarkable variation, with starvation, heat, and plant defensive stress altering the expression of 3985, 510, and 12 genes, respectively. Molecular responses showed little overlap across all three stressors. However, a common transcriptional stress response was identified under heat and starvation, involved with up-regulation of glycogen biosynthesis and molecular chaperones and down-regulation of bacterial endosymbiont cellular and insect cuticular components. Stressor-specific responses indicated heat affected expression of heat shock proteins and cuticular components, whereas starvation altered a diverse set of genes involved in primary metabolism, oxidative reductive processes, nucleosome and histone assembly, and the regulation of DNA repair and replication. Exposure to host plant defenses elicited the weakest response, of which half of the genes were of unknown function. This study highlights the need for standardizing stress levels when comparing across stress types and provides a basis for understanding the role of general vs. stressor specific molecular responses in aphids.

Soybean aphid (Hemiptera: Aphididae) response to soybean plant defense: stress levels, tradeoffs, and cross-virulence.

ABSTRACT A variety of management methods to control the soybean aphid (Aphis glycines Matsumura) have been investigated since its invasion into North America in 2000, among them plant resistance has emerged as a viable option for reducing aphid damage to soybeans and preventing outbreaks. Plant resistance methods often use natural soybean plant defenses that impose stress on aphids by reducing fitness and altering behavior. Research efforts have heavily focused on identification and development of aphid resistant soybean varieties, leaving much unknown about soybean aphid response to stressful host plant defenses. In this study, we aimed to 1) evaluate lifetime fitness consequences and phenotypic variation in response to host plant-induced stress and 2) investigate whether trade-offs involving fitness costs and/or cross-virulence to multiple antibiotic soybean varieties exists. We compared aphid survival and reproduction during and after a short period of exposure to soybeans with the Rag2 resistance gene and measured aphid clonal variation in response to Rag2 soybeans. In addition, we measured the performance of Rag2 virulent and avirulent aphids on five soybean varieties with various forms of antibiotic resistance. Our results indicate that plant defenses impose high levels of stress and have long-term fitness consequences, even after aphids are removed from resistant plants. We identified one aphid clone that was able to colonize Rag2 among the seven clones tested, suggesting that virulent genotypes may be prevalent in natural populations. Finally, although we did not find evidence of cross-virulence to multiple antibiotic soybean varieties, our results suggest independent mechanisms of aphid virulence to Ragl and Rag2 that may involve fitness costs.

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.

The pea aphid uses a version of the terminal system during oviparous, but not viviparous, development

BackgroundIn most species of aphid, female nymphs develop into either sexual or asexual adults depending on the length of the photoperiod to which their mothers were exposed. The progeny of these sexual and asexual females, in turn, develop in dramatically different ways. The fertilized oocytes of sexual females begin embryogenesis after being deposited on leaves (oviparous development) while the oocytes of asexual females complete embryogenesis within the mother (viviparous development). Compared with oviparous development, viviparous development involves a smaller transient oocyte surrounded by fewer somatic epithelial cells and a smaller early embryo that comprises fewer cells. To investigate whether patterning mechanisms differ between the earliest stages of the oviparous and viviparous modes of pea aphid development, we examined the expression of pea aphid orthologs of genes known to specify embryonic termini in other insects.ResultsHere we show that pea aphid oviparous ovaries express torso-like in somatic posterior follicle cells and activate ERK MAP kinase at the posterior of the oocyte. In addition to suggesting that some posterior features of the terminal system are evolutionarily conserved, our detection of activated ERK in the oocyte, rather than in the embryo, suggests that pea aphids may transduce the terminal signal using a mechanism distinct from the one used in Drosophila. In contrast with oviparous development, the pea aphid version of the terminal system does not appear to be used during viviparous development, since we did not detect expression of torso-like in the somatic epithelial cells that surround either the oocyte or the blastoderm embryo and we did not observe restricted activated ERK in the oocyte.ConclusionsWe suggest that while oviparous oocytes and embryos may specify posterior fate through an aphid terminal system, viviparous oocytes and embryos employ a different mechanism, perhaps one that does not rely on an interaction between the oocyte and surrounding somatic cells. Together, these observations provide a striking example of a difference in the fundamental events of early development that is both environmentally induced and encoded by the same genome.

Disentangling the Causes for Faster-X Evolution in Aphids

Abstract 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.