Search NASA⌕ Search

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.

At least 235 records · Page 13

Gut microbiome partially mediates and coordinates the effects of genetics on anxiety-like behavior in Collaborative Cross mice

Abstract Growing evidence suggests that the gut microbiome (GM) plays a critical role in health and disease. However, the contribution of GM to psychiatric disorders, especially anxiety, remains unclear. We used the Collaborative Cross (CC) mouse population-based model to identify anxiety associated host genetic and GM factors. Anxiety-like behavior of 445 mice across 30 CC strains was measured using the light/dark box assay and documented by video. A custom tracking system was developed to quantify seven anxiety-related phenotypes based on video. Mice were assigned to a low or high anxiety group by consensus clustering using seven anxiety-related phenotypes. Genome-wide association analysis (GWAS) identified 141 genes (264 SNPs) significantly enriched for anxiety and depression related functions. In the same CC cohort, we measured GM composition and identified five families that differ between high and low anxiety mice. Anxiety level was predicted with 79% accuracy and an AUC of 0.81. Mediation analyses revealed that the genetic contribution to anxiety was partially mediated by the GM. Our findings indicate that GM partially mediates and coordinates the effects of genetics on anxiety.

59 BASIC BIOLOGICAL SCIENCES↗

Toba volcano super eruption destroyed the ozone layer and caused a human population bottleneck

Volcanic eruptions trigger a broad spectrum of climatic responses. For example, the Mount Pinatubo eruption in 1991 forced an El Niño and global cooling, and the Tambora eruption in 1815 caused the "Year Without a Summer." Especially grand eruptions such as Toba around 74,000 years ago can push the Earth's climate into a volcanic winter state, significantly lowering the surface temperature and precipitation globally. Here we present a new, previously overlooked element of the volcanic effects spectrum: the radiative mechanism of stratospheric ozone depletion. We found that the volcanic plume of Toba enhanced the UV optical depth and suppressed the primary formation of stratospheric ozone from O2 photolysis. Sulfate aerosols additionally reflect the photons needed to break the O2 bond (λ < 242 nm), otherwise controlled by ozone absorption and Rayleigh scattering alone during volcanically quiescent conditions. Our NASA GISS ModelE simulations of the Toba eruption reveal up to 50% global ozone loss due to the overall photochemistry perturbations of the sulfate aerosols. We also consider and quantify the radiative effects of SO2, which partially compensated for the ozone loss by inhibiting the photolytic O3 sink. Our analysis shows that the magnitude of the ozone loss and UV-induced health-hazardous effects after the Toba eruption are similar to those in the aftermath of a potential nuclear conflict. These findings suggest a “Toba ozone catastrophe" as a likely contributor to the historic population decline in this period, consistent with a genetic bottleneck in human evolution.

Toba volcano↗

Dynamic nitrogen fixation in an aerobic endophyte of Populus

Biological nitrogen fixation by microbial diazotrophs can contribute significantly to nitrogen availability in non-nodulating plant species. In this study of molecular mechanisms and gene expression relating to biological nitrogen fixation, the aerobic nitrogen-fixing endophyte Burkholderia vietnamiensis, strain WPB, isolated from Populus trichocarpa served as a model for endophyte–poplar interactions. Nitrogen-fixing activity was observed to be dynamic on nitrogen-free medium with a subset of colonies growing to form robust, raised globular like structures. Secondary ion mass spectrometry (NanoSIMS) confirmed that N-fixation was uneven within the population. A fluorescent transcriptional reporter (GFP) revealed that the nitrogenase subunit nifH is not uniformly expressed across genetically identical colonies of WPB and that only ~11% of the population was actively expressing the nifH gene. Higher nifH gene expression was observed in clustered cells through monitoring individual bacterial cells using single-molecule fluorescence in situ hybridization. Through 15 N 2 enrichment, we identified key nitrogenous metabolites and proteins synthesized by WPB and employed targeted metabolomics in active and inactive populations. We cocultivated WPB Pnif-GFP with poplar within a RhizoChip, a synthetic soil habitat, which enabled direct imaging of microbial nifH expression within root epidermal cells. We observed that nifH expression is localized to the root elongation zone where the strain forms a unique physical interaction with the root cells. This work employed comprehensive experimentation to identify novel mechanisms regulating both biological nitrogen fixation and beneficial plant–endophyte interactions.

15 N-tracking metabolomics↗

Introgression and persistence of cultivar alleles in wild carrot ( Daucus carota ) populations in the United States

Abstract Premise Cultivated species and their wild relatives often hybridize in the wild, and the hybrids can survive and reproduce in some environments. However, it is unclear whether cultivar alleles are permanently incorporated into the wild genomes or whether they are purged by natural selection. This question is key to accurately assessing the risk of escape and spread of cultivar genes into wild populations. Methods We used genomic data and population genomic methods to study hybridization and introgression between cultivated and wild carrot (Daucus carota) in the United States. We used single nucleotide polymorphisms (SNPs) obtained via genotyping by sequencing for 450 wild individuals from 29 wild georeferenced populations in seven states and 144 cultivars from the United States, Europe, and Asia. Results Cultivated and wild carrot formed two genetically differentiated groups, and evidence of crop–wild admixture was detected in several but not all wild carrot populations in the United States. Two regions were identified where cultivar alleles were present in wild carrots: California and Nantucket Island (Massachusetts). Surprisingly, there was no evidence of introgression in some populations with a long‐known history of sympatry with the crop, suggesting that post‐hybridization barriers might prevent introgression in some areas. Conclusions Our results provide support for the introgression and long‐term persistence of cultivar alleles in wild carrots populations. We thus anticipate that the release of genetically engineered (GE) cultivars would lead to the introduction and spread of GE alleles in wild carrot populations.

Plant Sciences↗

Assessing electrogenetic activation via a network model of biological signal propagation

Introduction: Molecular communication is the transfer of information encoded by molecular structure and activity. We examine molecular communication within bacterial consortia as cells with diverse biosynthetic capabilities can be assembled for enhanced function. Their coordination, both in terms of engineered genetic circuits within individual cells as well as their population-scale functions, is needed to ensure robust performance. We have suggested that “electrogenetics,” the use of electronics to activate specific genetic circuits, is a means by which electronic devices can mediate molecular communication, ultimately enabling programmable control. Methods: Here, we have developed a graphical network model for dynamically assessing electronic and molecular signal propagation schemes wherein nodes represent individual cells, and their edges represent communication channels by which signaling molecules are transferred. We utilize graph properties such as edge dynamics and graph topology to interrogate the signaling dynamics of specific engineered bacterial consortia. Results: We were able to recapitulate previous experimental systems with our model. In addition, we found that networks with more distinct subpopulations (high network modularity) propagated signals more slowly than randomized networks, while strategic arrangement of subpopulations with respect to the inducer source (an electrode) can increase signal output and outperform otherwise homogeneous networks. Discussion: We developed this model to better understand our previous experimental results, but also to enable future designs wherein subpopulation composition, genetic circuits, and spatial configurations can be varied to tune performance. We suggest that this work may provide insight into the signaling which occurs in synthetically assembled systems as well as native microbial communities.

Chun, Kayla↗

Identification and Clinical Evaluation of Potential Biomarkers for Breast Cancer Resistance Protein ( BCRP / ABCG2 )

Clinical inhibition and genetic variation of the Breast Cancer Resistance Protein (BCRP/ABCG2) efflux transporter can significantly influence drug exposure, highlighting the need for reliable BCRP functional biomarkers. This study aimed to identify and evaluate biomarkers predictive of BCRP function in humans. A comprehensive analysis of metabolomic genome‐wide association studies (mGWAS) was conducted to discover potential BCRP biomarkers, followed by evaluation inin vitrotransporter assays and a clinical drug–drug interaction (DDI) study. Across multiple mGWAS datasets, plasma concentrations of three herbicide derivatives—4‐hydroxychlorothalonil (4HC), 3‐bromo‐5‐chloro‐2,6‐dihydroxybenzoic acid (BCDBA), and 3,5‐dichloro‐2,6‐dihydroxybenzoic acid (DCDBA)—were significantly elevated (P < 5E‐8) in individuals carrying reduced functionABCG2polymorphisms. These compounds were confirmed as novel BCRP substrates via transporter uptake assays and selected for clinical evaluation alongside riboflavin, a known BCRP substrate and potential BCRP biomarker. In a DDI study with 11 healthy subjects, eltrombopag, a BCRP inhibitor, increased rosuvastatin concentrations by approximately twofold (P = 0.002). No significant changes in the plasma concentrations of organic anion transporting polypeptide 1B (OATP1B) biomarkers (CP‐I and CP‐III) or potential BCRP biomarkers (4HC, BCDBA, DCDBA, or riboflavin) were observed. Notably, two subjects were heterozygous carriers for theABCG2p.Q141K variant and exhibited significantly higher baseline concentrations of 4HC (P = 0.004) and BCDBA (P = 0.0003), consistent with reduced BCRP function. These findings suggest that 4HC and BCDBA are promising biomarkers for baseline BCRP function in specific populations, such as those harboring reduced function genetic polymorphisms, but do not appear suitable for detecting acute BCRP inhibition.

Pharmacology & Pharmacy↗

Visualizing novel connections and genetic similarities across diseases using a network-medicine based approach

Understanding the genetic relationships between human disorders could lead to better treatment and prevention strategies, especially for individuals with multiple comorbidities. A common resource for studying genetic-disease relationships is the GWAS Catalog, a large and well curated repository of SNP-trait associations from various studies and populations. Some of these populations are contained within mega-biobanks such as the Million Veteran Program (MVP), which has enabled the genetic classification of several diseases in a large well-characterized and heterogeneous population. Here we aim to provide a network of the genetic relationships among diseases and to demonstrate the utility of quantifying the extent to which a given resource such as MVP has contributed to the discovery of such relations. We use a network-based approach to evaluate shared variants among thousands of traits in the GWAS Catalog repository. Our results indicate many more novel disease relationships that did not exist in early studies and demonstrate that the network can reveal clusters of diseases mechanistically related. Finally, we show novel disease connections that emerge when MVP data is included, highlighting methodology that can be used to indicate the contributions of a given biobank.

59 BASIC BIOLOGICAL SCIENCES↗

Plant microbiomes harbor potential to promote nutrient turnover in impoverished substrates of a Brazilian biodiversity hotspot

The substrates of the Brazilian campos rupestres, a grassland ecosystem, have extremely low concentrations of phosphorus and nitrogen, imposing restrictions to plant growth. Despite that, this ecosystem harbors almost 15% of the Brazilian plant diversity, raising the question of how plants acquire nutrients in such a harsh environment. Here, we set out to uncover the taxonomic profile, the compositional and functional differences and similarities, and the nutrient turnover potential of microbial communities associated with two plant species of the campos rupestres-dominant family Velloziaceae that grow over distinct substrates (soil and rock). Using amplicon sequencing data, we show that, despite the pronounced composition differentiation, the plant-associated soil and rock communities share a core of highly efficient colonizers that tend to be highly abundant and is enriched in 21 bacterial families. Functional investigation of metagenomes and 522 metagenome-assembled genomes revealed that the microorganisms found associated to plant roots are enriched in genes involved in organic compound intake, and phosphorus and nitrogen turnover. We show that potential for phosphorus transport, mineralization, and solubilization are mostly found within bacterial families of the shared microbiome, such as Xanthobacteraceae and Bryobacteraceae. We also detected the full repertoire of nitrogen cycle-related genes and discovered a lineage of Isosphaeraceae that acquired nitrogen-fixing potential via horizontal gene transfer and might be also involved in nitrification via a metabolic handoff association with Binataceae. We highlight that plant-associated microbial populations in the campos rupestres harbor a genetic repertoire with potential to increase nutrient availability and that the microbiomes of biodiversity hotspots can reveal novel mechanisms of nutrient turnover.

59 BASIC BIOLOGICAL SCIENCES↗

A large sequenced mutant library – valuable reverse genetic resource that covers 98% of sorghum genes

SUMMARY Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum , a drought and heat‐tolerant C4 species, requires a vast, large‐scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large‐scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)‐induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis ‐regulatory elements. A comparison of the distribution of mutations in the large‐scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS‐induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5′ UTR, gene body, and 3′‐UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the K a / K s value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB ( https://www.depts.ttu.edu/igcast/sorbmutdb.php ) and SorghumBase ( https://sorghumbase.org/ ). This current large‐scale sequence‐indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

59 BASIC BIOLOGICAL SCIENCES↗

Multiple Mutations Associated with Emergent Variants Can Be Detected as Low-Frequency Mutations in Early SARS-CoV-2 Pandemic Clinical Samples

Genetic analysis of intra-host viral populations provides unique insight into pre-emergent mutations that may contribute to the genotype of future variants. Clinical samples positive for SARS-CoV-2 collected in California during the first months of the pandemic were sequenced to define the dynamics of mutation emergence as the virus became established in the state. Deep sequencing of 90 nasopharyngeal samples showed that many mutations associated with the establishment of SARS-CoV-2 globally were present at varying frequencies in a majority of the samples, even those collected as the virus was first detected in the US. A subset of mutations that emerged months later in consensus sequences were detected as subconsensus members of intra-host populations. Spike mutations P681H, H655Y, and V1104L were detected prior to emergence in variant genotypes, mutations were detected at multiple positions within the furin cleavage site, and pre-emergent mutations were identified in the nucleocapsid and the envelope genes. Because many of the samples had a very high depth of coverage, a bioinformatics pipeline, “Mappgene”, was established that uses both iVar and LoFreq variant calling to enable identification of very low-frequency variants. This enabled detection of a spike protein deletion present in many samples at low frequency and associated with a variant of concern.

60 APPLIED LIFE SCIENCES↗

Genetic time traveling: sequencing old herbarium specimens, including the oldest herbarium specimen sequenced from kingdom Fungi, reveals the population structure of an agriculturally significant rust

Summary Sequencing herbarium specimens can be instrumental in answering ecological, evolutionary, and taxonomic inquiries. We developed a protocol for sequencing herbarium specimens of rust fungi (Pucciniales) and proceeded to sequence specimens ranging from 4 to 211 yr old from five different genera. We then obtained sequences from an economically important biological control agent, Puccinia suaveolens , to highlight the potential of sequencing herbarium specimens in an ecological sense and to evaluate the following hypotheses: (1) The population structure of a plant pathogen changes over time, and (2) introduced pathogens are more diverse in their native range. Our efforts resulted in sequences from 87 herbarium specimens that revealed a high level of diversity with a population structure that exhibited spatial–temporal patterns. The specimens sequenced from Europe showed more diversity than the ones from North America, uncovering an invasion pattern likely related to its European native host in North America. Additionally, to the best of our knowledge, the specimen from France collected in c. 1811 is the oldest herbarium specimen sequenced from kingdom Fungi. In conclusion, sequencing old herbarium specimens is an important tool that can be extrapolated to better understand plant–microbe evolution and to evaluate old type specimens to solidify the taxonomy of plant pathogenic fungi.

Plant Sciences↗

Genetic models in applied physiology: selected contribution: effects of spaceflight on immunity in the C57BL/6 mouse. I. Immune population distributions

There are several aspects of the spaceflight environment that may lead to changes in immunity: mission-related psychological stress, radiation, and changes in gravity. On December 5, 2001, the space shuttle Endeavor launched for a 12-day mission to examine these effects on C57BL/6 mice for the first time. On their return, assays were performed on the spleen, blood, and bone marrow. In response to flight, there were no significant differences in the general circulating leukocyte proportions. In contrast, there was an increase in splenic lymphocyte percentages, with a corresponding decrease in granulocytes. There was an overall shift in splenic lymphocytes away from T cells toward B cells, and a decrease in the CD4-to-CD8 ratios due to a decrease in T helpers. In contrast, there were proportional increases in bone marrow T cells, with decreases in B cells. Although the blast percentage and count were decreased in flight mice, the CD34(+) population was increased. The data were more consistent with a shift in bone marrow populations rather than a response to changes in the periphery. Many of the results are similar to those using other models. Clearly, spaceflight can influence immune parameters ranging from hematopoiesis to mature leukocyte mechanisms.

Non-NASA Center↗

Large-scale survey of prion protein genetic variability in scrapie disease-free goats from the United States

Scrapie is a slowly progressive neurodegenerative disease of small ruminants caused by an accumulation of an abnormal isoform of prion protein in the central nervous system. Polymorphisms of the prion protein gene (PRNP) strongly modulate scrapie resistance and incubation period in goats. The aim of this study was to identify PRNP genetic variability in goats across the United States. Blood from a total of 6,029 apparent scrapie disease-free goats from 654 operations and 19 breeds were analyzed. Sequencing of PRNP revealed 26 genotypes with different rates based on eight codons. The GG127, RR154, and QQ222 genotypes were predominant and showed a remarkably high rate across all goats. The QK222 and NS146 genotypes, known to be protective against scrapie, were found in 0.6% [with 95% CI = (0.3, 1.2)] and 22.0% [95% CI = (19.1, 25.2)] of goats, respectively. The QK222 genotype was found in 23.1% of Oberhasli goats tested, with 95%CI = (3.9, 68.7)] and 22.0% of Toggenburg goats tested with 95%CI = (9.7, 42.5)], while NS146 was found in 65.5% of Savannah goats tested, with 95%CI = (30.8, 89.9), 36.7% of Boer goats tested, with 95%CI = (33.1, 40.4), 36.3% of Nubian goats tested, with 95%CI = (27.0, 46.7)], and 35.6% of LaMancha goats tested, with 95%CI = (22.8, 50.8%). The MM142 and IM142 genotypes were found more frequently in goats on dairy operations, while the HR143, NS146, and ND146 genotypes were found more frequently in goats on meat operations. Goats in the east region had a higher percentage of goats with RH154, RQ211, and QK222 genotypes than goats in the west region. The results of this study showed high genetic variability of PRNP among the U.S. goat population, with differences by location and breed, and may serve as a rationale for development of goat breeding programs at the national level to mitigate the risk of scrapie.

59 BASIC BIOLOGICAL SCIENCES↗

Genetics and Genomics of Pathogen Resistance in Switchgrass (Final Report)

This project was funded by DOE under Grant no. DE-SC0016108. Originally approved for the 2016-2019 period, two no-cost extensions were solicited and approved, which prolonged the lifespan through July 2021. This final report informs on the results obtained so far from the research implemented. The research hinged on integrating genomics (genomic selection, RNAseq, virus-plant interactions) with classical genetics (conventional breeding) to incorporate durable resistance to fungal (rust) and viral (mosaic) diseases in switchgrass (Panicum virgatum) populations being bred for bioenergy. Higher biomass yield, higher quality (low lignin content), and durable disease resistance are key features to make lignocellulosic switchgrass feedstocks economically competitive and sustainable. Genomic selection is being applied on three generations of a switchgrass population derived from crossing two ecotypes (Kanlow as lowland female and Summer as upland male) with differential performance in terms of biomass yield and quality, disease resistance, and winter survivability. Target populations were screened for rust and mosaic in field and/or lab and phenotyped for biomass yield and quality traits. Genetic analyses were applied across generations to capture the joint inheritance of the targeted traits and predict breeding values for parents and progeny with greater accuracy. Parental and a panel of different switchgrass populations were genotyped with the DArTseq technology to develop SNP (0, 1, 2) and in-silico (presence/absence) DArT markers. Rust inoculations techniques were developed and applied successfully on switchgrass. The original populations (Kanlow and Summer) were sequenced with RNAseq to capture the gene expression profiles across sequential time-points and appraise the basis of greater resistance in the Kanlow vs the Summer ecotype. Constructs of PMV and sPMV mosaic virus were assembled and tested first on proso millet to find the best protocol to use later on switchgrass. Results from the preliminary analyses indicate that 1) ample additive genetic variation is available for selection and improving this inter-ecotypic population for yield, quality, and disease traits, 2) significant gains are to be expected with the genetic correlations being favorable between yield and lignin content and between yield and disease ratings, 3) substantial differences exist in the genetic regions controlling rust resistance in the two ecotypes, 4) co-infection with PMV isolates from Nebraska and its satellite from Kansas elicit severe mosaic symptoms, and 5) two different genetic systems are responsible for imparting resistance to rust and virus in switchgrass.

59 BASIC BIOLOGICAL SCIENCES↗

Association mapping across a multitude of traits collected in diverse environments in maize

Abstract Classical genetic studies have identified many cases of pleiotropy where mutations in individual genes alter many different phenotypes. Quantitative genetic studies of natural genetic variants frequently examine one or a few traits, limiting their potential to identify pleiotropic effects of natural genetic variants. Widely adopted community association panels have been employed by plant genetics communities to study the genetic basis of naturally occurring phenotypic variation in a wide range of traits. High-density genetic marker data—18M markers—from 2 partially overlapping maize association panels comprising 1,014 unique genotypes grown in field trials across at least 7 US states and scored for 162 distinct trait data sets enabled the identification of of 2,154 suggestive marker-trait associations and 697 confident associations in the maize genome using a resampling-based genome-wide association strategy. The precision of individual marker-trait associations was estimated to be 3 genes based on a reference set of genes with known phenotypes. Examples were observed of both genetic loci associated with variation in diverse traits (e.g., above-ground and below-ground traits), as well as individual loci associated with the same or similar traits across diverse environments. Many significant signals are located near genes whose functions were previously entirely unknown or estimated purely via functional data on homologs. This study demonstrates the potential of mining community association panel data using new higher-density genetic marker sets combined with resampling-based genome-wide association tests to develop testable hypotheses about gene functions, identify potential pleiotropic effects of natural genetic variants, and study genotype-by-environment interaction.

59 BASIC BIOLOGICAL SCIENCES↗

Limits on the evolutionary rates of biological traits

Abstract This paper focuses on the maximum speed at which biological evolution can occur. I derive inequalities that limit the rate of evolutionary processes driven by natural selection, mutations, or genetic drift. These rate limits link the variability in a population to evolutionary rates. In particular, high variances in the fitness of a population and of a quantitative trait allow for fast changes in the trait’s average. In contrast, low variability makes a trait less susceptible to random changes due to genetic drift. The results in this article generalize Fisher’s fundamental theorem of natural selection to dynamics that allow for mutations and genetic drift, via trade-off relations that constrain the evolutionary rates of arbitrary traits. The rate limits can be used to probe questions in various evolutionary biology and ecology settings. They apply, for instance, to trait dynamics within or across species or to the evolution of bacteria strains. They apply to any quantitative trait, e.g., from species’ weights to the lengths of DNA strands.

59 BASIC BIOLOGICAL SCIENCES↗

Limits on the Evolutionary Rates of Biological Traits

This paper focuses on the maximum speed at which biological evolution can occur. I derive inequalities that limit the rate of evolutionary processes driven by natural selection, mutations, or genetic drift. These rate limits link the variability in a population to evolutionary rates. In particular, high variances in the fitness of a population and of a quantitative trait allow for fast changes in the trait’s average. In contrast, low variability makes a trait less susceptible to random changes due to genetic drift. The results in this article generalize Fisher’s fundamental theorem of natural selection to dynamics that allow for mutations and genetic drift, via trade-of relations that constrain the evolutionary rates of arbitrary traits. The rate limits can be used to probe questions in various evolutionary biology and ecology settings. They apply, for instance, to trait dynamics within or across species or to the evolution of bacteria strains. They apply to any quantitative trait, e.g., from species’ weights to the lengths of DNA strands.

59 BASIC BIOLOGICAL SCIENCES↗

Automatic Molecular Design using Evolutionary Techniques

Molecular nanotechnology is the precise, three-dimensional control of materials and devices at the atomic scale. An important part of nanotechnology is the design of molecules for specific purposes. This paper describes early results using genetic software techniques to automatically design molecules under the control of a fitness function. The fitness function must be capable of determining which of two arbitrary molecules is better for a specific task. The software begins by generating a population of random molecules. The population is then evolved towards greater fitness by randomly combining parts of the better individuals to create new molecules. These new molecules then replace some of the worst molecules in the population. The unique aspect of our approach is that we apply genetic crossover to molecules represented by graphs, i.e., sets of atoms and the bonds that connect them. We present evidence suggesting that crossover alone, operating on graphs, can evolve any possible molecule given an appropriate fitness function and a population containing both rings and chains. Prior work evolved strings or trees that were subsequently processed to generate molecular graphs. In principle, genetic graph software should be able to evolve other graph representable systems such as circuits, transportation networks, metabolic pathways, computer networks, etc.

Globus, Al↗