Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence analysis”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The hydroxyl-water megamaser connection. I. Water emission toward OH megamaser hosts

Questions surround the connection of luminous extragalactic masers to galactic processes. The observation that water and hydroxyl megamasers rarely coexist in the same galaxy has given rise to a hypothesis that the two species appear in different phases of nuclear activity. The detection of simultaneous hydroxyl and water megamaser emission toward IC694 has called this hypothesis into question, but, because many megamasers have not been surveyed for emission in the other molecule, it remains unclear whether IC694 occupies a narrow phase of galaxy evolution or whether the relationship between megamaser species and galactic processes is more complicated than previously believed. In this paper, we present results of a systematic search for 22 GHz water maser emission among OH megamaser hosts to identify additional objects hosting both megamasers. Our work roughly doubles the number of galaxies searched for emission in both molecules, which host at least one confirmed maser. We confirm with a high degree of confidence ($\gt 8\sigma $) the detection of water emission toward IIZw96, firmly establishing it as the second object to cohost both water and hydroxyl megamasers after IC694. We find high luminosity, narrow features in the water feature in IIZw96. All dual megamaser candidates appear in merging galaxy systems suggestive that megamasers that coexistance may signal a brief phase along the merger sequence. In conclusion, a statistical analysis of the results of our observations provide possible evidence for an exclusion of H 2 O kilomasers among OH megamaser hosts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Whither the genus Caldicellulosiruptor and the order Thermoanaerobacterales: phylogeny, taxonomy, ecology, and phenotype

The order Thermoanaerobacterales currently consists of fermentative anaerobic bacteria, including the genus Caldicellulosiruptor. Caldicellulosiruptor are represented by thirteen species; all, but one, have closed genome sequences. Interest in these extreme thermophiles has been motivated not only by their high optimal growth temperatures (≥70°C), but also by their ability to hydrolyze polysaccharides including, for some species, both xylan and microcrystalline cellulose. Caldicellulosiruptor species have been isolated from geographically diverse thermal terrestrial environments located in New Zealand, China, Russia, Iceland and North America. Evidence of their presence in other terrestrial locations is apparent from metagenomic signatures, including volcanic ash in permafrost. Here, phylogeny and taxonomy of the genus Caldicellulosiruptor was re-examined in light of new genome sequences. Based on genome analysis of 15 strains, a new order, Caldicellulosiruptorales, is proposed containing the family Caldicellulosiruptoraceae, consisting of two genera, Caldicellulosiruptor and Anaerocellum. Furthermore, the order Thermoanaerobacterales also was re-assessed, using 91 genome-sequenced strains, and should now include the family Thermoanaerobacteraceae containing the genera Thermoanaerobacter, Thermoanaerobacterium, Caldanaerobacter, the family Caldanaerobiaceae containing the genus Caldanaerobius, and the family Calorimonaceae containing the genus Calorimonas. A main outcome of ANI/AAI analysis indicates the need to reclassify several previously designated species in the Thermoanaerobacterales and Caldicellulosiruptorales by condensing them into strains of single species. Comparative genomics of carbohydrate-active enzyme inventories suggested differentiating phenotypic features, even among strains of the same species, reflecting available nutrients and ecological roles in their native biotopes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Trimming and Decontamination of Metagenomic Data can Significantly Impact Assembly and Binning Metrics, Phylogenomic and Functional Analysis

Background: Investigators using metagenomic sequencing to study microbiomes often trim and decontaminate reads without knowing their effect on downstream analyses. Objective: This study was designed to evaluate the impacts JGI trimming and decontamination procedures have on assembly and binning metrics, placement of MAGs into species trees, and functional profiles of MAGs extracted from complex rhizosphere metagenomes, as well as how more aggressive trimming impacts these binning metrics. Methods: Twenty-three Miscanthus x giganteus rhizosphere metagenomes were subjected to different combinations and thresholds of force, kmer, and quality trimming and decontamination using BBDuk. Reads were assembled and binned in KBase. Phylogenomic and statistical analyses were applied to evaluate the effects of trimming and decontamination on downstream analyses. Results: We found that JGI trimmed and decontaminated reads had significant impacts on assembly and binning metrics compared to raw reads, including significantly higher total contig counts, more contigs greater than 10k bp in length, and larger total lengths of raw assemblies compared to QC assemblies, and 2.0% lower average contamination of QC MAGs compared to raw MAGs. We also found that differences in the placement of MAGs in species trees increased with decreasing completeness and contamination thresholds. Furthermore, aggressive trimming (Q20) was found to significantly reduce MAG counts. Conclusion: Trimming and decontamination of metagenomics reads prior to assembly can change an investigator’s answer to the questions, “Who is there and what are they doing?” However, mild trimming and decontamination of metagenomic reads with high-quality scores are recommended for removing sample processing and sequencing artifacts.

Whitham, Jason M.↗
DOE OSTI2024

Metagenome-assembled genomes provide insight into the metabolic potential during early production of Hydraulic Fracturing Test Site 2 in the Delaware Basin

Demand for natural gas continues to climb in the United States, having reached a record monthly high of 104.9 billion cubic feet per day (Bcf/d) in November 2023. Hydraulic fracturing, a technique used to extract natural gas and oil from deep underground reservoirs, involves injecting large volumes of fluid, proppant, and chemical additives into shale units. This is followed by a “shut-in” period, during which the fracture fluid remains pressurized in the well for several weeks. The microbial processes that occur within the reservoir during this shut-in period are not well understood; yet, these reactions may significantly impact the structural integrity and overall recovery of oil and gas from the well. To shed light on this critical phase, we conducted an analysis of both pre-shut-in material alongside production fluid collected throughout the initial production phase at the Hydraulic Fracturing Test Site 2 (HFTS 2) located in the prolific Wolfcamp formation within the Permian Delaware Basin of west Texas, USA. Specifically, we aimed to assess the microbial ecology and functional potential of the microbial community during this crucial time frame. Prior analysis of 16S rRNA sequencing data through the first 35 days of production revealed a strong selection for a Clostridia species corresponding to a significant decrease in microbial diversity. Here, we performed a metagenomic analysis of produced water sampled on Day 33 of production. This analysis yielded three high-quality metagenome-assembled genomes (MAGs), one of which was a Clostridia draft genome closely related to the recently classified Petromonas tenebris. This draft genome likely represents the dominant Clostridia species observed in our 16S rRNA profile. Annotation of the MAGs revealed the presence of genes involved in critical metabolic processes, including thiosulfate reduction, mixed acid fermentation, and biofilm formation. These findings suggest that this microbial community has the potential to contribute to well souring, biocorrosion, and biofouling within the reservoir. Our research provides unique insights into the early stages of production in one of the most prolific unconventional plays in the United States, with important implications for well management and energy recovery.

natural gas↗
DOE OSTI2026

A Chemoselective and Stereodivergent Platform of Heme‐Nitrene Transferases to Access Chiral Aryl‐β‐Amino Esters and An Investigation of the Sequence‐Activity Landscape

Engineered biocatalysts can utilize nitrene precursors to access enantioenriched amination products, yet they have not been applied to produce valuable, enantiomerically enriched noncanonical β-amino esters. Current approaches to synthesizing β-amino acids rely on pre-oxidized precursors and multistep synthetic approaches involving various protecting groups. We engineered a platform of heme enzymes for stereoselective C–H bond amination of readily available carboxylic ester derivatives to install primary amines. A directed evolution campaign coupled with sequencing of over 1000 variants enabled us to develop engineered variants that use either O-pivaloylhydroxylamine triflic acid (PONT) or hydroxylamine hydrochloride (H 2 NOH∙HCl) as aminating reagents. An analysis of the resulting sequence–activity dataset revealed additional improvements that could be made to the final variant, highlighting the utility of sequencing data to guide future steps in directed evolution campaigns. Furthermore, the evolved nitrene transferases expand the scope of accessible chiral β-amino acid building blocks for peptidomimetic applications and provide new starting points for the design and synthesis of enantioenriched β-amino acid motifs.

amino ester building blocks↗
DOE OSTI2016

The standard operating procedure of the DOE-JGI Metagenome Annotation Pipeline (MAP v.4)

Abstract The DOE-JGI Metagenome Annotation Pipeline (MAP v.4) performs structural and functional annotation for metagenomic sequences that are submitted to the Integrated Microbial Genomes with Microbiomes (IMG/M) system for comparative analysis. The pipeline runs on nucleotide sequences provided via the IMG submission site. Users must first define their analysis projects in GOLD and then submit the associated sequence datasets consisting of scaffolds/contigs with optional coverage information and/or unassembled reads in fasta and fastq file formats. The MAP processing consists of feature prediction including identification of protein-coding genes, non-coding RNAs and regulatory RNAs, as well as CRISPR elements. Structural annotation is followed by functional annotation including assignment of protein product names and connection to various protein family databases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

ChIP-seq analysis of the σ E regulon of Salmonella enterica serovar typhimurium reveals new genes implicated in heat shock and oxidative stress response

The alternative sigma factor σ E functions to maintain bacterial homeostasis and membrane integrity in response to extracytoplasmic stress by regulating thousands of genes both directly and indirectly. The transcriptional regulatory network governed by σ E in Salmonella and E. coli has been examined using microarray, however a genome-wide analysis of σ E –binding sites inSalmonella has not yet been reported. We infected macrophages with Salmonella Typhimurium over a select time course. Using chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-seq), 31 σ E –binding sites were identified. Seventeen sites were new, which included outer membrane proteins, a quorum-sensing protein, a cell division factor, and a signal transduction modulator. The consensus sequence identified for σ E in vivo binding was similar to the one previously reported, except for a conserved G and A between the -35 and -10 regions. One third of the σ E –binding sites did not contain the consensus sequence, suggesting there may be alternative mechanisms by which σ E modulates transcription. By dissecting direct and indirect modes of σ E -mediated regulation, we found that σ E activates gene expression through recognition of both canonical and reversed consensus sequence. Lastly, new σ E regulated genes ( greA, luxS, ompA and ompX ) are shown to be involved in heat shock and oxidative stress responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A Pilot Study of Whole-Blood Transcriptomic Analysis to Identify Genes Associated with Repetitive Low-Level Blast Exposure in Career Breachers

Repetitive low-level blast exposure is one of the major occupational health concerns among US military service members and law enforcement. This study seeks to identify gene expression using microRNA and RNA sequencing in whole-blood samples from experienced breachers and unexposed controls. We performed experimental RNA sequencing using Illumina’s HiSeq 2500 Sequencing System, and microRNA analysis using NanoString Technology nCounter miRNA expression panel in whole-blood total RNA samples from 15 experienced breachers and 14 age-, sex-, and race-matched unexposed controls. We identified 10 significantly dysregulated genes between experienced breachers and unexposed controls, with FDR corrected <0.05: One upregulated gene, LINC00996 (long intergenic non-protein coding RNA 996); and nine downregulated genes, IGLV3-16 (immunoglobulin lambda variable 3-16), CD200 (CD200 molecule), LILRB5 (leukocyte immunoglobulin-like receptor B5), ZNF667-AS1 (ZNF667 antisense RNA 1), LMOD1 (leiomodin 1), CNTNAP2 (contactin-associated protein 2), EVPL (envoplakin), DPF3 (double PHD fingers 3), and IGHV4-34 (immunoglobulin heavy variable 4-34). The dysregulated gene expressions reported here have been associated with chronic inflammation and immune response, suggesting that these pathways may relate to the risk of lasting neurological symptoms following high exposures to blast over a career.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Single sample resolution of rare microbial dark matter in a marine invertebrate metagenome

Direct, untargeted sequencing of environmental samples (metagenomics) and de novo genome assembly enable the study of uncultured and phylogenetically divergent organisms. However, separating individual genomes from a mixed community has often relied on the differential-coverage analysis of multiple, deeply sequenced samples. In the metagenomic investigation of the marine bryozoan Bugula neritina, we uncovered seven bacterial genomes associated with a single B. neritina individual that appeared to be transient associates, two of which were unique to one individual and undetectable using certain “universal” 16S rRNA primers and probes. We recovered high quality genome assemblies for several rare instances of “microbial dark matter,” or phylogenetically divergent bacteria lacking genomes in reference databases, from a single tissue sample that was not subjected to any physical or chemical pre-treatment. One of these rare, divergent organisms has a small (593 kbp), poorly annotated genome with low GC content (20.9%) and a 16S rRNA gene with just 65% sequence similarity to the closest reference sequence. Lastly, our findings illustrate the importance of sampling strategy and de novo assembly of metagenomic reads to understand the extent and function of bacterial biodiversity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

MAGNETIC ACTIVITY ANALYSIS FOR A SAMPLE OF G-TYPE MAIN SEQUENCE KEPLER TARGETS

The variation of a stellar light curve owing to rotational modulation by magnetic features (starspots and faculae) on the star’s surface can be used to investigate the magnetic properties of the host star. In this paper, we use the periodicity and magnitude of the light-curve variation as two proxies to study the stellar magnetic properties for a large sample of G-type main sequence Kepler targets, for which the rotation periods were recently determined. By analyzing the correlation between the two magnetic proxies, it is found that: (1) the two proxies are positively correlated for most of the stars in our sample, and the percentages of negative, zero, and positive correlations are 4.27%, 6.81%, and 88.91%, respectively; (2) negative correlation stars cannot have a large magnitude of light-curve variation; and (3) with the increase of rotation period, the relative number of positive correlation stars decreases and the negative correlation one increases. These results indicate that stars with shorter rotation period tend to have positive correlation between the two proxies, and a good portion of the positive correlation stars have a larger magnitude of light-curve variation (and hence more intense magnetic activities) than negative correlation stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A kernel machine method for detecting higher order interactions in multimodal datasets: Application to schizophrenia

Technological advances are enabling us to collect multimodal datasets at an increasing depth and resolution while with decreasing labors. Understanding complex interactions among multimodal datasets, however, is challenging. Here in this study, we tested the interaction effect of multimodal datasets using a novel method called the kernel machine for detecting higher order interactions among biologically relevant multimodal data. Using a semiparametric method on a reproducing kernel Hilbert space, we formulated the proposed method as a standard mixed-effects linear model and derived a score-based variance component statistic to test higher order interactions between multimodal datasets. The method was evaluated using extensive numerical simulation and real data from the Mind Clinical Imaging Consortium with both schizophrenia patients and healthy controls. Our method identified 13-triplets that included 6 gene-derived SNPs, 10 ROIs, and 6 gene-specific DNA methylations that are correlated with the changes in hippocampal volume, suggesting that these triplets may be important for explaining schizophrenia-related neurodegeneration. The performance of the proposed method is compared with the following methods: test based on only first and first few principal components followed by multiple regression, and full principal component analysis regression, and the sequence kernel association test. With strong evidence (p-value ≤0.000001), the triplet (MAGI2, CRBLCrus1.L, FBXO28) is a significant biomarker for schizophrenia patients. This novel method can be applicable to the study of other disease processes, where multimodal data analysis is a common task.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A mutagenic analysis of NahE, a hydratase-aldolase in the naphthalene degradative pathway

NahE is a hydratase-aldolase that converts o-substituted trans-benzylidenepyruvates (H, OH, or CO 2 - ) to benzaldehyde, salicylaldehyde, or 2-carboxybenzaldehyde, respectively, and pyruvate. The enzyme is in a bacterial degradative pathway for naphthalene, which is a toxic and persistent environmental contaminant. Sequence, crystallographic, and mutagenic analysis identified the enzyme as a member of the N-acetylneuraminate lyase (NAL) subgroup in the aldolase superfamily. As such, it has a conserved lysine (Lys183) and tyrosine (Tyr155), for Schiff base formation, as well as a GXXGE motif for binding of the pyruvoyl carboxylate group. A crystal structure of the selenomethionine derivative of NahE shows these active site elements along with nearby residues that might be involved in the mechanism and/or specificity. Mutations of five active site amino acids (Thr65, Trp128, Tyr155, Asn157, and Asn281) were constructed and kinetic parameters measured in order to assess the effect(s) on catalysis. Here, the results show that the two Trp128 mutants (Phe and Tyr) have the least effect on catalysis, whereas amino acids with bulky side chains at Thr65 (Val) and Asn281 (Leu) have the greatest effect. Changing Tyr155 to Phe and Asn157 to Ala also hinders catalysis, and the effects fall in between these extremes. These observations are put into a structural context using a crystal structure of the Schiff base of the reaction intermediate. Trapping experiments with substrate, Na(CN)BH 3 , and wild type enzyme and selected mutants mostly paralleled the kinetic analysis, and identified two salicylaldehyde-modified lysines: the active site lysine (Lys183) and one outside the active site (Lys279). The latter could be responsible for the observed inhibition of NahE by salicylaldehyde. Together, the results provide new insights into the NahE-catalyzed reaction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Pond Crash Forensics: Presumptive identification of pond crash agents by next generation sequencing in replicate raceway mass cultures of Nannochloropsis salina

Productivity of algal mass culture can be severely reduced by contaminating organisms. It is, therefore, important to identify contaminants, determine their effect on productivity and, ultimately, develop countermeasures against such contamination. In this paper, we utilized microbiome analysis by second-generation sequencing of small subunit rRNA genes to characterize the predator and pathogen burden of open raceway cultures of Nannochloropsis salina. Samples were analyzed from replicate raceways before and after crashes. In one culture cycle, we identified two algivorous species, the rotifer Brachionus and gastrotrich Chaetonotus, the presence of which may have contributed to the loss of algal biomass. In the second culture cycle, the raceways were treated with hypochlorite in an unsuccessful attempt to interdict the crash. Finally, our analyses were shown to be an effective strategy for the identification of the biological contaminants and the characterization of intervention strategies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass)

Background: MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass ( Panicum virgatum ), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass. Results: Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession “Alamo”, upland accession “Dacotah”, and two upland/lowland F 1 hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors. Conclusion: Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The relation between lignin sequence and its 3D structure

Background: Lignin, the second most abundant biopolymer on earth, plays a major structural role in plants, conferring mechanical strength and regulating water conduction. Understanding the three-dimensional structure of lignin is important for fundamental reasons as well as engineering plants towards lignin valorization. Lignin lacks a specific primary sequence, making its average chemical composition the focus of most recent studies. However, it remains unclear whether the 3D structure of lignin molecules depends on their sequence. Methods: In this work, we performed all-atom molecular dynamics simulation of three S/G-lignin molecules with the same average composition but different sequence. Results: A detailed statistical analysis of the radius of gyration and relative shape anisotropy reveals that the lignin sequence has no statistically significant effect on the global three-dimensional structure. We found however, that homopolymers of C-lignin with the same molecular weight have smaller radii of gyration than S/G-lignin. We attribute this to lower hydroxyl content of C-lignin, which makes it more compact and rigid. Conclusions: The 3D structure of lignin is influenced by the overall content of monomeric units and interunit linkages and not by its precise primary sequence. General Significance: Lignin is assumed to not have a well-defined primary structure. The results presented here demonstrate there are no significant differences in the global 3D structure of lignin molecules with the same average composition but different primary sequence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Fluorescent Amplification of Next Generation Sequencing Libraries Tools (FA-NGS tools) v1

This is a software utility for calculating pooling volumes and graphing quality control data of next generation sequencing libraries amplified with fluorescent intercalating agents and monitored by real-time quantitative PCR (qPCR). The tool is used to calculate pooling volumes of individual next generation sequencing libraries based on experimentally determined fluorescence values. To do this, fluorescent values are used as a proxy for relative concentration to determine the ratio of transfer volumes. The output is a "transfer file" that specifies the transfer volume and destination well of a given library. The tool can also be used to assess the quality of individual next generation sequencing libraries by melting curve analysis. The output are small-multiples graphs of the melting curves of individual sequencing libraries according to the specified plate layout.

Garber, Megan↗
DOE OSTI2017

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and noncoding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Vertical Structure of Phyllosphere Fungal Communities in a Tropical Forest in Thailand Uncovered by High-Throughput Sequencing

Phyllosphere fungi harbor a tremendous species diversity and play important ecological roles. However, little is known about their distribution patterns within forest ecosystems. We examined how species diversity and community composition of phyllosphere fungi change along a vertical structure in a tropical forest in Thailand. Fungal communities in 144 leaf samples from 19 vertical layers (1.28–34.4 m above ground) of 73 plant individuals (27 species) were investigated by metabarcoding analysis using Ion Torrent sequencing. In total, 1,524 fungal operational taxonomic units (OTUs) were detected among 890,710 reads obtained from the 144 leaf samples. Taxonomically diverse fungi belonging to as many as 24 orders of Ascomycota and 21 orders of Basidiomycota were detected, most of which inhabited limited parts of the lowest layers closest to the forest floor. Species diversity of phyllosphere fungi was the highest in the lowest layers closest to the forest floor, decreased with increasing height, and lowest in the canopy; 742 and 55 fungal OTUs were detected at the lowest and highest layer, respectively. On the layers close to the forest floor, phyllosphere fungal communities were mainly composed of low frequency OTUs and largely differentiated among plant individuals. Conversely, in the canopy, fungal communities consisted of similar OTUs across plant individuals, and as many as 86.1%–92.7% of the OTUs found in the canopy (≥22 m above ground) were also distributed in the lower layers. Overall, our study showed the variability of phyllosphere fungal communities along the vertical gradient of plant vegetation and environmental conditions, suggesting the significance of biotic and abiotic variation for the species diversity of phyllosphere fungi.

59 BASIC BIOLOGICAL SCIENCES↗