Search NASA⌕ Search

SEARCH · Search NASA

Results for “population genomics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Factors underlying a latitudinal gradient in the S/G lignin monomer ratio in natural poplar variants

The chemical composition of wood plays a pivotal role in the adaptability and structural integrity of trees. However, few studies have investigated the environmental factors that determine lignin composition and its biological significance in plants. Here, we examined the lignin syringyl-to-guaiacyl (S/G) ratio in members of a Populus trichocarpa population sourced from their native habitat and conducted a genome wide association study to identify genes linked to lignin formation. Our results revealed many significant associations, suggesting that lignin biosynthesis is a complex polygenic trait. Additionally, we found an increase in the S/G ratio from northern to southern geographic origin of the trees sampled, along with a corresponding metabolic and transcriptional reprogramming of xylem cell wall biosynthesis. Further molecular analysis identified a mutation in a cell wall laccase genetically associated with higher S/G ratios that predominate in trees from warmer lower latitudes. Collectively, our findings suggest that lignin heterogeneity arises from an evolutionary process enabling poplar adaptation to different climatic challenges.

adaptation↗
DOE OSTI2020

Bacterial Origin and Reductive Evolution of the CPR Group

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.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust

Biological soil crusts (biocrusts) are photosynthetic “hot spots” in deserts and cover ~12% of the Earth’s terrestrial surface, and yet they face an uncertain future given expected shifts in rainfall events. Laboratory wetting of biocrust communities is known to cause a bloom of Firmicutes which rapidly become dominant community members within 2 days after emerging from a sporulated state. We hypothesized that their bacteriophages (phages) would respond to such a dramatic increase in their host’s abundance. In our experiment, wetting caused Firmicutes to bloom and triggered a significant depletion of cyanobacterial diversity. We used genome-resolved metagenomics to link phage to their hosts and found that the bloom of the genus Bacillus correlated with a dramatic increase in the number of Caudovirales phages targeting these diverse spore-formers ( r = 0.762). After 2 days, we observed dramatic reductions in the relative abundances of Bacillus , while the number of Bacillus phages continued to increase, suggestive of a predator-prey relationship. We found predicted auxiliary metabolic genes (AMGs) associated with sporulation in several Caudovirales genomes, suggesting that phages may influence and even benefit from sporulation dynamics in biocrusts. Prophage elements and CRISPR-Cas repeats in Firmicutes metagenome-assembled genomes (MAGs) provide evidence of recent infection events by phages, which were corroborated by mapping viral contigs to their host MAGs. Combined, these findings suggest that the blooming Firmicutes become primary targets for biocrust Caudovirales phages, consistent with the classical “kill-the-winner” hypothesis. This work forms part of an overarching research theme studying the effects of a changing climate on biological soil crust (biocrust) in the Southwestern United States. To our knowledge, this study was the first to characterize bacteriophages in biocrust and offers a view into the ecology of phages in response to a laboratory wetting experiment. The phages identified here represent lineages of Caudovirales , and we found that the dynamics of their interactions with their Firmicutes hosts explain the collapse of a bacterial bloom that was induced by wetting. Moreover, we show that phages carried host-altering metabolic genes and found evidence of proviral infection and CRISPR-Cas repeats within host genomes. Our results suggest that phages exert controls on population density by lysing dominant bacterial hosts and that they further impact biocrust by acquiring host genes for sporulation. Future research should explore how dominant these phages are in other biocrust communities and quantify how much the control and lysis of blooming populations contributes to nutrient cycling in biocrusts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Syntrophic Co-Cultures of Clostridium Organisms to Produce Higher Alcohols & Other C6-C8 Metabolites (Final Report)

The goal of this project was to advance the systems biology understanding and predictive modeling of synthetic & syntrophic Clostridium microbial consortia, focusing on elucidation of metabolic networks and environmental signals in the consortia. The project has direct applicability to lignocellulosic-biomass based production of higher alcohols as advanced biofuels and C6-C8 metabolites, that can be used as chemicals or serve as biofuel precursors. This project studied and optimized three synthetic syntrophic systems. It examined the population dynamics using flow-cytometry, time-lapse microscopy and PCR analysis. 13C-based tracer analyses was used to examine the metabolite exchange between the syntrophic cell populations and the impact of those interactions on the transcriptome of the individual populations. To enhance our analytical and predictive capabilities, genome-scale models (GSMs) for these syntrophies was developed. RNAseq data for these syntrophic coculture systems were acquired to enable a molecular level understanding of the syntrophies aiming to identify the genetic networks of each organism in the co-culture and compare those against the networks of pure cultures.

09 BIOMASS FUELS↗
DOE OSTI2020

Steelhead ( Oncorhynchus mykiss ) lineages and sexes show variable patterns of association of adult migration timing and age‐at‐maturity traits with two genomic regions

Abstract As life history diversity plays a critical role in supporting the resilience of exploited populations, understanding the genetic basis of those life history variations is important for conservation management. However, effective application requires a robust understanding of the strength and universality of genetic associations. Here, we examine genetic variation of single nucleotide polymorphisms in genomic regions previously associated with migration phenology and age‐at‐maturity in steelhead ( Oncorhynchus mykiss ) from the Columbia River. We found chromosome 28 markers (GREB1L, ROCK1 genes) explained significant variance in migration timing in both coastal and inland steelhead. However, strength of association was much greater in coastal than inland steelhead ( R 2 0.51 vs. 0.08), suggesting that genomic background and challenging inland migration pathways may act to moderate effects of this region. Further, we found that chromosome 25 candidate markers (SIX6 gene) were significantly associated with age and size at first return migration for inland steelhead, and this pattern was mediated by sex in a predictable pattern (males R 2 = 0.139–0.170; females R 2 = 0.096–0.111). While this encourages using these candidate regions in predicting life history characteristics, we suggest that stock‐specific associations and haplotype frequencies will be useful in guiding implementation of genetic assays to inform management.

Willis, Stuart C.↗
DOE OSTI2020

Reduced meiotic recombination in rhesus macaques and the origin of the human recombination landscape

Characterizing meiotic recombination rates across the genomes of nonhuman primates is important for understanding the genetics of primate populations, performing genetic analyses of phenotypic variation and reconstructing the evolution of human recombination. Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primates in biomedical research. We constructed a high-resolution genetic map of the rhesus genome based on whole genome sequence data from Indian-origin rhesus macaques. The genetic markers used were approximately 18 million SNPs, with marker density 6.93 per kb across the autosomes. We report that the genome-wide recombination rate in rhesus macaques is significantly lower than rates observed in apes or humans, while the distribution of recombination across the macaque genome is more uniform. These observations provide new comparative information regarding the evolution of recombination in primates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Building a FAIR data ecosystem for incorporating single-cell transcriptomics data into agricultural genome to phenome research

Introduction The agriculture genomics community has numerous data submission standards available, but the standards for describing and storing single-cell (SC, e.g., scRNA- seq) data are comparatively underdeveloped. Methods To bridge this gap, we leveraged recent advancements in human genomics infrastructure, such as the integration of the Human Cell Atlas Data Portal with Terra, a secure, scalable, open-source platform for biomedical researchers to access data, run analysis tools, and collaborate. In parallel, the Single Cell Expression Atlas at EMBL-EBI offers a comprehensive data ingestion portal for high-throughput sequencing datasets, including plants, protists, and animals (including humans). Developing data tools connecting these resources would offer significant advantages to the agricultural genomics community. The FAANG data portal at EMBL-EBI emphasizes delivering rich metadata and highly accurate and reliable annotation of farmed animals but is not computationally linked to either of these resources. Results Herein, we describe a pilot-scale project that determines whether the current FAANG metadata standards for livestock can be used to ingest scRNA-seq datasets into Terra in a manner consistent with HCA Data Portal standards. Importantly, rich scRNA-seq metadata can now be brokered through the FAANG data portal using a semi-automated process, thereby avoiding the need for substantial expert curation. We have further extended the functionality of this tool so that validated and ingested SC files within the HCA Data Portal are transferred to Terra for further analysis. In addition, we verified data ingestion into Terra, hosted on Azure, and demonstrated the use of a workflow to analyze the first ingested porcine scRNA-seq dataset. Additionally, we have also developed prototype tools to visualize the output of scRNA-seq analyses on genome browsers to compare gene expression patterns across tissues and cell populations. This JBrowse tool now features distinct tracks, showcasing PBMC scRNA-seq alongside two bulk RNA-seq experiments. Discussion We intend to further build upon these existing tools to construct a scientist-friendly data resource and analytical ecosystem based on Findable, Accessible, Interoperable, and Reusable (FAIR) SC principles to facilitate SC-level genomic analysis through data ingestion, storage, retrieval, re-use, visualization, and comparative annotation across agricultural species.

Genetics & Heredity↗
DOE OSTI2022

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species’ geographic ranges, increasing extinction risk and losses of species’ genetic diversity. Although preserving genetic diversity is key to maintaining species’ adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations’ post-2020 targets for genetic preservation.

Science & Technology - Other Topics↗
DOE OSTI2025

Mondo: integrating disease terminology across communities

Precision medicine aims to enhance diagnosis, treatment, and prognosis by integrating multimodal data at the point of care. However, challenges arise due to the vast number of diseases, differing methods of classification, and conflicting terminological coding systems and practices used to represent molecular definitions of disease. This lack of interoperability artificially constrains the potential for diagnosis, clinical decision support, care outcome analysis, as well as data linkage across research domains to support the development or repurposing of therapeutics. There is a clear and pressing need for a unified system for managing disease entities⁠—including identifiers, synonyms, and definitions. To address these issues, we created the Mondo disease ontology—a community-driven, open-source, unified disease classification system that harmonizes diverse terminologies into a consistent, computable framework. Mondo integrates key medical and biomedical terminologies, including Online Mendelian Inheritance in Man (OMIM), Orphanet, Medical Subject Headings (MeSH), National Cancer Institute Thesaurus (NCIt), and more, to provide a comprehensive and accurate representation of disease concepts with fully provenanced and attributed links back to the sources. Mondo can be used as the handle for curation of gene–disease associations utilized in diagnostic applications, research applications such as computational phenotyping, and in clinical coding systems in clinical decision support by pointing the clinician to the numerous knowledge resources linked to the Mondo identifier. Mondo's community-centric approach, stewarded by the Monarch Initiative's expertise in ontologies, ensures that the ontology remains adaptable to the evolving needs of biomedical research and clinical communities, as well as the knowledge providers.

biomedical informatics↗
DOE OSTI2025

Nitrogen fixation in Pampinifervens , a new species-rich genus of Aquificaceae that inhabits a wide pH range in terrestrial hot springs

Members of the Aquificota play important roles in hydrogen, sulfur, and carbon cycles in geothermal systems as ubiquitous and abundant chemolithoautotrophs; however, their roles in the nitrogen cycle are poorly defined. Here, we show that Aquificaceae isolate T-2 T from Tengchong, China, fixes nitrogen based on growth with dinitrogen as the sole nitrogen source and incorporation of 15 N atoms when grown with 15 N 2 . We further show evidence that suggests its nitrogen fixation (nif) genes are transcribed in situ in a nearby spring. We propose the name Pampinifervens diazotrophicum gen. et sp. nov., with strain T-2 T (= JCM 35475 T = DSM 116324 T ) as the nomenclatural type. Its close relative designated Pampinifervens florentissimum sp. nov. T-8 T (= CGMC 1.5214 T = JCM 33569 T ) does not fix nitrogen but respires nitrate. Analysis of 104 metagenome-assembled genomes from 84 metagenomes allowed us to circumscribe the genus Pampinifervens and identify at least 16 sympatric Pampinifervens species that dominate the Aquificaceae populations in diverse hot springs in southwest China, with distinct species globally. Analysis of these Pampinifervens genomes revealed variable presence of nif and denitrification genes and frequent gene loss and horizontal transfer. Genes for pH homeostasis were also present in Pampinifervens genomes, with a proliferation in species from Tengchong, consistent with the wide pH growth range of both Pampinifervens isolates and the prevalence of Pampinifervens in both alkaline and acidic springs. The dynamic evolution of nif and denitrification genes, along with pH homeostasis genes, may facilitate the diversification of Pampinifervens into diverse springs with different nitrogen availability.

Aquificae↗
DOE OSTI2022

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↗
DOE OSTI2026

A haplotype-resolved, chromosome-scale genome assembly for the southern live oak, Quercus virginiana

Hybridization is a major force driving diversification, migration, and adaptation in Quercus species. While population genetics and phylogenetics have traditionally been used for studying these processes, advances in sequencing technology now enable us to incorporate comparative and pan-genomic approaches as well. Here, we present a highly contiguous, chromosome-scale and haplotype-resolved genome assembly for the southern live oak, Quercus virginiana, the first reference genome for section Virentes, as part of the American Campus Tree Genomes program. Originating from a clone of Auburn University's historic “Toomer's Oak,” this assembly contributes to the pool of genomic resources for investigating recombination, haplotype variation, and structural genomic changes influencing hybridization potential in this clade and across Quercus. It also provides insights into the architecture of the putative centromeric regions within the genus. Alongside other oak references, the Q. virginiana genome will support research into the evolution and adaptation of the Quercus genus.

Quercus virginiana↗
DOE OSTI2024

Elemental profiling and genomewide association studies reveal genomic variants modulating ionomic composition in Populus trichocarpa leaves

Samples were collected from a population of 1,089 black cottonwood genotypes (P. trichocarpa) assembled from native stands to encompass the central portion of the natural range of the species, stretching from 38.8° to 54.3° N Q13 latitude from California, USA, to British Columbia, Canada. Establishment of the common garden, growth conditions, and site maintenance have been described by Muchero et al (2015). In this study, leaf samples for ionomic profiling were collected from 4-year-old trees, during the growing season, in July 2012, from a field located in Clatskanie, Oregon, USA (46°6′11″N 123°12′13″W). The field site was located in a protected alluvial floodplain containing a uniform Wauna-Locoda silt loam soil area characterized by an acidic pH, in Columbia County, Oregon. A subset of 584 out of the 1,089 P. trichocarpa genotypes were represented in this sampling. These genotypes were randomly selected to represent the geographical distribution of the population. A single fully mature (LPI 7-9) leaf on the south side of the tree exposed to full sunlight conditions was removed from the tree within a 6-hour window centering on solar noon and immediately frozen under dry ice before processing. Leaf samples of 584 P. trichocarpa genotypes were finely ground to 40 mm particle size using a mortar and pestle, and ionomic composition was analyzed using ICP-MS. In total, 20 elements were profiled, including aluminum (Al27), arsenic (As75), boron (B11), cadmium (Cd111), calcium (Ca43), cobalt (Co), copper (Cu), iron (Fe57), magnesium (Mg25), manganese (Mn55), molybdenum (Mo), nickel (Ni60), phosphorus (P31), potassium (K39), rubidium (Rb85), selenium (Se82), sodium (Na23), strontium (Sr88), sulfur (S34), and zinc (Zn66), following a protocol established by Ziegler et al. (2013). For each sample, 75mg of powder was digested overnight in 2.5 mL HNO3 containing 20 parts per billion (ppb) indium as an internal standard, following the protocol described in Ziegler et al. (2013). Following a dilution, concentration of the 20 elements was measured using an Elan 6000 DRC-e mass spectrometer (Perkin-Elmer SCIEX) connected to a PFA microflow nebulizer (Elemental Scientific) and Apex HF desolvator (Elemental Scientific). One measurement per sample per genotype was done. For subsequent analyses, the quantifications were converted to total element concentration.

CBI ionomics, GWAS, plasma-mass spectrometry, neut↗
DOE OSTI2023

Phylogeography of the blacklegged tick ( Ixodes scapularis ) throughout the USA identifies candidate loci for differences in vectorial capacity

Abstract The blacklegged tick ( Ixodes scapularis ( Journal of the Academy of Natural Sciences of Philadelphia , 1821, 2 , 59)) is a vector of Borrelia burgdorferi sensu stricto ( s.s .) ( International Journal of Systematic Bacteriology , 1984, 34 , 496), the causative bacterial agent of Lyme disease, part of a slow‐moving epidemic of Lyme borreliosis spreading across the northern hemisphere. Well‐known geographical differences in the vectorial capacity of these ticks are associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographical studies of these ticks have been restricted to relatively few populations or few genetic loci. Here we present the most comprehensive phylogeographical study of genome‐wide markers in I. scapularis , conducted by using 3RAD (triple‐enzyme restriction‐site associated sequencing) and surveying 353 ticks from 33 counties throughout the species' range. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify five spatially associated genetic clusters of I. scapularis . In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme‐endemic regions. Finally, we identify 10 variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome‐scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis .

3RAD↗
DOE OSTI2023

The first two chromosome‐scale genome assemblies of American hazelnut enable comparative genomic analysis of the genus Corylus

Summary The native, perennial shrub American hazelnut ( Corylus americana ) is cultivated in the Midwestern United States for its significant ecological benefits, as well as its high‐value nut crop. Implementation of modern breeding methods and quantitative genetic analyses of C. americana requires high‐quality reference genomes, a resource that is currently lacking. We therefore developed the first chromosome‐scale assemblies for this species using the accessions ‘Rush’ and ‘Winkler’. Genomes were assembled using HiFi PacBio reads and Arima Hi‐C data, and Oxford Nanopore reads and a high‐density genetic map were used to perform error correction. N50 scores are 31.9 Mb and 35.3 Mb, with 90.2% and 97.1% of the total genome assembled into the 11 pseudomolecules, for ‘Rush’ and ‘Winkler’, respectively. Gene prediction was performed using custom RNAseq libraries and protein homology data. ‘Rush’ has a BUSCO score of 99.0 for its assembly and 99.0 for its annotation, while ‘Winkler’ had corresponding scores of 96.9 and 96.5, indicating high‐quality assemblies. These two independent assemblies enable unbiased assessment of structural variation within C. americana , as well as patterns of syntenic relationships across the Corylus genus. Furthermore, we identified high‐density SNP marker sets from genotyping‐by‐sequencing data using 1343 C. americana , C. avellana and C. americana × C. avellana hybrids, in order to assess population structure in natural and breeding populations. Finally, the transcriptomes of these assemblies, as well as several other recently published Corylus genomes, were utilized to perform phylogenetic analysis of sporophytic self‐incompatibility (SSI) in hazelnut, providing evidence of unique molecular pathways governing self‐incompatibility in Corylus .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Metagenomics-resolved genomics provides novel insights into chitin turnover, metabolic specialization, and niche partitioning in the octocoral microbiome

Abstract Background The role of bacterial symbionts that populate octocorals (Cnidaria, Octocorallia) is still poorly understood. To shed light on their metabolic capacities, we examined 66 high-quality metagenome-assembled genomes (MAGs) spanning 30 prokaryotic species, retrieved from microbial metagenomes of three octocoral species and seawater. Results Symbionts of healthy octocorals were affiliated with the taxa Endozoicomonadaceae , Candidatus Thioglobaceae , Metamycoplasmataceae , unclassified Pseudomonadales , Rhodobacteraceae , unclassified Alphaproteobacteria and Ca. Rhabdochlamydiaceae . Phylogenomics inference revealed that the Endozoicomonadaceae symbionts uncovered here represent two species of a novel genus unique to temperate octocorals, here denoted Ca. Gorgonimonas eunicellae and Ca. Gorgonimonas leptogorgiae . Their genomes revealed metabolic capacities to thrive under suboxic conditions and high gene copy numbers of serine-threonine protein kinases, type 3-secretion system, type-4 pili, and ankyrin-repeat proteins, suggesting excellent capabilities to colonize, aggregate, and persist inside their host. Contrarily, MAGs obtained from seawater frequently lacked symbiosis-related genes. All Endozoicomonadaceae symbionts harbored endo-chitinase and chitin-binging protein-encoding genes, indicating that they can hydrolyze the most abundant polysaccharide in the oceans. Other symbionts, including Metamycoplasmataceae and Ca. Thioglobaceae , may assimilate the smaller chitin oligosaccharides resulting from chitin breakdown and engage in chitin deacetylation, respectively, suggesting possibilities for substrate cross-feeding and a role for the coral microbiome in overall chitin turnover. We also observed sharp differences in secondary metabolite production potential between symbiotic lineages. Specific Proteobacteria taxa may specialize in chemical defense and guard other symbionts, including Endozoicomonadaceae , which lack such capacity. Conclusion This is the first study to recover MAGs from dominant symbionts of octocorals, including those of so-far unculturable Endozoicomonadaceae , Ca. Thioglobaceae and Metamycoplasmataceae symbionts. We identify a thus-far unanticipated, global role for Endozoicomonadaceae symbionts of corals in the processing of chitin, the most abundant natural polysaccharide in the oceans and major component of the natural zoo- and phytoplankton feed of octocorals. We conclude that niche partitioning, metabolic specialization, and adaptation to low oxygen conditions among prokaryotic symbionts likely contribute to the plasticity and adaptability of the octocoral holobiont in changing marine environments. These findings bear implications not only for our understanding of symbiotic relationships in the marine realm but also for the functioning of benthic ecosystems at large.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Harnessing the predicted maize pan-interactome for putative gene function prediction and prioritization of candidate genes for important traits

Abstract The recent assembly and annotation of the 26 maize nested association mapping population founder inbreds have enabled large-scale pan-genomic comparative studies. These studies have expanded our understanding of agronomically important traits by integrating pan-transcriptomic data with trait-specific gene candidates from previous association mapping results. In contrast to the availability of pan-transcriptomic data, obtaining reliable protein–protein interaction (PPI) data has remained a challenge due to its high cost and complexity. We generated predicted PPI networks for each of the 26 genomes using the established STRING database. The individual genome-interactomes were then integrated to generate core- and pan-interactomes. We deployed the PPI clustering algorithm ClusterONE to identify numerous PPI clusters that were functionally annotated using gene ontology (GO) functional enrichment, demonstrating a diverse range of enriched GO terms across different clusters. Additional cluster annotations were generated by integrating gene coexpression data and gene description annotations, providing additional useful information. We show that the functionally annotated PPI clusters establish a useful framework for protein function prediction and prioritization of candidate genes of interest. Our study not only provides a comprehensive resource of predicted PPI networks for 26 maize genomes but also offers annotated interactome clusters for predicting protein functions and prioritizing gene candidates. The source code for the Python implementation of the analysis workflow and a standalone web application for accessing the analysis results are available at https://github.com/eporetsky/PanPPI.

Genetics & Heredity↗