Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species
Not provided.
Biochemistry & Molecular Biology↗
SEARCH · Search NASA
Results for “population size”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Getting more from heterogeneous HIV-1 surveillance data in a high immigration country: estimation of incidence and undiagnosed population size using multiple biomarkers
Most HIV infections originate from individuals who are undiagnosed and unaware of their infection. Estimation of this quantity from surveillance data is hard because there is incomplete knowledge about (i) the time between infection and diagnosis (TI) for the general population, and (ii) the time between immigration and diagnosis for foreign-born persons. We developed a new statistical method for estimating the incidence of HIV-1 and the number of undiagnosed people living with HIV (PLHIV), based on dynamic modelling of heterogeneous HIV-1 surveillance data. The methods consist of a Bayesian non-linear mixed effects model using multiple biomarkers to estimate TI of HIV-1-positive individuals, and a novel incidence estimator which distinguishes between endogenous and exogenous infections by modelling explicitly the probability that a foreign-born person was infected either before or after immigration. The incidence estimator allows for direct calculation of the number of undiagnosed persons. The new methodology is illustrated combining heterogeneous surveillance data from Sweden between 2003 and 2015. A leave-one-out cross-validation study showed that the multiple-biomarker model was more accurate than single biomarkers (mean absolute error 1.01 vs ≥1.95). We estimate that 816 [95% credible interval (CI) 775-865] PLHIV were undiagnosed in 2015, representing a proportion of 10.8% (95% CI 10.3-11.4%) of all PLHIV. Here, the proposed methodology will enhance the utility of standard surveillance data streams and will be useful to monitor progress towards and compliance with the 90–90-90 UNAIDS target.
60 APPLIED LIFE SCIENCES↗
Impact of the HIV-1 genetic background and HIV-1 population size on the evolution of raltegravir resistance
Background Emergence of resistance against integrase inhibitor raltegravir in human immunodeficiency virus type 1 (HIV-1) patients is typically associated with selection of one of three signature mutations: Y143C/R, Q148K/H/R or N155H, representing three distinct resistance pathways. The mechanisms that drive selection of a specific pathway are still poorly understood. We reviewed the effect of the HIV-1 genetic background and population dynamics on the emergence of raltegravir resistance. Using deep sequencing we analyzed the integrase coding sequence (CDS) in longitudinal samples from five patients who initiated raltegravir plus optimized background therapy at viral loads > 5000 copies/ml. To investigate the role of the HIV-1 genetic background we created recombinant viruses containing the viral integrase coding region from pre-raltegravir samples from two patients in whom raltegravir resistance developed through different pathways. The in vitro selections performed with these recombinant viruses were designed to mimic natural population bottlenecks. Results Deep sequencing analysis of the viral integrase CDS revealed that the virological response to raltegravir containing therapy inversely correlated with the relative amount of unique sequence variants that emerged suggesting diversifying selection during drug pressure. In 4/5 patients multiple signature mutations representing different resistance pathways were observed. Interestingly, the resistant population can consist of a single resistant variant that completely dominates the population but also of multiple variants from different resistance pathways that coexist in the viral population. We also found evidence for increased diversification after stronger bottlenecks. In vitro selections with low viral titers, mimicking population bottlenecks, revealed that both recombinant viruses and HXB2 reference virus were able to select mutations from different resistance pathways, although typically only one resistance pathway emerged in each individual culture. Conclusions The generation of a specific raltegravir resistant variant is not predisposed in the genetic background of the viral integrase CDS. Typically, in the early phases of therapy failure the sequence space is explored and multiple resistance pathways emerge and then compete for dominance which frequently results in a switch of the dominant population over time towards the fittest variant or even multiple variants of similar fitness that can coexist in the viral population.
60 APPLIED LIFE SCIENCES↗
Multi-Objective Hyperparameter Optimization for Spiking Neural Network Neuroevolution
Neuroevolution has had significant success over recent years, but there has been relatively little work applying neuroevolution approaches to spiking neural networks (SNNs). SNNs are a type of neural network that includes temporal processing component, are not easily trained using other methods, and can be deployed into energy-efficient neuromorphic hardware. In this work, we investigate two evolutionary approaches for training SNNs. We explore the impact of the hyperparameters of the evolutionary approaches, including tournament size, population size, and representation type, on the performance of the algorithms. We present a multi-objective Bayesian-based hyperparameter optimization approach to tune the hyperparameters to produce the most accurate and smallest SNNs. We show that the hyperparameters can significantly affect the performance of these algorithms. We also perform sensitivity analysis and demonstrate that every hyperparameter value has the potential to perform well, assuming other hyperparameter values are set correctly.
Parsa, Maryam↗
Do Genetic Drift and Gene Flow Affect the Geographic Distribution of Female Plants in Gynodioecious Lobelia siphilitica?
Variation in population sex ratio is particularly pronounced in gynodioecious angiosperms. Extremely high female frequencies in gynodioecious populations cannot be readily explained by selective forces alone. To assess the contributions of drift and gene flow to extreme sex-ratio variation, we documented sex ratio and population size in 92 populations of Lobelia siphilitica across its range and genotyped plants using plastid and nuclear genetic markers. Similarity in spatial patterns of genetic and demographic variables may suggest that drift and/or gene flow have contributed to population sex-ratio variation in L. siphilitica. We found strong spatial structuring of extremely high female frequencies: populations with >50% female plants are restricted to the south–central portion of the range. However, we did not detect any spatial structuring in population size nor metrics of genetic diversity, suggesting that extreme variation in female frequency is not strongly affected by drift or gene flow. Extreme sex-ratio variation is frequently observed in gynodioecious plants, but its causes are difficult to identify. Further investigation into mechanisms that create or maintain the spatial structure of sex ratios in gynodioecious species will provide much needed insight.
54 ENVIRONMENTAL SCIENCES↗
Investigating the Particle Growth in Bimodal Pt/C Catalysts by In-Situ Small-Angle X-ray Scattering: Challenges in the Evaluation of Stress Test Protocol-Dependent Degradation Mechanisms
The influence of different combinations of accelerated stress test (AST) protocols simulating load-cycle and start/stop conditions of a proton exchange membrane fuel cell (PEMFC) vehicle is investigated on a bimodal Pt/C catalyst. The bimodal Pt/C catalyst, prepared by mixing two commercial catalysts, serves as a model system and consists of two distinguishable size populations. The change in mean particle size was investigated by in situ small-angle X-ray scattering (SAXS). The comparison to the reference catalysts, i.e., the two single-size population catalysts, uncovers the presence of electrochemical Ostwald ripening as a degradation mechanism in the bimodal catalyst. Increasing the harshness of the applied AST protocol combinations by faster changing between load-cycle or start/stop conditions, the particle size of the larger population of the bimodal catalyst increases faster than expected. Surprisingly, the change in mean particle size of the smaller size population indicates a smaller increase for harsher AST protocols, which might be explained by a substantial electrochemical Ostwald ripening.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Dynamic regimes of local homogeneous population model with time lag
We investigated Moran - Ricker model with time lag 1. It is made analytical and numerical study of the model. It is shown there is co-existence of various dynamic regimes under the same values of parameters. The model simultaneously possesses several different limit regimes: stable state, periodic fluctuations, and chaotic attractor. The research results show if present population size substantially depends on population number of previous year then it is observed quasi-periodic oscillations. Fluctuations with period 2 occur when the growth of population size is regulated by density dependence in the current year.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Footprints of Worldwide Adaptation in Structured Populations of Drosophila melanogaster Through the Expanded DEST 2.0 Genomic Resource
Abstract Large-scale genomic resources can place genetic variation into an ecologically informed context. To advance our understanding of the population genetics of the fruit fly Drosophila melanogaster, we present an expanded release of the community-generated population genomics resource Drosophila Evolution over Space and Time (DEST 2.0; https://dest.bio/). This release includes 530 high-quality pooled libraries from flies collected across six continents over more than a decade (2009 to 2021), most at multiple time points per year; 211 of these libraries are sequenced and shared here for the first time. We used this enhanced resource to elucidate several aspects of the species' demographic history and identify novel signs of adaptation across spatial and temporal dimensions. For example, we showed that the spatial genetic structure of populations is stable over time, but that drift due to seasonal contractions of population size causes populations to diverge over time. We identified signals of adaptation that vary between continents in genomic regions associated with xenobiotic resistance, consistent with independent adaptation to common pesticides. Moreover, by analyzing samples collected during spring and fall across Europe, we provide new evidence for seasonal adaptation related to loci associated with pathogen response. Furthermore, we have also released an updated version of the DEST genome browser. This is a useful tool for studying spatiotemporal patterns of genetic variation in this classic model system.
Biochemistry & Molecular Biology↗
The Impact of Variability in Patient Exposure During Premarket Clinical Development on Postmarket Safety Outcomes
We characterized the size of the premarket safety population for 278 small-molecule new molecular entities (NMEs) and 61 new therapeutic biologics (NTBs) approved by the US Food and Drug Administration (FDA) between October 1, 2002, and December 31, 2014, evaluating the relationship of premarket safety population size to regulatory characteristics and postmarket safety outcomes. The median size of the safety population was 1,044, and was lower for NTBs than NMEs (median: 920 vs. 1,138, P = 0.04), orphan products than nonorphan products (393 vs. 1,606, P < 0.001), and for products with fast-track designation (617 vs. 1,455, P < 0.001), priority review (630 vs. 1,735, P < 0.001), and accelerated approval (475 vs. 1,164, P < 0.001), than products without that designation. The median number of postmarket safety label updates and issues added to the label were higher with larger premarket exposure among nonorphan products, but not among orphan products. Products with accelerated approval using a surrogate end point had a higher median number of safety issues added to the label than those with full approval, but this did not vary with the size of the safety population; fast-track and priority review were not associated with the number of safety issues added to the label. Here, a smaller safety population size was associated with a longer time to first safety outcome for nonorphan products but not orphan products. For orphan and nonorphan products combined, smaller premarket safety population size is not associated with the number or timing of postmarket safety outcomes, regardless of expedited program participation.
60 APPLIED LIFE SCIENCES↗
COMPILE: a GWAS computational pipeline for gene discovery in complex genomes
Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Acoustic and Genetic Data Can Reduce Uncertainty Regarding Populations of Migratory Tree-Roosting Bats Impacted by Wind Energy
Wind turbine-related mortality may pose a population-level threat for migratory tree-roosting bats, such as the hoary bat (Lasiurus cinereus) in North America. These species are dispersed within their range, making it impractical to estimate census populations size using traditional survey methods. Nonetheless, understanding population size and trends is essential for evaluating and mitigating risk from wind turbine mortality. Using various sampling techniques, including systematic acoustic sampling and genetic analyses, we argue that building a weight of evidence regarding bat population status and trends is possible to (1) assess the sustainability of mortality associated with wind turbines; (2) determine the level of mitigation required; and (3) evaluate the effectiveness of mitigation measures to ensure population viability for these species. Long-term, systematic data collection remains the most viable option for reducing uncertainty regarding population trends for migratory tree-roosting bats. We recommend collecting acoustic data using the statistically robust North American Bat Monitoring Program (NABat) protocols and that genetic diversity is monitored at repeated time intervals to show species trends. There are no short-term actions to resolve these population-level questions; however, we discuss opportunities for relatively short-term investments that will lead to long-term success in reducing uncertainty.
17 WIND ENERGY↗
Asynchronous demographic responses to Pleistocene climate change in Eastern Nearctic vertebrates
Abstract Pleistocene climatic cycles altered species distributions in the Eastern Nearctic of North America, yet the degree of congruent demographic response to the Pleistocene among codistributed taxa remains unknown. We use a hierarchical approximate Bayesian computational approach to test if population sizes across lineages of snakes, lizards, turtles, mammals, birds, salamanders and frogs in this region expanded synchronously to Late Pleistocene climate changes. Expansion occurred in 75% of 74 lineages, and of these, population size trajectories across the community were partially synchronous, with coexpansion found in at least 50% of lineages in each taxonomic group. For those taxa expanding outside of these synchronous pulses, factors related to when they entered the community, ecological thresholds or biotic interactions likely condition their timing of response to Pleistocene climate change. Unified timing of population size change across communities in response to Pleistocene climate cycles is likely rare in North America.
Burbrink, Frank T.↗
Ringed seal demography in a changing climate
Abstract Climate change is affecting species’ distributions and abundances worldwide. Baseline population estimates, against which future observations may be compared, are necessary if we are to detect ecological change. Arctic sea ice ecosystems are changing rapidly and we lack baseline population estimates for many ice‐associated species. Provided we can detect them, changes in Arctic marine ecosystems may be signaled by changes in indicator species such as ringed seals ( Pusa hispida ). Ringed seal monitoring has provided estimates of survival and fertility rates, but these have not been used for population‐level inference. Using matrix population models, we synthesized existing demographic parameters to obtain estimates of historical ringed seal population growth and structure in Amundsen Gulf and Prince Albert Sound, Canada. We then formalized existing hypotheses about the effects of emerging environmental stressors (i.e., earlier spring ice breakup and reduced snow depth) on ringed seal pup survival. Coupling the demographic model to ice and snow forecasts available from the Coupled Model Intercomparison Project resulted in projections of ringed seal population size and structure up to the year 2100. These projections showed median declines in population size ranging from 50% to 99%. Corresponding to these projected declines were substantial changes in population structure, with increasing proportions of ringed seal pups and adults and declining proportions of juveniles. We explored if currently collected, harvest‐based data could be used to detect the projected changes in population stage structure. Our model suggests that at a present sample size of 100 seals per year, the projected changes in stage structure would only be reliably detected by mid‐century, even for the most extreme climate models. This modeling process revealed inconsistencies in existing estimates of ringed seal demographic rates. Mathematical population models such as these can contribute both to understanding past population trends as well as predicting future ones, both of which are necessary if we are to detect and interpret future observations.
Reimer, Jody R.↗
Scaling methods for accelerating kinetic Monte Carlo simulations of chemical reaction networks
Various kinetic Monte Carlo algorithms become inefficient when some of the population sizes in a system are large, which gives rise to a large number of reaction events per unit time. In this work, we present a new acceleration algorithm based on adaptive and heterogeneous scaling of reaction rates and stoichiometric coefficients. The algorithm is conceptually related to the commonly used idea of accelerating a stochastic simulation by considering a subvolume λΩ (0 < λ < 1) within a system of interest, which reduces the number of reaction events per unit time occurring in a simulation by a factor 1/λ at the cost of greater error in unbiased estimates of first moments and biased overestimates of second moments. Our new approach offers two unique benefits. First, scaling is adaptive and heterogeneous, which eliminates the pitfall of overaggressive scaling. Second, there is no need for an a priori classification of populations as discrete or continuous (as in a hybrid method), which is problematic when discreteness of a chemical species changes during a simulation. The method requires specification of only a single algorithmic parameter, N c , a global critical population size above which populations are effectively scaled down to increase simulation efficiency. The method, which we term partial scaling, is implemented in the open-source BioNetGen software package. We demonstrate that partial scaling can significantly accelerate simulations without significant loss of accuracy for several published models of biological systems. Finally, these models characterize activation of the mitogen-activated protein kinase ERK, prion protein aggregation, and T-cell receptor signaling.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The waiting time problem in a model hominin population
Functional information is normally communicated using specific, context-dependent strings of symbolic characters. This is true within the human realm (texts and computer programs), and also within the biological realm (nucleic acids and proteins). In biology, strings of nucleotides encode much of the information within living cells. How do such information-bearing nucleotide strings arise and become established? This paper uses comprehensive numerical simulation to understand what types of nucleotide strings can realistically be established via the mutation/selection process, given a reasonable timeframe. The program Mendel’s Accountant realistically simulates the mutation/selection process, and was modified so that a starting string of nucleotides could be specified, and a corresponding target string of nucleotides could be specified. We simulated a classic pre-human hominin population of at least 10,000 individuals, with a generation time of 20 years, and with very strong selection (50 % selective elimination). Random point mutations were generated within the starting string. Whenever an instance of the target string arose, all individuals carrying the target string were assigned a specified reproductive advantage. When natural selection had successfully amplified an instance of the target string to the point of fixation, the experiment was halted, and the waiting time statistics were tabulated. Using this methodology we tested the effect of mutation rate, string length, fitness benefit, and population size on waiting time to fixation. Biologically realistic numerical simulations revealed that a population of this type required inordinately long waiting times to establish even the shortest nucleotide strings. To establish a string of two nucleotides required on average 84 million years. To establish a string of five nucleotides required on average 2 billion years. We found that waiting times were reduced by higher mutation rates, stronger fitness benefits, and larger population sizes. However, even using the most generous feasible parameters settings, the waiting time required to establish any specific nucleotide string within this type of population was consistently prohibitive. We show that the waiting time problem is a significant constraint on the macroevolution of the classic hominin population. Routine establishment of specific beneficial strings of two or more nucleotides becomes very problematic.
54 ENVIRONMENTAL SCIENCES↗
Surveying Grassland Islands: the genetics and performance of Appalachian switchgrass (Panicum virgatum L.) collections
The interior Southeastern United States could contain novel germplasm for the bioenergy crop switchgrass due to its diverse habitats and geographic location between genetic subpopulations (Atlantic, Midwest, and Gulf). Collections from this region could accelerate breeding progress, contribute to conservation efforts, and improve understanding of isolated grasslands in the region. This study located 22 sites in the Midsouth region and obtained 1,521,210 single nucleotide polymorphism markers of 202 individuals through genotype-by-sequencing. Here, individuals were evaluated for flowering time, winter survival and tiller number. Comparison to a national diversity panel revealed that branches of two major subpopulations occur in the region with two levels of polyploidy: Atlantic tetraploids and Midwest octoploids. Two locations contained admixed octoploid individuals with Midwest and Gulf genetics. Field performance of the Midwest octoploids conformed with prior reported performance of the Midwest subpopulation, although three sites contained promising late flowering traits. The Atlantic tetraploids had moderate winter survival, short stature, and anomalously early flowering. Atlantic populations mostly occurred in marginal sites and their morphological and flowering time adaptations may be a resource conservation strategy. Demographic inference of historical effective population size variation in a subset of tetraploid locations indicated a widespread recent decline in effective population size. This pattern is consistent with isolation of these switchgrass communities from larger populations and is further supported by evidence of inbreeding within the populations (F I = 0.18). The populations documented in this study contain novel genetic diversity and adaptations to a range of marginal habitats. Therefore, this study provides a new source of germplasm for future breeding and conservation programs.
59 BASIC BIOLOGICAL SCIENCES↗
Invasion and rapid adaptation of guppies ( Poecilia reticulata ) across the Hawaiian Archipelago
How much does natural selection, as opposed to genetic drift, admixture, and gene flow, contribute to the evolution of invasive species following introduction to a new environment? Here we assess how evolution can shape biological invasions by examining population genomic variation in non-native guppies (Poecilia reticulata) introduced to the Hawaiian Islands approximately a century ago. By examining 18 invasive populations from four Hawaiian islands and four populations from the native range in northern South America, we reconstructed the history of introductions and evaluated population structure as well as the extent of ongoing gene flow across watersheds and among islands. Patterns of differentiation indicate that guppies have developed significant population structure, with little natural or human-mediated gene flow having occurred among populations following introduction. Demographic modeling and admixture graph analyses together suggest that guppies were initially introduced to O‘ahu and Maui and then translocated to Hawai‘i and Kaua‘i. We detected evidence for only one introduction event from the native range, implying that any adaptive evolution in introduced populations likely utilized the genetic variation present in the founding population. Environmental association tests accounting for population structure identified loci exhibiting signatures of adaptive variation related to predators and landscape characteristics but not nutrient regimes. When paired with high estimates of effective population sizes and detectable population structure, the presence of environment-associated loci supports the role of natural selection in shaping contemporary evolution of Hawaiian guppy populations. Our findings indicate that local adaptation may engender invasion success, particularly in species with life histories that facilitate rapid evolution. Finally, evidence of low gene flow between populations suggests that removal could be an effective approach to control invasive guppies across the Hawaiian archipelago.
59 BASIC BIOLOGICAL SCIENCES↗
Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta
Genomic variation is widespread, and both neutral and selective processes can generate similar patterns in the genome. These processes are not mutually exclusive, so it is difficult to infer the evolutionary mechanisms that govern population and species divergence. Boechera stricta is a perennial relative of Arabidopsis thaliana native to largely undisturbed habitats with two geographic and ecologically divergent subspecies. Here, we delineate the evolutionary processes driving the genetic diversity and population differentiation in this species. Using whole-genome re-sequencing data from 517 B. stricta accessions, we identify four genetic groups that diverged around 30-180 thousand years ago, with long-term small effective population sizes and recent population expansion after the Last Glacial Maximum. We find three genomic regions with elevated nucleotide diversity, totaling about 10% of the genome. These three regions of elevated nucleotide diversity show excess of intermediate-frequency alleles, higher absolute divergence ( d XY ), and lower relative divergence ( F ST ) than genomic background, and significant enrichment in immune-related genes, reflecting long-term balancing selection. Scattered across the genome, we also find regions with both high F ST and d XY among the groups, termed F ST -islands. Population genetic signatures indicate that F ST -islands with elevated divergence, which have experienced directional selection, are derived from divergent sorting of ancient polymorphisms. Our results suggest that long-term balancing selection on disease resistance genes may have maintained ancestral haplotypes across different geographical lineages, and unequal sorting of balanced polymorphisms may have generated genomic regions with elevated divergence. This study highlights the importance of ancestral balanced polymorphisms as crucial components of genome-wide variation.
59 BASIC BIOLOGICAL SCIENCES↗