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 487 records · Page 27

Guide-bound structures of an RNA-targeting A-cleaving CRISPR–Cas13a enzyme

CRISPR adaptive immune systems protect bacteria from infections by deploying CRISPR RNA (crRNA)-guided enzymes to recognize and cut foreign nucleic acids. Type VI-A CRISPR–Cas systems include the Cas13a enzyme, an RNA-activated RNase capable of crRNA processing and single-stranded RNA degradation upon target-transcript binding. Here we present the 2.0-Å resolution crystal structure of a crRNA-bound Lachnospiraceae bacterium Cas13a (LbaCas13a), representing a recently discovered Cas13a enzyme subtype. This structure and accompanying biochemical experiments define the Cas13a catalytic residues that are directly responsible for crRNA maturation. In addition, the orientation of the foreign-derived target-RNA-specifying sequence in the protein interior explains the conformational gating of Cas13a nuclease activation. These results describe how Cas13a enzymes generate functional crRNAs and how catalytic activity is blocked before target-RNA recognition, with implications for both bacterial immunity and diagnostic applications.

59 BASIC BIOLOGICAL SCIENCES↗

Sequence polymorphism in an insect RNA virus field population: A snapshot from a single point in space and time reveals stochastic differences among and within individual hosts

Population structure of Homalodisca coagulata Virus-1 (HoCV-1) among and within field-collected insects sampled from a single point in space and time was examined. Polymorphism in complete consensus sequences among single-insect isolates was dominated by synonymous substitutions. The mutant spectrum of the C2 helicase region within each single-insect isolate was unique and dominated by nonsynonymous singletons. Bootstrapping was used to correct the within-isolate nonsynonymous:synonymous arithmetic ratio (N:S) for RT-PCR error, yielding an N:S value ~one log-unit greater than that of consensus sequences. Probability of all possible single-base substitutions for the C2 region predicted N:S values within 95% confidence limits of the corrected within-isolate N:S when the only constraint imposed was viral polymerase error bias for transitions over transversions. These results indicate that bottlenecks coupled with strong negative/purifying selection drive consensus sequences toward neutral sequence space, and that most polymorphism within single-insect isolates is composed of newly-minted mutations sampled prior to selection. -- Highlights: •Sampling protocol minimized differential selection/history among isolates. •Polymorphism among consensus sequences dominated by negative/purifying selection. •Within-isolate N:S ratio corrected for RT-PCR error by bootstrapping. •Within-isolate mutant spectrum dominated by new mutations yet to undergo selection.

60 APPLIED LIFE SCIENCES↗

CRISPR–Cas9 bends and twists DNA to read its sequence

In bacterial defense and genome editing applications, the CRISPR-associated protein Cas9 searches millions of DNA base pairs to locate a 20-nucleotide, guide RNA-complementary target sequence that abuts a protospacer-adjacent motif (PAM). Target capture requires Cas9 to unwind DNA at candidate sequences using an unknown ATP-independent mechanism. Here we show that Cas9 sharply bends and undertwists DNA on PAM binding, thereby flipping DNA nucleotides out of the duplex and toward the guide RNA for sequence interrogation. Cryogenic-electron microscopy (cryo-EM) structures of Cas9–RNA–DNA complexes trapped at different states of the interrogation pathway, together with solution conformational probing, reveal that global protein rearrangement accompanies formation of an unstacked DNA hinge. Bend-induced base flipping explains how Cas9 ‘reads’ snippets of DNA to locate target sites within a vast excess of nontarget DNA, a process crucial to both bacterial antiviral immunity and genome editing. This mechanism establishes a physical solution to the problem of complementarity-guided DNA search and shows how interrogation speed and local DNA geometry may influence genome editing efficiency.

59 BASIC BIOLOGICAL SCIENCES↗

The western red cedar (Thuja plicata) 8-8' DIRIGENT family displays diverse expression patterns and conserved monolignol coupling specificity

The isolation and characterization of a multigene family of the first class of dirigent proteins (namely that mainly involved in 8-8' coupling leading to (+)-pinoresinol in this case) is reported, this comprising of nine western red cedar (Thuja plicata) DIRIGENT genes (DIR1-9) of 72-99.5% identity to each other. Their corresponding cDNA clones had coding regions for 180-183 amino acids with each having a predicted molecular mass of ca. 20 kDa including the signal peptide. Real time-PCR established that the DIRIGENT isovariants were differentially expressed during growth and development of T. plicata (P < 0.05). The phylogenetic relationships and the rates and patterns of nucleotide substitution suggest that the DIRIGENT gene may have evolved via paralogous expansion at an early stage of vascular plant diversification. Thereafter, western red cedar paralogues have maintained an high homogeneity presumably via a concerted evolutionary mode. This, in turn, is assumed to be the driving force for the differential formation of 8-8'-linked pinoresinol derived (poly)lignans in the needles, stems, bark and branches, as well as for massive accumulation of 8-8'-linked plicatic acid-derived (poly)lignans in heartwood.

Non-NASA Center↗

Strand-Specific RNA-Seq Analyses of Fruiting Body Development in Coprinopsis cinerea

We report that the basidiomycete fungus Coprinopsis cinerea is an important model system for multicellular development. Fruiting bodies of C. cinerea are typical mushrooms, which can be produced synchronously on defined media in the laboratory. To investigate the transcriptome in detail during fruiting body development, high-throughput sequencing (RNA-seq) was performed using cDNA libraries strand-specifically constructed from 13 points (stages/tissues) with two biological replicates. The reads were aligned to 14,245 predicted transcripts, and counted for forward and reverse transcripts. Differentially expressed genes (DEGs) between two adjacent points and between vegetative mycelium and each point were detected by Tag Count Comparison (TCC). To validate RNA-seq data, expression levels of selected genes were compared using RPKM values in RNA-seq data and qRT-PCR data, and DEGs detected in microarray data were examined in MA plots of RNA-seq data by TCC. We discuss events deduced from GO analysis of DEGs. In addition, we uncovered both transcription factor candidates and antisense transcripts that are likely to be involved in developmental regulation for fruiting.

59 BASIC BIOLOGICAL SCIENCES↗

Identification and Removal of Contaminant Sequences From Ribosomal Gene Databases: Lessons From the Census of Deep Life

Earth’s subsurface environment is one of the largest, yet least studied, biomes on Earth, and many questions remain regarding what microorganisms are indigenous to the subsurface. Through the activity of the Census of Deep Life (CoDL) and the Deep Carbon Observatory, an open access 16S ribosomal RNA gene sequence database from diverse subsurface environments has been compiled. However, due to low quantities of biomass in the deep subsurface, the potential for incorporation of contaminants from reagents used during sample collection, processing, and/or sequencing is high. Thus, to understand the ecology of subsurface microorganisms (i.e., the distribution, richness, or survival), it is necessary to minimize, identify, and remove contaminant sequences that will skew the relative abundances of all taxa in the sample. In this meta-analysis, we identify putative contaminants associated with the CoDL dataset, recommend best practices for removing contaminants from samples, and propose a series of best practices for subsurface microbiology sampling. The most abundant putative contaminant genera observed, independent of evenness across samples, were Propionibacterium, Aquabacterium, Ralstonia, and Acinetobacter. While the top five most frequently observed genera were Pseudomonas, Propionibacterium, Acinetobacter, Ralstonia, and Sphingomonas. The majority of the most frequently observed genera (high evenness) were associated with reagent or potential human contamination. Additionally, in DNA extraction blanks, we observed potential archaeal contaminants, including methanogens, which have not been discussed in previous contamination studies. Such contaminants would directly affect the interpretation of subsurface molecular studies, as methanogenesis is an important subsurface biogeochemical process. Utilizing previously identified contaminant genera, we found that ~27% of the total dataset were identified as contaminant sequences that likely originate from DNA extraction and DNA cleanup methods. Thus, controls must be taken at every step of the collection and processing procedure when working with low biomass environments such as, but not limited to, portions of Earth’s deep subsurface. Taken together, we stress that the CoDL dataset is an incredible resource for the broader research community interested in subsurface life, and steps to remove contamination derived sequences must be taken prior to using this dataset.

16S rRNA↗

Bacterial Diversity in Two Neonatal Intensive Care Units (NICUs)

Infants in Neonatal Intensive Care Units (NICUs) are particularly susceptible to opportunistic infection. Infected infants have high mortality rates, and survivors often suffer life-long neurological disorders. The causes of many NICU infections go undiagnosed, and there is debate as to the importance of inanimate hospital environments (IHEs) in the spread of infections. Here, we used culture-independent next-generation sequencing to survey bacterial diversity in two San Diego NICUs and to track the sources of microbes in these environments. Thirty IHE samples were collected from two Level-Three NICU facilities. We extracted DNA from these samples and amplified the bacterial small subunit (16S) ribosomal RNA gene sequence using ‘universal’ barcoded primers. The purified PCR products were pooled into a single reaction for pyrosequencing, and the data were analyzed using QIIME. On average, we detected 93+/-39 (mean +/- standard deviation) bacterial genera per sample in NICU IHEs. Many of the bacterial genera included known opportunistic pathogens, and many were skin-associated (e.g., Propionibacterium). In one NICU, we also detected fecal coliform bacteria (Enterobacteriales) in a high proportion of the surface samples. Comparison of these NICU-derived sequences to previously published high-throughput 16S rRNA amplicon studies of other indoor environments (offices, restrooms and healthcare facilities), as well as human- and soil-associated environments, found the majority of the NICU samples to be similar to typical building surface and air samples, with the notable exception of the IHEs which were dominated by Enterobacteriaceae. Our findings provide evidence that NICU IHEs harbor a high diversity of human-associated bacteria and demonstrate the potential utility of molecular methods for identifying and tracking bacterial diversity in NICUs.

59 BASIC BIOLOGICAL SCIENCES↗

Circularization of 23S rRNA but not 16S rRNA within archaeal ribosomes

Background Processing of archaeal 16S and 23S rRNAs is believed to involve excision of individual rRNAs from polycistronic precursors, circularization of excised rRNAs, and re-linearization before the incorporation into ribosomes. However, all the knowledge is derived from several isolated species, leaving open the possibility that different processes may occur in other archaeal groups. Results Here, we investigate rRNAs from diverse and mostly uncultivated archaea. Sequencing of total cellular RNA from eight phylum-level lineages indicates that archaeal circular 23S rRNA transcript abundances vastly exceed those of linear counterparts, and linear versions are often undetectable. As the majority of rRNAs derive from mature ribosomes, the data suggest that ribosomes contain circular 23S rRNAs. Thus, we directly sequence RNA extracted from isolated ribosomes of a model archaeon, Methanosarcina acetivorans, and confirm that the 23S rRNAs in the ribosomes are circular. Structural modeling places the 5′ and 3′ ends of the linear precursors of archaeal 23S rRNAs in close proximity to form a GNRA tetraloop (in which N is A, C, G, or U and R is A or G), consistent with their existence as circular molecules. We also confirm the existence of circular 16S rRNA intermediates in transcriptomes of most archaea, yet a circular form is not evident in some distinct archaeal groups, suggesting that certain archaea do not circularize 16S rRNA during processing. Conclusions Our findings uncover unexpected variations in the processing required to generate mature rRNAs and the conformation of functional molecules in archaeal ribosomes.

Archaea↗

Estimating the Timing of Early Simian-Human Immunodeficiency Virus Infections: a Comparison between Poisson Fitter and BEAST

Many HIV prevention strategies are currently under consideration where it is highly informative to know the study participants’ times of infection. These can be estimated using viral sequence data sampled early in infection. However, there are several scenarios that, if not addressed, can skew timing estimates. These include multiple transmitted/founder (TF) viruses, APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like)-mediated mutational enrichment, and recombination. Here, we suggest a pipeline to identify these problems and resolve the biases that they introduce. We then compare two modeling strategies to obtain timing estimates from sequence data. The first, Poisson Fitter (PF), is based on a Poisson model of random accumulation of mutations relative to the TF virus (or viruses) that established the infection. The second uses a coalescence-based phylogenetic strategy as implemented in BEAST. The comparison is based on timing predictions using plasma viral RNA (cDNA) sequence data from 28 simian-human immunodeficiency virus (SHIV)-infected animals for which the exact day of infection is known. In this particular setting, based on nucleotide sequences from samples obtained in early infection, the Poisson method yielded more accurate, more precise, and unbiased estimates for the time of infection than did the explored implementations of BEAST. The inference of the time of infection is a critical parameter in testing the efficacy of clinical interventions in protecting against HIV-1 infection. For example, in clinical trials evaluating the efficacy of passively delivered antibodies (Abs) for preventing infections, accurate time of infection data are essential for discerning levels of the Abs required to confer protection, given the natural Ab decay rate in the human body. In such trials, genetic sequences from early in the infection are regularly sampled from study participants, generally prior to immune selection, when the viral population is still expanding and genetic diversity is low. In this particular setting of early viral growth, the Poisson method is superior to the alternative approach based on coalescent methods. This approach can also be applied in human vaccine trials, where accurate estimates of infection times help ascertain if vaccine-elicited immune protection wanes over time.

59 BASIC BIOLOGICAL SCIENCES↗

RNA Splicing Events in Circulation Distinguish Individuals With and Without New-onset Type 1 Diabetes

Context: Alterations in RNA splicing may influence protein isoform diversity that contributes to or reflects the pathophysiology of certain diseases. Whereas specific RNA splicing events in pancreatic islets have been investigated in models of inflammation in vitro, how RNA splicing in the circulation correlates with or is reflective of type 1 diabetes (T1D) disease pathophysiology in humans remains unexplored. Objective: To use machine learning to investigate if alternative RNA splicing events differ between individuals with and without new-onset T1D and to determine if these splicing events provide insight into T1D pathophysiology. Methods: RNA deep sequencing was performed on whole blood samples from 2 independent cohorts: a training cohort consisting of 12 individuals with new-onset T1D and 12 age- and sex-matched nondiabetic controls and a validation cohort of the same size and demographics. Machine learning analysis was used to identify specific isoforms that could distinguish individuals with T1D from controls. Results: Distinct patterns of RNA splicing differentiated participants with T1D from unaffected controls. Notably, certain splicing events, particularly involving retained introns, showed significant association with T1D. Machine learning analysis using these splicing events as features from the training cohort demonstrated high accuracy in distinguishing between T1D subjects and controls in the validation cohort. Gene Ontology pathway enrichment analysis of the retained intron category showed evidence for a systemic viral response in T1D subjects. Conclusion: Alternative RNA splicing events in whole blood are significantly enriched in individuals with new-onset T1D and can effectively distinguish these individuals from unaffected controls. Further, our findings also suggest that RNA splicing profiles offer the potential to provide insights into disease pathogenesis.

60 APPLIED LIFE SCIENCES↗

Use of a Capture-Based Pathogen Transcript Enrichment Strategy for RNA-Seq Analysis of the Francisella Tularensis LVS Transcriptome during Infection of Murine Macrophages

Francisella tularensis is a zoonotic intracellular pathogen that is capable of causing potentially fatal human infections. Like all successful bacterial pathogens, F. tularensis rapidly responds to changes in its environment during infection of host cells, and upon encountering different microenvironments within those cells. This ability to appropriately respond to the challenges of infection requires rapid and global shifts in gene expression patterns. In this study, we use a novel pathogen transcript enrichment strategy and whole transcriptome sequencing (RNA-Seq) to perform a detailed characterization of the rapid and global shifts in F. tularensis LVS gene expression during infection of murine macrophages. We performed differential gene expression analysis on all bacterial genes at two key stages of infection: phagosomal escape, and cytosolic replication. By comparing the F. tularensis transcriptome at these two stages of infection to that of the bacteria grown in culture, we were able to identify sets of genes that are differentially expressed over the course of infection. This analysis revealed the temporally dynamic expression of a number of known and putative transcriptional regulators and virulence factors, providing insight into their role during infection. In addition, we identified several F. tularensis genes that are significantly up-regulated during infection but had not been previously identified as virulence factors. These unknown genes may make attractive therapeutic or vaccine targets.

59 BASIC BIOLOGICAL SCIENCES↗

Circular RNA profiling reveals that circCOL3A1-859267 regulate type I collagen expression in photoaged human dermal fibroblasts

Production of type I collagen declines is a main characteristic during photoaging, but the mechanism is still not fully understood. Circular RNAs (circRNAs) are a class of newly identified non-coding RNAs with regulatory potency by sequestering miRNAs like a sponge. It's more stable than linear RNAs, and would be a useful tool for regulation of gene expression. However, the role of circRNAs in collagen expression during photoaging is still unclear. Here we performed deep sequencing of RNA generated from UVA irradiated and no irradiated human dermal fibroblasts (HDFs) and identified 29 significantly differentially expressed circRNAs (fold change ≥ 1.5, P < 0.05), 12 circRNAs were up-regulated and 17 circRNAs were down-regulated.3 most differentially expressed circRNAs were verified by qRT-PCR and the down-regulated circCOL3A1-859267 exhibited the most significantly altered in photoaged HDFs. Overexpression of circCOL3A1-859267 inhibited UVA-induced decrease of type I collagen expression and silencing of it reduced type I collagen intensity. Via a bioinformatic method, 44 miRNAs were predicted to binding with circCOL3A1-859267, 5 of them have been confirmed or predicted to interact with type I collagen. This study show that circCOL3A1-859267 regulate type I collagen expression in photoaged HDFs, suggesting it may be a novel target for interfering photoaging.

60 APPLIED LIFE SCIENCES↗

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net).

59 BASIC BIOLOGICAL SCIENCES↗

An anti-silencer– and SATB1-dependent chromatin hub regulates Rag1 and Rag2 gene expression during thymocyte development

Rag1 and Rag2 gene expression in CD4 + CD8 + double-positive (DP) thymocytes depends on the activity of a distant anti-silencer element (ASE) that counteracts the activity of an intergenic silencer. However, the mechanistic basis for ASE activity is unknown. Here, we show that the ASE physically interacts with the distant Rag1 and Rag2 gene promoters in DP thymocytes, bringing the two promoters together to form an active chromatin hub. Moreover, we show that the ASE functions as a classical enhancer that can potently activate these promoters in the absence of the silencer or other locus elements. In thymocytes lacking the chromatin organizer SATB1, we identified a partial defect in Tcra gene rearrangement that was associated with reduced expression of Rag1 and Rag2 at the DP stage. SATB1 binds to the ASE and Rag promoters, facilitating inclusion of Rag2 in the chromatin hub and the loading of RNA polymerase II to both the Rag1 and Rag2 promoters. Our results provide a novel framework for understanding ASE function and demonstrate a novel role for SATB1 as a regulator of Rag locus organization and gene expression in DP thymocytes.

histone H3 lysine 27 trimethylation↗

High expression of diffuse panbronchiolitis critical region 1 gene promotes cell proliferation, migration and invasion in pancreatic ductal adenocarcinoma

Highlights: • DPCR1 is overexpressed in PDAC and correlates with prognosis of PDAC patients. • Inhibition of DPCR1 suppresses PDAC cell proliferation in vitro and in vivo. • Inhibition of DPCR1 suppresses PDAC cell migration and invasion. • DPCR1 may promote tumor progression by regulating NF-κB pathway. Diffuse panbronchiolitis critical region 1 (DPCR1) is located in the major histocompatibility complex (MHC) class I. It was reported to be downregulated in invasive pituitary adenoma compared with that in non-invasive tumors, but upregulated in the precursor of gastric carcinogenesis. However, the direct effect of DPCR1 on cancer cells has rarely been reported, and the role DPCR1 in pancreatic ductal adenocarcinoma (PDAC) remains unclear. The clinical sample validation and public data analysis of the present study demonstrated that DPCR1 was upregulated markedly in PDAC and this high expression was negatively correlated with the patient prognosis. Functionally, knocking down DPCR1 in PDAC cell lines inhibited cell proliferation, migration and invasion in vitro. Tumor xenograft experiments further showed that suppression of DPCR1 inhibited tumor growth in vivo. In addition, the results of RNA deep sequencing and qRT-PCR assay showed that DPCR1 participated in PADC progression by regulating nuclear factor-kappa B signaling pathway, suggesting that it might be a novel oncogene in tumor progression and a potential therapeutic target in PDAC as well.

60 APPLIED LIFE SCIENCES↗

Identification of candidate host-specificity genes in Exserohilum turcicum using comparative genomics and transcriptomics

Abstract Exserohilum turcicum causes northern corn leaf blight and sorghum leaf blight. While the same species cause disease in both crops, the strains are host-specific. Here, we report the sequence and de novo annotated assemblies of one sorghum- and one maize-specific E. turcicum strain. The strains were sequenced using the PacBio Sequel II system. The total genome length for both assemblies was between 44 and 45 Mb with N50 of ∼2.5 Mb. Ninety-eight percent of the Benchmarking Universal Single-Copy Orthologs (BUSCO) for both assemblies had complete status. The estimated number of genes was 11,762 and 12,029 in the sorghum- and maize-specific isolates, respectively. Funannotate, EffectorP, SignalP, and transcriptome data were used to create functional annotation of each genome. The whole-genome comparison identified ten large-scale inversions and three translocations between the maize- and sorghum-specific strains, along with homologous genes and gene duplications. RNA was sequenced from the maize- and sorghum-specific isolate 10 days post-inoculation in maize and sorghum and from axenic cultures. Gene expression data from planta and axenic growth experiments were compared for each strain. Candidate host-specificity genes were identified by combining results from whole-genome comparison, synteny analysis, gene annotations, and transcriptome data. Overall, this study identified several candidate host-specificity genes that provide insights into E. turcicum interaction with its hosts.

Krone, Mara J. (ORCID:0000000159006624)↗

High Bacterial Diversity of Biological Soil Crusts in Water Tracks over Permafrost in the High Arctic Polar Desert

In this paper we report the bacterial diversity of biological soil crusts (biocrusts) inhabiting polar desert soils at the northern land limit of the Arctic polar region (83° 05 N). Employing pyrosequencing of bacterial 16S rRNA genes this study demonstrated that these biocrusts harbor diverse bacterial communities, often as diverse as temperate latitude communities. The effect of wetting pulses on the composition of communities was also determined by collecting samples from soils outside and inside of permafrost water tracks, hill slope flow paths that drain permafrost-affected soils. The intermittent flow regime in the water tracks was correlated with altered relative abundance of phylum level taxonomic bins in the bacterial communities, but the alterations varied between individual sampling sites. Bacteria related to the Cyanobacteria and Acidobacteria demonstrated shifts in relative abundance based on their location either inside or outside of the water tracks. Among cyanobacterial sequences, the proportion of sequences belonging to the family Oscillatoriales consistently increased in relative abundance in the samples from inside the water tracks compared to those outside. Acidobacteria showed responses to wetting pulses in the water tracks, increasing in abundance at one site and decreasing at the other two sites. Subdivision 4 acidobacterial sequences tended to follow the trends in the total Acidobacteria relative abundance, suggesting these organisms were largely responsible for the changes observed in the Acidobacteria. Finally, taken together, these data suggest that the bacterial communities of these high latitude polar biocrusts are diverse but do not show a consensus response to intermittent flow in water tracks over high Arctic permafrost.

59 BASIC BIOLOGICAL SCIENCES↗