ASU Scholarship Showcase

The molecular mechanisms that allow generalist parasitoids to exploit many, often very distinct hosts are practically unknown. The wasp Aphidius ervi, a generalist koinobiont parasitoid of aphids, was introduced from Europe into Chile in the late 1970s to control agriculturally important aphid species. A recent study showed significant differences in host preference and host acceptance (infectivity) depending on the host A. ervi were reared on. In contrast, no genetic differentiation between A. ervi populations parasitizing different aphid species and aphids of the same species reared on different host plants was found in Chile. Additionally, the same study did not find any fitness effects in A. ervi if offspring were reared on a different host as their mothers. Here, we determined the effect of aphid host species (Sitobion avenae versus Acyrthosiphon pisum reared on two different host plants alfalfa and pea) on the transcriptome of adult A. ervi females.

We found a large number of differentially expressed genes (between host species: head: 2,765; body: 1,216; within the same aphid host species reared on different host plants: alfalfa versus pea: head 593; body 222). As expected, the transcriptomes from parasitoids reared on the same host species (pea aphid) but originating from different host plants (pea versus alfalfa) were more similar to each other than the transcriptomes of parasitoids reared on a different aphid host and host plant (head: 648 and 1,524 transcripts; body: 566 and 428 transcripts). We found several differentially expressed odorant binding proteins and olfactory receptor proteins in particular, when we compared parasitoids from different host species. Additionally, we found differentially expressed genes involved in neuronal growth and development as well as signaling pathways.

These results point towards a significant rewiring of the transcriptome of A. ervi depending on aphid-plant complex where parasitoids develop, even if different biotypes of a certain aphid host species (A. pisum) are reared on the same host plant. This difference seems to persist even after the different wasp populations were reared on the same aphid host in the laboratory for more than 50 generations. This indicates that either the imprinting process is very persistent or there is enough genetic/allelic variation between A. ervi populations. The role of distinct molecular mechanisms is discussed in terms of the formation of host fidelity.

Community associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has become a major cause of skin and soft tissue infections (SSTIs) in the US. We developed an age-structured compartmental model to study the spread of CA-MRSA at the population level and assess the effect of control intervention strategies. We used Monte-Carlo Markov Chain (MCMC) techniques to parameterize our model using monthly time series data on SSTIs incidence in children (≤19 years) during January 2004 -December 2006 in Maricopa County, Arizona. Our model-based forecast for the period January 2007–December 2008 also provided a good fit to data. We also carried out an uncertainty and sensitivity analysis on the control reproduction number, R_c which we estimated at 1.3 (95% CI [1.2,1.4]) based on the model fit to data. Using our calibrated model, we evaluated the effect of typical intervention strategies namely reducing the contact rate of infected individuals owing to awareness of infection and decolonization strategies targeting symptomatic infected individuals on both and the long-term disease dynamics. We also evaluated the impact of hypothetical decolonization strategies targeting asymptomatic colonized individuals. We found that strategies focused on infected individuals were not capable of achieving disease control when implemented alone or in combination. In contrast, our results suggest that decolonization strategies targeting the pediatric population colonized with CA-MRSA have the potential of achieving disease elimination.

Nasonia, a genus of four closely related parasitoid insect species, is a model system for genetic research. Their haplodiploid genetics (haploid males and diploid females) and interfertile species are advantageous for the genetic analysis of complex traits and the genetic basis of species differences. A fine-scale genomic map is an important tool for advancing genetic studies in this system. We developed and used a hybrid genotyping microarray to generate a high-resolution genetic map that covers 79% of the sequenced genome of Nasonia vitripennis. The microarray is based on differential hybridization of species-specific oligos between N. vitripennis and Nasonia giraulti at more than 20,000 markers spanning the Nasonia genome. The map places 729 scaffolds onto the five linkage groups of Nasonia, including locating many smaller scaffolds that would be difficult to map by other means. The microarray was used to characterize 26 segmental introgression lines containing chromosomal regions from one species in the genetic background of another. These segmental introgression lines have been used for rapid screening and mapping of quantitative trait loci involved in species differences. Finally, the microarray is extended to bulk-segregant analysis and genotyping of other Nasonia species combinations. These resources should further expand the usefulness of Nasonia for studies of the genetic basis and architecture of complex traits and speciation.

mRNA expression dynamics promote and maintain the identity of somatic tissues in living organisms; however, their impact in post-transcriptional gene regulation in these processes is not fully understood. Here, we applied the PAT-Seq approach to systematically isolate, sequence, and map tissue-specific mRNA from five highly studied Caenorhabditis elegans somatic tissues: GABAergic and NMDA neurons, arcade and intestinal valve cells, seam cells, and hypodermal tissues, and studied their mRNA expression dynamics. The integration of these datasets with previously profiled transcriptomes of intestine, pharynx, and body muscle tissues, precisely assigns tissue-specific expression dynamics for 60% of all annotated C. elegans protein-coding genes, providing an important resource for the scientific community. The mapping of 15,956 unique high-quality tissue-specific polyA sites in all eight somatic tissues reveals extensive tissue-specific 3′untranslated region (3′UTR) isoform switching through alternative polyadenylation (APA) . Almost all ubiquitously transcribed genes use APA and harbor miRNA targets in their 3′UTRs, which are commonly lost in a tissue-specific manner, suggesting widespread usage of post-transcriptional gene regulation modulated through APA to fine tune tissue-specific protein expression. Within this pool, the human disease gene C. elegans orthologs rack-1 and tct-1 use APA to switch to shorter 3′UTR isoforms in order to evade miRNA regulation in the body muscle tissue, resulting in increased protein expression needed for proper body muscle function. Our results highlight a major positive regulatory role for APA, allowing genes to counteract miRNA regulation on a tissue-specific basis.

Biomarkers encompass a wide range of different measurable indicators, representing a tangible link to physiological changes occurring within the body. Accessibility, sensitivity, and specificity are significant factors in biomarker suitability. New biomarkers continue to be discovered, and questions over appropriate selection and assessment of their usefulness remain. If traditional markers of inflammation are not sufficiently robust in their specificity, then perhaps alternative means of detection may provide more information. Epigenetic drift (epigenetic modifications as they occur as a direct function with age), and its ancillary elements, including platelets, secreted microvesicles (MVs), and microRNA (miRNA), may hold enormous predictive potential. The majority of epigenetic drift observed in blood is independent of variations in blood cell composition, addressing concerns affecting traditional blood-based biomarker efficacy. MVs are found in plasma and other biological fluids in healthy individuals. Altered MV/miRNA profiles may also be found in individuals with various diseases. Platelets are also highly reflective of physiological and lifestyle changes, making them extremely sensitive biomarkers of human health. Platelets release increased levels of MVs in response to various stimuli and under a plethora of disease states, which demonstrate a functional effect on other cell types.

Carefully calibrated transmission models have the potential to guide public health officials on the nature and scale of the interventions required to control epidemics. In the context of the ongoing Ebola virus disease (EVD) epidemic in Liberia, Drake and colleagues, in this issue of PLOS Biology, employed an elegant modeling approach to capture the distributions of the number of secondary cases that arise in the community and health care settings in the context of changing population behaviors and increasing hospital capacity. Their findings underscore the role of increasing the rate of safe burials and the fractions of infectious individuals who seek hospitalization together with hospital capacity to achieve epidemic control. However, further modeling efforts of EVD transmission and control in West Africa should utilize the spatial-temporal patterns of spread in the region by incorporating spatial heterogeneity in the transmission process. Detailed datasets are urgently needed to characterize temporal changes in population behaviors, contact networks at different spatial scales, population mobility patterns, adherence to infection control measures in hospital settings, and hospitalization and reporting rates.

Background: Ebola virus disease (EVD) has generated a large epidemic in West Africa since December 2013. This mini-review is aimed to clarify and illustrate different theoretical concepts of infectiousness in order to compare the infectiousness across different communicable diseases including EVD.

Methods: We employed a transmission model that rests on the renewal process in order to clarify theoretical concepts on infectiousness, namely the basic reproduction number, R₀, which measures the infectiousness per generation of cases, the force of infection (i.e. the hazard rate of infection), the intrinsic growth rate (i.e. infectiousness per unit time) and the per-contact probability of infection (i.e. infectiousness per effective contact).

Results: Whereas R₀ of EVD is similar to that of influenza, the growth rate (i.e. the measure of infectiousness per unit time) for EVD was shown to be comparatively lower than that for influenza. Moreover, EVD and influenza differ in mode of transmission whereby the probability of transmission per contact is lower for EVD compared to that of influenza.

Conclusions: The slow spread of EVD associated with the need for physical contact with body fluids supports social distancing measures including contact tracing and case isolation. Descriptions and interpretations of different variables quantifying infectiousness need to be used clearly and objectively in the scientific community and for risk communication.

A central goal of biology is to uncover the genetic basis for the origin of new phenotypes. A particularly effective approach is to examine the genomic architecture of species that have secondarily lost a phenotype with respect to their close relatives. In the eusocial Hymenoptera, queens and workers have divergent phenotypes that may be produced via either expression of alternative sets of caste-specific genes and pathways or differences in expression patterns of a shared set of multifunctional genes. To distinguish between these two hypotheses, we investigated how secondary loss of the worker phenotype in workerless ant social parasites impacted genome evolution across two independent origins of social parasitism in the ant genera Pogonomyrmex and Vollenhovia. We sequenced the genomes of three social parasites and their most-closely related eusocial host species and compared gene losses in social parasites with gene expression differences between host queens and workers. Virtually all annotated genes were expressed to some degree in both castes of the host, with most shifting in queen-worker bias across developmental stages. As a result, despite >1 My of divergence from the last common ancestor that had workers, the social parasites showed strikingly little evidence of gene loss, damaging mutations, or shifts in selection regime resulting from loss of the worker caste. This suggests that regulatory changes within a multifunctional genome, rather than sequence differences, have played a predominant role in the evolution of social parasitism, and perhaps also in the many gains and losses of phenotypes in the social insects.

MicroRNAs (miRNAs) regulate gene output by targeting degenerate elements in mRNAs and have undergone drastic expansions in higher metazoan genomes. The evolutionary advantage of maintaining copies of highly similar miRNAs is not well understood, nor is it clear what unique functions, if any, miRNA family members possess. Here, we study evolutionary patterns of metazoan miRNAs, focusing on the targeting preferences of the let-7 and miR-10 families. These studies reveal hotspots for sequence evolution with implications for targeting and secondary structure. High-throughput screening for functional targets reveals that each miRNA represses sites with distinct features and regulates a large number of genes with cooperative function in regulatory networks. Unexpectedly, given the high degree of similarity, single-nucleotide changes grant miRNA family members with distinct targeting preferences. Together, our data suggest complex functional relationships among miRNA duplications, novel expression patterns, sequence change, and the acquisition of new targets.

The ongoing Zika virus (ZIKV) epidemic in the Americas poses a major global public health emergency. While ZIKV is transmitted from human to human by bites of Aedes mosquitoes, recent evidence indicates that ZIKV can also be transmitted via sexual contact with cases of sexually transmitted ZIKV reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA. Yet, the role of sexual transmission on the spread and control of ZIKV infection is not well-understood. We introduce a mathematical model to investigate the impact of mosquito-borne and sexual transmission on the spread and control of ZIKV and calibrate the model to ZIKV epidemic data from Brazil, Colombia, and El Salvador. Parameter estimates yielded a basic reproduction number R_0 = 2.055 (95% CI: 0.523–6.300), in which the percentage contribution of sexual transmission is 3.044% (95% CI: 0.123–45.73). Our sensitivity analyses indicate that R_0 is most sensitive to the biting rate and mortality rate of mosquitoes while sexual transmission increases the risk of infection and epidemic size and prolongs the outbreak. Prevention and control efforts against ZIKV should target both the mosquito-borne and sexual transmission routes.