Faculty and Staff

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Inbreeding in hermaphroditic plants can occur through two different mechanisms: biparental inbreeding, when a plant mates with a related individual, or self-fertilization, when a plant mates with itself. To avoid inbreeding, many hermaphroditic plants have evolved self-incompatibility (SI) systems which prevent or limit self-fertilization. One particular SI system—homomorphic SI—can also reduce biparental inbreeding. Homomorphic SI is found in many angiosperm species, and it is often assumed that the additional benefit of reduced biparental inbreeding may be a factor in the success of this SI system. To test this assumption, we developed a spatially-explicit, individual-based simulation of plant populations that displayed three different types of homomorphic SI. We measured the total level of inbreeding avoidance by comparing each population to a self-compatible population (NSI), and we measured biparental inbreeding avoidance by comparing to a population of self-incompatible plants that were free to mate with any other individual (PSI).

Because biparental inbreeding is more common when offspring dispersal is limited, we examined the levels of biparental inbreeding over a range of dispersal distances. We also tested whether the introduction of inbreeding depression affected the level of biparental inbreeding avoidance. We found that there was a statistically significant decrease in autozygosity in each of the homomorphic SI populations compared to the PSI population and, as expected, this was more pronounced when seed and pollen dispersal was limited. However, levels of homozygosity and inbreeding depression were not reduced. At low dispersal, homomorphic SI populations also suffered reduced female fecundity and had smaller census population sizes. Overall, our simulations showed that the homomorphic SI systems had little impact on the amount of biparental inbreeding in the population especially when compared to the overall reduction in inbreeding compared to the NSI population. With further study, this observation may have important consequences for research into the origin and evolution of homomorphic self-incompatibility systems.

Background: Blindness has evolved repeatedly in cave-dwelling organisms, and many hypotheses have been proposed to explain this observation, including both accumulation of neutral loss-of-function mutations and adaptation to darkness. Investigating the loss of sight in cave dwellers presents an opportunity to understand the operation of fundamental evolutionary processes, including drift, selection, mutation, and migration.

Results: Here we model the evolution of blindness in caves. This model captures the interaction of three forces: (1) selection favoring alleles causing blindness, (2) immigration of sightedness alleles from a surface population, and (3) mutations creating blindness alleles. We investigated the dynamics of this model and determined selection-strength thresholds that result in blindness evolving in caves despite immigration of sightedness alleles from the surface. We estimate that the selection coefficient for blindness would need to be at least 0.005 (and maybe as high as 0.5) for blindness to evolve in the model cave-organism, Astyanax mexicanus.

Conclusions: Our results indicate that strong selection is required for the evolution of blindness in cave-dwelling organisms, which is consistent with recent work suggesting a high metabolic cost of eye development.

The Dawn Framing Camera (FC) has imaged the northern hemisphere of the Asteroid (4) Vesta at high spatial resolution and coverage. This study represents the first investigation of the overall geology of the northern hemisphere (22–90°N, quadrangles Av-1, 2, 3, 4 and 5) using these unique Dawn mission observations. We have compiled a morphologic map and performed crater size–frequency distribution (CSFD) measurements to date the geologic units. The hemisphere is characterized by a heavily cratered surface with a few highly subdued basins up to ∼200 km in diameter. The most widespread unit is a plateau (cratered highland unit), similar to, although of lower elevation than the equatorial Vestalia Terra plateau. Large-scale troughs and ridges have regionally affected the surface. Between ∼180°E and ∼270°E, these tectonic features are well developed and related to the south pole Veneneia impact (Saturnalia Fossae trough unit), elsewhere on the hemisphere they are rare and subdued (Saturnalia Fossae cratered unit). In these pre-Rheasilvia units we observed an unexpectedly high frequency of impact craters up to ∼10 km in diameter, whose formation could in part be related to the Rheasilvia basin-forming event. The Rheasilvia impact has potentially affected the northern hemisphere also with S–N small-scale lineations, but without covering it with an ejecta blanket. Post-Rheasilvia impact craters are small (<60 km in diameter) and show a wide range of degradation states due to impact gardening and mass wasting processes. Where fresh, they display an ejecta blanket, bright rays and slope movements on walls. In places, crater rims have dark material ejecta and some crater floors are covered by ponded material interpreted as impact melt.

Oppia Quadrangle Av-10 (288–360°E, ±22°) is a junction of key geologic features that preserve a rough history of Asteroid (4) Vesta and serves as a case study of using geologic mapping to define a relative geologic timescale. Clear filter images, stereo-derived topography, slope maps, and multispectral color-ratio images from the Framing Camera on NASA’s Dawn spacecraft served as basemaps to create a geologic map and investigate the spatial and temporal relationships of the local stratigraphy. Geologic mapping reveals the oldest map unit within Av-10 is the cratered highlands terrain which possibly represents original crustal material on Vesta that was then excavated by one or more impacts to form the basin Feralia Planitia. Saturnalia Fossae and Divalia Fossae ridge and trough terrains intersect the wall of Feralia Planitia indicating that this impact basin is older than both the Veneneia and Rheasilvia impact structures, representing Pre-Veneneian crustal material. Two of the youngest geologic features in Av-10 are Lepida (∼45 km diameter) and Oppia (∼40 km diameter) impact craters that formed on the northern and southern wall of Feralia Planitia and each cross-cuts a trough terrain. The ejecta blanket of Oppia is mapped as ‘dark mantle’ material because it appears dark orange in the Framing Camera ‘Clementine-type’ color-ratio image and has a diffuse, gradational contact distributed to the south across the rim of Rheasilvia. Mapping of surface material that appears light orange in color in the Framing Camera ‘Clementine-type’ color-ratio image as ‘light mantle material’ supports previous interpretations of an impact ejecta origin. Some light mantle deposits are easily traced to nearby source craters, but other deposits may represent distal ejecta deposits (emplaced >5 crater radii away) in a microgravity environment.

We used Dawn spacecraft data to identify and delineate geological units and landforms in the Marcia quadrangle of Vesta as a means to assess the role of the large, relatively young impact craters Marcia (∼63 km diam.) and Calpurnia (∼53 km diam.) and their surrounding ejecta field on the local geology. We also investigated a local topographic high with a dark-rayed crater named Aricia Tholus, and the impact crater Octavia that is surrounded by a distinctive diffuse mantle. Crater counts and stratigraphic relations suggest that Marcia is the youngest large crater on Vesta, in which a putative impact melt on the crater floor ranges in age between ∼40 and 60 Ma (depending upon choice of chronology system), and Marcia’s ejecta blanket ranges in age between ∼120 and 390 Ma (depending upon choice of chronology system).

We interpret the geologic units in and around Marcia crater to mark a major vestan time-stratigraphic event, and that the Marcia Formation is one of the geologically youngest formations on Vesta. Marcia crater reveals pristine bright and dark material in its walls and smooth and pitted terrains on its floor. The smooth unit we interpret as evidence of flow of impact melts and (for the pitted terrain) release of volatiles during or after the impact process. The distinctive dark ejecta surrounding craters Marcia and Calpurnia is enriched in OH- or H-bearing phases and has a variable morphology, suggestive of a complex mixture of impact ejecta and impact melts including dark materials possibly derived from carbonaceous chondrite-rich material. Aricia Tholus, which was originally interpreted as a putative vestan volcanic edifice based on lower resolution observations, appears to be a fragment of an ancient impact basin rim topped by a dark-rayed impact crater. Octavia crater has a cratering model formation age of ∼280–990 Ma based on counts of its ejecta field (depending upon choice of chronology system), and its ejecta field is the second oldest unit in this quadrangle. The relatively young craters and their related ejecta materials in this quadrangle are in stark contrast to the surrounding heavily cratered units that are related to the billion years old or older Rheasilvia and Veneneia impact basins and Vesta’s ancient crust preserved on Vestalia Terra.

We studied high-resolution images of asteroid Vesta's surface (~70 and 20–25 m/pixel) obtained during the High- and Low-Altitude Mapping Orbits (HAMO, LAMO) of NASA's Dawn mission to assess the formation mechanisms responsible for a variety of lobate, flow-like features observed across the surface. We searched for evidence of volcanic flows, based on prior mathematical modeling and the well-known basaltic nature of Vesta's crust, but no unequivocal morphologic evidence of ancient volcanic activity has thus far been identified. Rather, we find that all lobate, flow-like features on Vesta appear to be related either to impact or erosional processes. Morphologically distinct lobate features occur in and around impact craters, and most of these are interpreted as impact ejecta flows, or possibly flows of impact melt. Estimates of melt production from numerical models and scaling laws suggests that large craters like Marcia (~60 km diameter) could have potentially produced impact melt volumes ranging from tens of millions of cubic meters to a few tens of cubic kilometers, which are relatively small volumes compared to similar-sized lunar craters, but which are consistent with putative impact melt features observed in Dawn images. There are also examples of lobate flows that trend downhill both inside and outside of crater rims and basin scarps, which are interpreted as the result of gravity-driven mass movements (slumps and landslides).

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.