Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phylogenetics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

New species and records of the symbiotic shrimp genus Leptalpheus Williams, 1965, with notes on Fenneralpheus Felder & Manning, 1986, and preliminary molecular analysis of phylogenetic relationships (Crustacea: Decapoda: Alpheidae)

The shrimp genera Leptalpheus Williams, 1965 and Fenneralpheus Felder & Manning, 1986 are composed entirely of symbiotic species that co-inhabit burrows of infaunal macrocrustaceans. We report extensive collections of these genera from western Atlantic, eastern Pacific and Indo-West Pacific regions. Integrative taxonomy methods, including morphological comparisons and analysis of three mitochondrial genetic markers, are used to test species hypotheses and evolutionary relationships among members of these genera. Our molecular analysis failed to recover Leptalpheus or Fenneralpheus as monophyletic groups. Our results strongly supported the monophyly of three clades composed of species of Leptalpheus, loosely corresponding to previously proposed species groups. Three new species closely related to Leptalpheus forceps Williams, 1965, L. marginalis Anker, 2011, and L. mexicanus Ríos & Carvacho, 1983 are described. Leptalpheus ankeri n. sp., from the Caribbean Sea, Atlantic coast of Florida, and Gulf of Mexico, is a polymorphic species that exhibits two major cheliped morphotypes. Leptalpheus sibo n. sp., from the Pacific coast of Nicaragua, is morphologically very similar to L. ankeri n. sp., likely its transisthmian sister species, and shares its cheliped polymorphism. A reassessment of L. forceps concluded that records of this species from the Caribbean Sea and Brazil are not conspecific with L. forceps sensu stricto from the Atlantic coast of the USA and the Gulf of Mexico, and they are herein described as Leptalpheus degravei n. sp. Based on both molecular and morphological evidence, we found Leptalpheus bicristatus Anker, 2011 to be a junior synonym of L. mexicanus and Leptalpheus canterakintzi Anker & Lazarus, 2015 to be a junior synonym of Leptalpheus azuero Anker, 2011. First reports of Leptalpheus axianassae Dworschak & Coelho, 1999 in Texas and Mexico, Leptalpheus denticulatus Anker & Marin, 2009 in the Mariana Islands, Leptalpheus felderi Anker, Vera Caripe & Lira, 2006 and Leptalpheus lirai Vera Caripe, Pereda & Anker, 2021 in the USA, and Leptalpheus pereirai Anker & Vera Caripe, 2016 in Cuba are included.

Zoology↗
DOE OSTI2025

Data for Phylogenetic diversity of light-dependent phosphorylation of Thr78 in Rubisco activase

Rubisco activase is an ATP-dependent chaperone that facilitates dissociation of inhibitory sugar phosphates from the catalytic sites of Rubisco during photosynthesis. In Arabidopsis, Rubisco activase is negatively regulated by dark-dependent phosphorylation of Thr78. The prevalence of Thr78 in Rubisco activase was investigated across sequences from 91 plant species, finding that 29 (∼32%) species shared a threonine in the same position. Analysis of seven C3 species with an antibody raised against a Thr78 phospho-peptide demonstrated that this position is phosphorylated in multiple genera. However, light-dependent dephosphorylation of Thr78 was observed only in Arabidopsis. Further, phosphorylation of Thr78 could not be detected in any of the four C4 grass species examined. The results suggest that despite conservation of Thr78 in Rubisco activase from a wide range of species, a regulatory role for phosphorylation at this site is more limited. This provides a case study for how variation in post-translational regulation can amplify functional divergence across the phylogeny of plants beyond what is explained by sequence variation in a metabolically important protein.

photosynthesis↗
DOE OSTI2024

Severe and mild drought cause distinct phylogenetically linked shifts in the blue grama (Bouteloua gracilis) rhizobiome

Plants rely on a diverse rhizobiome to regulate nutrient acquisition and plant health. With increasing severity and frequency of droughts worldwide due to climate change, untangling the relationships between plants and their rhizobiomes is vital to maintaining agricultural productivity and protecting ecosystem diversity. While some plant physiological responses to drought are generally conserved, patterns of root exudation (release of small metabolites shown to influence microbes) and the consequential effects on the plant rhizobiome can differ widely across plant species under drought. To address this knowledge gap, we conducted a greenhouse study using blue grama ( Bouteloua gracilis ), a drought-tolerant C4 grass native to shortgrass prairie across North American plains, as a model organism to study the effect of increasing drought severity (ambient, mild drought, severe drought) on root exudation and the rhizobiome. Our previous results demonstrated physiological effects of increasing drought severity including an increase in belowground carbon allocation through root exudation and shifts in root exudate composition concurrent with the gradient of drought severity. This work is focused on the rhizobiome community structure using targeted sequencing and found that mild and severe drought resulted in unique shifts in the bacterial + archaeal and fungal communities relative to ambient, non-droughted controls. Specifically, using the change in relative abundance between ambient and drought conditions for each ZOTU as a surrogate for population-scale drought tolerance (e.g., as a response trait), we found that rhizobiome response to drought was non-randomly distributed across the phylogenies of both communities, suggesting that Planctomycetota , Thermoproteota (formerly Thaumarchaeota ), and the Glomeromycota were the primary clades driving these changes. Correlation analyses indicated weak correlations between droughted community composition and a select few root exudate compounds previously implicated in plant drought responses including pyruvic acid, D-glucose, and myoinositol. This study demonstrates the variable impacts of drought severity on the composition of the blue grama rhizobiome and provides a platform for hypothesis generation for targeted functional studies of specific taxa involved in plant-microbe drought responses.

Goemann, Hannah M.↗
DOE OSTI2022

Description and Genome Characterization of Three Novel Fungal Strains Isolated from Mars 2020 Mission-Associated Spacecraft Assembly Facility Surfaces—Recommendations for Two New Genera and One Species

National Aeronautics and Space Administration’s (NASA) spacecraft assembly facilities are monitored for the presence of any bacteria or fungi that might conceivably survive a transfer to an extraterrestrial environment. Fungi present a broad and diverse range of phenotypic and functional traits to adapt to extreme conditions, hence the detection of fungi and subsequent eradication of them are needed to prevent forward contamination for future NASA missions. During the construction and assembly for the Mars 2020 mission, three fungal strains with unique morphological and phylogenetic properties were isolated from spacecraft assembly facilities. The reconstruction of phylogenetic trees based on several gene loci (ITS, LSU, SSU, RPB, TUB, TEF1) using multi-locus sequence typing (MLST) and whole genome sequencing (WGS) analyses supported the hypothesis that these were novel species. Here we report the genus or species-level classification of these three novel strains via a polyphasic approach using phylogenetic analysis, colony and cell morphology, and comparative analysis of WGS. The strain FJI-L9-BK-P1 isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility (JPL-SAF) exhibited a putative phylogenetic relationship with the strain Aaosphaeria arxii CBS175.79 but showed distinct morphology and microscopic features. Another JPL-SAF strain, FJII-L3-CM-DR1, was phylogenetically distinct from members of the family Trichomeriaceae and exhibited morphologically different features from the genera Lithohypha and Strelitziana. The strain FKI-L1-BK-DR1 isolated from the Kennedy Space Center facility was identified as a member of Dothideomycetes incertae sedis and is closely related to the family Kirschsteiniotheliaceae according to a phylogenetic analysis. The polyphasic taxonomic approach supported the recommendation for establishing two novel genera and one novel species. The names Aaosphaeria pasadenensis (FJI-L9-BK-P1 = NRRL 64424 = DSM 114621), Pasadenomyces melaninifex (FJII-L3-CM-DR1 = NRRL 64433 = DSM 114623), and Floridaphiala radiotolerans (FKI-L1-BK-DR1 = NRRL 64434 = DSM 114624) are proposed as type species. Furthermore, resistance to ultraviolet-C and presence of specific biosynthetic gene cluster(s) coding for metabolically active compounds are unique to these strains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. In this work, we demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Description of Pegethrix niliensis sp. nov., a Novel Cyanobacterium from the Nile River Basin, Egypt: A Polyphasic Analysis and Comparative Study of Related Genera in the Oculatellales Order

In this paper, we examine the filamentous cyanobacterial strain NILCB16 and describe it as a new species within the genus Pegethrix. The original population was sampled from a mat growing in an irrigation canal in the Nile River, Egypt. Initially classified under Plectonema or Planktolyngbya, the strain is a potential producer of the toxins microcystin and β-N-Methylamino-L-Alanine (BMAA). Additionally, we reviewed the taxonomic relationships between the Oculatellales genera. To describe the new species, we conducted a polyphasic study, encompassing 16S rRNA gene phylogenetic analyses performed using both Maximum Likelihood and Bayesian methods, sequence identity (p-distance) analysis, 16S-23S ITS secondary structures, and morphological and habitat comparisons. The phylogenetic analysis revealed that strain NILCB16 clustered within the Pegethrix clade with strong phylogenetic support, but in a distinct position from other species in the genus. The strain shared a maximum 16S rRNA gene identity of 97.3% with P. qiandaoensis and 96.1% with the type species, P. bostrychoides. Morphologically, NILCB16 can be differentiated from other species in the genus by its lack of false branching. Our phylogenetic analyses also show that Pegethrix, Cartusia, Elainella, and Maricoleus are clustered with strong phylogenetic support. They exhibit high 16S rRNA gene identity and are morphologically indistinguishable, suggesting they could potentially be merged into a single genus in the future.

Hentschke, Guilherme Scotta (ORCID:000000034396024↗
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 OSTI2021

New Data Define the Molecular Phylogeny and Taxonomy of Four Freshwater Suctorian Ciliates With Redefinition of Two Families Heliophryidae and Cyclophryidae (Ciliophora, Phyllopharyngea, Suctoria)

Four suctorian ciliates, Cyclophrya magna Gönnert, 1935, Peridiscophrya florea (Kormos & Kormos, 1958) Dovgal, 2002, Heliophrya rotunda (Hentschel, 1916) Matthes, 1954 and Dendrosoma radians Ehrenberg, 1838, were collected from a freshwater lake in Ningbo, China. The morphological redescription and molecular phylogenetic analyses of these ciliates were investigated. Phylogenetic analyses inferred from SSU rDNA sequences show that all three suctorian orders, Endogenida, Evaginogenida, and Exogenida, are monophyletic and that the latter two clusters as sister clades. The newly sequenced P. florea forms sister branches with C. magna , while sequences of D. radians group with those from H. rotunda within Endogenida. The family Heliophryidae, which is comprised of only two genera, Heliophrya and Cyclophrya , was previously assigned to Evaginogenida. There is now sufficient evidence, however, that the type genus Heliophrya reproduces by endogenous budding, which corresponds to the definitive feature of Endogenida. In line with this and with the support of molecular phylogenetic analyses, we therefore transfer the family Heliophryidae with the type genus Heliophrya to Endogenida. The other genus, Cyclophrya , still remains in Evaginogenida because of its evaginative budding. Therefore, combined with morphological and phylogenetic analysis, Cyclophyidae are reactivated, and it belongs to Evaginogenida.

Ma, Mingzhen↗
DOE OSTI2022

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development

The transcription factor WRINKLED1 ( WRI1 ) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI 1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase ( FRK 3) and of phosphoglucose isomerase ( PGI 1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Fire and Drought Affect Multiple Aspects of Diversity in a Migratory Bird Stopover Community

Drought and high-severity, stand-replacing wildfires can have substantial impacts on the composition of avian communities, including stop-over communities during migration. An inextricable link exists between drought and wildfire, each operating and impacting across different timescales. Many studies have found nonlinear avian abundance trends in breeding community time series data that include pre- and post-fire observations, describing an initial decrease in abundance followed by rapid increases that can attenuate over time. Here, we use a fall bird-banding dataset to evaluate shifts in a drought-impacted avian community following wildfire from taxonomic, functional, and phylogenetic perspectives. We looked at the community as a whole and also categorized birds as residents, migrants, and breeders to assess potential varying responses at the study site. We observed post-fire shifts in functional and phylogenetic diversity that corresponded to changes in vegetation. An influx of migratory insectivores post-fire drove much of the variation between pre- and post-fire avian communities and toward a more related, less phylogenetically dispersed community. A concurrent monsoon season drought was also associated with functional and phylogenetic diversity, highlighting the intertwined pulse press effects on avian communities. Overall, our results suggest that, although bird communities are immediately impacted by fire-driven resource changes, they can rebound over time, it is unclear how long-term drought may continue to shape the composition of these avian communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Origin and evolution of HIV-1 subtype A6

Background: HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes. Methods: For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used. Results: Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1. Conclusion: The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Modest functional diversity decline and pronounced composition shifts of microbial communities in a mixed waste-contaminated aquifer

Background: Microbial taxonomic diversity declines with increased environmental stress. Yet, few studies have explored whether phylogenetic and functional diversities track taxonomic diversity along the stress gradient. Here, we investigated microbial communities within an aquifer in Oak Ridge, Tennessee, USA, which is characterized by a broad spectrum of stressors, including extremely high levels of nitrate, heavy metals like cadmium and chromium, radionuclides such as uranium, and extremely low pH (< 3). Results: Both taxonomic and phylogenetic α-diversities were reduced in the most impacted wells, while the decline in functional α-diversity was modest and statistically insignificant, indicating a more robust buffering capacity to environmental stress. Differences in functional gene composition (i.e., functional β-diversity) were pronounced in highly contaminated wells, while convergent functional gene composition was observed in uncontaminated wells. The relative abundances of most carbon degradation genes were decreased in contaminated wells, but genes associated with denitrification, adenylylsulfate reduction, and sulfite reduction were increased. Compared to taxonomic and phylogenetic compositions, environmental variables played a more significant role in shaping functional gene composition, suggesting that niche selection could be more closely related to microbial functionality than taxonomy. Conclusions: Overall, we demonstrated that despite a reduced taxonomic α-diversity, microbial communities under stress maintained functionality underpinned by environmental selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Convergent reductive evolution and host adaptation in Mycoavidus bacterial endosymbionts of Mortierellaceae fungi

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Codon Optimization Improves the Prediction of Xylose Metabolism from Gene Content in Budding Yeasts

Xylose is the second most abundant monomeric sugar in plant biomass. Consequently, xylose catabolism is an ecologically important trait for saprotrophic organisms, as well as a fundamentally important trait for industries that hope to convert plant mass to renewable fuels and other bioproducts using microbial metabolism. Although common across fungi, xylose catabolism is rare within Saccharomycotina, the subphylum that contains most industrially relevant fermentative yeast species. The genomes of several yeasts unable to consume xylose have been previously reported to contain the full set of genes in the XYL pathway, suggesting the absence of a gene–trait correlation for xylose metabolism. Here, we measured growth on xylose and systematically identified XYL pathway orthologs across the genomes of 332 budding yeast species. Although the XYL pathway coevolved with xylose metabolism, we found that pathway presence only predicted xylose catabolism about half of the time, demonstrating that a complete XYL pathway is necessary, but not sufficient, for xylose catabolism. We also found that XYL1 copy number was positively correlated, after phylogenetic correction, with xylose utilization. We then quantified codon usage bias of XYL genes and found that XYL3 codon optimization was significantly higher, after phylogenetic correction, in species able to consume xylose. Finally, we showed that codon optimization of XYL2 was positively correlated, after phylogenetic correction, with growth rates in xylose medium. We conclude that gene content alone is a weak predictor of xylose metabolism and that using codon optimization enhances the prediction of xylose metabolism from yeast genome sequence data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Substitution Models of Protein Evolution with Selection on Enzymatic Activity

Abstract Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme–substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

Ferreiro, David↗