Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolutionary biology”

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.

750 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evolutionary Diversification of Host-Targeted Bartonella Effectors Proteins Derived from a Conserved FicTA Toxin-Antitoxin Module

Proteins containing a FIC domain catalyze AMPylation and other post-translational modifications (PTMs). In bacteria, they are typically part of FicTA toxin-antitoxin modules that control conserved biochemical processes such as topoisomerase activity, but they have also repeatedly diversified into host-targeted virulence factors. Among these, Bartonella effector proteins (Beps) comprise a particularly diverse ensemble of FIC domains that subvert various host cellular functions. However, no comprehensive comparative analysis has been performed to infer molecular mechanisms underlying the biochemical and functional diversification of FIC domains in the vast Bep family. Here, we used X-ray crystallography, structural modelling, and phylogenetic analyses to unravel the expansion and diversification of Bep repertoires that evolved in parallel in three Bartonella lineages from a single ancestral FicTA toxin-antitoxin module. Our analysis is based on 99 non-redundant Bep sequences and nine crystal structures. Inferred from the conservation of the FIC signature motif that comprises the catalytic histidine and residues involved in substrate binding, about half of them represent AMP transferases. A quarter of Beps show a glutamate in a strategic position in the putative substrate binding pocket that would interfere with triphosphate-nucleotide binding but may allow binding of an AMPylated target for deAMPylation or another substrate to catalyze a distinct PTM. The β-hairpin flap that registers the modifiable target segment to the active site exhibits remarkable structural variability. The corresponding sequences form few well-defined groups that may recognize distinct target proteins. The binding of Beps to promiscuous FicA antitoxins is well conserved, indicating a role of the antitoxin to inhibit enzymatic activity or to serve as a chaperone for the FIC domain before translocation of the Bep into host cells. Taken together, our analysis indicates a remarkable functional plasticity of Beps that is mostly brought about by structural changes in the substrate pocket and the target dock. These findings may guide future structure–function analyses of the highly versatile FIC domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Improving deep learning-based protein distance prediction in CASP14

Abstract Motivation Accurate prediction of residue–residue distances is important for protein structure prediction. We developed several protein distance predictors based on a deep learning distance prediction method and blindly tested them in the 14th Critical Assessment of Protein Structure Prediction (CASP14). The prediction method uses deep residual neural networks with the channel-wise attention mechanism to classify the distance between every two residues into multiple distance intervals. The input features for the deep learning method include co-evolutionary features as well as other sequence-based features derived from multiple sequence alignments (MSAs). Three alignment methods are used with multiple protein sequence/profile databases to generate MSAs for input feature generation. Based on different configurations and training strategies of the deep learning method, five MULTICOM distance predictors were created to participate in the CASP14 experiment. Results Benchmarked on 37 hard CASP14 domains, the best performing MULTICOM predictor is ranked 5th out of 30 automated CASP14 distance prediction servers in terms of precision of top L/5 long-range contact predictions [i.e. classifying distances between two residues into two categories: in contact (<8 Angstrom) and not in contact otherwise] and performs better than the best CASP13 distance prediction method. The best performing MULTICOM predictor is also ranked 6th among automated server predictors in classifying inter-residue distances into 10 distance intervals defined by CASP14 according to the precision of distance classification. The results show that the quality and depth of MSAs depend on alignment methods and sequence databases and have a significant impact on the accuracy of distance prediction. Using larger training datasets and multiple complementary features improves prediction accuracy. However, the number of effective sequences in MSAs is only a weak indicator of the quality of MSAs and the accuracy of predicted distance maps. In contrast, there is a strong correlation between the accuracy of contact/distance predictions and the average probability of the predicted contacts, which can therefore be more effectively used to estimate the confidence of distance predictions and select predicted distance maps. Availability and implementation The software package, source code and data of DeepDist2 are freely available at https://github.com/multicom-toolbox/deepdist and https://zenodo.org/record/4712084#.YIIM13VKhQM. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Transcriptome Dynamics of the Inflorescence in Reciprocally Formed Allopolyploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary mechanism and is prevalent among land plants. Most polyploid species examined have multiple origins, which provide genetic diversity and may enhance the success of polyploids. In some polyploids, recurrent origins can result from reciprocal crosses between the same diploid progenitors. Although great progress has been made in understanding the genetic consequences of polyploidy, the genetic implications of reciprocal polyploidization remain poorly understood, especially in natural polyploids. Tragopogon (Asteraceae) has become an evolutionary model system for studies of recent and recurrent polyploidy. Allotetraploid T. miscellus has formed reciprocally in nature with resultant distinctive floral and inflorescence morphologies (i.e., short- vs. long-liguled forms). In this study, we performed comparative inflorescence transcriptome analyses of reciprocally formed T. miscellus and its diploid parents, T. dubius and T. pratensis. In both forms of T. miscellus, homeolog expression of ~70% of the loci showed vertical transmission of the parental expression patterns (i.e., parental legacy), and ~20% of the loci showed biased homeolog expression, which was unbalanced toward T. pratensis. However, 17.9% of orthologous pairs showed different homeolog expression patterns between the two forms of T. miscellus. No clear effect of cytonuclear interaction on biased expression of the maternal homeolog was found. In terms of the total expression level of the homeologs studied, 22.6% and 16.2% of the loci displayed non-additive expression in short- and long-liguled T. miscellus, respectively. Unbalanced expression level dominance toward T. pratensis was observed in both forms of T. miscellus. Significantly, genes annotated as being involved in pectin catabolic processes were highly expressed in long-liguled T. miscellus relative to the short-liguled form, and the majority of these differentially expressed genes were transgressively down-regulated in short-liguled T. miscellus. Given the known role of these genes in cell expansion, they may play a role in the differing floral and inflorescence morphologies of the two forms. In summary, the overall inflorescence transcriptome profiles are highly similar between reciprocal origins of T. miscellus. However, the dynamic homeolog-specific expression and non-additive expression patterns observed in T. miscellus emphasize the importance of reciprocal origins in promoting the genetic diversity of polyploids.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO 2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Adaptive modification of antiviral defense systems in microbial community under Cr-induced stress

Background The prokaryotic antiviral defense systems are crucial for mediating prokaryote-virus interactions that influence microbiome functioning and evolutionary dynamics. Despite the prevalence and significance of prokaryotic antiviral defense systems, their responses to abiotic stress and ecological consequences remain poorly understood in soil ecosystems. We established microcosm systems with varying concentrations of hexavalent chromium (Cr(VI)) to investigate the adaptive modifications of prokaryotic antiviral defense systems under abiotic stress. Results Utilizing hybrid metagenomic assembly with long-read and short-read sequencing, we discovered that anti- viral defense systems were more diverse and prevalent in heavily polluted soils, which was corroborated by meta-analyses of public datasets from various heavy metal-contaminated sites. As the Cr(VI) concentration increased, prokaryotes with defense systems favoring prokaryote-virus mutualism gradually supplanted those with defense systems incurring high adaptive costs. Additionally, as Cr(VI) concentrations increased, enriched antiviral defense systems exhibited synchronization with microbial heavy metal resistance genes. Furthermore, the proportion of antiviral defense systems carried by mobile genetic elements (MGEs), including plasmids and viruses, increased by approximately 43% and 39%, respectively, with rising Cr concentrations. This trend is conducive to strengthening the dissemination and sharing of defense resources within microbial communities. Conclusions Overall, our study reveals the adaptive modification of prokaryotic antiviral defense systems in soil ecosystems under abiotic stress, as well as their positive contributions to establishing prokaryote-virus mutualism and the evolution of microbial heavy metal resistance. These findings advance our understanding of microbial adaptation in stressful environments and may inspire novel approaches for microbiome manipulation and bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Comparing Life Histories of the Shortest-Lived Turtle Known (Chicken Turtles, Deirochelys reticularia ) with Long-Lived Blanding's Turtles ( Emydoidea blandingii )

Evolutionary theories predict major differences in life-history trait values of long- and short-lived organisms. Such comparisons have not been possible for chelonians because no shortlived turtle was known until research revealed that chicken turtles (Deirochelys reticularia; DR) have a maximum longevity of 21 yrs. Life-history trait values of DR females are 1) age at maturity of females = 5–6 yrs; 2) clutches per season = 1.6; 3) annual fecundity = 6–8 female eggs per female; 4) average juvenile survivorship from age 1 to maturity = 0.60; and 5) low average annual adult survivorship = 0.66. We compared DR with the very long-lived Blanding’s turtles (Emydoidea blandingii; EB) in Michigan. Over 14 yrs with no mortality (the minimum age at maturity of EB), the maximum potential fecundity produced by a single female embryo and her mature female offspring was 5 female eggs for EB and 1040 eggs for DR. Comparisons of life table output for approximately stable populations of DR and EB resulted in cohort generation times of 7 and 37 yrs, respectively. The life-history prediction that short-lived organisms should produce smaller offspring was not supported. Average wet mass of eggs is 10 g (8.4–11.3 g) for DR and 12 g (10–14 g) for EB; and average wet mass of hatchlings is 7.3 g (6–9 g) for DR and 9.3 g (6–13 g) for EB. Both differences are smaller than expected based on the difference in longevity. Here, short-lived female DR have an unusual tactic of investing in high fecundity and making substantial body sizespecific investment in large eggs, which may reflect why juvenile survivorship had greater influence on population change rates than did adult survivorship. In contrast, adult survivorship had the greatest influence on population change rates of EB. Comparison of cohorts of 1000 female DR and EB hatchlings highlights the differences in life histories of short- and long-lived turtles: all DR would be dead by the time the last female EB had reached maturity at 21 yrs of age.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Long-reads↗
DOE OSTI2025

Convergent expansions of keystone gene families drive metabolic innovation in Saccharomycotina yeasts

Many remarkable phenotypes have repeatedly occurred across vast evolutionary distances. When convergent traits emerge on the tree of life, they are sometimes driven by the same underlying gene families, while other times, many different gene families are involved. Conversely, a gene family may be repeatedly recruited for a single trait or many different traits. To understand the general rules governing convergence at both genomic and phenotypic levels, we systematically tested associations between 56 binary metabolic traits and gene count in 14,785 gene families from 993 Saccharomycotina yeasts. Using a recently developed phylogenetic approach that reduces spurious correlations, we found that gene family expansion and contraction were significantly linked to trait gain and loss in 45/56 (80%) traits. While 595/739 (81%) significant gene families were associated with only one trait, we also identified several “keystone” gene families that were significantly associated with up to 13/56 (23%) of all traits. Strikingly, most of these families are known to encode metabolic enzymes and transporters, including all members of the industrially relevant MAL tose fermentation loci in the baker’s yeast Saccharomyces cerevisiae. These results indicate that convergent evolution on the gene family level may be more widespread across deeper timescales than previously believed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Incorporating the Acclimation of Photosynthesis and Leaf Respiration in the Noah‐MP Land Surface Model: Model Development and Evaluation

Abstract Realistic simulation of leaf photosynthetic and respiratory processes is needed for accurate prediction of the global carbon cycle. These two processes systematically acclimate to long‐term environmental changes by adjusting photosynthetic and respiratory traits (e.g., the maximum photosynthetic capacity at 25°C ( V cmax,25 ) and the leaf respiration rate at 25°C ( R 25 )) following increasingly well‐understood principles. While some land surface models (LSMs) now account for thermal acclimation, they do so by assigning empirical parameterizations for individual plant functional types (PFTs). Here, we have implemented an Eco‐Evolutionary Optimality (EEO)‐based scheme to represent the universal acclimation of photosynthesis and leaf respiration to multiple environmental effects, and that therefore requires no PFT‐specific parameterizations, in a standard version of the widely used LSM, Noah MP. We evaluated model performance with plant trait data from a 5‐year experiment and extensive global field measurements, and carbon flux measurements from FLUXNET2015. We show that observed R 25 and V cmax,25 vary substantially both temporally and spatially within the same PFT ( C.V. >20%). Our EEO‐based scheme captures 62% of the temporal and 70% of the spatial variations in V cmax,25 (73% and 54% of the variations in R 25 ). The standard scheme underestimates gross primary production by 10% versus 2% for the EEO‐based scheme and generates a larger spread in r (correlation coefficient) across flux sites (0.79 ± 0.16 vs. 0.84 ± 0.1, mean ± S.D.). The standard scheme greatly overestimates canopy respiration (bias: ∼200% vs. 8% for the EEO scheme), resulting in less CO 2 uptake by terrestrial ecosystems. Our approach thus simulates climate‐carbon coupling more realistically, with fewer parameters.

Ren, Yanghang [Department of Earth System Science ↗
DOE OSTI2023

Transcriptional Control of hgcAB by an ArsR-Like Regulator in Pseudodesulfovibrio mercurii ND132

The hgcAB gene pair encodes mercury (Hg) methylation capability in a diverse group of microorganisms, but its evolution and transcriptional regulation remain unknown. Working from the possibility that the evolutionary function of HgcAB may not be Hg methylation, we test a possible link to arsenic resistance. Using model Hg methylator Pseudodesulfovibrio mercurii ND132, we evaluated transcriptional control of hgcAB by a putative ArsR encoded upstream and cotranscribed with hgcAB. This regulator shares homology with ArsR repressors of arsenic resistance and S-adenosylhomocysteine (SAH)- responsive regulators of methionine biosynthesis but is distinct from other ArsR/SahR proteins in P. mercurii. Using quantitative PCR (qPCR) and RNA sequencing (RNA-seq) transcriptome analyses, we confirmed this ArsR regulates hgcAB transcription and is responsive to arsenic and SAH. Additionally, RNA-seq indicated a possible link between hgcAB activity and arsenic transformations, with significant upregulation of other ArsRregulated arsenic resistance operons alongside hgcAB. Interestingly, wild-type ND132 was less sensitive to As(V) (but not As(III)) than an hgcAB knockout strain, supporting the idea that hgcAB may be linked to arsenic resistance. Arsenic significantly impacted rates of Hg methylation by ND132; however, responses varied with culture conditions. Differences in growth and metabolic activity did not account for arsenic impacts on methylation. While arsenic significantly increased hgcAB expression, hgcAB gene and transcript abundance was not a good predictor of Hg methylation rates. Taken together, these results support the idea that Hg and As cycling are linked in P. mercurii ND132. Our results may hold clues to the evolution of hgcAB and the controls on Hg methylation in nature. This work reveals a link between microbial mercury methylation and arsenic resistance and may hold clues to the evolution of mercury methylation genes (hgcAB). Microbes with hgcAB produce methylmercury, a strong neurotoxin that readily accumulates in the food web. This study addresses a critical gap in our understanding about the environmental factors that control hgcAB expression. We show that hgcAB expression is controlled by an ArsR-like regulator responsive to both arsenic and S-adenosylhomocysteine in our model organism, Pseudodesulfovibrio mercurii ND132. Exposure to arsenic also significantly impacted Pseudodesulfovibrio mercurii ND132 mercury methylation rates. However, expression of hgcAB was not always a good predictor of Hg methylation rates, highlighting the roles of Hg bioavailability and other biochemical mechanisms in methylmercury production. This study improves our understanding of the controls on hgcAB expression, which is needed to better predict environmental methylmercury production

RNA-seq↗
DOE OSTI2024

Catabolic pathway acquisition by rhizosphere bacteria readily enables growth with a root exudate component but does not affect root colonization

Horizontal gene transfer (HGT) is a fundamental evolutionary process that plays a key role in bacterial evolution. The likelihood of a successful transfer event is expected to depend on the precise balance of costs and benefits resulting from pathway acquisition. Most experimental analyses of HGT have focused on phenotypes that have large fitness benefits under appropriate selective conditions, such as antibiotic resistance. However, many examples of HGT involve phenotypes that are predicted to provide smaller benefits, such as the ability to catabolize additional carbon sources. We have experimentally simulated the consequences of one such HGT event in the laboratory, studying the effects of transferring a pathway for catabolism of the plant-derived aromatic compound salicyl alcohol between rhizosphere isolates from the Pseudomonas genus. We find that pathway acquisition enables rapid catabolism of salicyl alcohol with only minor disruptions to the existing metabolic and regulatory networks of the new host. However, this new catabolic potential does not confer a measurable fitness advantage during competitive growth in the rhizosphere. We conclude that the phenotype of salicyl alcohol catabolism is readily transferable but is selectively neutral under environmentally relevant conditions. We propose that this condition is common and that HGT of many pathways will be self-limiting because the selective benefits are small.

59 BASIC BIOLOGICAL SCIENCES↗