A large genomic insertion containing a duplicated follistatin gene is linked to the pea aphid male wing dimorphism.

Wing dimorphisms have long served as models for examining the ecological and evolutionary tradeoffs associated with alternative phenotypes. Here, we investigated the genetic cause of the pea aphid (Acyrthosiphon pisum) male wing dimorphism, wherein males exhibit one of two morphologies that differ in correlated traits that include the presence or absence of wings. We mapped this trait difference to a single genomic region and, using third generation, long-read sequencing, we identified a 120 kb insertion in the wingless allele. This insertion includes a duplicated follistatin gene, which is a strong candidate gene in the minimal mapped interval to cause the dimorphism. We found that both alleles were present prior to pea aphid biotype lineage diversification, we estimated that the insertion occurred millions of years ago, and we propose that both alleles have been maintained in the species, likely due to balancing selection.

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.

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.

Moving the Pharmacy Practice Model Initiative forward within a state affiliate.

PURPOSE: Efforts to advance the ASHP Pharmacy Practice Model Initiative (PPMI) in the Michigan Society of Health-System Pharmacists (MSHP) are described. SUMMARY: After the Pharmacy Practice Model Summit in November 2010, the board of directors of MSHP began to strategize ways to help health-system pharmacists in Michigan achieve the vision and concepts envisioned by the PPMI. The ultimate goal was to develop a process for acting on recommendations developed by the PPMI to advance the practice of health-system pharmacy in Michigan. A task force was formed and reviewed the 147 national recommendations from the ASHP Pharmacy Practice Model Summit and grouped them into related areas of focus. Six focus areas were identified: acute care, ambulatory care, education and training, organizational affairs and leadership, pharmacy technicians, and technology and information systems. A PPMI Michigan conference was arranged in which focus groups would assess these six areas. Each focus group was limited to six or seven participants, with a member of the task force serving as the facilitator for the group. Individual focus groups then formulated recommendations MSHP could develop into actionable strategies to address the key issues identified during the morning session. A total of 56 recommendations were submitted by the focus groups for consideration by all conference participants. Over 80% of the recommendations were deemed to be high impact/high feasibility. CONCLUSION: A process for acting on recommendations of the ASHP PPMI to advance the practice of health-system pharmacy within the state of Michigan was developed.

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.

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

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

BACKGROUND: In 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. RESULTS: Here 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. CONCLUSIONS: We 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.

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.

BACKGROUND: Allele-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. RESULTS: We 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. CONCLUSION: We 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.

Contrasting patterns of sequence evolution at the functionally redundant bric à brac paralogs in Drosophila melanogaster.

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.

Inpatient medication history verification by pharmacy students.

PURPOSE: The ability of pharmacy students to improve the accuracy of patients' medication histories was studied. METHODS: This prospective study was conducted between January and April 2007 at a 424-bed community teaching hospital. Pharmacy students were assigned one or two patients daily admitted to the inpatient internal medicine service and were required to perform a thorough medication history for each. Patients were included in the study if a medication history was performed and recorded on the medication reconciliation form. Students were instructed to obtain medication histories by interviewing the patient, a family member, or both and calling the patient's community pharmacy to verify all medications. If there were any discrepancies between these sources of information and the initial medication reconciliation form, the information was reconfirmed with the patient. Any information obtained by the students that had not already been documented on the medication reconciliation record was updated in the patient's chart. RESULTS: A total of 326 charts were included in this analysis. Student-obtained medication histories resulted in the addition of previously undocumented prescription medications to 175 charts (53.7%) and nonprescription medications or natural products to 167 charts (51.2%). Calling the patients' community pharmacy helped to identify omissions or discrepancies approximately 75% of the time. Overall, the students improved the accuracy of medication histories for 220 (67%) of the 326 patients. CONCLUSION: Pharmacy students' participation in obtaining medication histories improved the completeness and accuracy of patient medication records.

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.

Reovirus-induced apoptosis requires mitochondrial release of Smac/DIABLO and involves reduction of cellular inhibitor of apoptosis protein levels.

Many viruses belonging to diverse viral families with differing structure and replication strategies induce apoptosis both in cultured cells in vitro and in tissues in vivo. Despite this fact, little is known about the specific cellular apoptotic pathways induced during viral infection. We have previously shown that reovirus-induced apoptosis of HEK cells is initiated by death receptor activation but requires augmentation by mitochondrial apoptotic pathways for its maximal expression. We now show that reovirus infection of HEK cells is associated with selective cytosolic release of the mitochondrial proapoptotic factors cytochrome c and Smac/DIABLO, but not the release of apoptosis-inducing factor. Release of these factors is not associated with loss of mitochondrial transmembrane potential and is blocked by overexpression of Bcl-2. Stable expression of caspase-9b, a dominant-negative form of caspase-9, blocks reovirus-induced caspase-9 activation but fails to significantly reduce activation of the key effector caspase, caspase-3. Smac/DIABLO enhances apoptosis through its action on cellular inhibitor of apoptosis proteins (IAPs). Reovirus infection is associated with selective down-regulation of cellular IAPs, including c-IAP1, XIAP, and survivin, effects that are blocked by Bcl-2 expression, establishing the dependence of IAP down-regulation on mitochondrial events. Taken together, these results are consistent with a model in which Smac/DIABLO-mediated inhibition of IAPs, rather than cytochrome c-mediated activation of caspase-9, is the key event responsible for mitochondrial augmentation of reovirus-induced apoptosis. These studies provide the first evidence for the association of Smac/DIABLO with virus-induced apoptosis.

Reovirus-induced alterations in gene expression related to cell cycle regulation.

Mammalian reovirus infection results in perturbation of host cell cycle progression. Since reovirus infection is known to activate cellular transcription factors, we investigated alterations in cell cycle-related gene expression following HEK293 cell infection by using the Affymetrix U95A microarray. Serotype 3 reovirus infection results in differential expression of 10 genes classified as encoding proteins that function at the G(1)-to-S transition, 11 genes classified as encoding proteins that function at G(2)-to-M transition, and 4 genes classified as encoding proteins that function at the mitotic spindle checkpoint. Serotype 1 reovirus infection results in differential expression of four genes classified as encoding proteins that function at the G(1)-to-S transition and three genes classified as encoding proteins that function at G(2)-to-M transition but does not alter any genes classified as encoding proteins that function at the mitotic spindle checkpoint. We have previously shown that serotype 3, but not serotype 1, reovirus infection induces a G(2)-to-M transition arrest resulting from an inhibition of cdc2 kinase activity. Of the differentially expressed genes encoding proteins regulating the G(2)-to-M transition, chk1, wee1, and GADD45 are known to inhibit cdc2 kinase activity. A hypothetical model describing serotype 3 reovirus-induced inhibition of cdc2 kinase is presented, and reovirus-induced perturbations of the G(1)-to-S, G(2)-to-M, and mitotic spindle checkpoints are discussed.