Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RNA sequencing”

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 199 records · Page 11

Models of the Protocellular Structures, Functions and Evolution

In the absence of extinct or extant record of protocells, the most direct way to test our understanding of the origin of cellular life is to construct laboratory models that capture important features of protocellular systems. Such efforts are currently underway in a collaborative project between NASA-Ames, Harvard medical School and University of California. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures. The centerpiece of this project is a method for the in vitro evolution of protein enzymes toward arbitrary catalytic targets. A similar approach has already been developed for nucleic acids: First, a very large population of candidate molecules is generated using a random synthetic approach. Next, the small numbers of molecules that can accomplish the desired task are selected. These molecules are next vastly multiplied using the polymerase chain reaction. A mutagenic approach, in which the sequences of selected molecules are randomly altered, can yield further improvements in performance or alterations of specificities. Unfortunately, the catalytic potential of nucleic acids is rather limited. Proteins are more catalytically capable but cannot be directly amplified. In the new technique, this problem is circumvented by covalently linking each protein of the initial, diverse, pool to the RNA sequence that codes for it. Then, selection is performed on the proteins, but the nucleic acids are replicated. To date, we have obtained "a proof of concept" by evolving simple, novel proteins capable of selectively binding adenosine tri-phosphate (ATP). Our next goal is to create an enzyme that can phosphorylate amino acids and another to catalyze the formation of peptide bonds in the absence of nucleic acid templates. This latter reaction does not take place in contemporary cells. once developed, these enzymes will be encapsulated in liposomes so that they will function in a simulated cellular environment. To provide a continuous energy supply, usually needed to activate the substrates, an energy transduction complex which generates ATP from adenosine diphosphate, inorganic phosphate and light will be used. This system, consisting of two modern proteins, ATP synthase and bacteriorhodopsin, has already been built and shown to work efficiently. By coupling chemical synthesis to such a system, it will be possible to drive chemical reactions by light if only the substrates for these reactions are supplied.

Pohorille, Andrew↗

Models of Protocellular Structure, Function and Evolution

In the absence of any record of protocells, the most direct way to test our understanding of the origin of cellular life is to construct laboratory models that capture important features of protocellular systems. Such efforts are currently underway in a collaborative project between NASA-Ames, Harvard Medical School and University of California. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures. The centerpiece of this project is a method for the in vitro evolution of protein enzymes toward arbitrary catalytic targets. A similar approach has already been developed for nucleic acids in which a small number of functional molecules are selected from a large, random population of candidates. The selected molecules are next vastly multiplied using the polymerase chain reaction. A mutagenic approach, in which the sequences of selected molecules are randomly altered, can yield further improvements in performance or alterations of specificities. Unfortunately, the catalytic potential of nucleic acids is rather limited. Proteins are more catalytically capable but cannot be directly amplified. In the new technique, this problem is circumvented by covalently linking each protein of the initial, diverse, pool to the RNA sequence that codes for it. Then, selection is performed on the proteins, but the nucleic acids are replicated. Additional information is contained in the original extended abstract.

New, Michael H.↗

Development and validation of a next-generation sequencing-based multigene assay to predict the prognosis of estrogen receptor-positive, HER2-negative breast cancer

Purpose: Multigene assays provide useful prognostic information regarding hormone receptor (HR)-positive breast cancer. Next-generation sequencing (NGS)-based platforms have numerous advantages including reproducibility and adaptability in local laboratories. Here we aimed to develop and validate an NGS-based multigene assay to predict the distant recurrence risk. Experimental Design: In total, 179 genes including 30 reference genes highly correlated with the 21-gene recurrence score (RS) algorithm were selected from public databases. Targeted RNA-sequencing was performed using 250 and 93 archived breast cancer samples with a known RS in the training and verification sets, respectively, to develop the algorithm and NGS–Prognostic Score (NGS-PS). The assay was validated in 413 independent samples with long-term follow-up data on distant metastasis. Results: In the verification set, the NGS-PS and 21-gene RS displayed 91.4% concurrence (85/93 samples). In the validation cohort of 413 samples, area under the receiver operating characteristic curve plotted using NGS-PS values classified for distant recurrence was 0.76. The best NGS-PS cut-off value predicting distant metastasis was 20. Furthermore, 269 and 144 patients were classified as low- and high-risk patients in accordance with the cut-off. Five- and 10-year estimates of distant metastasis–free survival (DMFS) for low- versus high-risk groups were 97.0% versus 77.8% and 93.2% versus 64.4%, respectively. The age-related HR for distant recurrence without chemotherapy was 9.73 (95% CI, 3.59–26.40) and 3.19 (95% CI, 1.40–7.29) for patients aged ≤50 and >50 years, respectively. Conclusions: The newly developed and validated NGS-based multigene assay can predict the distant recurrence risk in ER-positive, HER2-negative breast cancer.

60 APPLIED LIFE SCIENCES↗

Host genotype controls ecological change in the leaf fungal microbiome

Leaf fungal microbiomes can be fundamental drivers of host plant success, as they contain pathogens that devastate crop plants and taxa that enhance nutrient uptake, discourage herbivory, and antagonize pathogens. We measured leaf fungal diversity with amplicon sequencing across an entire growing season in a diversity panel of switchgrass ( Panicum virgatum ). We also sampled a replicated subset of genotypes across 3 additional sites to compare the importance of time, space, ecology, and genetics. We found a strong successional pattern in the microbiome shaped both by host genetics and environmental factors. Further, we used genome-wide association (GWA) mapping and RNA sequencing to show that 3 cysteine-rich receptor-like kinases (crRLKs) were linked to a genetic locus associated with microbiome structure. We confirmed GWAS results in an independent set of genotypes for both the internal transcribed spacer (ITS) and large subunit (LSU) ribosomal DNA markers. Fungal pathogens were central to microbial covariance networks, and genotypes susceptible to pathogens differed in their expression of the 3 crRLKs, suggesting that host immune genes are a principal means of controlling the entire leaf microbiome.

59 BASIC BIOLOGICAL SCIENCES↗

Structures of MERS1, the 5' processing enzyme of mitochondrial mRNAs in Trypanosoma brucei

Most mitochondrial mRNAs are transcribed as polycistronic precursors that are cleaved by endonucleases to produce mature mRNA transcripts. However, recent studies have shown that mitochondrial transcripts in the kinetoplastid protozoan, Trypanosoma brucei, are transcribed individually. Also unlike most mitochondrial mRNAs, the 5' end of these transcripts harbor a triphosphate that is hydrolyzed. This modification is carried out by a putative Nudix hydrolase called MERS1. The Nudix motif in MERS1 is degenerate, lacking a conserved glutamic acid, thus it is unclear how it may bind its substrates and whether it contains a Nudix fold. To obtain insight into this unusual hydrolase, we determined structures of apo, GTP-bound and RNA-bound T. brucei MERS1 to 2.30 Å, 2.45 Å, and 2.60 Å, respectively. The MERS1 structure has a unique fold that indeed contains a Nudix motif. The nucleotide bound structures combined with binding studies reveal that MERS1 shows preference for RNA sequences with a central guanine repeat which it binds in a single-stranded conformation. The apo MERS1 structure indicates that a significant portion of its nucleotide binding site folds upon substrate binding. Finally, a potential interaction region for a binding partner, MERS2, that activates MERS1 was identified. The MERS2-like peptide inserts a glutamate near the missing Nudix acidic residue in the RNA binding pocket, suggesting how the enzyme may be activated. Thus, the combined studies reveal insight into the structure and enzyme properties of MERS1 and its substrate-binding activities.

59 BASIC BIOLOGICAL SCIENCES↗

Egr-1 mediated cardiac miR-99 family expression diverges physiological hypertrophy from pathological hypertrophy

The physiological cardiac hypertrophy is an adaptive condition without myocyte cell death, while pathological hypertrophy is a maladaptive condition associated with myocyte cell death. This study explores the miRNome of α-2M-induced physiologically hypertrophied cardiomyocytes and the role of miRNA-99 family during cardiac hypertrophy. Physiological and pathological cardiac hypertrophy was induced in H9c2 cardiomyoblast cell lines using α-2M and isoproterenol respectively. Total RNA isolation and small RNA sequencing were executed for physiological hypertrophy model. The differentially expressed miRNAs and its target mRNAs were validated in animal models. Transcription factor binding sites were predicted in the promoter of specific miRNAs and validated by ChIP-PCR. Subsequently, the selected miRNA was functionally characterized by overexpression and silencing. The effects of silencing of upstream regulator and downstream target gene were studied. Analysis of small RNA reads revealed the differential expression of a large set of miRNAs during hypertrophy, of which miR-99 family was highly downregulated upon α-2M treatment. However, this miR-99 family expression was upregulated during pathological hypertrophy and confirmed in animal models. ChIP-PCR confirms the binding of Egr-1 transcription factor to the miR-99 promoter. Further, silencing of Egr-1 decreased the expression of miR-99. The overexpression or silencing of miR-99 diverges the physiological hypertrophy to pathological hypertrophy and vice versa by regulating Akt-1 pathway. Silencing of Akt-1 replicates the effect of overexpression of miR-99.

60 APPLIED LIFE SCIENCES↗

Genomic insights into the origin, domestication and diversification of Brassica juncea

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000–14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

59 BASIC BIOLOGICAL SCIENCES↗

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis -element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

59 BASIC BIOLOGICAL SCIENCES↗

Identification and characterization of a class of MALAT1 -like genomic loci

The MALAT1 (Metastasis-Associated Lung Adenocarcinoma Transcript 1) gene encodes a noncoding RNA that is processed into a long nuclear retained transcript ( MALAT1 ) and a small cytoplasmic tRNA-like transcript (mascRNA). Using an RNA sequence- and structure-based covariance model, we identified more than 130 genomic loci in vertebrate genomes containing the MALAT1 3' end triple-helix structure and its immediate downstream tRNA-like structure, including 44 in the green lizard Anolis carolinensis . Structural and computational analyses revealed a co-occurrence of components of the 3' end module. MALAT1 -like genes in Anolis carolinensis are highly expressed in adult testis, thus we named them testis-abundant long noncoding RNAs (tancRNAs). MALAT1 -like loci also produce multiple small RNA species, including PIWI-interacting RNAs (piRNAs), from the antisense strand. The 3' ends of tancRNAs serve as potential targets for the PIWI-piRNA complex. Furthermore, we have identified an evolutionarily conserved class of long noncoding RNAs (lncRNAs) with similar structural constraints, post-transcriptional processing, and subcellular localization and a distinct function in spermatocytes.

59 BASIC BIOLOGICAL SCIENCES↗

scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators

Abstract Marek’s disease (MD), a T cell lymphoma disease in chickens, is caused by the Marek’s disease virus (MDV) found ubiquitously in the poultry industry. Genetically resistant Line 6 3 (L6) and susceptible Line 7 2 (L7) chickens have been instrumental to research on avian immune system response to MDV infection. In this study we characterized molecular signatures unique to splenic immune cell types across different genetic backgrounds 6 days after infection. Using three populations, L6, L7, and an F1 cross between L6xL7, we evaluated the immune cell transcriptome of responding cell types using single cell RNA sequencing. Several MDV genes were found expressed mainly in cytotoxic T cells while ICP4 and MEQ MDV genes were expressed across infected cell types. Using the F1 we quantified allele specific expression (ASE) of biallelic SNPs and found biased expression of parental alleles specific to immune cell subtypes. We identified 22 SNPs with ASE in response to MDV infection mapped to gene rich regions surrounding 59 genes of critical importance for chromatin remodeling and transcriptional regulation. Histone deacetylase genes (HDAC1 and HDAC8) had increased expression of L6 alleles, while small nuclear RNA genes (SNORA68 and SNORA72) expressed higher levels of L7 alleles with infection in T cell subsets. SNPs with ASE also mapped genes important for an adequate immune response including GNLY (cytotoxic activity) and PDIA3 (component of MHC class I peptide loading complex), and genes known to promote viral replication (MCM5 and EIF3M). These results show that functional variants associated with susceptibility to MD may have a bigger impact in subsets of immune cell types, and by characterizing the transcriptomes of these subtypes we can unravel molecular signatures specific to MD genomic resistance.

Science & Technology - Other Topics↗

SATB2 expression increased anchorage-independent growth and cell migration in human bronchial epithelial cells

The special AT-rich sequence-binding protein 2 (SATB2) is a protein that binds to the nuclear matrix attachment region of the cell and regulates gene expression by altering chromatin structure. In our previous study, we reported that SATB2 gene expression was induced in human bronchial epithelial BEAS-2B cells transformed by arsenic, chromium, nickel and vanadium. In this study, we show that ectopic expression of SATB2 in the normal human bronchial epithelial cell-line BEAS-2B increased anchorage-independent growth and cell migration, meanwhile, shRNA-mediated knockdown of SATB2 significantly decreased anchorage-independent growth in Ni transformed BEAS-2B cells. RNA sequencing analyses of SATB2 regulated genes revealed the enrichment of those involved in cytoskeleton, cell adhesion and cell-movement pathways. Our evidence supports the hypothesis that SATB2 plays an important role in BEAS-2B cell transformation. - Highlights: • We performed SATB2 overexpression in the BEAS-2B cell line. • We performed SATB2 knockdown in a Ni transformed BEAS-2B cell line. • SATB2 induced anchorage-independent growth and increased cell migration. • SATB2 knockdown significantly decreased anchorage-independent growth. • We identified alterations in gene involved in cytoskeleton, cell adhesion.

60 APPLIED LIFE SCIENCES↗

MacoNPV baculovirus midgut-specific gene expression during infection of the bertha armyworm, Mamestra configurata

Baculoviruses have two forms, occlusion derived virus (ODV) which is responsible for primary infection in host midgut tissue and budded virus (BV), which infects all other host tissues during secondary infection. This study examined the primary infection by ODV of midgut cells of bertha armyworm Mamestra configurata fourth instar larvae and measured the expression of viral genes over a time course of infection. Both digital PCR and RNA sequencing methods showed the profile of transcription to be different from those produced by AcMNPV BV infection of in vitro cell cultures. This included having unique collections of genes expressed early, as well as much greater late gene expression of p6.9 and much reduced expression of polh and p10. These differences likely reflect characteristics unique to the critical step of in vivo midgut cell infection, and provide insights into the processes that regulate viral gene expression in different host tissues. -- Highlights: •The transcriptome of MacoNPV ODV in larval midgut was measured by RNA-seq and digital PCR. •The earliest genes expressed included fusion protein, hoar, and me53. •p6.9 was highly expressed late but polH and p10 were less so. •These patterns are unique from BV of other baculoviruses in tissue culture cells.

60 APPLIED LIFE SCIENCES↗

Stacking interactions in PUF-RNA complexes

Stacking interactions between amino acids and bases are common in RNA-protein interactions. Many proteins that regulate mRNAs interact with single-stranded RNA elements in the 3' UTR (3'-untranslated region) of their targets. PUF proteins are exemplary. Here we focus on complexes formed between a Caenorhabditis elegans PUF protein, FBF, and its cognate RNAs. Stacking interactions are particularly prominent and involve every RNA base in the recognition element. To assess the contribution of stacking interactions to formation of the RNA-protein complex, we combine in vivo selection experiments with site-directed mutagenesis, biochemistry, and structural analysis. Our results reveal that the identities of stacking amino acids in FBF affect both the affinity and specificity of the RNA-protein interaction. Substitutions in amino acid side chains can restrict or broaden RNA specificity. We conclude that the identities of stacking residues are important in achieving the natural specificities of PUF proteins. Similarly, in PUF proteins engineered to bind new RNA sequences, the identity of stacking residues may contribute to 'target' versus 'off-target' interactions, and thus be an important consideration in the design of proteins with new specificities.

59 BASIC BIOLOGICAL SCIENCES↗

A human cell type similar to murine central nervous system perivascular fibroblasts

Highlights: • Human diploid cell lines were established from tumor and non-tumor brain tissue. • Human cells share a common transcriptomic profile with murine perivascular fibroblasts. • Human cells resemble commercially available cell of presumed pericyte origin. The brain vasculature has several specific features, one of them being the blood-brain barrier (BBB), which supports and protects the brain by allowing for the passage of oxygen and nutrients, while at the same time preventing passage of pathogens and toxins. The BBB also prevents efficient delivery of drugs to the brain, e.g. for treatment of brain tumors. In the murine brain, perivascular fibroblasts were recently identified as a novel potential constituent of the BBB. Here we present the existence of human cells that could be the equivalent to the murine brain perivascular fibroblasts. Using RNA sequencing, we show a similar transcriptomic profile of cultured human brain cells and murine perivascular fibroblasts. These data open up a window for new hypotheses on cell types involved in human CNS diseases.

60 APPLIED LIFE SCIENCES↗

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of sgRNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Lastly, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

59 BASIC BIOLOGICAL SCIENCES↗

A-674563 increases chondrocyte marker expression in cultured chondrocytes by inhibiting Sox9 degradation

Highlights: • We established a high-throughput screening system for induction of chondrocytic characters. • We screened 5822 compounds and found that A-674563 enhanced transcription of chondrocyte markers. • A-674563 inhibited the degradation of Sox9 protein. • A-674563 increased expression of ubiquitin-specific peptidase 29, which is known to induce deubiquitination. The implantation of autologous chondrocytes is a therapeutic treatment for articular cartilage damage. However, the benefits are limited due to the expansion of chondrocytes in monolayer culture, which causes loss of chondrocytic characters. Therefore, culture conditions that enhance chondrocytic characters are needed. We screened 5822 compounds and found that A-674563 enhanced the transcription of several chondrocyte marker genes, including Col2a1, Acan and Col11a2, in mouse primary chondrocytes. Experiments using cycloheximide, MG132 and bafilomycin A1 have revealed that Sox9 is degraded through the ubiquitin-proteasome pathway and that A-674563 inhibits this degradation, resulting in larger amount of Sox9 protein. RNA sequencing transcriptome analysis showed that A-674563 increases the expression of the gene that encodes ubiquitin-specific peptidase 29, which is known to induce the deubiquitination of proteins. Although the precise mechanism remains to be determined, our findings indicated that A-674563 could contribute to culture conditions that expand chondrocytes without losing chondrocytic characters.

60 APPLIED LIFE SCIENCES↗

Transcriptome and Proteome Dynamics of a Light-Dark Synchronized Bacterial Cell Cycle

Growth of the ocean’s most abundant primary producer, the cyanobacterium Prochlorococcus, is tightly synchronized to the natural 24-hour light-dark cycle. We sought to quantify the relationship between transcriptome and proteome dynamics that underlie this obligate photoautotroph’s highly choreographed response to the daily oscillation in energy supply. Using RNA-sequencing transcriptomics and mass spectrometry-based quantitative proteomics, we measured timecourses of paired mRNA-protein abundances for 312 genes every 2 hours over a light-dark cycle. These temporal expression patterns reveal strong oscillations in transcript abundance that are broadly damped at the protein level, with mRNA levels varying on average 2.3 times more than the corresponding protein. The single strongest observed protein-level oscillation is in a ribonucleotide reductase, which may reflect a defense strategy against phage infection. The peak in abundance of most proteins also lags that of their transcript by 2–8 hours, and the two are completely antiphase for some genes. While abundant antisense RNA was detected, it apparently does not account for the observed divergences between expression levels. The redirection of flux through central carbon metabolism from daytime carbon fixation to nighttime respiration is associated with quite small changes in relative enzyme abundances. Our results indicate that expression responses to periodic stimuli that are common in natural ecosystems (such as the diel cycle) can diverge significantly between the mRNA and protein levels. Protein expression patterns that are distinct from those of cognate mRNA have implications for the interpretation of transcriptome and metatranscriptome data in terms of cellular metabolism and its biogeochemical impact.

54 ENVIRONMENTAL SCIENCES↗

Using next generation transcriptome sequencing to predict an ectomycorrhizal metablome.

Mycorrhizae, symbiotic interactions between soil fungi and tree roots, are ubiquitous in terrestrial ecosystems. The fungi contribute phosphorous, nitrogen and mobilized nutrients from organic matter in the soil and in return the fungus receives photosynthetically-derived carbohydrates. This union of plant and fungal metabolisms is the mycorrhizal metabolome. Understanding this symbiotic relationship at a molecular level provides important contributions to the understanding of forest ecosystems and global carbon cycling. We generated next generation short-read transcriptomic sequencing data from fully-formed ectomycorrhizae between Laccaria bicolor and aspen (Populus tremuloides) roots. The transcriptomic data was used to identify statistically significantly expressed gene models using a bootstrap-style approach, and these expressed genes were mapped to specific metabolic pathways. Integration of expressed genes that code for metabolic enzymes and the set of expressed membrane transporters generates a predictive model of the ectomycorrhizal metabolome. The generated model of mycorrhizal metabolome predicts that the specific compounds glycine, glutamate, and allantoin are synthesized by L. bicolor and that these compounds or their metabolites may be used for the benefit of aspen in exchange for the photosynthetically-derived sugars fructose and glucose. The analysis illustrates an approach to generate testable biological hypotheses to investigate the complex molecular interactions that drive ectomycorrhizal symbiosis. These models are consistent with experimental environmental data and provide insight into the molecular exchange processes for organisms in this complex ecosystem. The method used here for predicting metabolomic models of mycorrhizal systems from deep RNA sequencing data can be generalized and is broadly applicable to transcriptomic data derived from complex systems.

54 ENVIRONMENTAL SCIENCES↗