Metatranscriptomics of N2-fixing cyanobacteria in the Amazon River plume

Global marine microbial diversity and its potential in bioprospecting

The past two decades has witnessed a remarkable increase in the number of microbial genomes retrieved from marine systems1,2. However, it has remained challenging to translate this marine genomic diversity into biotechnological and biomedical applications3,4. Here we recovered 43,191 bacterial and archaeal genomes from publicly available marine metagenomes, encompassing a wide range of diversity with 138 distinct phyla, redefining the upper limit of marine bacterial genome size and revealing complex trade-offs between the occurrence of CRISPR-Cas systems and antibiotic resistance genes. In silico bioprospecting of these marine genomes led to the discovery of a novel CRISPR-Cas9 system, ten antimicrobial peptides, and three enzymes that degrade polyethylene terephthalate. In vitro experiments confirmed their effectiveness and efficacy. This work provides evidence that global-scale sequencing initiatives advance our understanding of how microbial diversity has evolved in the oceans and is maintained, and demonstrates how such initiatives can be sustainably exploited to advance biotechnology and biomedicine.

Cascade damming impacts on microbial mediated nitrogen cycling in rivers

Rivers display essential role in nitrogen (N) cycling in terrestrial and aquatic ecosystems, but now they are suffering from damming worldwide, especially from cascade damming. Despite of the importance of microorganisms in biogeochemical nutrient cycling, little attention has been paid to microbial functional biogeography under damming disturbances. Here, the Geochip microarray was applied to investigate the microbial mediated N cycling across the single-dammed Yarlung Tsangpo-Brahmaputra River and the cascade-dammed Lancang-Mekong River in southwest China. Our results showed that the N cycling processes (nitrogen fixation, ammonification, denitrification, nitrification and anammox) were stimulated in reservoirs in both rivers and the enhancement was inversely coupled with hydraulic retention time, but the recovery of N-cycling gene abundance in downstream of dam was intervened by cascade damming. Moreover, N-cycling gene composition was significantly altered in the single-dammed river, while no remarkable change was found in the cascade-dammed reaches. However, different from the unvaried gene composition, cascade damming intervened the recovery of N-cycling gene flow connectivity and resulted in the continuous decrease of connectivity in cascade damming reaches. In addition, in the single-dammed river, nutrients were the important drivers for variation in gene abundance, while they did not influence gene composition. Meanwhile, the abundance and composition of N-cycling genes in the cascade-dammed river were both significantly correlated to geographical parameters and water physical characteristics. Therefore, our study has vital implications for anticipating microbial functional response and biogeochemical feedback to ongoing cascade damming, contributing to the protection of river ecosystems under river regulation.

Genome-Resolved Metagenomics and Metatranscriptomics Reveal that Aquificae Dominates Arsenate Reduction in Tengchong Geothermal Springs

Elevated arsenic (As) is common in geothermal springs, shaping the evolution of As metabolism genes and As transforming microbes. Herein, genome-level microbial metabolisms and As cycling strategies in Tengchong geothermal springs were demonstrated for the first time based on metagenomic and metatranscriptomic analyses. Sulfur cycling was dominated by Aquificae oxidizing thiosulfate via the sox system, fueling the respiration and carbon dioxide fixation processes. Arsenate reduction via arsC [488.63 ± 271.60 transcripts per million (TPM)] and arsenite efflux via arsB (442.98 ± 284.81 TPM) were the primary detoxification pathway, with most genes and transcripts contributed by the members in phylum Aquificae. A complete arsenotrophic cycle was also transcriptionally active as evidenced by the detection of aioA transcripts and arrA transcript reads mapped onto metagenome-assembled genomes (MAGs) affiliated with Crenarchaeota. MAGs affiliated with Aquificae had great potential of reducing arsenate via arsC and fixing nitrogen and carbon dioxide via nifDHK and reductive tricarboxylic acid (rTCA) cycle, respectively. Aquificae's arsenate reduction potential via arsC was observed for the first time at the transcriptional level. This study expands the diversity of the arsC-based arsenate-reducing community and highlights the importance of Aquificae to As biogeochemistry.

Depth-dependent variability of biological nitrogen fixation and diazotrophic communities in mangrove sediments

BACKGROUND

Nitrogen-fixing prokaryotes (diazotrophs) contribute substantially to nitrogen input in mangrove sediments, and their structure and nitrogen fixation rate (NFR) are significantly controlled by environmental conditions. Despite the well-known studies on diazotrophs in surficial sediments, the diversity, structure, and ecological functions of diazotrophic communities along environmental gradients of mangrove sediment across different depths are largely unknown. Here, we investigated how biological nitrogen fixation varied with the depth of mangrove sediments from the perspectives of both NFR and diazotrophic communities.

RESULTS

Through acetylene reduction assay, nifH gene amplicon and metagenomic sequencing, we found that the NFR increased but the diversity of diazotrophic communities decreased with the depth of mangrove sediments. The structure of diazotrophic communities at different depths was largely driven by salinity and exhibited a clear divergence at the partitioning depth of 50 cm. Among diazotrophic genera correlated with NFR, Agrobacterium and Azotobacter were specifically enriched at 50-100 cm sediments, while Anaeromyxobacter, Rubrivivax, Methylocystis, Dickeya, and Methylomonas were more abundant at 0-50 cm. Consistent with the higher NFR, metagenomic analysis demonstrated the elevated abundance of nitrogen fixation genes (nifH/D/K) in deep sediments, where nitrification genes (amoA/B/C) and denitrification genes (nirK and norB) became less abundant. Three metagenome-assembled genomes (MAGs) of diazotrophs from deep mangrove sediments indicated their facultatively anaerobic and mixotrophic lifestyles as they contained genes for low-oxygen-dependent metabolism, hydrogenotrophic respiration, carbon fixation, and pyruvate fermentation.

CONCLUSIONS

This study demonstrates the depth-dependent variability of biological nitrogen fixation in terms of NFR and diazotrophic communities, which to a certain extent relieves the degree of nitrogen limitation in deep mangrove sediments. Video Abstract.

Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters

Rivers connect the carbon cycle in land with that in aquatic ecosystems by transporting and transforming terrestrial organic matter (TeOM). The Amazon River receives huge loads of TeOM from the surrounding rainforest, promoting a substantial microbial heterotrophic activity and consequently, CO₂ outgassing. In the Amazon River, microbes degrade up to 55% of the lignin present in the TeOM. Yet, the main microbial genomes involved in TeOM degradation were unknown. Here, we characterize 51 population genomes (PGs) representing some of the most abundant microbes in the Amazon River deriving from 106 metagenomes. The 51 reconstructed PGs are among the most abundant microbes in the Amazon River, and 53% of them are not able to degrade TeOM. Among the PGs capable of degrading TeOM, 20% were exclusively cellulolytic, while the others could also oxidize lignin. The transport and consumption of lignin oxidation byproducts seemed to be decoupled from the oxidation process, being apparently performed by different groups of microorganisms. By connecting the genomic features of abundant microbes in the Amazon River with the degradation machinery of TeOM, we suggest that a complex microbial consortium could explain the quick turnover of TeOM previously observed in this ecosystem.

Transcriptional patterns identify resource controls on the diazotroph Trichodesmium in the Atlantic and Pacific oceans

The N₂-fixing cyanobacterium Trichodesmium is intensely studied because of the control this organism exerts over the cycling of carbon and nitrogen in the low nutrient ocean gyres. Although iron (Fe) and phosphorus (P) bioavailability are thought to be major drivers of Trichodesmium distributions and activities, identifying resource controls on Trichodesmium is challenging, as Fe and P are often organically complexed and their bioavailability to a single species in a mixed community is difficult to constrain. Further, Fe and P geochemistries are linked through the activities of metalloenzymes, such as the alkaline phosphatases (APs) PhoX and PhoA, which are used by microbes to access dissolved organic P (DOP). Here we identified significant correlations between Trichodesmium-specific transcriptional patterns in the North Atlantic (NASG) and North Pacific Subtropical Gyres (NPSG) and patterns in Fe and P biogeochemistry, with the relative enrichment of Fe stress markers in the NPSG, and P stress markers in the NASG. We also observed the differential enrichment of Fe-requiring PhoX transcripts in the NASG and Fe-insensitive PhoA transcripts in the NPSG, suggesting that metalloenzyme switching may be used to mitigate Fe limitation of DOP metabolism in Trichodesmium. This trait may underpin Trichodesmium success across disparate ecosystems.

Capturing the genetic makeup of the active microbiome in situ

More than any other technology, nucleic acid sequencing has enabled microbial ecology studies to be complemented with the data volumes necessary to capture the extent of microbial diversity and dynamics in a wide range of environments. In order to truly understand and predict environmental processes, however, the distinction between active, inactive and dead microbial cells is critical. Also, experimental designs need to be sensitive toward varying population complexity and activity, and temporal as well as spatial scales of process rates. There are a number of approaches, including single-cell techniques, which were designed to study in situ microbial activity and that have been successively coupled to nucleic acid sequencing. The exciting new discoveries regarding in situ microbial activity provide evidence that future microbial ecology studies will indispensably rely on techniques that specifically capture members of the microbiome active in the environment. Herein, we review those currently used activity-based approaches that can be directly linked to shotgun nucleic acid sequencing, evaluate their relevance to ecology studies, and discuss future directions.

Metatranscriptomic analysis of prokaryotic communities active in sulfur and arsenic cycling in Mono Lake, California, USA

This study evaluates the transcriptionally active, dissimilatory sulfur- and arsenic-cycling components of the microbial community in alkaline, hypersaline Mono Lake, CA, USA. We sampled five depths spanning the redox gradient (10, 15, 18, 25 and 31 m) during maximum thermal stratification. We used custom databases to identify transcripts of genes encoding complex iron-sulfur molybdoenzyme (CISM) proteins, with a focus on arsenic (arrA, aioA and arxA) and sulfur cycling (dsrA, aprA and soxB), and assigned them to taxonomic bins. We also report on the distribution of transcripts related to the ars arsenic detoxification pathway. Transcripts from detoxification pathways were not abundant in oxic surface waters (10 m). Arsenic cycling in the suboxic and microaerophilic zones of the water column (15 and 18 m) was dominated by arsenite-oxidizing members of the Gammaproteobacteria most closely affiliated with Thioalkalivibrio and Halomonas, transcribing arxA. We observed a transition to arsenate-reducing bacteria belonging to the Deltaproteobacteria and Firmicutes transcribing arsenate reductase (arrA) in anoxic bottom waters of the lake (25 and 31 m). Sulfur cycling at 15 and 18 m was dominated by Gammaproteobacteria (Thioalkalivibrio and Thioalkalimicrobium) oxidizing reduced S species, with a transition to sulfate-reducing Deltaproteobacteria at 25 and 31 m. Genes related to arsenic and sulfur oxidation from Thioalkalivibrio were more highly transcribed at 15 m relative to other depths. Our data highlight the importance of Thioalkalivibrio to arsenic and sulfur biogeochemistry in Mono Lake and identify new taxa that appear capable of transforming arsenic.

Bacterial Biogeography across the Amazon River-Ocean Continuum

Spatial and temporal patterns in microbial biodiversity across the Amazon river-ocean continuum were investigated along ∼675 km of the lower Amazon River mainstem, in the Tapajós River tributary, and in the plume and coastal ocean during low and high river discharge using amplicon sequencing of 16S rRNA genes in whole water and size-fractionated samples (0.2–2.0 μm and >2.0 μm). River communities varied among tributaries, but mainstem communities were spatially homogeneous and tracked seasonal changes in river discharge and co-varying factors. Co-occurrence network analysis identified strongly interconnected river assemblages during high (May) and low (December) discharge periods, and weakly interconnected transitional assemblages in September, suggesting that this system supports two seasonal microbial communities linked to river discharge. In contrast, plume communities showed little seasonal differences and instead varied spatially tracking salinity. However, salinity explained only a small fraction of community variability, and plume communities in blooms of diatom-diazotroph assemblages were strikingly different than those in other high salinity plume samples. This suggests that while salinity physically structures plumes through buoyancy and mixing, the composition of plume-specific communities is controlled by other factors including nutrients, phytoplankton community composition, and dissolved organic matter chemistry. Co-occurrence networks identified interconnected assemblages associated with the highly productive low salinity near-shore region, diatom-diazotroph blooms, and the plume edge region, and weakly interconnected assemblages in high salinity regions. This suggests that the plume supports a transitional community influenced by immigration of ocean bacteria from the plume edge, and by species sorting as these communities adapt to local environmental conditions. Few studies have explored patterns of microbial diversity in tropical rivers and coastal oceans. Comparison of Amazon continuum microbial communities to those from temperate and arctic systems suggest that river discharge and salinity are master variables structuring a range of environmental conditions that control bacterial communities across the river-ocean continuum.

Piecewise Structural Equation Model (SEM) Disentangles the Environmental Conditions Favoring Diatom Diazotroph Associations (DDAs) in the Western Tropical North Atlantic (WTNA)

Diatom diazotroph associations (DDAs) are important components in the world’s oceans, especially in the western tropical north Atlantic (WTNA), where blooms have a significant impact on carbon and nitrogen cycling. However, drivers of their abundances and distribution patterns remain unknown. Here, we examined abundance and distribution patterns for two DDA populations in relation to the Amazon River (AR) plume in the WTNA. Quantitative PCR assays, targeting two DDAs (het-1 and het-2) by their symbiont’s nifH gene, served as input in a piecewise structural equation model (SEM). Collections were made during high (spring 2010) and low (fall 2011) flow discharges of the AR. The distributions of dissolved nutrients, chlorophyll-a, and DDAs showed coherent patterns indicative of areas influenced by the AR. A symbiotic Hemiaulus hauckii-Richelia (het-2) bloom (>10⁶ cells L^-1) occurred during higher discharge of the AR and was coincident with mesohaline to oceanic (30–35) sea surface salinities (SSS), and regions devoid of dissolved inorganic nitrogen (DIN), low concentrations of both DIP (>0.1 μmol L^-1) and Si (>1.0 μmol L^-1). The Richelia (het-1) associated with Rhizosolenia was only present in 2010 and at lower densities (10-1.76 × 10⁵nifH copies L^-1) than het-2 and limited to regions of oceanic SSS (>36). The het-2 symbiont detected in 2011 was associated with H. membranaceus (>10³nifH copies L^-1) and were restricted to regions with mesohaline SSS (31.8–34.3), immeasurable DIN, moderate DIP (0.1–0.60 μmol L^-1) and higher Si (4.19–22.1 μmol L^-1). The piecewise SEM identified a profound direct negative effect of turbidity on the het-2 abundance in spring 2010, while DIP and water turbidity had a more positive influence in fall 2011, corroborating our observations of DDAs at subsurface maximas. We also found a striking difference in the influence of salinity on DDA symbionts suggesting a niche differentiation and preferences in oceanic and mesohaline salinities by het-1 and het-2, respectively. The use of the piecewise SEM to disentangle the complex and concomitant hydrography of the WTNA acting on two biogeochemically relevant populations was novel and underscores its use to predict conditions favoring abundance and distributions of microbial populations.

Expression patterns of elemental cycling genes in the Amazon River Plume

Metatranscriptomics and metagenomics data sets benchmarked with internal standards were used to characterize the expression patterns for biogeochemically relevant bacterial and archaeal genes mediating carbon, nitrogen, phosphorus and sulfur uptake and metabolism through the salinity gradient of the Amazon River Plume. The genes were identified in 48 metatranscriptomic and metagenomic data sets summing to >500 million quality-controlled reads from six locations in the plume ecosystem. The ratio of transcripts per gene copy (a direct measure of expression made possible by internal standard additions) showed that the free-living bacteria and archaea exhibited only small changes in the expression levels of biogeochemically relevant genes through the salinity and nutrient zones of the plume. In contrast, the expression levels of genes in particle-associated cells varied over orders of magnitude among the stations, with the largest differences measured for genes mediating aspects of nitrogen cycling (nifH, amtB and amoA) and phosphorus acquisition (pstC, phoX and phoU). Taxa varied in their baseline gene expression levels and extent of regulation, and most of the spatial variation in the expression level could be attributed to changes in gene regulation after removing the effect of shifting taxonomic composition. We hypothesize that changes in microbial element cycling along the Amazon River Plume are largely driven by shifting activities of particle-associated cells, with most activities peaking in the mesohaline regions where N₂ fixation rates are elevated.

Variable depth distribution of Trichodesmium clades in the North Pacific Ocean

Populations of nitrogen-fixing cyanobacteria in the genus Trichodesmium are critical to ocean ecosystems, yet predicting patterns of Trichodesmium distribution and their role in ocean biogeochemistry is an ongoing challenge. This may, in part, be due to differences in the physiological ecology of Trichodesmium species, which are not typically considered independently in field studies. In this study, the abundance of the two dominant Trichodesmium clades (Clade I and Clade III) was investigated during a survey at Station ALOHA in the North Pacific Subtropical Gyre (NPSG) using a clade-specific qPCR approach. While Clade I dominated the Trichodesmium community, Clade III abundance was >50% in some NPSG samples, in contrast to the western North Atlantic where Clade III abundance was always <10%. Clade I populations were distributed down to depths >80 m, while Clade III populations were only observed in the mixed layer and found to be significantly correlated with depth and temperature. These data suggest active niche partitioning of Trichodesmium species from different clades, as has been observed in other cyanobacteria. Tracking the distribution and physiology of Trichodesmium spp. would contribute to better predictions of the physiological ecology of this biogeochemically important genus in the present and future ocean.

Use of a Filter Cartridge for Filtration of Water Samples and Extraction of Environmental DNA

Recent studies demonstrated the use of environmental DNA (eDNA) from fishes to be appropriate as a non-invasive monitoring tool. Most of these studies employed disk fiber filters to collect eDNA from water samples, although a number of microbial studies in aquatic environments have employed filter cartridges, because the cartridge has the advantage of accommodating large water volumes and of overall ease of use. Here we provide a protocol for filtration of water samples using the filter cartridge and extraction of eDNA from the filter without having to cut open the housing. The main portions of this protocol consists of 1) filtration of water samples (water volumes ≤4 L or >4 L); (2) extraction of DNA on the filter using a roller shaker placed in a preheated incubator; and (3) purification of DNA using a commercial kit. With the use of this and previously-used protocols, we perform metabarcoding analysis of eDNA taken from a huge aquarium tank (7,500 m3) with known species composition, and show the number of detected species per library from the two protocols as the representative results. This protocol has been developed for metabarcoding eDNA from fishes, but is also applicable to eDNA from other organisms.

Methodologies for probing the metatranscriptome of grassland soil

Metatranscriptomics provides an opportunity to identify active microbes and expressed genes in complex soil communities in response to particular conditions. Currently, there are a limited number of soil metatranscriptome studies to provide guidance for using this approach in this challenging matrix. Hence, we evaluated the technical challenges of applying soil metatranscriptomics to a highly diverse, low activity natural system. We used a non-targeted rRNA removal approach, duplex nuclease specific (DSN) normalization, to generate a metatranscriptomic library from field collected soil supporting a perennial grass, Miscanthus x giganteus (a biofuel crop), and evaluated its ability to provide insight into its active community members and their expressed protein-coding genes. We also evaluated various bioinformatics approaches for analyzing our soil metatranscriptome, including annotation of unassembled transcripts, de novo assembly, and aligning reads to known genomes. Further, we evaluated various databases for their ability to provide annotations for our metatranscriptome. Overall, our results emphasize that low activity, highly genetically diverse and relatively stable microbiomes, like soil, requires very deep sequencing to sample the transcriptome beyond the common core functions. We identified several key areas that metatranscriptomic analyses will benefit from including increased rRNA removal, assembly of short read transcripts, and more relevant reference bases while providing a priority set of expressed genes for functional assessment.

Patterns of Transcript Abundance of Eukaryotic Biogeochemically-Relevant Genes in the Amazon River Plume

The Amazon River has the largest discharge of all rivers on Earth, and its complex plume system fuels a wide array of biogeochemical processes, across a large area of the western tropical North Atlantic. The plume thus stimulates microbial processes affecting carbon sequestration and nutrient cycles at a global scale. Chromosomal gene expression patterns of the 2.0 to 156 μm size-fraction eukaryotic microbial community were investigated in the Amazon River Plume, generating a robust dataset (more than 100 million mRNA sequences) that depicts the metabolic capabilities and interactions among the eukaryotic microbes. Combining classical oceanographic field measurements with metatranscriptomics yielded characterization of the hydrographic conditions simultaneous with a quantification of transcriptional activity and identity of the community. We highlight the patterns of eukaryotic gene expression for 31 biogeochemically significant gene targets hypothesized to be valuable within forecasting models. An advantage to this targeted approach is that the database of reference sequences used to identify the target genes was selectively constructed and highly curated optimizing taxonomic coverage, throughput, and the accuracy of annotations. A coastal diatom bloom highly expressed nitrate transporters and carbonic anhydrase presumably to support high growth rates and enhance uptake of low levels of dissolved nitrate and CO₂. Diatom-diazotroph association (DDA: diatoms with nitrogen fixing symbionts) blooms were common when surface salinity was mesohaline and dissolved nitrate concentrations were below detection, and hence did not show evidence of nitrate utilization, suggesting they relied on ammonium transporters to aquire recently fixed nitrogen. These DDA blooms in the outer plume had rapid turnover of the photosystem D1 protein presumably caused by photodegradation under increased light penetration in clearer waters, and increased expression of silicon transporters as silicon became limiting. Expression of these genes, including carbonic anhydrase and transporters for nitrate and phosphate, were found to reflect the physiological status and biogeochemistry of river plume environments. These relatively stable patterns of eukaryotic transcript abundance occurred over modest spatiotemporal scales, with similarity observed in sample duplicates collected up to 2.45 km in space and 120 minutes in time. These results confirm the use of metatranscriptomics as a valuable tool to understand and predict microbial community function.

Experimental Identification of Small Non-Coding RNAs in the Model Marine Bacterium Ruegeria pomeroyi DSS-3

In oligotrophic ocean waters where bacteria are often subjected to chronic nutrient limitation, community transcriptome sequencing has pointed to the presence of highly abundant small RNAs (sRNAs). The role of sRNAs in regulating response to nutrient stress was investigated in a model heterotrophic marine bacterium Ruegeria pomeroyi grown in continuous culture under carbon (C) and nitrogen (N) limitation. RNAseq analysis identified 99 putative sRNAs. Sixty-nine were cis-encoded and located antisense to a presumed target gene. Thirty were trans-encoded and initial target prediction was performed computationally. The most prevalent functional roles of genes anti-sense to the cis-sRNAs were transport, cell-cell interactions, signal transduction, and transcriptional regulation. Most sRNAs were transcribed equally under both C and N limitation, and may be involved in a general stress response. However, 14 were regulated differentially between the C and N treatments and may respond to specific nutrient limitations. A network analysis of the predicted target genes of the R. pomeroyi cis-sRNAs indicated that they average fewer connections than typical protein-encoding genes, and appear to be more important in peripheral or niche-defining functions encoded in the pan genome.

Phosphite utilization by the globally important marine diazotroph Trichodesmium

Species belonging to the filamentous cyanobacterial genus Trichodesmium are responsible for a significant fraction of oceanic nitrogen fixation. The availability of phosphorus (P) likely constrains the growth of Trichodesmium in certain regions of the ocean. Moreover, Trichodesmium species have recently been shown to play a role in an emerging oceanic phosphorus redox cycle, further highlighting the key role these microbes play in many biogeochemical processes in the contemporary ocean. Here, we show that Trichodesmium erythraeum IMS101 can grow on the reduced inorganic compound phosphite as its sole source of P. The components responsible for phosphite utilization are identified through heterologous expression of the T. erythraeum IMS101 Tery_0365-0368 genes, encoding a putative adenosine triphosphate (ATP)-binding cassette transporter and nicotinamide adenine dinucleotide (NAD)-dependent dehydrogenase, in the model cyanobacteria Synechocystis sp. PCC6803. We demonstrate that only combined expression of both the transporter and the dehydrogenase enables Synechocystis to utilize phosphite, confirming the function of Tery_0365-0367 as a phosphite uptake system (PtxABC) and Tery_0368 as a phosphite dehydrogenase (PtxD). Our findings suggest that reported uptake of phosphite by Trichodesmium consortia in the field likely reflects an active biological process by Trichodesmium. These results highlight the diversity of phosphorus sources available to Trichodesmium in a resource-limited ocean.

Metagenomic and metatranscriptomic inventories of the lower Amazon River, May 2011

BACKGROUND

The Amazon River runs nearly 6500 km across the South American continent before emptying into the western tropical North Atlantic Ocean. In terms of both volume and watershed area, it is the world's largest riverine system, affecting elemental cycling on a global scale.

RESULTS

A quantitative inventory of genes and transcripts benchmarked with internal standards was obtained at five stations in the lower Amazon River during May 2011. At each station, metagenomes and metatranscriptomes were obtained in duplicate for two microbial size fractions (free-living, 0.2 to 2.0 μm; particle-associated, 2.0 to 297 μm) using 150 × 150 paired-end Illumina sequencing. Forty eight sample datasets were obtained, averaging 15 × 10(6) potential protein-encoding reads each (730 × 10(6) total). Prokaryotic metagenomes and metatranscriptomes were dominated by members of the phyla Actinobacteria, Planctomycetes, Betaproteobacteria, Verrucomicrobia, Nitrospirae, and Acidobacteria. The actinobacterium SCGC AAA027-L06 reference genome recruited the greatest number of reads overall, with this single bin contributing an average of 50 billion genes and 500 million transcripts per liter of river water. Several dominant taxa were unevenly distributed between the free-living and particle-associated size fractions, such as a particle-associated bias for reads binning to planctomycete Schlesneria paludicola and a free-living bias for actinobacterium SCGC AAA027-L06. Gene expression ratios (transcripts to gene copy ratio) increased downstream from Óbidos to Macapá and Belém, indicating higher per cell activity of Amazon River bacteria and archaea as river water approached the ocean.

CONCLUSION

This inventory of riverine microbial genes and transcripts, benchmarked with internal standards for full quantitation, provides an unparalleled window into microbial taxa and functions in the globally important Amazon River ecosystem.

Highly Iterated Palindromic Sequences (HIPs) and Their Relationship to DNA Methyltransferases

The sequence GCGATCGC (Highly Iterated Palindrome, HIP1) is commonly found in high frequency in cyanobacterial genomes. An important clue to its function may be the presence of two orphan DNA methyltransferases that recognize internal sequences GATC and CGATCG. An examination of genomes from 97 cyanobacteria, both free-living and obligate symbionts, showed that there are exceptional cases in which HIP1 is at a low frequency or nearly absent. In some of these cases, it appears to have been replaced by a different GC-rich palindromic sequence, alternate HIPs. When HIP1 is at a high frequency, GATC- and CGATCG-specific methyltransferases are generally present in the genome. When an alternate HIP is at high frequency, a methyltransferase specific for that sequence is present. The pattern of 1-nt deviations from HIP1 sequences is biased towards the first and last nucleotides, i.e., those distinguish CGATCG from HIP1. Taken together, the results point to a role of DNA methylation in the creation or functioning of HIP sites. A model is presented that postulates the existence of a GmeC-dependent mismatch repair system whose activity creates and maintains HIP sequences.