SEARCH · Search NASA
Results for “population genetics”
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.
Population Genetics of Sugar Kelp Throughout the Northeastern United States Using Genome-Wide Markers
An assessment of genetic diversity of marine populations is critical not only for the understanding and preserving natural biodiversity but also for its commercial potential. As commercial demand rises for marine resources, it is critical to generate baseline information for monitoring wild populations. Furthermore, anthropogenic stressors on the coastal environment, such as warming sea temperatures and overharvesting of wild populations, are leading to the destruction of keystone marine species such as kelps. In this study, we conducted a fine-scale genetic analysis using genome-wide high-density markers on Northwest Atlantic sugar kelp. The population structure for a total of 149 samples from the Gulf of Maine (GOM) and Southern New England (SNE) was investigated using AMOVA, FST, admixture, and PCoA. Genome-wide association analyses were conducted for six morphological traits, and the extended Lewontin and Krakauer (FLK) test was used to detect selection signatures. Our results indicate that the GOM region is more heterogeneous than SNE. These two regions have large genetic difference (between-location FST ranged from 0.21 to 0.32) and were separated by Cape Cod, which is known to be the biogeographic barrier for other taxa. We detected one significant SNP (P = 2.03 × 10 -7 ) associated with stipe length, and 248 SNPs with higher-than-neutral differentiation. The findings of this study provide baseline knowledge on sugar kelp population genetics for future monitoring, managing and potentially restoring wild populations, as well as assisting in selective breeding to improve desirable traits for future commercialization opportunities.
Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach
Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.
Social network structure and the spread of complex contagions from a population genetics perspective
Ideas, behaviors, and opinions spread through social networks. If the probability of spreading to a new individual is a non-linear function of the fraction of the individuals’ affected neighbors, such a spreading process becomes a “complex contagion”. This non-linearity does not typically appear with physically spreading infections, but instead can emerge when the concept that is spreading is subject to game theoretical considerations (e.g. for choices of strategy or behavior) or psychological effects such as social reinforcement and other forms of peer influence (e.g. for ideas, preferences, or opinions). Here we study how the stochastic dynamics of such complex contagions are affected by the underlying network structure. Motivated by simulations of complex contagions on real social networks, we present a framework for analyzing the statistics of contagions with arbitrary non-linear adoption probabilities based on the mathematical tools of population genetics. The central idea is to use an effective lower-dimensional diffusion process to approximate the statistics of the contagion. This leads to a tradeoff between the effects of ”selection” (microscopic tendencies for an idea to spread or die out), random drift, and network structure. Our framework illustrates intuitively several key properties of complex contagions: stronger community structure and network sparsity can significantly enhance the spread, while broad degree distributions dampen the effect of selection compared to random drift. Finally, we show that some structural features can exhibit critical values that demarcate regimes where global contagions become possible for networks of arbitrary size. Our results draw parallels between the competition of genes in a population and memes in a world of minds and ideas. Our tools provide insight into the spread of information, behaviors, and ideas via social influence, and highlight the role of macroscopic network structure in determining their fate.
Reply to: Population genetic considerations regarding the interpretation of within-patient SARS-CoV-2 polymorphism data
Not Available
An integrated model of population genetics and community ecology
Explore the source record for details and available documents.
CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations
CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of sgRNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Lastly, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.
Exploring genetic diversity, population structure, and subgenome differences in the allopolyploid Camelina sativa : implications for future breeding and research studies
Abstract Camelina (Camelina sativa), an allohexaploid species, is an emerging aviation biofuel crop that has been the focus of resurgent interest in recent decades. To guide future breeding and crop improvement efforts, the community requires a deeper comprehension of subgenome dominance, often noted in allopolyploid species, “alongside an understanding of the genetic diversity” and population structure of material present within breeding programs. We conducted population genetic analyses of a C. sativa diversity panel, leveraging a new genome, to estimate nucleotide diversity and population structure, and analyzed for patterns of subgenome expression dominance among different organs. Our analyses confirm that C. sativa has relatively low genetic diversity and show that the SG3 subgenome has substantially lower genetic diversity compared to the other two subgenomes. Despite the low genetic diversity, our analyses identified 13 distinct subpopulations including two distinct wild populations and others putatively representing founders in existing breeding populations. When analyzing for subgenome composition of long non-coding RNAs, which are known to play important roles in (a)biotic stress tolerance, we found that the SG3 subgenome contained significantly more lincRNAs compared to other subgenomes. Similarly, transcriptome analyses revealed that expression dominance of SG3 is not as strong as previously reported and may not be universal across all organ types. From a global analysis, SG3 “was only significant higher expressed” in flower, flower bud, and fruit organs, which is an important discovery given that the crop yield is associated with these organs. Collectively, these results will be valuable for guiding future breeding efforts in camelina.
Increasing temperature weakens the positive effect of genetic diversity on population growth
Abstract Genetic diversity and temperature increases associated with global climate change are known to independently influence population growth and extinction risk. Whether increasing temperature may influence the effect of genetic diversity on population growth, however, is not known. We address this issue in the model protist system Tetrahymena thermophila . We test the hypothesis that at temperatures closer to the species’ thermal optimum (i.e., the temperature at which population growth is maximal, or T opt ), genetic diversity should have a weaker effect on population growth compared to temperatures away from the thermal optimum. To do so, we grew populations of T . thermophila with varying levels of genetic diversity at increasingly warmer temperatures and quantified their intrinsic population growth rate, r . We found that genetic diversity increases population growth at cooler temperatures, but that as temperature increases, this effect weakens. We also show that a combination of changes in the amount of expressed genetic diversity (G) in r , plastic changes in population growth across temperatures (E), and strong G × E interactions underlie this temperature effect. Our results uncover important but largely overlooked temperature effects that have implications for the management of small populations with depauperate genetic stocks in an increasingly warming world.
Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes
Abstract Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.
BLOOD-BASED MULTI-SCALE MODEL FOR CANCER RISK FROM GCR IN GENETICALLY DIVERSE POPULATIONS
OBJECTIVES AND METHODS This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2). QUANTIFICATION OF 53BP1+ FOCI IN VITRO AND ASSOCIATIONS TO IN VIVO RADIATION SUSCEPTIBILITY IN 15 MOUSE STRAINS We reported in vitro repair kinetic and repairable fractions of RIF for the 15 mouse strains and introduced a mathematical model for RIF as a function of time, dose and LET. We noted that the metabolic activity of cells modulates the RIF response, and we introduced the open access tool terRIFic (Tool for Enhanced Results of RIF In Cells, https://radbiolab.shinyapps.io/terrific/) to correct for such bias using confluence level. Notably, at 4h post-irradiation, RIF/Gy decreased with dose or LET: as the dose or LET increases, so does the proximity of DNA double-strand-breaks (DSB) and our data suggest that proximal DSBs are brought together inside isolated RIF for repair. The RIF/Gy trend was inverted at 24h, suggesting RIF with high DSB content are more difficult to repair. We showed that in vitro metrics correlate with in vivo measurements in the same 15 mouse strains, such as survival levels of immune cells or spontaneous cancer incidence, suggesting a relationship between the efficiency of DSB repair and cancer risk or radiation toxicity. In addition to the efficiency of repair and persistent RIF, the amount of spontaneous foci before irradiation was also found to be strain dependent. Finally, we performed genome-wide association study in the same 15 mouse strains using all RIF phenotypes measured in vitro, identifying genes of interest and validating RIF as an ideal biomarker for individual radiation sensitivity. BASELINE 53BP1+ FOCI PREDICTS INDIVIDUAL HUMAN RESPONSE TO GCR COMPONENTS Based on the analysis of radiation responses of 576 donor PBMCs (using quantification of 53BP1+ foci, oxidative stress and cell death), we observed a wide variability of subject- and LET-dependent radiation responses, with radiation-induced DNA repair foci increasing with LET, though oxidative stress being notably reduced by high-LET irradiation, potentially due to a switch between hydrogen peroxide and oxygen radical-based mechanisms. We identified a relationship between few spontaneous DNA foci at baseline and increased DNA repair after irradiation, accompanied by an alteration in immunoregulatory cytokine secretion, which might be adapted as biomarkers to predict ionizing radiation sensitivity. Among demographic variables, only latent cytomegalovirus infection and age were predictive of high baseline foci formation. Finally, we have performed low-throughput whole genome sequencing of all samples and are currently in the process of identifying the genes and pathways associated with low and high-LET ionizing radiation sensitivity in humans.
An Open-Science Approach to Address Individual Response to Simulated GCR In Genetically Diverse Populations of Mice and Humans
This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm). A genome-wide association study (GWAS) was performed on the mouse strains between DNA damage responses to space radiation and single nucleotide polymorphisms (SNPs). We found SNPs, which were significantly associated to the RIF phenotype, mapped to genes and pathways that are functionally linked to health hazards for deep space exploration (e.g. carcinogenesis, nervous system damage and immune dysfunction). Some of these SNPs were located within protein coding regions, potentially interfering with protein functions and providing promising genetic targets for countermeasures. We also found correlations between both spontaneous and radiation-induced DNA damage and SNPs mapped to pathways associated with cellular metabolism. GWAS is undergoing for the human data. All data have been made available via the NASA Space Biology Open-Science database (genelab.nasa.gov) and we will discuss how various genomic and transcriptomic datasets can be accessed for modeling and integrated using machine learning methods for discovering new radiation biology.
Genetic diversity, population structure and anthracnose resistance response in a novel sweet sorghum diversity panel
Sweet sorghum is an attractive feedstock for the production of renewable chemicals and fuels due to the readily available fermentable sugars that can be extracted from the juice, and the additional stream of fermentable sugars that can be obtained from the cell wall polysaccharides in the bagasse. An important selection criterion for new sweet sorghum germplasm is resistance to anthracnose, a disease caused by the fungal pathogen Colletotrichum sublineolum. The identification of novel anthracnose-resistance sources present in sweet sorghum germplasm offers a fast track towards the development of new resistant sweet sorghum germplasm. We established a sweet sorghum diversity panel (SWDP) of 272 accessions from the USDA-ARS National Plant Germplasm (NPGS) collection that includes landraces from 22 countries and advanced breeding material, and that represents ~15% of the NPGS sweet sorghum collection. Genomic characterization of the SWDP identified 171,954 single nucleotide polymorphisms (SNPs) with an average of one SNP per 4,071 kb. Population structure analysis revealed that the SWDP could be stratified into four populations and one admixed group, and that this population structure could be aligned to sorghum’s racial classification. Results from a two-year replicated trial of the SWDP for anthracnose resistance response in Texas, Georgia, Florida, and Puerto Rico showed 27 accessions to be resistant across locations, while 145 accessions showed variable resistance response against local pathotypes. A genome-wide association study identified 16 novel genomic regions associated with anthracnose resistance. Four resistance loci on chromosomes 3, 6, 8 and 9 were identified against pathotypes from Puerto Rico, and two resistance loci on chromosomes 3 and 8 against pathotypes from Texas. In Georgia and Florida, three resistance loci were detected on chromosomes 4, 5, 6 and four on chromosomes 4, 5 (two loci) and 7, respectively. One resistance locus on chromosome 2 was effective against pathotypes from Texas and Puerto Rico and a genomic region of 41.6 kb at the tip of chromosome 8 was associated with resistance response observed in Georgia, Texas, and Puerto Rico. This publicly available SWDP and the extensive evaluation of anthracnose resistance represent a valuable genomic resource for the improvement of sorghum.
KBase Narrative - Supplemental Data for "Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes"
Here we provide the complete set of circular elements produced by SCAPP from metagenomic data that may represent mobile genetic elements (MGEs). For the publication, these sequences were de-replicated within the high [U] and low [U] sample sets. Additionally, 98 of the circular elements were removed for downstream analysis due to high suspicion of being (1) chloroplast or mitochondrial sequences or (2) erroneously circularized long repeat regions. The complete list of circular elements used for the analysis along with metadata can be found in Tables S5 and S6.
Skin-Based DNA Repair Phenotype for Cancer Risk from GCR in Genetically Diverse Populations
Predicting cancer risk associated with cosmic radiation remains a mission-critical challenge for NASA radiation health scientists and mission planners. Epidemiological data are lacking and risk methods do not take individual radiation sensitivity into account. In our approach we hypothesize that genetic factors strongly influence risk of cancer from space radiation and that biomarkers reflecting DNA damage and cell death are ideal tools to predict risk and monitor potential health effects post-flight. At this workshop, we will be reporting the work we have done over the first 9 months of this proposal. Skin cells from 15 different strains of mice already characterized for radiation-induced cancer sensitivity (B6C3F; BALB/cByJ, C57BL/6J, CBA/CaJ, C3H/HeMsNrsf), and 10 strains from the DOE collaborative cross-mouse model were expanded from ear biopsy and cultivated until Passage 3. On average, 3 males and 3 females for each strain were expanded and frozen for further characterization at the NSRL beam line during the NSRL16C run for three LET (350 MeV/n Si, 350 MeV/n Ar and 600 MeV/n Fe) and two ion fluences (1 and 3 particles per cell). The mice work has established new metrics for the usage of Radiation Induced Foci as a marker for various aspect of DNA repair deficiencies. In year 2, we propose to continue characterization of the mouse lines with low LET to identify loci specific to high- versus low- LET and establish genetic linkage for the various DNA repair biomarkers. Correlation with cancer risk from each animals strain and gender will also be investigated. On the human side, we will start characterizing the DNA damage response induced ex-vivo in 200 human's blood donors for radiation sensitivity with a tentative 500 donors by the end of this project. All ex-vivo phenotypic data will be correlated to genetic characterization of each individual human donors using SNP arrays characterization as done for mice. Similarly, ex-vivo phenotypic features from mice will be associated to cancer risk, to identify which biomarkers correlate the most with cancer risk. Genetic traits across humans will also be associated to radiation phenotypic features as a function of age and gender.
PyPop: a mature open-source software pipeline for population genomics
Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations.
Exploring the Genetic Basis of Wild Boar ( Sus scrofa ) and Its Connection to Classical Swine Fever Spread
Classical swine fever (CSF) is the one of the most devastating contagious diseases in domestic swine and wild boar/pigs (Sus scrofa). Population genetics is often used to estimate animal dispersal and can also help evaluate host population connectivity, which is crucial for understanding pathogen dispersal. We surveyed genetic population structure of boars using MIG-seq analysis to clarify the geographic barriers that influence boar dispersal in north-central Japan and to demonstrate the relationship between the spread of CSF infection among boars and their population structure. We obtained 382 single-nucleotide polymorphisms from 348 wild boar samples, and the results of STRUCTURE analysis indicated that the highest ΔK value was at K = 2, followed by K = 4. Based on these results, it is evident that the Abukuma river, a major river within north-central Japan, does not act as a barrier to the gene flow of boars, but rather that human infrastructure hinders their dispersal. Further, according to the time series change in the capture site of CSF-infected wild boar and the sum of the probability of belonging to each of the four clades in individual CSF-infected wild boar, our results indicated that the genetic structure of boar populations was correlated with the outbreak pathway of CSF across our study region. Our study suggests that predictions of disease spread, especially for widely distributed host species, is challenging because of the risk of cryptic breaks and changes in wide range connectivity; however, understanding the genetic population structure of wild boar can be a useful tool for predicting the spread of CSF. We concluded that genetic analysis of host population structure may have the possibility to improve predictions of the future dynamics of disease spread.
Uncovering novel sources of anthracnose resistance in populations of genetically diverse sorghums Sorghum bicolor (L.) Moench
Anthracnose (Colletotrichum sublineolum P. Henn., in Kabat and Bubák) is the major disease of sorghum [ Sorghum bicolor (L.) Moench ] in the Southeastern United States and the Caribbean, where the acreage of sorghum has been expanding rapidly in the past decade. The causal pathogen of anthracnose is a hemibiotrophic fungus: after an initial invasion of live plant cells, the fungus kills the plant tissue and obtains its nutrients from the necrotic tissue. Anthracnose resistance is the most effective and sustainable solution to ensuring sorghum biomass yield and biomass quality. During the past few years, multiple research efforts have been directed to identify new sources of resistance in the Sorghum Association Panel (SAP) and tropical germplasm, and to delimit genomic regions associated with the observed anthracnose resistant response. Nevertheless, these resistance loci only explained a limited portion of the observed variation indicating the existence of other resistance loci. This research project was directed to uncover anthracnose resistance loci present in the SAP and to validate previously identified resistance loci via introgression into susceptible sweet sorghums. In addition, U.S. National Plant Germplasm System (NPGS) tropical sorghum germplasm was evaluated to identify new sources of anthracnose resistance.