Search NASA⌕ Search

SEARCH · Search NASA

Results for “transcriptomic profiling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Zn deficiency disrupts Cu and S homeostasis in Chlamydomonas resulting in over accumulation of Cu and Cysteine

Growth of Chlamydomonas reinhardtii in zinc (Zn) limited medium leads to disruption of copper (Cu) homeostasis, resulting in up to 40-fold Cu over-accumulation relative to its typical Cu quota. Here, we show that Chlamydomonas controls its Cu quota by balancing Cu import and export, which is disrupted in a Zn deficient cell, thus establishing a mechanistic connection between Cu and Zn homeostasis. Transcriptomics, proteomics and elemental profiling revealed that Zn-limited Chlamydomonas cells up-regulate a subset of genes encoding “first responder” proteins involved in sulfur (S) assimilation and consequently accumulate more intracellular S, which is incorporated into L-cysteine, γ-glutamylcysteine, and homocysteine. Most prominently, in the absence of Zn, free L-cysteine is increased ~80-fold, corresponding to ~2.8 × 10 9 molecules/cell. Interestingly, classic S-containing metal binding ligands like glutathione and phytochelatins do not increase. X-ray fluorescence microscopy showed foci of S accumulation in Zn-limited cells that co-localize with Cu, phosphorus and calcium, consistent with Cu-thiol complexes in the acidocalcisome, the site of Cu(I) accumulation. Notably, cells that have been previously starved for Cu do not accumulate S or Cys, causally connecting cysteine synthesis with Cu accumulation. We suggest that cysteine is an in vivo Cu(I) ligand, perhaps ancestral, that buffers cytosolic Cu.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Climate adaptation in Populus trichocarpa : key adaptive loci identified for stomata and leaf traits

We investigated adaptive genetic variation in Populus trichocarpa, a potential biofuel feedstock crop, to better understand how physiological traits may influence tolerance to water limitation. Our study focused on leaf and stomatal traits, given their roles in plant–water relations and adaptation. Using a diversity panel of over 1300 genotypes, we measured 14 leaf and stomatal traits under control (well-watered) and drought (water-limited) conditions. We conducted genome-wide association studies (GWAS), climate association analyses, and transcriptome (RNA-seq) profiling to identify genetic loci associated with phenotypic variation and adaptation. Stomatal traits, including size and density, were correlated with the climate of origin, with genotypes from more arid regions tending to have smaller but denser stomata. GWAS identified multiple loci associated with trait variation, including a major-effect region on chromosome 10 linked to stomatal size and abaxial contact angle. This locus overlapped with a tandem array of 3-ketoacyl-CoA synthase (KCS) genes and showed strong allele–climate and gene expression associations. Our findings reveal genetic and phenotypic variation consistent with local adaptation and suggest that future climates may favor alleles associated with smaller stomata, particularly under increasing aridity. This work provides insights into climate adaptation and breeding strategies for resilience in perennial crops.

Populus trichocarpa↗
DOE OSTI2022

Unified epigenomic, transcriptomic, proteomic, and metabolomic taxonomy of Alzheimer’s disease progression and heterogeneity

Alzheimer’s disease (AD) is a heterogeneous disorder with abnormalities in multiple biological domains. In an advanced machine learning analysis of postmortem brain and in vivo blood multi-omics molecular data ( N = 1863), we integrated epigenomic, transcriptomic, proteomic, and metabolomic profiles into a multilevel biological AD taxonomy. We obtained a personalized multilevel molecular index of AD dementia progression that predicts severity of neuropathologies, and identified three robust molecular-based subtypes that explain much of the pathologic and clinical heterogeneity of AD. These subtypes present distinct patterns of alteration in DNA methylation, RNA, proteins, and metabolites, identifiable in the brain and subsequently in blood. In addition, the genetic variations that predispose to the various AD subtypes in brain predict distinct spatial patterns of alteration in cell types, suggesting a unique influence of each putative AD variant on neuropathological mechanisms. These observations support that an individually tailored multi-omics molecular taxonomy of AD may represent distinct targets for preventive or treatment interventions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Cis-regulatory evolution shapes facial diversity in birds and mammals

Birds and mammals exhibit extraordinary facial diversity, reflecting adaptations to distinct ecological niches and feeding strategies. While core face-building developmental programs are conserved and orchestrated by interactions between ectodermal organizers and the underlying mesenchyme, mechanisms driving facial shape variation remain poorly understood. Here, we integrate single-cell transcriptomic and chromatin accessibility profiling of mouse and chicken developing face to construct a comparative regulatory map. Although both ectodermal and mesenchymal populations display distinct regulatory features in each species, the mesenchyme exhibits markedly greater divergence, pointing to its central role in shaping facial morphology. We further reveal unexpected molecular complexity in the main face-shaping organizer, including a mouse-specific Shh/Wnt5a expression domain. At key morphogen loci (Bmp4, Fgf8, and Wnt5a), conserved and lineage-specific enhancers exhibit spatially restricted activity patterns that mirror divergent signaling domains. These findings demonstrate how cis-regulatory evolution modulates conserved developmental programs to generate morphological novelty, providing a valuable resource for studying vertebrate facial evolution.

Kyomen, Stella [Max Planck Institute for Evolution↗
NASA NTRS2007

Comparative MicroRNA Expression Patterns in Fibroblasts after Low and High Doses of Low-LET Radiation Exposure

Exposure to ionizing radiation causes DNA damage to cells, and provokes a plethora of cellular responses controlled by unique gene-directed signaling pathways. MicroRNAs (miRNAs) are small (~22-nucleotide), non-coding RNAs which functionally silence gene expression by either degrading the messages or inhibiting translation. Here we investigate radiation-dependent changes in these negative regulators by comparing the expression patterns of all 462 known human miRNAs in fibroblasts, after exposure to low (0.1 Gy) or high (2 Gy) doses of X-rays at 30 min, 2, 6 and 24 hrs post-treatment. The expression patterns of microRNAs after low and high doses of radiation show a similar qualitative down-regulation trend at early (0.5 hr) and late (24 hr) time points, with a quantitatively steeper slope following the 2 Gy exposures. Interestingly, an interruption of this downward trend is observed after the 2 Gy exposure, i.e. a significant up-regulation of microRNAs at 2 hrs, then reverting to the downward trend by 6 hrs; this interruption at the intermediate time point was not observed with the 0.1 Gy exposure. At the early time point (0.5 hr), candidate gene targets of selected down-regulated microRNAs, common to both 0.1 and 2 Gy exposures, were those functioning in chromatin remodeling. Candidate target genes of unique up-regulated microRNAs seen at a 2 hr intermediate time point, after the 2 Gy exposure only, are those involved in cell death signaling. Finally, putative target genes of down-regulated microRNAs seen at the late (24 hr) time point after either doses of radiation are those involved in the up-regulation of DNA repair, cell signaling and homeostasis. Thus we hypothesize that after radiation exposure, microRNAs acting as hub negative regulators for unique signaling pathways needed to be down-regulated so as to de-repress their target genes for the proper cellular responses, including DNA repair and cell maintenance. The unique microRNAs up-regulated at 2 hr after 2 Gy suggest the cellular response to functionally suppress the apoptotic death signaling reflex after exposure to high dose radiation. Further analyses with transcriptome and global proteomic profiling will validate the reciprocal expression of signature microRNAs selected in our radiation-exposed cells, and their candidate target gene families, and test our hypothesis that unique radiation-specific microRNAs are keys in governing signaling responses for damage control of this environmental hazard.

Maes, Olivier C.↗
NASA NTRS2016

Advances in Rodent Research Missions on the International Space Station

A research platform for rodent experiment on the ISS is a valuable tool for advancing biomedical research in space. Capabilities offered by the Rodent Research project developed at NASA Ames Research Center can support experiments of much longer duration on the ISS than previous experiments performed on the Space Shuttle. NASAs Rodent Research (RR)-1 mission was completed successfully and achieved a number of objectives, including validation of flight hardware, on-orbit operations, and science capabilities as well as support of a CASIS-sponsored experiment (Novartis) on muscle atrophy. Twenty C57BL6J adult female mice were launched on the Space-X (SpX) 4 Dragon vehicle, and thrived for up to 37 days in microgravity. Daily health checks of the mice were performed during the mission via downlinked video; all flight animals were healthy and displayed normal behavior, and higher levels of physical activity compared to ground controls. Behavioral analysis demonstrated that Flight and Ground Control mice exhibited the same range of behaviors, including eating, drinking, exploratory behavior, self- and allo-grooming, and social interactions indicative of healthy animals. The animals were euthanized on-orbit and select tissues were collected from some of the mice on orbit to assess the long-term sample storage capabilities of the ISS. In general, the data obtained from the flight mice were comparable to those from the three groups of control mice (baseline, vivarium and ground controls, which were housed in flight hardware), showing that the ISS has adequate capability to support long-duration rodent experiments. The team recovered 35 tissues from 40 RR-1 frozen carcasses, yielding 3300 aliquots of tissues to distribute to the scientific community in the U.S., including NASAs GeneLab project and scientists via Space Biology's Biospecimen Sharing Program Ames Life Science Data Archive. Tissues also were distributed to Russian research colleagues at the Institute for Biomedical Problems. The expression levels of select genes including albumin, catalase, GAPDH, HMGCoA Reductase, and IGF1 were determined using RNA isolated from the livers by qPCR and no significant differences by one factor ANOVA were found between flight and ground control groups. In addition, some of the liver samples were analyzed for transcriptomic, epigenomic and proteomic profiles; some of the data sets are now available to the scientific community through GeneLabs open science data website. A second long duration mission, Rodent Research-2 (RR-2) was completed on the ISS in 2015; 20 female C57BL6J mice were successfully maintained on the ISS for various durations, with the last group of 5 animals living on-orbit for 54 days. Furthermore, we continue to expand the ISSs capabilities by introducing new on-orbit technologies including blood collection and separation, bone densitometry scanning, muscle grip strength and anesthesia with recovery. In addition, series of ground-based verification testing to fly male mice and increase the total number of mice on-orbit from 20 to 40. Subsequent missions will provide the capability to return live mice from the ISS animals to evaluate recovery on Earth, further expanding operational and science capabilities of the RR project on the ISS.

Choi, S. Y.↗
DOE OSTI2026

Single-cell and spatial omics in plants: from cellular atlases to regulatory mechanisms

Single-cell RNA sequencing (scRNA-seq) has transformed transcriptomic studies by enabling gene expression profiling at the resolution of individual cells within and across a broad range of tissue types, revealing cellular heterogeneity that is obscured in bulk tissue transcriptomes. Over the past decade, improvements in microfluidics and library preparation have drastically increased throughput, allowing tens of thousands of cells to be assayed in a single experiment. Although initially developed in animal systems, scRNA-seq has rapidly emerged as a powerful and widely adopted approach in plant biology. Beyond transcriptomics, the integration of single-cell data with chromatin accessibility, proteomics, metabolomics, and spatial omics is enabling a system-level understanding of plant gene regulation and cellular organization. Network-based analytical frameworks further support the reconstruction of gene regulatory networks and the interpretation of complex single-cell data. In this review, we summarize the current technological landscape of plant single-cell studies, discuss key experimental and analytical challenges, and review emerging strategies for validating single-cell discoveries. We also discuss future directions in applying single-cell technologies to woody perennials plants and bioenergy-relevant crops, emphasizing their potential to accelerate the discovery of cell type-specific regulatory mechanisms underlying growth, stress resilience, and biomass production.

Li, Miaomiao [ORNL] (ORCID:0000000321326168)↗
DOE OSTI2026

Probing the limits of genetic recoding using multi-omics-guided evolution

Engineering the genetic code—by reassigning multiple of the 64 natural codons—enables making organisms resistant to all viruses, preventing genetic information exchange, and allowing the biosynthesis of genetically encoded unnatural polymers. However, synonymous codon replacement—recoding—is frequently lethal, and how recoding impacts fitness remains poorly explored. Here, we explore these effects using genome synthesis, directed evolution, and genome-transcriptome-translatome-proteome co-profiling on multiple synthetic Escherichia coli genomes. We construct six partially recoded E. coli strains bearing up to 45.8% of a synthetic genome with a deleterious 57-codon genetic code. As our analyses revealed widespread defects—including unassigned codons in Syn61 and Syn57—we apply multi-omics to revise our genome design and mitigate defects. Using multi-omics, we show that recoding induces transcriptional and translational changes leading to fitness defects under hundreds of conditions. Finally, we develop a multi-omics-guided evolution strategy that rapidly restores fitness, enabling genome synthesis with radical changes.

Nyerges, Akos [Harvard Medical School, Boston, MA ↗
DOE OSTI2024

Signature analysis of high-throughput transcriptomics screening data for mechanistic inference and chemical grouping

Abstract High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1,751 unique chemicals from the EPA’s ToxCast collection in MCF7 cells using 8 concentrations and an exposure duration of 6 h. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms of action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation, and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points of departure, predicting chemicals’ MeOA and grouping chemicals with similar bioactivity profiles.

Toxicology↗
DOE OSTI2024

Mechanical forces orchestrate the metabolism of the developing oilseed rape embryo

The initial free expansion of the embryo within a seed is at some point inhibited by its contact with the testa, resulting in its formation of folds and borders. Although less obvious, mechanical forces appear to trigger and accelerate seed maturation. However, the mechanistic basis for this effect remains unclear. Manipulation of the mechanical constraints affecting either the in vivo or in vitro growth of oilseed rape embryos was combined with analytical approaches, including magnetic resonance imaging and computer graphic reconstruction, immunolabelling, flow cytometry, transcriptomic, proteomic, lipidomic and metabolomic profiling. Our data implied that, in vivo, the imposition of mechanical restraints impeded the expansion of testa and endosperm, resulting in the embryo's deformation. An acceleration in embryonic development was implied by the cessation of cell proliferation and the stimulation of lipid and protein storage, characteristic of embryo maturation. The underlying molecular signature included elements of cell cycle control, reactive oxygen species metabolism and transcriptional reprogramming, along with allosteric control of glycolytic flux. Constricting the space allowed for the expansion of in vitro grown embryos induced a similar response. The conclusion is that the imposition of mechanical constraints over the growth of the developing oilseed rape embryo provides an important trigger for its maturation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Bioinformatic and experimental evidence for suicidal and catalytic plant THI4s

Like fungi and some prokaryotes, plants use a thiazole synthase (THI4) to make the thiazole precursor of thiamin. Fungal THI4s are suicide enzymes that destroy an essential active-site Cys residue to obtain the sulfur atom needed for thiazole formation. In contrast, certain prokaryotic THI4s have no active-site Cys, use sulfide as sulfur donor, and are truly catalytic. The presence of a conserved active-site Cys in plant THI4s and other indirect evidence implies that they are suicidal. To confirm this, we complemented the Arabidopsistz-1 mutant, which lacks THI4 activity, with a His-tagged Arabidopsis THI4 construct. LC–MS analysis of tryptic peptides of the THI4 extracted from leaves showed that the active-site Cys was predominantly in desulfurated form, consistent with THI4 having a suicide mechanism in planta. Unexpectedly, transcriptome data mining and deep proteome profiling showed that barley, wheat, and oat have both a widely expressed canonical THI4 with an active-site Cys, and a THI4-like paralog (non-Cys THI4) that has no active-site Cys and is the major type of THI4 in developing grains. Transcriptomic evidence also indicated that barley, wheat, and oat grains synthesize thiamin de novo, implying that their non-Cys THI4s synthesize thiazole. Structure modeling supported this inference, as did demonstration that non-Cys THI4s have significant capacity to complement thiazole auxotrophy in Escherichia coli. There is thus a prima facie case that non-Cys cereal THI4s, like their prokaryotic counterparts, are catalytic thiazole synthases. Bioenergetic calculations show that, relative to suicide THI4s, such enzymes could save substantial energy during the grain-filling period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Transitions of foliar mycobiota community and transcriptome in response to pathogenic conifer needle interactions

Profiling the host–mycobiota interactions in healthy vs. diseased forest ecosystems helps understand the dynamics of understudied yet increasingly important threats to forest health that are emerging due to climate change. We analyzed the structural and functional changes of the mycobiota and the responses of Pinus contorta in the Lophodermella needle cast pathosystem through metabarcoding and metatranscriptomics. When needles transitioned from asymptomatic to symptomatic, dysbiosis of the mycobiota occurred, but with an enrichment of Lophodermella pathogens. Many pathogenicity-related genes were highly expressed by the mycobiota at the necrotrophic phase, showing an active pathogen response that are absent in asymptomatic needles. This study also revealed that Lophodermella spp. are members of a healthy needle mycobiota that have latent lifestyles suggesting that other pine needle pathogens may have similar biology. Interestingly, Pinus contorta upregulated defense genes in healthy needles, indicating response to fungal recognition, while a variety of biotic and abiotic stresses genes were activated in diseased needles. Further investigation to elucidate the possible antagonistic interplay of other biotic members leading to disease progression and/or suppression is warranted. This study provides insights into microbial interactions in non-model pathosystems and contributes to the development of new forest management strategies against emerging latent pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genome evolution and transcriptome plasticity is associated with adaptation to monocot and dicot plants in Colletotrichum fungi

Colletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes. Here we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation. These results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Three pairs of fungal Trametes strains isolated from distinct geographic origins show conserved genomic features and adaptive response to plant biomass

The genomes of white-rot fungi hold extended repertoires of enzymes active on virtually all the chemical bonds that intertwine lignocellulose polymers, and several Trametes species have been identified as powerful tools for biorefinery or bioremediation. However, only few studies have addressed the intra-species polymorphism one would expect from fungal strains collected in contrasted environments. We compared the genome sequence of pairs of strains collected in different geographic areas, for each of three fungal species. Using an updated list of the predicted functions for fungal ligno- and cellulolytic enzymes (CAZymes), we observed a high conservation of the gene repertoires among the six strains. We compared the adaptative response of the fungi grown on crystalline cellulose, wheat straw, aspen or pine sawdust by transcriptomics and secretomics. The gene regulation profiles were determined by the species and the substrates, rather than the strain. The secretomes did not show marked differences in the sets of secreted CAZymes after 3 day-growth on the substrates. We identified five transcription factor genes and two sesquiterpenoid synthesis genes induced during growth on lignocellulose. Wider studies using larger sets of strains will be necessary to evaluate the genericity of our findings, and to assess the phenotype diversity one could expect from geographic diversity as compared to taxonomic diversity in Trametes fungi.

Drula, E. [French National Research Institute for ↗
DOE OSTI2021

CreA-mediated repression of gene expression occurs at low monosaccharide levels during fungal plant biomass conversion in a time and substrate dependent manner

Carbon catabolite repression enables fungi to utilize the most favourable carbon source in the environment, and is mediated by a key regulator, CreA, in most fungi. CreA-mediated regulation has mainly been studied at high monosaccharide concentrations, an uncommon situation in most natural biotopes. In nature, many fungi rely on plant biomass as their major carbon source by producing enzymes to degrade plant cell wall polysaccharides into metabolizable sugars. To determine the role of CreA when fungi grow in more natural conditions and in particular with respect to degradation and conversion of plant cell walls, we compared transcriptomes of a creA deletion and reference strain of the ascomycete Aspergillus niger during growth on sugar beet pulp and wheat bran. Transcriptomics, extracellular sugar concentrations and growth profiling of A. niger on a variety of carbon sources, revealed that also under conditions with low concentrations of free monosaccharides, CreA has a major effect on gene expression in a strong time and substrate composition dependent manner. In addition, we compared the CreA regulon from five fungi during their growth on crude plant biomass or cellulose. It showed that CreA commonly regulated genes related to carbon metabolism, sugar transport and plant cell wall degrading enzymes across different species. We therefore conclude that CreA has a crucial role for fungi also in adapting to low sugar concentrations as occurring in their natural biotopes, which is supported by the presence of CreA orthologs in nearly all fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula

Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

59 BASIC BIOLOGICAL SCIENCES↗