SEARCH · Search NASA
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.
Here, we utilized metagenomic and size-filtered virome sequencing to recover 4,157 medium, high, or complete quality viral genomes from soils taken from three high elevation sites in New Mexico, USA. Among recovered viral genomes, 90% were from size-filtered samples, indicating the importance of this enrichment in assessments of complex viromes.
Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth–resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into theAtabeyarchaeum deiterrae-1chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.
The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.
Pangenome analyses are increasingly being utilized to study the evolution of eukaryotic organisms. While pangenomes can provide insight into polymorphic gene content, inferences about the ecological and adaptive potential of such organisms also need to be accompanied by additional supportive genomic analyses. In this study we constructed a pangenome of Claviceps purpurea from 24 genomes and examined the positive selection and recombination landscape of an economically important fungal organism for pharmacology and agricultural research. Together, these analyses revealed that C. purpurea has a relatively large accessory genome (~ 38%), high recombination rates (ρ = 0.044), and transposon mediated gene duplication. However, due to observations of relatively low transposable element (TE) content (8.8%) and a lack of variability in genome sizes, prolific TE expansion may be controlled by frequent recombination. We additionally identified that within the ergoline biosynthetic cluster the lpsA1 and lpsA2 were the result of a recombination event. However, the high recombination rates observed in C. purpurea may be influencing an overall trend of purifying selection across the genome. These results showcase the use of selection and recombination landscapes to identify mechanisms contributing to pangenome structure and primary factors influencing the evolution of an organism.
Comparative analyses of mycobacteriophage genomes reveals extensive genetic diversity in genome organization and gene content, contributing to widespread mosaicism. We previously reported that the prophage of mycobacteriophage Butters (cluster N) provides defense against infection by Island3 (subcluster I1). To explore the anti-Island3 defense mechanism, we attempted to isolate Island3 defense escape mutants on a Butters lysogen, but only uncovered phages with recombinant genomes comprised of regions of Butters and Island3 arranged from left arm to right arm as Butters-Island3-Butters (BIBs). Recombination occurs within two distinct homologous regions that encompass lysin A, lysin B, and holin genes in one segment, and RecE and RecT genes in the other. Structural genes of mosaic BIB genomes are contributed by Butters while the immunity cassette is derived from Island3. Consequently, BIBs are morphologically identical to Butters (as shown by transmission electron microscopy) but are homoimmune with Island3. Recombinant phages overcome antiphage defense and silencing of the lytic cycle. We leverage this observation to propose a stratagem to generate novel phages for potential therapeutic use.
The phylogeny and timescale of life are becoming better understood as the analysis of genomic data from model organisms continues to grow. As a result, discoveries are being made about the early history of life and the origin and development of complex multicellular life. This emerging comparative framework and the emphasis on historical patterns is helping to bridge barriers among organism-based research communities.
Abstract Recent genomic analyses have revealed that microbial communities are predominantly composed of persistent, sequence-discrete species and intraspecies units (genomovars), but the mechanisms that create and maintain these units remain unclear. By analyzing closely-related isolate genomes from the same or related samples and identifying recent recombination events using a novel bioinformatics methodology, we show that high ecological cohesiveness coupled to frequent-enough and unbiased (i.e., not selection-driven) horizontal gene flow, mediated by homologous recombination, often underlie these diversity patterns. Ecological cohesiveness was inferred based on greater similarity in temporal abundance patterns of genomes of the same vs. different units, and recombination was shown to affect all sizable segments of the genome (i.e., be genome-wide) and have two times or greater impact on sequence evolution than point mutations. These results were observed in bothSalinibacter ruber, an environmental halophilic organism, andEscherichia coli, the model gut-associated organism and an opportunistic pathogen, indicating that they may be more broadly applicable to the microbial world. Therefore, our results represent a departure compared to previous models of microbial speciation that invoke either ecology or recombination, but not necessarily their synergistic effect, and answer an important question for microbiology: what a species and a subspecies are.
Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15–21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. In conclusion, these results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.
This award allowed us to develop novel hardware for phylogenetics, collect genomic data and produce several phylogenies of deuterostome organisms, communicate the results publicly, release software into the public domain, publish textbooks and papers, and prepare for the next research projects. There are no resulting subject inventions to report. We review these activities in three sections: 1) Hardware and software and development; 2) Evolutionary biology research; 3) Our proposed future direction, predictive analysis of pathogens in support of the NASA mission.
Conversion of CO 2 to value-added products presents an opportunity to reduce GHG emissions while generating revenue. Formate, which can be generated by the electrochemical reduction of CO 2 , has been proposed as a promising intermediate compound for microbial upgrading. Here we present progress towards improving the soil bacterium Cupriavidus necator H16, which is capable of growing on formate as its sole source of carbon and energy using the Calvin-Benson-Bassham (CBB) cycle, as a host for formate utilization. Using adaptive laboratory evolution, we generated several isolates that exhibited faster growth rates on formate. The genomes of these isolates were sequenced, and resulting mutations were systematically reintroduced by metabolic engineering, to identify those that improved growth. Further, the metabolic impact of several mutations was investigated further using RNA-seq transcriptomics. We found that deletion of a transcriptional regulator implicated in quorum sensing, PhcA, reduced expression of several operons and led to improved growth on formate. Growth was also improved by deleting large genomic regions present on the extrachromosomal megaplasmid pHG1, particularly two hydrogenase operons and the megaplasmid CBB operon, one of two copies present in the genome. Based on these findings, we generated a rationally engineered ..delta..phcA and megaplasmid-deficient strain that exhibited a 24% faster maximum growth rate on formate. Moreover, this strain achieved a 7% growth rate improvement on succinate and a 19% increase on fructose, demonstrating the broad utility of microbial genome reduction. This strain has the potential to serve as an improved microbial chassis for biological conversion of formate to value-added products.
From a micro-biology perspective, directed evolution is a technique that uses controlled environmental pressures to select for a desired phenotype. Directed evolution has the distinct advantage over rational design of not needing extensive knowledge of the genome or pathways associated with a microorganism to induce phenotypes. However, there are currently limitations to the applicability of this technique including being time-consuming, error-prone, and dependent on existing assays that may lack selectivity for the given phenotype. The AADEC (Autonomous Adaptive Directed Evolution Chamber) system is a proof-of-concept instrument to automate and improve the technique such that directed evolution can be used more effectively as a general bioengineering tool. A series of tests using the automated system and comparable by-hand survival assay measurements have been carried out using UV-C radiation and Escherichia coli cultures in order to demonstrate the advantages of the AADEC versus traditional implementations of directed evolution such as random mutagenesis. AADEC uses UV-C exposure as both a source of environmental stress and mutagenesis, so in order to evaluate the UV-C tolerance obtained from the cultures, a manual UV-C exposure survival assay was developed alongside the device to compare the survival fractions at a fixed dosage. This survival assay involves exposing E. coli to UV-C radiation using a custom-designed exposure hood to control the flux and dose. Surviving cells are counted then transferred to the next iteration and so on for several iterations to calculate the survival fractions for each exposure iteration.This survival assay primarily serves as a baseline for the AADEC device, allowing quantification of the differences between the AADEC system over the manual approach. The primary data of comparison is survival fractions; this is obtained by optical density and plate counts in the manual assay and by optical density growth curve fits pre- and post-exposure in the automated case. This data can then be compiled to calculate trends over the iterations to characterize increasing UV-C resistance of the E.coli strains. The observed trends are statistically indistinguishable through several iterations from both sources.
Expression profiling by high throughput sequencing The model yeast Saccharomyces cerevisiae has evolved aerobic fermentation in part through a whole genome duplication and subsequent retention of glycolytic genes. To identify independent yeast lineages with a high glycolytic rate phenotype we developed an assay measuring extracellular acidification rates (ECAR) to evaluate the glycolytic rates across diverse yeast species. This assay identified a novel group of yeasts in the genus Saturnispora with rapid glycolytic rates. Through comparative approaches, we found that several glycolytic genes encoding hexokinases and enzymes in lower glycolysis had higher expression and modifications in promoter sequences in rapid ECAR species. Intriguingly, many of the upregulated genes are the same genes that are duplicated in S. cerevisiae. When the transcription factor required for their activation was deleted, the mutants had a slow glycolytic rate and more respiratory phenotype. This work shows how the aerobic fermentation phenotype has convergently evolved in Saturnispora through transcriptional rewiring. This divergent mechanism nevertheless impacts the same glycolytic genes, which suggests that there are evolutionary constraints on how aerobic fermentation can arise.
Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.
The candidate phyla radiation (CPR) is a proposed subdivision within the bacterial domain comprising several candidate phyla. CPR organisms are united by small genome and physical sizes, lack several metabolic enzymes, and populate deep branches within the bacterial subtree of life. Here, these features raise intriguing questions regarding their origin and mode of evolution. In this study, we performed a comparative and phylogenomic analysis to investigate CPR origin and evolution. Unlike previous gene/protein sequence-based reports of CPR evolution, we used protein domain superfamilies classified by protein structure databases to resolve the evolutionary relationships of CPR with non-CPR bacteria, Archaea, Eukarya, and viruses. Across all supergroups, CPR shared maximum superfamilies with non-CPR bacteria and were placed as deep branching bacteria in most phylogenomic trees. CPR contributed 1.22% of new superfamilies to bacteria including the ribosomal protein L19e and encoded four core superfamilies that are likely involved in cell-to-cell interaction and establishing episymbiotic lifestyles. Although CPR and non-CPR bacterial proteomes gained common superfamilies over the course of evolution, CPR and Archaea had more common losses. These losses mostly involved metabolic superfamilies. In fact, phylogenies built from only metabolic protein superfamilies separated CPR and non-CPR bacteria. These findings indicate that CPR are bacterial organisms that have probably evolved in an Archaea-like manner via the early loss of metabolic functions. We also discovered that phylogenies built from metabolic and informational superfamilies gave contrasting views of the groupings among Archaea, Bacteria, and Eukarya, which add to the current debate on the evolutionary relationships among superkingdoms.
No abstract available
Organisms orchestrate cellular functions through transcription factor (TF) interactions with their target genes, although these regulatory relationships are largely unknown in most species. Here we report a high-throughput approach for characterizing TF-target gene interactions across species and its application to 354 TFs across 48 bacteria, generating 17,000 genome-wide binding maps. This dataset revealed themes of ancient conservation and rapid evolution of regulatory modules. We observed rewiring, where the TF sensing and regulatory role is maintained while the arrangement and identity of target genes diverges, in some cases encoding entirely new functions. We further integrated phenotypic information to define new functional regulatory modules and pathways. So, we identified 242 new TF DNA binding motifs, including a 70% increase of known Escherichia coli motifs and the first annotation in Pseudomonas simiae, revealing deep conservation in bacterial promoter architecture. Our method provides a versatile tool for functional characterization of genetic pathways in prokaryotes and eukaryotes.
Intracellular biomacromolecules, such as industrial enzymes and biopolymers, represent an important class of bio-derived products obtained from bacterial hosts. A common key step in the downstream separation of these biomolecules is lysis of the bacterial cell wall to effect release of cytoplasmic contents. Cell lysis is typically achieved either through mechanical disruption or reagent-based methods, which introduce issues of energy demand, material needs, high costs, and scaling problems. Osmolysis, a cell lysis method that relies on hypoosmotic downshock upon resuspension of cells in distilled water, has been applied for bioseparation of intracellular products from extreme halophiles and mammalian cells. However, most industrial bacterial strains are non-halotolerant and relatively resistant to hypoosmotic cell lysis. To overcome this limitation, we developed two strategies to increase the susceptibility of non-halotolerant hosts to osmolysis using Cupriavidus necator, a strain often used in electromicrobial production, as a prototypical strain. In one strategy, C. necator was evolved to increase its halotolerance from 1.5% to 3.25% (w/v) NaCl through adaptive laboratory evolution, and genes potentially responsible for this phenotypic change were identified by whole genome sequencing. The evolved halotolerant strain experienced an osmolytic efficiency of 47% in distilled water following growth in 3% (w/v) NaCl. In a second strategy, the cells were made susceptible to osmolysis by knocking out the large-conductance mechanosensitive channel (mscL) gene in C. necator. When these strategies were combined by knocking out the mscL gene from the evolved halotolerant strain, greater than 90% osmolytic efficiency was observed upon osmotic downshock. A modified version of this strategy was applied to E. coli BL21 by deleting the mscL and mscS (small-conductance mechanosensitive channel) genes. When grown in medium with 4% NaCl and subsequently resuspended in distilled water, this engineered strain experienced 75% cell lysis, although decreases in cell growth rate due to higher salt concentrations were observed. Our strategy is shown to be a simple and effective way to lyse cells for the purification of intracellular biomacromolecules and may be applicable in many bacteria used for bioproduction.
Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.