ASU Scholarship Showcase

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.

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.

Megacities are major sources of anthropogenic fossil fuel CO₂(FFCO₂) emissions. The spatial extents of these large urban systems cover areas of 10 000 km² or more with complex topography and changing landscapes. We present a high-resolution land–atmosphere modelling system for urban CO₂ emissions over the Los Angeles (LA) megacity area. The Weather Research and Forecasting (WRF)-Chem model was coupled to a very high-resolution FFCO₂ emission product, Hestia-LA, to simulate atmospheric CO₂ concentrations across the LA megacity at spatial resolutions as fine as  ∼  1 km. We evaluated multiple WRF configurations, selecting one that minimized errors in wind speed, wind direction, and boundary layer height as evaluated by its performance against meteorological data collected during the CalNex-LA campaign (May–June 2010). Our results show no significant difference between moderate-resolution (4 km) and high-resolution (1.3 km) simulations when evaluated against surface meteorological data, but the high-resolution configurations better resolved planetary boundary layer heights and vertical gradients in the horizontal mean winds. We coupled our WRF configuration with the Vulcan 2.2 (10 km resolution) and Hestia-LA (1.3 km resolution) fossil fuel CO₂ emission products to evaluate the impact of the spatial resolution of the CO₂ emission products and the meteorological transport model on the representation of spatiotemporal variability in simulated atmospheric CO₂ concentrations. We find that high spatial resolution in the fossil fuel CO₂ emissions is more important than in the atmospheric model to capture CO₂ concentration variability across the LA megacity. Finally, we present a novel approach that employs simultaneous correlations of the simulated atmospheric CO₂ fields to qualitatively evaluate the greenhouse gas measurement network over the LA megacity. Spatial correlations in the atmospheric CO₂ fields reflect the coverage of individual measurement sites when a statistically significant number of sites observe emissions from a specific source or location. We conclude that elevated atmospheric CO₂ concentrations over the LA megacity are composed of multiple fine-scale plumes rather than a single homogenous urban dome. Furthermore, we conclude that FFCO₂ emissions monitoring in the LA megacity requires FFCO₂ emissions modelling with  ∼1 km resolution because coarser-resolution emissions modelling tends to overestimate the observational constraints on the emissions estimates.

We present a high-resolution atmospheric inversion system combining a Lagrangian Particle Dispersion Model (LPDM) and the Weather Research and Forecasting model (WRF), and test the impact of assimilating meteorological observation on transport accuracy. A Four Dimensional Data Assimilation (FDDA) technique continuously assimilates meteorological observations from various observing systems into the transport modeling system, and is coupled to the high resolution CO₂ emission product Hestia to simulate the atmospheric mole fractions of CO₂. For the Indianapolis Flux Experiment (INFLUX) project, we evaluated the impact of assimilating different meteorological observation systems on the linearized adjoint solutions and the CO₂ inverse fluxes estimated using observed CO₂ mole fractions from 11 out of 12 communications towers over Indianapolis for the Sep.-Nov. 2013 period. While assimilating WMO surface measurements improved the simulated wind speed and direction, their impact on the planetary boundary layer (PBL) was limited. Simulated PBL wind statistics improved significantly when assimilating upper-air observations from the commercial airline program Aircraft Communications Addressing and Reporting System (ACARS) and continuous ground-based Doppler lidar wind observations. Wind direction mean absolute error (MAE) decreased from 26 to 14 degrees and the wind speed MAE decreased from 2.0 to 1.2 m s_-1, while the bias remains small in all configurations (< 6 degrees and 0.2 m s_-1). Wind speed MAE and ME are larger in daytime than in nighttime. PBL depth MAE is reduced by ~10%, with little bias reduction. The inverse results indicate that the spatial distribution of CO₂ inverse fluxes were affected by the model performance while the overall flux estimates changed little across WRF simulations when aggregated over the entire domain. Our results show that PBL wind observations are a potent tool for increasing the precision of urban meteorological reanalyses, but that the impact on inverse flux estimates is dependent on the specific urban environment.

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.

The INFLUX experiment has taken multiple approaches to estimate the carbon dioxide (CO₂) flux in a domain centered on the city of Indianapolis, Indiana. One approach, Hestia, uses a bottom-up technique relying on a mixture of activity data, fuel statistics, direct flux measurement and modeling algorithms. A second uses a Bayesian atmospheric inverse approach constrained by atmospheric CO₂ measurements and the Hestia emissions estimate as a prior CO₂ flux. The difference in the central estimate of the two approaches comes to 0.94 MtC (an 18.7% difference) over the eight-month period between September 1, 2012 and April 30, 2013, a statistically significant difference at the 2-sigma level. Here we explore possible explanations for this apparent discrepancy in an attempt to reconcile the flux estimates. We focus on two broad categories: 1) biases in the largest of bottom-up flux contributions and 2) missing CO₂ sources. Though there is some evidence for small biases in the Hestia fossil fuel carbon dioxide (FFCO2) flux estimate as an explanation for the calculated difference, we find more support for missing CO₂ fluxes, with biological respiration the largest of these. Incorporation of these differences bring the Hestia bottom-up and the INFLUX inversion flux estimates into statistical agreement and are additionally consistent with wintertime measurements of atmospheric ¹⁴CO₂. We conclude that comparison of bottom-up and top-down approaches must consider all flux contributions and highlight the important contribution to urban carbon budgets of animal and biotic respiration. Incorporation of missing CO₂ fluxes reconciles the bottom-up and inverse-based approach in the INFLUX domain.

We assess the detectability of city emissions via a tower-based greenhouse gas (GHG) network, as part of the Indianapolis Flux (INFLUX) experiment. By examining afternoon-averaged results from a network of carbon dioxide (CO₂), methane (CH₄), and carbon monoxide (CO) mole fraction measurements in Indianapolis, Indiana for 2011–2013, we quantify spatial and temporal patterns in urban atmospheric GHG dry mole fractions. The platform for these measurements is twelve communications towers spread across the metropolitan region, ranging in height from 39 to 136 m above ground level, and instrumented with cavity ring-down spectrometers. Nine of the sites were deployed as of January 2013 and data from these sites are the focus of this paper. A background site, chosen such that it is on the predominantly upwind side of the city, is utilized to quantify enhancements caused by urban emissions. Afternoon averaged mole fractions are studied because this is the time of day during which the height of the boundary layer is most steady in time and the area that influences the tower measurements is likely to be largest. Additionally, atmospheric transport models have better performance in simulating the daytime convective boundary layer compared to the nighttime boundary layer. Averaged from January through April of 2013, the mean urban dormant-season enhancements range from 0.3 ppm CO₂ at the site 24 km typically downwind of the edge of the city (Site 09) to 1.4 ppm at the site at the downwind edge of the city (Site 02) to 2.9 ppm at the downtown site (Site 03). When the wind is aligned such that the sites are downwind of the urban area, the enhancements are increased, to 1.6 ppm at Site 09, and 3.3 ppm at Site 02. Differences in sampling height affect the reported urban enhancement by up to 50%, but the overall spatial pattern remains similar. The time interval over which the afternoon data are averaged alters the calculated urban enhancement by an average of 0.4 ppm. The CO₂ observations are compared to CO₂ mole fractions simulated using a mesoscale atmospheric model and an emissions inventory for Indianapolis. The observed and modeled CO₂ enhancements are highly correlated (r² = 0.94), but the modeled enhancements prior to inversion average 53% of those measured at the towers. Following the inversion, the enhancements follow the observations closely, as expected. The CH₄ urban enhancement ranges from 5 ppb at the site 10 km predominantly downwind of the city (Site 13) to 21 ppb at the site near the landfill (Site 10), and for CO ranges from 6 ppb at the site 24 km downwind of the edge of the city (Site 09) to 29 ppb at the downtown site (Site 03). Overall, these observations show that a dense network of urban GHG measurements yield a detectable urban signal, well-suited as input to an urban inversion system given appropriate attention to sampling time, sampling altitude and quantification of background conditions.

The ‘Hestia Project’ uses a bottom-up approach to quantify fossil fuel CO₂(FFCO₂) emissions spatially at the building/street level and temporally at the hourly level. Hestia FFCO₂ emissions are provided in the form of a group of sector-specific vector layers with point, line, and polygon sources to support carbon cycle science and climate policy. Application to carbon cycle science, in particular, requires regular gridded data in order to link surface carbon fluxes to atmospheric transport models. However, the heterogeneity and complexity of FFCO₂ sources within regular grids is sensitive to spatial resolution. From the perspective of a data provider, we need to find a balance between resolution and data volume so that the gridded data product retains the maximum amount of information content while maintaining an efficient data volume.

The Shannon entropy determines the minimum bits that are needed to encode an information source and can serve as a metric for the effective information content. In this paper, we present an analysis of the Shannon entropy of gridded FFCO₂ emissions with varying resolutions in four Hestia study areas, and find: (1) the Shannon entropy increases with smaller grid resolution until it reaches a maximum value (the max-entropy resolution); (2) total emissions (the sum of several sector-specific emission fields) show a finer max-entropy resolution than each of the sector-specific fields; (3) the residential emissions show a finer max-entropy resolution than the commercial emissions; (4) the max-entropy resolution of the onroad emissions grid is closely correlated to the density of the road network. These findings suggest that the Shannon entropy can detect the information effectiveness of the spatial resolution of gridded FFCO₂ emissions. Hence, the resolution-entropy relationship can be used to assist in determining an appropriate spatial resolution for urban CO₂ flux studies. We conclude that the optimal spatial resolution for providing Hestia total FFCO₂ emissions products is centered around 100 m, at which the FFCO₂ emissions data can not only fully meet the requirement of urban flux integration, but also be effectively used in understanding the relationships between FFCO₂ emissions and various social-economic variables at the U.S. census block group level.

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.

Atmospheric CO₂ inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO₂ measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr^-1 (±0.5 Pg C yr^-1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr^-1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr^-1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr[superscript −1] for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr^-1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr^-1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr^-1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.