Search NASA⌕ Search

SEARCH · Search NASA

Results for “genome evolution”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Leveraging computational genomics to understand the molecular basis of metal homeostasis

Genome-based data is helping to reveal the diverse strategies plants and algae use to maintain metal homeostasis. In addition to acquisition, distribution and storage of metals, acclimating to feast or famine can involve a wealth of genes that we are just now starting to understand. The fast-paced acquisition of genome-based data, however, is far outpacing our ability to experimentally characterize protein function. Computational genomic approaches are needed to fill the gap between what is known and unknown. To avoid misconstruing bioinformatically derived data, which is the root cause of the inaccurate functional annotations that plague databases, functional inferences from diverse sources and contextualization of that evidence with a robust understanding of protein family evolution is needed. Phylogenomic- and comparative-genomic-based studies can aid in the interpretation of experimental data or provide a spark for the discovery of a new function. These analyses not only lead to novel insight into a target protein's function but can generate thought-provoking insights across protein families.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker’s yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Abstract Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole‐genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri , a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency‐based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%–75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests c erk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.

Chen, Zhongqi↗
DOE OSTI2022

Cross-Kingdom Comparative Genomics of Aromatic Catabolic Pathways in Fungi and Bacteria

The overall goal of this project is to test the hypothesis that white-rot fungi can simultaneously depolymerize lignin extracellularly and catabolize depolymerization products intracellularly as carbon and energy sources. The results from this project will lead to improved understanding of lignin utilization by white-rot fungi, and enable identification of promising fungal strains for lignin catabolism and valorization. As part of this effort, we have conducted a comparative genomic study to investigate a diversity of aromatic catabolic pathways in the fungal and bacterial kingdoms. Furthermore, we have combined genomic and phylogenetic approaches to decipher the evolution of certain aromatic catabolic enzymes and discover new classes of enzymes, which represents a foundation for future biochemical and molecular genetic studies.

aromatic compounds↗
DOE OSTI2025

Comparative genomics reveals the high diversity and adaptation strategies of Polaromonas from polar environments

Abstract Background Bacteria from the genus Polaromonas are dominant phylotypes found in a variety of low-temperature environments in polar regions. The diversity and biogeographic distribution of Polaromonas have been largely expanded on the basis of 16 S rRNA gene amplicon sequencing. However, the evolution and cold adaptation mechanisms of Polaromonas from polar regions are poorly understood at the genomic level. Results A total of 202 genomes of the genus Polaromonas were analyzed, and 121 different species were delineated on the basis of average nucleotide identity (ANI) and phylogenomic placements. Remarkably, 8 genomes recovered from polar environments clustered into a separate clade (‘polar group’ hereafter). The genome size, coding density and coding sequences (CDSs) of the polar group were significantly different from those of other nonpolar Polaromonas . Furthermore, the enrichment of genes involved in carbohydrate and peptide metabolism was evident in the polar group. In addition, genes encoding proteins related to betaine synthesis and transport were increased in the genomes from the polar group. Phylogenomic analysis revealed that two different evolutionary scenarios may explain the adaptation of Polaromonas to cold environments in polar regions. Conclusions The global distribution of the genus Polaromonas highlights its strong adaptability in both polar and nonpolar environments. Species delineation significantly expands our understanding of the diversity of the Polaromonas genus on a global scale. In this study, a polar-specific clade was found, which may represent a specific ecotype well adapted to polar environments. Collectively, genomic insight into the metabolic diversity, evolution and adaptation of the genus Polaromonas at the genome level provides a genetic basis for understanding the potential response mechanisms of Polaromonas to global warming in polar regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tracking ebolavirus genomic drift with a resequencing microarray

Filoviruses are emerging pathogens that cause acute fever with high fatality rate and present a global public health threat. During the 2013–2016 Ebola virus outbreak, genome sequencing allowed the study of virus evolution, mutations affecting pathogenicity and infectivity, and tracing the viral spread. In 2018, early sequence identification of the Ebolavirus as EBOV in the Democratic Republic of the Congo supported the use of an Ebola virus vaccine. However, field-deployable sequencing methods are needed to enable a rapid public health response. Resequencing microarrays (RMA) are a targeted method to obtain genomic sequence on clinical specimens rapidly, and sensitively, overcoming the need for extensive bioinformatic analysis. This study presents the design and initial evaluation of an ebolavirus resequencing microarray (Ebolavirus-RMA) system for sequencing the major genomic regions of four Ebolaviruses that cause disease in humans. The design of the Ebolavirus-RMA system is described and evaluated by sequencing repository samples of three Ebolaviruses and two EBOV variants. The ability of the system to identify genetic drift in a replicating virus was achieved by sequencing the ebolavirus glycoprotein gene in a recombinant virus cultured under pressure from a neutralizing antibody. Comparison of the Ebolavirus-RMA results to the Genbank database sequence file with the accession number given for the source RNA and Ebolavirus-RMA results compared to Next Generation Sequence results of the same RNA samples showed up to 99% agreement.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genome organization and botanical diversity

Abstract The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes—both in content and in dynamics—has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE ‘nests’, or dense aggregations of TEs. Here, some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Biogeographic patterns in populations of marine Pseudoalteromonas atlantica isolates

Intra-specific genomic diversity is well documented in microbes. The question, however, remains whether natural selection or neutral evolution is the major contributor to this diversity. We undertook this study to estimate genomic diversity in Pseudoalteromonas atlantica populations and whether the diversity, if present, could be attributed to environmental factors or distance effects. Here, we isolated and sequenced twenty-three strains of P. atlantica from three geographically distant deep marine basins and performed comparative genomic analyses to study the genomic diversity of populations among these basins. Average nucleotide identity followed a strictly geographical pattern. In two out of three locations, the strains within the location exhibited >99.5% identity, whereas, among locations, the strains showed <98.11% identity. Phylogenetic and pan-genome analysis also reflected the biogeographical separation of the strains. Strains from the same location shared many accessory genes and clustered closely on the phylogenetic tree. Phenotypic diversity between populations was studied in ten out of twenty-three strains testing carbon and nitrogen source utilization and osmotolerance. A genetic basis for phenotypic diversity could be established in most cases but was apparently not influenced by local environmental conditions. Our study suggests that neutral evolution may have a substantial role in the biodiversity of P. atlantica.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

The landscape of regulatory element evolution in a C4 perennial grass

Gene regulatory evolution is a well-known source of phenotypic diversity and adaptive evolution. Although cis-regulatory elements (CREs) play a vital role in gene expression evolution, the molecular evolution of CREs remains mostly unknown due to the difficulty in identifying and characterizing these functional elements. Comparative genomic analyses of noncoding DNA can be leveraged to identify conserved noncoding sequences (CNS), many of which may harbor functional CREs conserved by purifying selection. However, purely computational inference of CREs from putative CNS can be erroneous due to the complex genomic architecture in plants. One promising experimental approach to identify CREs is by profiling accessible chromatin regions (ACRs) that are often associated with the location of CREs. In this study, we use comparative genomics along with the profiling of ACRs to study the molecular evolution of putative functional noncoding regulatory regions in Panicoid grasses. We identified sets of CNS that varied in relationship to the degree of evolutionary divergence among the studied taxa, including identifying core-Panicoid-CNS. We augmented this analysis by profiling ACRs in Panicum hallii ecotypes using ATAC-seq. ACRs had low SNP density at the summit, harbored a high frequency of core-Panicoid-CNS, and were enriched with expression QTL. These data help to annotate the P. hallii genome for putative functional elements and suggest that a large proportion of these ACRs are evolving under purifying selection. Turnover in CNS and ACR between ecotypes of P. hallii identifies a small set of putatively divergent CREs that may underlie differences in gene regulation between genotypes from inland and coastal habitats. In summary, we profiled ACRs in Panicoid grasses and integrated this data with our putative CNS prediction framework, which provides unique insight into patterns of polymorphism and divergence in CREs in C4 perennial grasses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Postglacial migration shaped the genomic diversity and global distribution of the wild ancestor of lager-brewing hybrids

The wild, cold-adapted parent of hybrid lager-brewing yeasts, Saccharomyces eubayanus, has a complex and understudied natural history. The exploration of this diversity can be used both to develop new brewing applications and to enlighten our understanding of the dynamics of yeast evolution in the wild. Here, we integrate whole genome sequence and phenotypic data of 200 S. eubayanus strains, the largest collection known to date. S. eubayanus has a multilayered population structure, consisting of two major populations that are further structured into six subpopulations. Four of these subpopulations are found exclusively in the Patagonian region of South America; one is found predominantly in Patagonia and sparsely in Oceania and North America; and one is specific to the Holarctic ecozone. Plant host associations differed between subpopulations and between S. eubayanus and its sister species, Saccharomyces uvarum. S. eubayanus is most abundant and genetically diverse in northern Patagonia, where some locations harbor more genetic diversity than is found outside of South America, suggesting that northern Patagonia east of the Andes was a glacial refugium for this species. All but one subpopulation shows isolation-by-distance, and gene flow between subpopulations is low. However, there are strong signals of ancient and recent outcrossing, including two admixed lineages, one that is sympatric with and one that is mostly isolated from its parental populations. Using our extensive biogeographical data, we build a robust model that predicts all known and a handful of additional regions of the globe that are climatically suitable for S. eubayanus, including Europe where host accessibility and competitive exclusion by other Saccharomyces species may explain its continued elusiveness. We conclude that this industrially relevant species has rich natural diversity with many factors contributing to its complex distribution and natural history.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Genomic analysis of Klebsiella aerogenes circulating in New Mexico

Klebsiella aerogenes is an opportunistic pathogen and a growing cause of healthcare-associated infections, characterized by multidrug resistance and the emergence of global high-risk clones. However, regional genomic surveillance data remain limited. Here, we sought to characterize the population structure, transmission dynamics and resistance mechanisms of clinical K. aerogenes in Albuquerque, New Mexico. We sequenced 177 clinical isolates collected between 2021 and 2023. We also developed a novel, species-specific PopPUNK database to facilitate rapid, high-resolution typing. The New Mexico K. aerogenes population was diverse but dominated by two global pandemic lineages, ST93 (47.5%) and ST4 (7.9%), which were significantly enriched for the virulence factors yersiniabactin and colibactin. Genomic evidence for recent local transmission was rare, with only four putative transmission pairs identified. The resistome was characterized by intrinsic and adaptive mutations. Nearly all isolates possessed gyrA mutations associated with decreased fluoroquinolone susceptibility. Mutations in the AmpC regulator AmpD and the outer membrane porin Omp36 were common, particularly within the dominant ST93 lineage. These mutations have been associated with increased AmpC-mediated carbapenem resistance. Our findings underscore the critical importance of genomic surveillance to monitor the transmission and evolution of adaptive resistance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Plastome evolution in annual Brachypodium species reveals widespread heteroplasmy and chloroplast capture, lineage-specific codon usage bias, and low positive selection

Comparative genomics and plastome phylogenomics have advanced significantly in recent years, highlighting the diversity, possible admixture, and non-neutral evolution of the predominantly considered non-recombinant chloroplast genomes in angiosperms. The grass genus Brachypodium serves as a powerful model for studying evolutionary processes in monocots. We analyzed 287 plastomes across the native circum-Mediterranean range of the three annual Brachypodium species ( B. distachyon, B. stacei, B. hybridum ), focusing on their structural variation, selection patterns and phylogenomic relationships. Our analyses confirmed the differentiation of the S and D plastomes, inherited respectively from the diploid progenitor species B. stacei and B. distachyon . We identified novel structural rearrangements and indels, and unique repeat motifs, along with widespread heteroplasmy, particularly in ancestral B. hybridum -D plastotypes. SNP diversity varied among plastotypes, reflecting population dynamics and evolutionary histories, with B. hybridum -D plastotypes showing the highest normalized diversity and B. hybridum -S the lowest. Positive selection was detected in 29 plastid genes by Tajima’s neutrality test, and in nine genes by site and branch-site evolutionary models, including matK, ndhF, rbcL, and rpoC2. Phylogenomic analyses revealed well-supported clades corresponding to the S and D plastome lineages, with frequent chloroplast capture events and long-distance dispersals shaping their evolutionary trajectories.

allopolyploidy↗
DOE OSTI2023

Large-scale genomic analyses with machine learning uncover predictive patterns associated with fungal phytopathogenic lifestyles and traits

Abstract Invasive plant pathogenic fungi have a global impact, with devastating economic and environmental effects on crops and forests. Biosurveillance, a critical component of threat mitigation, requires risk prediction based on fungal lifestyles and traits. Recent studies have revealed distinct genomic patterns associated with specific groups of plant pathogenic fungi. We sought to establish whether these phytopathogenic genomic patterns hold across diverse taxonomic and ecological groups from the Ascomycota and Basidiomycota, and furthermore, if those patterns can be used in a predictive capacity for biosurveillance. Using a supervised machine learning approach that integrates phylogenetic and genomic data, we analyzed 387 fungal genomes to test a proof-of-concept for the use of genomic signatures in predicting fungal phytopathogenic lifestyles and traits during biosurveillance activities. Our machine learning feature sets were derived from genome annotation data of carbohydrate-active enzymes (CAZymes), peptidases, secondary metabolite clusters (SMCs), transporters, and transcription factors. We found that machine learning could successfully predict fungal lifestyles and traits across taxonomic groups, with the best predictive performance coming from feature sets comprising CAZyme, peptidase, and SMC data. While phylogeny was an important component in most predictions, the inclusion of genomic data improved prediction performance for every lifestyle and trait tested. Plant pathogenicity was one of the best-predicted traits, showing the promise of predictive genomics for biosurveillance applications. Furthermore, our machine learning approach revealed expansions in the number of genes from specific CAZyme and peptidase families in the genomes of plant pathogens compared to non-phytopathogenic genomes (saprotrophs, endo- and ectomycorrhizal fungi). Such genomic feature profiles give insight into the evolution of fungal phytopathogenicity and could be useful to predict the risks of unknown fungi in future biosurveillance activities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High Accuracy Indicators of Androgen Suppression Therapy Failure for Prostate Cancer—A Modeling Study

Prostate cancer is a serious public health concern in the United States. The primary obstacle to effective long-term management for prostate cancer patients is the eventual development of treatment resistance. Due to the uniquely chaotic nature of the neoplastic genome, it is difficult to determine the evolution of tumor composition over the course of treatment. Hence, a drug is often applied continuously past the point of effectiveness, thereby losing any potential treatment combination with that drug permanently to resistance. If a clinician is aware of the timing of resistance to a particular drug, then they may have a crucial opportunity to adjust the treatment to retain the drug’s usefulness in a potential treatment combination or strategy. In this study, we investigate new methods of predicting treatment failure due to treatment resistance using a novel mechanistic model built on an evolutionary interpretation of Droop cell quota theory. We analyze our proposed methods using patient PSA and androgen data from a clinical trial of intermittent treatment with androgen deprivation therapy. Our results produce two indicators of treatment failure. The first indicator, proposed from the evolutionary nature of the cancer population, is calculated using our mathematical model with a predictive accuracy of 87.3% (sensitivity: 96.1%, specificity: 65%). The second indicator, conjectured from the implication of the first indicator, is calculated directly from serum androgen and PSA data with a predictive accuracy of 88.7% (sensitivity: 90.2%, specificity: 85%). Our results demonstrate the potential and feasibility of using an evolutionary tumor dynamics model in combination with the appropriate data to aid in the adaptive management of prostate cancer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes

Abstract Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Engineering improved ethylene production: Leveraging systems biology and adaptive laboratory evolution

Ethylene is a small hydrocarbon gas widely used in the chemical industry. Annual worldwide production currently exceeds 150 million tons, producing considerable amounts of CO 2 contributing to climate change. The need for a sustainable alternative is therefore imperative. Ethylene is natively produced by several different microorganisms, including Pseudomonas syringae pv. phaseolicola via a process catalyzed by the ethylene forming enzyme (EFE), subsequent heterologous expression of EFE has led to ethylene production in non-native bacterial hosts including E. coli and cyanobacteria. However, solubility of EFE and substrate availability remain rate limiting steps in biological ethylene production. We employed a combination of genome scale metabolic modelling, continuous fermentation, and protein evolution to enable the accelerated development of a high efficiency ethylene producing E. coli strain, yielding a 49-fold increase in production, the most significant improvement reported to date. Furthermore, we have clearly demonstrated that this increased yield resulted from metabolic adaptations that were uniquely linked to the EFE enzyme (WT vs mutant). Our findings provide a novel solution to deregulate metabolic bottlenecks in key pathways, which can be readily applied to address other engineering challenges.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Functional and genetic markers of niche partitioning among enigmatic members of the human oral microbiome

Abstract Introduction Microbial residents of the human oral cavity have long been a major focus of microbiology due to their influence on host health and intriguing patterns of site specificity amidst the lack of dispersal limitation. However, the determinants of niche partitioning in this habitat are yet to be fully understood, especially among taxa that belong to recently discovered branches of microbial life. Results Here, we assemble metagenomes from tongue and dental plaque samples from multiple individuals and reconstruct 790 non-redundant genomes, 43 of which resolve to TM7, a member of the Candidate Phyla Radiation, forming six monophyletic clades that distinctly associate with either plaque or tongue. Both pangenomic and phylogenomic analyses group tongue-specific clades with other host-associated TM7 genomes. In contrast, plaque-specific TM7 group with environmental TM7 genomes. Besides offering deeper insights into the ecology, evolution, and mobilome of cryptic members of the oral microbiome, our study reveals an intriguing resemblance between dental plaque and non-host environments indicated by the TM7 evolution, suggesting that plaque may have served as a stepping stone for environmental microbes to adapt to host environments for some clades of microbes. Additionally, we report that prophages are widespread among oral-associated TM7, while absent from environmental TM7, suggesting that prophages may have played a role in adaptation of TM7 to the host environment. Conclusions Our data illuminate niche partitioning of enigmatic members of the oral cavity, including TM7, SR1, and GN02, and provide genomes for poorly characterized yet prevalent members of this biome, such as uncultivated Flavobacteriaceae.

60 APPLIED LIFE SCIENCES↗