Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Analysis, RNA”

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

Decoding co-/post-transcriptional complexities of plant transcriptomes and epitranscriptome using next-generation sequencing technologies

Next-generation sequencing (NGS) technologies - Illumina RNA-seq, Pacific Biosciences isoform sequencing (PacBio Iso-seq), and Oxford Nanopore direct RNA sequencing (DRS) - have revealed the complexity of plant transcriptomes and their regulation at the co-/post-transcriptional level. Global analysis of mature mRNAs, transcripts from nuclear run-on assays, and nascent chromatin-bound mRNAs using short as well as full-length and single-molecule DRS reads have uncovered potential roles of different forms of RNA polymerase II during the transcription process, and the extent of co-transcriptional pre-mRNA splicing and polyadenylation. These tools have also allowed mapping of transcriptome-wide start sites in cap-containing RNAs, poly(A) site choice, poly(A) tail length, and RNA base modifications. The emerging theme from recent studies is that reprogramming of gene expression in response to developmental cues and stresses at the co-/post-transcriptional level likely plays a crucial role in eliciting appropriate responses for optimal growth and plant survival under adverse conditions. Although the mechanisms by which developmental cues and different stresses regulate co-/post-transcriptional splicing are largely unknown, a few recent studies indicate that the external cues target spliceosomal and splicing regulatory proteins to modulate alternative splicing. In this review, we provide an overview of recent discoveries on the dynamics and complexities of plant transcriptomes, mechanistic insights into splicing regulation, and discuss critical gaps in co-/post-transcriptional research that need to be addressed using diverse genomic and biochemical approaches.

Biochemistry & Molecular Biology↗

Differential expression of members of the annexin multigene family in Arabidopsis

Although in most plant species no more than two annexin genes have been reported to date, seven annexin homologs have been identified in Arabidopsis, Annexin Arabidopsis 1-7 (AnnAt1--AnnAt7). This establishes that annexins can be a diverse, multigene protein family in a single plant species. Here we compare and analyze these seven annexin gene sequences and present the in situ RNA localization patterns of two of these genes, AnnAt1 and AnnAt2, during different stages of Arabidopsis development. Sequence analysis of AnnAt1--AnnAt7 reveals that they contain the characteristic four structural repeats including the more highly conserved 17-amino acid endonexin fold region found in vertebrate annexins. Alignment comparisons show that there are differences within the repeat regions that may have functional importance. To assess the relative level of expression in various tissues, reverse transcription-PCR was carried out using gene-specific primers for each of the Arabidopsis annexin genes. In addition, northern blot analysis using gene-specific probes indicates differences in AnnAt1 and AnnAt2 expression levels in different tissues. AnnAt1 is expressed in all tissues examined and is most abundant in stems, whereas AnnAt2 is expressed mainly in root tissue and to a lesser extent in stems and flowers. In situ RNA localization demonstrates that these two annexin genes display developmentally regulated tissue-specific and cell-specific expression patterns. These patterns are both distinct and overlapping. The developmental expression patterns for both annexins provide further support for the hypothesis that annexins are involved in the Golgi-mediated secretion of polysaccharides.

Non-NASA Center↗

Derivatization of Mirror-Image L -Nucleic Acids with 2'-OMe Modification for Thermal and Structural Stabilization

To further expand the functionality and enhance the stability of mirror-image nucleic acids as advanced agents for basic research and therapeutic design, we have synthesized 2'-deoxy-2'-methoxy- L -uridine phosphoramidite and incorporated it into L -DNA and L -RNA by solid-phase synthesis quantitatively. We found that the thermostability of L -nucleic acids is dramatically improved after introducing the modifications. Moreover, we successfully crystallized both L -DNA and L -RNA duplexes containing the 2'-OMe modifications and sharing identical sequences. Crystal structure determination and analysis revealed the overall structures of the mirror-image nucleic acids, and for the first time it was possible to interpret the structural deviations caused by 2'-OMe and 2'-OH groups in the oligonucleotides, which are very similar. In conclusion, this novel chemical nucleic acid modification has the potential to be used to design nucleic acid-based therapeutics and materials in the future.

2'-OMe-modifications↗

Molecular Regulatory Mechanisms Drive Emergent Pathogenetic Properties of Neisseria gonorrhoeae

Neisseria gonorrhoeae is the causative agent of the sexually transmitted infection (STI) gonorrhea, with an estimated 87 million annual cases worldwide. N. gonorrhoeae predominantly colonizes the male and female genital tract (FGT). In the FGT, N. gonorrhoeae confronts fluctuating levels of nutrients and oxidative and non-oxidative antimicrobial defenses of the immune system, as well as the resident microbiome. One mechanism utilized by N. gonorrhoeae to adapt to this dynamic FGT niche is to modulate gene expression primarily through DNA-binding transcriptional regulators. Here, we describe the major N. gonorrhoeae transcriptional regulators, genes under their control, and how these regulatory processes lead to pathogenic properties of N. gonorrhoeae during natural infection. We also discuss the current knowledge of the structure, function, and diversity of the FGT microbiome and its influence on gonococcal survival and transcriptional responses orchestrated by its DNA-binding regulators. We conclude with recent multi-omics data and modeling tools and their application to FGT microbiome dynamics. Understanding the strategies utilized by N. gonorrhoeae to regulate gene expression and their impact on the emergent characteristics of this pathogen during infection has the potential to identify new effective strategies to both treat and prevent gonorrhea.

59 BASIC BIOLOGICAL SCIENCES↗

Single‐Cell Nanodroplet Processing Proteomics Pipeline for Analysis of Human‐Derived Microglia

Single-cell omics tools provide unique insights into heterogeneous cell populations and their responses to stimuli. For example, single-cell RNA sequencing has identified several transcriptionally distinct populations of microglia, which are resident immune cells of the central nervous system (CNS) that are responsive to CNS injury, infection, and neurodegeneration. To date, single-cell studies of microglia have focused on RNA-sequencing or cytometry by time of flight (CyTOF), which provide indirect readouts of protein abundance or quantification of a limited number of targets. Herein, we present a workflow based on FACS-assisted isolation, cryopreservation, and nanodroplet-based processing for single-cell mass spectrometry proteomics analysis of the postmortem human brain cortex-derived microglia. From a single microglial cell, 1039 proteins could be identified on average. As a proof-of-principle, we applied single-cell proteomics for exploring the heterogeneity of brain microglia at the cellular level. This pilot proteomics data partially recapitulates the prior microglia subtypes. Specifically, we determined that mitochondrial proteins, in particular members of NADH dehydrogenase (Complex I), cytochrome b-c1 (Complex III), cytochrome c oxidase (Complex IV), F1-ATPase (Complex V), and Na+/K+-ATPase complex, drive variation across microglia. This pipeline offers the potential for identifying functionally and analytically relevant protein targets for microglia in Alzheimer's disease and other neurological disorders.

59 BASIC BIOLOGICAL SCIENCES↗

De novo transcriptome sequencing of Capsicum frutescens . L and comprehensive analysis of salt stress alleviating mechanism by Bacillus atrophaeus WU ‐9

Abstract Salt stress, as one of the most severe environmental stresses, can cause a series of changes in plants. However, the explanation of plant salt stress alleviating mechanism of plant growth–promoting rhizobacteria (PGPR) was hindered by the limited availability of transcriptomic information for salt stress‐treated plants grown in a microorganism‐controlled environment. Our previous reports have selected Bacillus atrophaeus WU‐9 as PGPR significantly alleviating pepper ( Capsicum frutescens . L) salt stress. In this work, the RNA‐seq analysis of salt stress‐treated and untreated plants, grown with and without WU‐9 in a microorganism‐controlled environment, was used to reveal the plant salt stress alleviating mechanisms of WU‐9. Twelve sequencing libraries, prepared by treating with WU‐9 and salt (150 mM NaCl for 36 h), were constructed by RNA‐Seq technique. Non‐inoculated seedlings mainly respond to salt stress through regulation of signal transduction, such as ethylene‐activated signaling pathway, signaling and cell communication, etc. And ethylene signal participated in salt stress response in pepper through regulating defense responses, fruit ripening and senescence. WU‐9 inoculation under salt stress mainly improves salt tolerance and plant growth by regulating salt stress‐responding ethylene and auxin signal transduction, utilization of proline, photosynthesis, antioxidant enzyme activities and cell enlargement. Furthermore, 86 differentially expressed genes and 20 transcription factors were identified as associated with salt stress response and tolerance. Thus, this innovative transcriptomic study identified the salt stress response and alleviation in C. frutescens . L with PGPR inoculation. This result provided novel insights into the salinity alleviation in pepper regulated by PGPR.

Wang, Wenfei↗

Impact of Soil Salinity on the Cowpea Nodule-Microbiome and the Isolation of Halotolerant PGPR Strains to Promote Plant Growth under Salinity Stress

Cowpea is one of the major legumes cultivated in arid and semiarid regions of the world. Four soil-microbial samples (SS-1 through SS-4) collected from semiarid soils in Punjab, Pakistan were planted with cowpea (Vigna unguiculata) crops, which were grown under salinity stress to analyze bacterial composition in the rhizosphere and within nodules using cultivation-dependent and -independent methods. Two varieties, 603 and the salt-tolerant CB 46, were each inoculated with or without the four different native soil samples or grown in medium either N-deficient (-N) or supplemented with N (+N). Plants inoculated with soil samples SS-2 and SS-4 grew better than plants inoculated with SS-1- and SS-3 and grew comparably with the +N controls. Environmental DNA (eDNA) was isolated from SS-1 and SS-4, and, by 16S ribosomal RNA sequencing, the soil microbiomes consisted mainly of Actinobacteria, Firmicutes, Proteobacteria, and other nonproteobacterial genera. However, analysis of eDNA isolated from cowpea nodules established by the trap plants showed that the nodule microbiome consisted almost exclusively of proteobacterial sequences, particularly species of Bradyrhizobium. Bacteria were isolated from both soils and nodules, and 34 of the 51 isolates tested positive for plant-growth-promoting rhizobacteria traits in plate assays. Many could serve as future inocula for crops in arid soils. The discrepancy between the types of bacteria isolated by culturing bacteria isolated from surface-sterilized cowpea nodules (proteobacteria and nonproteobacteria) versus those detected by sequencing DNA isolated from the nodules (proteobacteria) from cowpea nodules (proteobacteria and nonproteobacteria) versus those detected in the nodule microbiome (proteobacteria) needs further study.

59 BASIC BIOLOGICAL SCIENCES↗

Development of Solid-State Nanopore Technology for Life Detection

Biomarkers for life on Earth are an important starting point to guide the search for life elsewhere. However, the search for life beyond Earth should incorporate technologies capable of recognizing an array of potential biomarkers beyond what we see on Earth, in order to minimize the risk of false negatives from life detection missions. With this in mind, charged linear polymers may be a universal signature for life, due to their ability to store information while also inherently reducing the tendency of complex tertiary structure formation that significantly inhibit replication. Thus, these molecules are attractive targets for biosignature detection as potential "self-sustaining chemical signatures." Examples of charged linear polymers, or polyelectrolytes, include deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) as well as synthetic polyelectrolytes that could potentially support life, including threose nucleic acid (TNA) and other xenonucleic acids (XNAs). Nanopore analysis is a novel technology that has been developed for singlemolecule sequencing with exquisite single nucleotide resolution which is also well-suited for analysis of polyelectrolyte molecules. Nanopore analysis has the ability to detect repeating sequences of electrical charges in organic linear polymers, and it is not molecule- specific (i.e. it is not restricted to only DNA or RNA). In this sense, it is a better life detection technique than approaches that are based on specific molecules, such as the polymerase chain reaction (PCR), which requires that the molecule being detected be composed of DNA.

Bywaters, K. B.↗

Space Flown Rodent Liver RNA Sequencing Data for Machine Learning in Space Biology Research

High-throughput nucleic acid sequencing (DNA-seq, RNA-seq) has become widespread in biomedical research due to the growing availability and affordability of these assays. Data analysis has been accelerated in recent years by the adoption of artificial intelligence (AI) and machine learning (ML) techniques by biomedical researchers. In space biology research, RNAseq datasets from space-flown experimental samples are critical for characterizing the gene expression aberrations associated with exposure to spaceflight stressors. However, space biological experiments tend to be very low sample size, so identifying proper AI/ML algorithms for sequencing data analysis is an ongoing challenge since these algorithms typically require large sample size. The NASA Science Mission Directorate (SMD) has started the “Benchmark Initiative for AI/ML”, focused on creating datasets meant for three main applications: 1) scientific benchmarking, which finds the best algorithm for a specific problem; 2) application benchmarking, which measures algorithm performance against a set of parameters; and 3) system benchmarking, which evaluates performance of hardware and software architecture. These scientific benchmarks consist of an AI-ready dataset and a reference implementation on a specific scientific question. In this work, we focused on generating standardized datasets to allow the scientific community to benchmark AI/ML algorithms in the domain of space biology. We present here a standardized, AI-ready, publicly available benchmark dataset for space biology RNA-seq data as a collaboration between the NASA AI4LS (Artificial Intelligence for Life Sciences) working group. and NASA’s SMD. This dataset consists of space-flown and ground control mouse liver found in the NASA GeneLab omics database. However, to amplify the small sample number (n=112 samples) for ML purposes, we employ Gaussian noise and a generative adversarial network to extend this dataset to 6,000 synthetic samples, matching the original gene expression characteristics.

James Casaletto↗

Sequence, overproduction and purification of Vibrio proteolyticus ribosomal protein L18 for in vitro and in vivo studies

A strategy suggested by comparative genomic studies was used to amplify the entire Vibrio proteolyticus (Vp) gene for ribosomal protein L18. Vp L18 and its flanking regions were sequenced and compared with the deduced amino acid (aa) sequences of other known L18 proteins. A 26-aa residue segment at the carboxy terminus contains many strongly conserved residues and may be critical for the L18 interaction with 5S rRNA. This approach should allow rapid characterization of L18 from large numbers of bacteria. Both Vp L18 and Escherichia coli (Ec) L18 were overproduced and purified using a T7 expression vector which fuses an N-terminal peptide segment (His-tag) containing 6 histidine residues to the recombinant protein. The purified fusion proteins, Vp His::L18 and Ec His::L18, were both found to bind to either the Vp 5S or Ec 5S rRNAs in vitro. Vp His::L18 protein was also shown to incorporate into Ec ribosomes in vivo. This His-tag strategy likely will have general applicability for the study of ribosomal proteins in vitro and in vivo.

Non-NASA Center↗

Polyomavirus JCV excretion and genotype analysis in HIV-infected patients receiving highly active antiretroviral therapy

OBJECTIVE: To assess the frequency of shedding of polyomavirus JC virus (JCV) genotypes in urine of HIV-infected patients receiving highly active antiretroviral therapy (HAART). METHODS: Single samples of urine and blood were collected prospectively from 70 adult HIV-infected patients and 68 uninfected volunteers. Inclusion criteria for HIV-infected patients included an HIV RNA viral load < 1000 copies, CD4 cell count of 200-700 x 106 cells/l, and stable HAART regimen. PCR assays and sequence analysis were carried out using JCV-specific primers against different regions of the virus genome. RESULTS: JCV excretion in urine was more common in HIV-positive patients but not significantly different from that of the HIV-negative group [22/70 (31%) versus 13/68 (19%); P = 0.09]. HIV-positive patients lost the age-related pattern of JCV shedding (P = 0.13) displayed by uninfected subjects (P = 0.01). Among HIV-infected patients significant differences in JCV shedding were related to CD4 cell counts (P = 0.03). Sequence analysis of the JCV regulatory region from both HIV-infected patients and uninfected volunteers revealed all to be JCV archetypal strains. JCV genotypes 1 (36%) and 4 (36%) were the most common among HIV-infected patients, whereas type 2 (77%) was the most frequently detected among HIV-uninfected volunteers. CONCLUSION: These results suggest that JCV shedding is enhanced by modest depressions in immune function during HIV infection. JCV shedding occurred in younger HIV-positive persons than in the healthy controls. As the common types of JCV excreted varied among ethnic groups, JCV genotypes associated with progressive multifocal leukoencephalopathy may reflect demographics of those infected patient populations.

Non-NASA Center↗

HNRNPUL1 inhibits cisplatin sensitivity of esophageal squamous cell carcinoma through regulating the formation of circMAN1A2

Highlights: • HNRNPUL1 is an independent prognostic factor of ESCC patients receiving platinum-based chemotherapy. • Knockdown of HNRNPUL1 enhances cisplatin sensitivity of ESCC cells. • HNRNPUL1 inhibits the cisplatin sensitivity of ESCC cells by regulating the formation of circMAN1A2. Cisplatin (CDDP) is widely used for chemotherapy of esophageal squamous cell carcinoma (ESCC) but the drug resistance limits its therapeutic benefit. Heterogeneous nuclear ribonucleoprotein U-like 1 (HNRNPUL1) belongs to the family of RNA-binding proteins (RBPs) and is involved in DNA damage repair. To investigate whether and how HNRNPUL1 affects CDDP resistance of ESCC, we evaluated the expression of HNRNPUL1 and found that it was associated with recurrence in ESCC patients receiving postoperative platinum-based chemotherapy and was an independent prognostic factor for disease-free survival (DFS). Besides, we showed that the reduced expression of HNRNPUL1 enhanced the CDDP sensitivity of ESCC cells. Furthermore, RNA immunoprecipitation coupled with high-throughput sequencing (RIP-seq) were performed and a range of HNRNPUL1-binding RNAs influenced by CDDP treatment were identified followed by bioinformatics analysis. In terms of mechanism, we found that HNRNPUL1 inhibited CDDP sensitivity of ESCC cells by regulating the CDDP sensitivity-inhibited circular RNA (circRNA) MAN1A2 formation. Taken together, our results first demonstrated the role of HNRNPUL1 in CDDP resistance of ESCC and suggested that HNRNPUL1 may be a potential target of ESCC chemotherapy.

60 APPLIED LIFE SCIENCES↗

Exploring Connectivity in Sequence Space of Functional RNA

Emergence of replicable genetic molecules was one of the marking points in the origin of life, evolution of which can be conceptualized as a walk through the space of all possible sequences. A theoretical concept of fitness landscape helps to understand evolutionary processes through assigning a value of fitness to each genotype. Then, evolution of a phenotype is viewed as a series of consecutive, single-point mutations. Natural selection biases evolution toward peaks of high fitness and away from valleys of low fitness. whereas neutral drift occurs in the sequence space without direction as mutations are introduced at random. Large networks of neutral or near-neutral mutations on a fitness landscape, especially for sufficiently long genomes, are possible or even inevitable. Their detection in experiments, however, has been elusive. Although a few near-neutral evolutionary pathways have been found, recent experimental evidence indicates landscapes consist of largely isolated islands. The generality of these results, however, is not clear, as the genome length or the fraction of functional molecules in the genotypic space might have been insufficient for the emergence of large, neutral networks. Thorough investigation on the structure of the fitness landscape is essential to understand the mechanisms of evolution of early genomes. RNA molecules are commonly assumed to play the pivotal role in the origin of genetic systems. They are widely believed to be early, if not the earliest, genetic and catalytic molecules, with abundant biochemical activities as aptamers and ribozymes, i.e. RNA molecules capable, respectively, to bind small molecules or catalyze chemical reactions. Here, we present results of our recent studies on the structure of the sequence space of RNA ligase ribozymes selected through in vitro evolution. Several hundred thousands of sequences active to a different degree were obtained by way of deep sequencing. Analysis of these sequences revealed several large clusters defined such that every sequence in a cluster can be reached from any other sequence in the same cluster through a series of single point mutations. Sequences in a single cluster appear to adopt more than one secondary structure. The mechanism of refolding within a single cluster was examined. To shed light on possible evolutionary paths in the space of ribozymes, the connectivity between clusters was investigated. The effect of length of RNA molecules on the structure of the fitness landscape and possible evolutionary paths was examined by way of comparing functional sequences of 20 and 80 nucleobases in length. It was found that sequences of different lengths shared secondary structure motifs that were presumed responsible for catalytic activity, with increasing complexity and global structural rearrangements emerging in longer molecules.

Wei, Chenyu↗

Somatic Mutation in Mice on the International Space Station (ISS): Guanine Substitution Suggests Link to Cancer Risk

We conducted comprehensive analysis of single nucleotide somatic mutations in mice exposed to microgravity and other factors aboard the International Space Station (ISS), using data archived in GeneLab. Animals in the experimental cohort consisted of mice that spent 37 days on the ISS within the Rodent Habitat. Ground control animals consisted of mice of identical age, sex, strain, in a terrestrial Rodent Habitat controlled for temperature, humidity and carbon dioxide levels, to match ISS conditions as closely as possible. RNA extracted from eye, liver, skeletal muscle, and kidney tissue specimens was subjected to next-generation sequencing to acquire primary data. Our analysis employed cutting-edge software developed at NASA Ames Research Center, executed on the NASA Ames Supercomputer and on another high-performance computer, for accurate variant calling of single point mutations. ISS-flown mice exhibited a notably heightened level of somatic mutation compared to control mice. The degree of somatic mutation correlated with the degree of gene expression across the four tissue types, i.e., the greatest rate of mutation accumulation was seen in highly expressed genes. We discovered that guanine substitutions were the most common type of somatic mutation. This observation is consistent with the hypothesis that DNA mutation events stem from reactive oxygen/nitrogen/chlorine species-mediated guanine oxidation induced by the spaceflight environment. Since guanine oxidation is a prominent feature of the DNA mutation landscape that accompanies malignant transformation, our findings suggest a possible link between the spaceflight environment and cancer risk that is independent of radiation carcinogenesis.

ISS↗

Human Host Cellular Response to HCoV-229E Infection Proteomics (ACS-JM-DP2)

The purpose of this experiment was to evaluate the human host cellular response to wild-type Human coronavirus strain 229E (HCoV-229E) infection. Sample data was obtained for mock and infected immortalized human lung epithelial cells (A549) (MOI 5) nuclear extracts, immortalized human lung fibroblasts cells (MRC5) (MOI5) nuclear extracts, and primary human airway epithelial (HAE) (MOI 3) cells from lung tissue and processed for proteome analysis. Processed datasets are openly accessible from the download button and contain secondary processed proteomic results files and supporting metadata materials. Experimental proteomics samples were prepared using Limited Proteolysis (LiP) methods for Label-free quantification (LFQ) and global proteomic evaluation. Sample data was acquired using a Q-Exactive HF-X mass spectrometer and was processed and compiled using MaxQuant software (v.1.6.17.0). Processed proteomic data downloads include a sample naming key, processed MaxQuant results/parameters, and protein annotated relative abundance files. See corresponding primary data accessions below and Viral Experiment LiP Analysis source code supporting data transparency and reuse. Experimental transcriptomics samples were collected in parallel and processed for RNA sequencing (RNA-Seq) as summarized under ACS-DP1 (https://data.pnnl.gov/group/nodes/dataset/34069).

59 BASIC BIOLOGICAL SCIENCES↗

CrossMP: Enabling Cross-Modality Translation between Single-Cell RNA-Seq and Single-Cell ATAC-Seq through Web-Based Portal

In recent years, there has been a growing interest in profiling multiomic modalities within individual cells simultaneously. One such example is integrating combined single-cell RNA sequencing (scRNA-seq) data and single-cell transposase-accessible chromatin sequencing (scATAC-seq) data. Integrated analysis of diverse modalities has helped researchers make more accurate predictions and gain a more comprehensive understanding than with single-modality analysis. However, generating such multimodal data is technically challenging and expensive, leading to limited availability of single-cell co-assay data. Here, we propose a model for cross-modal prediction between the transcriptome and chromatin profiles in single cells. Our model is based on a deep neural network architecture that learns the latent representations from the source modality and then predicts the target modality. It demonstrates reliable performance in accurately translating between these modalities across multiple paired human scATAC-seq and scRNA-seq datasets. Additionally, we developed CrossMP, a web-based portal allowing researchers to upload their single-cell modality data through an interactive web interface and predict the other type of modality data, using high-performance computing resources plugged at the backend.

59 BASIC BIOLOGICAL SCIENCES↗

A 5.8S nuclear ribosomal RNA gene sequence database: applications to ecology and evolution

We complied a 5.8S nuclear ribosomal gene sequence database for animals, plants, and fungi using both newly generated and GenBank sequences. We demonstrate the utility of this database as an internal check to determine whether the target organism and not a contaminant has been sequenced, as a diagnostic tool for ecologists and evolutionary biologists to determine the placement of asexual fungi within larger taxonomic groups, and as a tool to help identify fungi that form ectomycorrhizae.

Databases, Factual↗