ASU Scholarship Showcase

Under models of isolation-by-distance, population structure is determined by the probability of identity-by-descent between pairs of genes according to the geographic distance between them. Well established analytical results indicate that the relationship between geographical and genetic distance depends mostly on the neighborhood size of the population which represents a standardized measure of gene flow. To test this prediction, we model local dispersal of haploid individuals on a two-dimensional landscape using seven dispersal kernels: Rayleigh, exponential, half-normal, triangular, gamma, Lomax and Pareto. When neighborhood size is held constant, the distributions produce similar patterns of isolation-by-distance, confirming predictions. Considering this, we propose that the triangular distribution is the appropriate null distribution for isolation-by-distance studies. Under the triangular distribution, dispersal is uniform over the neighborhood area which suggests that the common description of neighborhood size as a measure of an effective, local panmictic population is valid for popular families of dispersal distributions. We further show how to draw random variables from the triangular distribution efficiently and argue that it should be utilized in other studies in which computational efficiency is important.

In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.

Mathematical models of infectious diseases are a valuable tool in understanding the mechanisms and patterns of disease transmission. It is, however, a difficult subject to teach, requiring both mathematical expertise and extensive subject-matter knowledge of a variety of disease systems. In this article, we explore several uses of zombie epidemics to make mathematical modeling and infectious disease epidemiology more accessible to public health professionals, students, and the general public. We further introduce a web-based simulation, White Zed (http://cartwrig.ht/apps/whitezed/), that can be deployed in classrooms to allow students to explore models before implementing them. In our experience, zombie epidemics are familiar, approachable, flexible, and an ideal way to introduce basic concepts of infectious disease epidemiology.

Mutation is the ultimate source of all genetic variation and is, therefore, central to evolutionary change. Previous work on Paramecium tetraurelia found an unusually low germline base-substitution mutation rate in this ciliate. Here, we tested the generality of this result among ciliates using Tetrahymena thermophila. We sequenced the genomes of 10 lines of T. thermophila that had each undergone approximately 1,000 generations of mutation accumulation (MA). We applied an existing mutation-calling pipeline and developed a new probabilistic mutation detection approach that directly models the design of an MA experiment and accommodates the noise introduced by mismapped reads. Our probabilistic mutation-calling method provides a straightforward way of estimating the number of sites at which a mutation could have been called if one was present, providing the denominator for our mutation rate calculations. From these methods, we find that T. thermophila has a germline base-substitution mutation rate of 7.61 × 10 ^-12 per-site, per cell division, which is consistent with the low base-substitution mutation rate in P. tetraurelia. Over the course of the evolution experiment, genomic exclusion lines derived from the MA lines experienced a fitness decline that cannot be accounted for by germline base-substitution mutations alone, suggesting that other genetic or epigenetic factors must be involved. Because selection can only operate to reduce mutation rates based upon the "visible" mutational load, asexual reproduction with a transcriptionally silent germline may allow ciliates to evolve extremely low germline mutation rates.

Salmonella enterica serovar Typhimurium genome encodes 13 fimbrial operons. Most of the fimbriae encoded by these operons are not produced under laboratory conditions but are likely to be synthesized in vivo. We used an in vivo expression technology (IVET) strategy to identify four fimbrial operons, agf, saf, sti, and stc that are expressed in the spleen. When any three of these operons were deleted, the strain retained wild-type virulence. However, when all four operons were deleted, the resulting strain was completely attenuated, indicating that these four fimbriae play functionally redundant roles critical for virulence. In mice, oral doses of as low as 1 × 10[superscript 5] CFU of the strain with four fimbrial operons deleted provided 100% protection against challenge with 1 × 10[superscript 9] CFU of wild-type S. Typhimurium. We also examined the possible effect of these fimbriae on the ability of a Salmonella vaccine strain to deliver a guest antigen. We modified one of our established attenuated vaccine strains, χ9088, to delete three fimbrial operons while the fourth operon was constitutively expressed. Each derivative was modified to express the Streptococcus pneumoniae antigen PspA. Strains that constitutively expressed saf or stc elicited a strong Th1 response with significantly greater levels of anti-PspA serum IgG and greater protective efficacy than strains carrying saf or stc deletions. The isogenic strain in which all four operons were deleted generated the lowest anti-PspA levels and did not protect against challenge with virulent S. pneumoniae. Our results indicate that these fimbriae play important roles, as yet not understood, in Salmonella virulence and immunogenicity.

The most common evolutionary events at the molecular level are single-base substitutions, as well as insertions and deletions (indels) of short DNA segments. A large body of research has been devoted to develop probabilistic substitution models and to infer their parameters using likelihood and Bayesian approaches. In contrast, relatively little has been done to model indel dynamics, probably due to the difficulty in writing explicit likelihood functions. Here, we contribute to the effort of modeling indel dynamics by presenting SpartaABC, an approximate Bayesian computation (ABC) approach to infer indel parameters from sequence data (either aligned or unaligned). SpartaABC circumvents the need to use an explicit likelihood function by extracting summary statistics from simulated sequences. First, summary statistics are extracted from the input sequence data. Second, SpartaABC samples indel parameters from a prior distribution and uses them to simulate sequences. Third, it computes summary statistics from the simulated sets of sequences. By computing a distance between the summary statistics extracted from the input and each simulation, SpartaABC can provide an approximation to the posterior distribution of indel parameters as well as point estimates. We study the performance of our methodology and show that it provides accurate estimates of indel parameters in simulations. We next demonstrate the utility of SpartaABC by studying the impact of alignment errors on the inference of positive selection. A C ++ program implementing SpartaABC is freely available in http://spartaabc.tau.ac.il.

Three-dimensional models of human intestinal epithelium mimic the differentiated form and function of parental tissues often not exhibited by two-dimensional monolayers and respond to Salmonella in key ways that reflect in vivo infections. To further enhance the physiological relevance of three-dimensional models to more closely approximate in vivo intestinal microenvironments encountered by Salmonella, we developed and validated a novel three-dimensional co-culture infection model of colonic epithelial cells and macrophages using the NASA Rotating Wall Vessel bioreactor. First, U937 cells were activated upon collagen-coated scaffolds. HT-29 epithelial cells were then added and the three-dimensional model was cultured in the bioreactor until optimal differentiation was reached, as assessed by immunohistochemical profiling and bead uptake assays. The new co-culture model exhibited in vivo-like structural and phenotypic characteristics, including three-dimensional architecture, apical-basolateral polarity, well-formed tight/adherens junctions, mucin, multiple epithelial cell types, and functional macrophages. Phagocytic activity of macrophages was confirmed by uptake of inert, bacteria-sized beads. Contribution of macrophages to infection was assessed by colonization studies of Salmonella pathovars with different host adaptations and disease phenotypes (Typhimurium ST19 strain SL1344 and ST313 strain D23580; Typhi Ty2). In addition, Salmonella were cultured aerobically or microaerobically, recapitulating environments encountered prior to and during intestinal infection, respectively. All Salmonella strains exhibited decreased colonization in co-culture (HT-29-U937) relative to epithelial (HT-29) models, indicating antimicrobial function of macrophages. Interestingly, D23580 exhibited enhanced replication/survival in both models following invasion. Pathovar-specific differences in colonization and intracellular co-localization patterns were observed. These findings emphasize the power of incorporating a series of related three-dimensional models within a study to identify microenvironmental factors important for regulating infection.

Salmonella enterica serovar Typhimurium strains belonging to sequence type ST313 are a major cause of fatal bacteremia among HIV-infected adults and children in sub-Saharan Africa. Unlike “classical” non-typhoidal Salmonella (NTS), gastroenteritis is often absent during ST313 infections and isolates are most commonly recovered from blood, rather than from stool. This is consistent with observations in animals, in which ST313 strains displayed lower levels of intestinal colonization and higher recovery from deeper tissues relative to classic NTS isolates. A better understanding of the key environmental factors regulating these systemic infections is urgently needed. Our previous studies using dynamic Rotating Wall Vessel (RWV) bioreactor technology demonstrated that physiological levels of fluid shear regulate virulence, gene expression, and stress response profiles of classic S. Typhimurium. Here we provide the first demonstration that fluid shear alters the virulence potential and pathogenesis-related stress responses of ST313 strain D23580 in a manner that differs from classic NTS.