Sequencing our way to more accurate community abundance

Over the last two decades, there has been a huge increase in our understanding of microbial diversity, structure and composition enabled by high-throughput sequencing technologies. Yet, it is unclear how the number of sequences translates to the number of cells or species within the community. In some cases, additional observational data may be required to ensure relative abundance patterns from sequence reads are biologically meaningful. The goal of DNA-based methods for biodiversity assessments is to obtain robust community abundance data, simultaneously, from environmental samples. In this issue of Molecular Ecology Resources, Pierella Karlusich et al. (2022) describe a new method for quantifying phytoplankton cell abundance. Using Tara Oceans data sets, the authors propose the photosynthetic gene psbO for reporting accurate relative abundance of the entire phytoplankton community from metagenomic data. The authors demonstrate higher correlations with traditional optical methods (including microscopy and flow cytometry), using their new method, improving upon molecular abundance assessments using multicopy marker genes. Furthermore, to facilitate application of their approach, the authors curated a psbO gene database for accessible taxonomic queries. This is an important step towards improving species abundance estimates from molecular data and eventually reporting of absolute species abundance, enhancing our understanding of community dynamics.     

Quantification of microbial communities in near-surface and deeply buried marine sediments on the Peru continental margin using real-time PCR

Deeply buried marine sediments harbour a large fraction of all prokaryotes on Earth but it is still unknown which phylogenetic and physiological microbial groups dominate the deep biosphere. In this study real-time PCR allowed a comparative quantitative microbial community analysis in near-surface and deeply buried marine sediments from the Peru continental margin. The 16S rRNA gene copy numbers of prokaryotes and Bacteria were almost identical with a maximum of 10(8)-10(10) copies cm(-3) in the near-surface sediments. Archaea exhibited one to three orders of magnitude lower 16S rRNA gene copy numbers. The 18S rRNA gene of Eukarya was always at least three orders of magnitude less abundant than the 16S rRNA gene of prokaryotes. The 16S rRNA gene of the Fe(III)- and Mn(IV)-reducing bacterial family Geobacteraceae and the dissimilatory (bi)sulfite reductase gene (dsrA) of sulfate-reducing prokaryotes were abundant with 10(6)-10(8) copies cm(-3) in near-surface sediments but showed lower numbers and an irregular distribution in the deep sediments. The copy numbers of all genes decreased with sediment depth exponentially. The depth gradients were steeper for the gene copy numbers than for numbers of total prokaryotes (acridine orange direct counts), which reflects the ongoing degradation of the high-molecular-weight DNA with sediment age and depth. The occurrence of eukaryotic DNA also suggests DNA preservation in the deeply buried sediments.     

Community Structure and Function of Planktonic Crenarchaeota: Changes with Depth in the South China Sea

Marine Crenarchaeota represent a widespread and abundant microbial group in marine ecosystems. Here, we investigated the abundance, diversity, and distribution of planktonic Crenarchaeota in the epi-, meso-, and bathypelagic zones at three stations in the South China Sea (SCS) by analysis of crenarchaeal 16S rRNA gene, ammonia monooxygenase gene amoA involved in ammonia oxidation, and biotin carboxylase gene accA putatively involved in archaeal CO₂ fixation. Quantitative PCR analyses indicated that crenarchaeal amoA and accA gene abundances varied similarly with archaeal and crenarchaeal 16S rRNA gene abundances at all stations, except that crenarchaeal accA genes were almost absent in the epipelagic zone. Ratios of the crenarchaeal amoA gene to 16S rRNA gene abundances decreased ~2.6 times from the epi- to bathypelagic zones, whereas the ratios of crenarchaeal accA gene to marine group I crenarchaeal 16S rRNA gene or to crenarchaeal amoA gene abundances increased with depth, suggesting that the metabolism of Crenarchaeota may change from the epi- to meso- or bathypelagic zones. Denaturing gradient gel electrophoresis profiling of the 16S rRNA genes revealed depth partitioning in archaeal community structures. Clone libraries of crenarchaeal amoA and accA genes showed two clusters: the “shallow” cluster was exclusively derived from epipelagic water and the “deep” cluster was from meso- and/or bathypelagic waters, suggesting that niche partitioning may take place between the shallow and deep marine Crenarchaeota. Overall, our results show strong depth partitioning of crenarchaeal populations in the SCS and suggest a shift in their community structure and ecological function with increasing depth.     

Quantification of mcrA by fluorescent PCR in methanogenic and methanotrophic microbial communities

A quantitative fluorogenic PCR method for detecting methanogenic and methanotrophic orders was established using a refined primer set for the methyl coenzyme M reductase subunit A gene (mcrA). The method developed was applied to several microbial communities in which diversity and abundance of methanogens or anaerobic methanotrophs (ANMEs) was identified by 16S rRNA gene clone analysis, and strong correlations between the copy numbers of mcrA with those of archaeal 16S rRNA genes in the communities were observed. The assay can be applied to detecting and assessing the abundance of methanogens and/or ANMEs in anoxic environments that could not be detected by 16S rRNA gene sequence analyses.     

Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: Accelerated sequence evolution,altered promoter structure,and tRNA pseudogenes

The nucleotide sequence of a 7.4 kb region containing the entire plastid ribosomal RNA operon of the nongreen parasitic plant Epifagus virginiana has been determined. Analysis of the sequence indicates that all four rRNA genes are intact and almost certainly functional. In contrast, the split genes for tRNA^Ile and tRNA^Ala present in the 16S-23S rRNA spacer region have become pseudogenes, and deletion upstream of the 16S rRNA gene has removed a tRNA^Val gene and most of the promoter region for the rRNA operon. The rate of nucleotide substitution in 16S and 23S rRNAs is several times higher in Epifagus than in tobacco, a related photosynthetic plant. Possible reasons for this, including relaxed translational constraints, are discussed.     

Microbial diversity of active layer and permafrost in an acidic wetland from the Canadian High Arctic

The abundance and structure of archaeal and bacterial communities from the active layer and the associated permafrost of a moderately acidic (pH < 5.0) High Arctic wetland (Axel Heiberg Island, Nunavut, Canada) were investigated using culture- and molecular-based methods. Aerobic viable cell counts from the active layer were ～100-fold greater than those from the permafrost (2.5 × 10(5) CFU·(g soil dry mass)(-1)); however, a greater diversity of isolates were cultured from permafrost, as determined by 16S rRNA gene sequencing. Isolates from both layers demonstrated growth characteristics of a psychrotolerant, halotolerant, and acidotolerant community. Archaea constituted 0.1% of the total 16S rRNA gene copy number and, in the 16S rRNA gene clone library, predominantly (71% and 95%) consisted of Crenarchaeota related to Group I. 1b. In contrast, bacterial communities were diverse (Shannon's diversity index, H = ～4), with Acidobacteria constituting the largest division of active layer clones (30%) and Actinobacteria most abundant in permafrost (28%). Direct comparisons of 16S rRNA gene sequence data highlighted significant differences between the bacterial communities of each layer, with the greatest differences occurring within Actinobacteria. Comparisons of 16S rRNA gene sequences with those from other Arctic permafrost and cold-temperature wetlands revealed commonly occurring taxa within the phyla Chloroflexi, Acidobacteria, and Actinobacteria (families Intrasporangiaceae and Rubrobacteraceae).     

Novel groups of Gammaproteobacteria catalyse sulfur oxidation and carbon fixation in a coastal, intertidal sediment

The oxidation of hydrogen sulfide is essential to sulfur cycling in marine habitats. However, the role of microbial sulfur oxidation in marine sediments and the microorganisms involved are largely unknown, except for the filamentous, mat‐forming bacteria. In this study we explored the diversity, abundance and activity of sulfur‐oxidizing prokaryotes (SOP) in sulfidic intertidal sediments using 16S rRNA and functional gene sequence analyses, fluorescence in situ hybridization (FISH) and microautoradiography. The 16S rRNA gene analysis revealed that distinct clades of uncultured Gammaproteobacteria are important SOP in the tidal sediments. This was supported by the dominance of gammaproteobacterial sequences in clone libraries of genes encoding the reverse dissimilatory sulfite reductase (rDSR) and the adenosine phosphosulfate reductase (APR). Numerous sequences of all three genes grouped with uncultured autotrophic SOP. Accordingly, Gammaproteobacteria accounted for 40–70% of all ¹⁴CO₂‐incorporating cells in surface sediments as shown by microautoradiography. Furthermore, phylogenetic analysis of all three genes consistently suggested a discrete population of SOP that was most closely related to the sulfur‐oxidizing endosymbionts of the tubeworm Oligobrachia spp. FISH showed that members of this population (WS‐Gam209 group) were abundant, reaching up to 1.3 × 10⁸ cells ml^?1 (4.6% of all cells). Approximately 25% of this population incorporated CO₂, consistent with a chemolithoautotrophic metabolism most likely based on sulfur oxidation. Thus, we hypothesize that novel, gammaproteobacterial SOP attached to sediment particles may play a more important role for sulfide removal and primary production in marine sediments than previously assumed.     

Salinity controls soil microbial community structure and function in coastal estuarine wetlands

Soil salinity acts as a critical environmental filter on microbial communities, but the consequences for microbial diversity and biogeochemical processes are poorly understood. Here, we characterized soil bacterial communities and microbial functional genes in a coastal estuarine wetland ecosystem across a gradient (~5 km) ranging from oligohaline to hypersaline habitats by applying the PCR-amplified 16S rRNA (rRNA) genes sequencing and microarray-based GeoChip 5.0 respectively. Results showed that saline soils in marine intertidal and supratidal zone exhibited higher bacterial richness and Faith's phylogenetic diversity than that in the freshwater-affected habitats. The relative abundance of taxa assigned to Gammaproteobacteria, Bacteroidetes and Firmicutes was higher with increasing salinity, while those affiliated with Acidobacteria, Chloroflexi and Cyanobacteria were more prevalent in wetland soils with low salinity. The phylogenetic inferences demonstrated the deterministic role of salinity filtering on the bacterial community assembly processes. The abundance of most functional genes involved in carbon degradation and nitrogen cycling correlated negatively with salinity, except for the hzo gene, suggesting a critical role of the anammox process in tidal affected zones. Overall, the salinity filtering effect shapes the soil bacterial community composition, and soil salinity act as a critical inhibitor in the soil biogeochemical processes in estuary ecosystems.     

A comparison of primer sets for detecting 16S rRNA and hydrazine oxidoreductase genes of anaerobic ammonium-oxidizing bacteria in marine sediments

Published polymerase chain reaction primer sets for detecting the genes encoding 16S rRNA gene and hydrazine oxidoreductase (hzo) in anammox bacteria were compared by using the same coastal marine sediment samples. While four previously reported primer sets developed to detect the 16S rRNA gene showed varying specificities between 12% and 77%, an optimized primer combination resulted in up to 98% specificity, and the recovered anammox 16S rRNA gene sequences were >95% sequence identical to published sequences from anammox bacteria in the Candidatus “Scalindua” group. Furthermore, four primer sets used in detecting the hzo gene of anammox bacteria were highly specific (up to 92%) and efficient, and the newly designed primer set in this study amplified longer hzo gene segments suitable for phylogenetic analysis. The optimized primer set for the 16S rRNA gene and the newly designed primer set for the hzo gene were successfully applied to identify anammox bacteria from marine sediments of aquaculture zone, coastal wetland, and deep ocean where the three ecosystems form a gradient of anthropogenic impact. Results indicated a broad distribution of anammox bacteria with high niche-specific community structure within each marine ecosystem.     

rRNA Operon Copy Number Reflects Ecological Strategies of Bacteria

Although natural selection appears to favor the elimination of gene redundancy in prokaryotes, multiple copies of each rRNA-encoding gene are common on bacterial chromosomes. Despite this conspicuous deviation from single-copy genes, no phenotype has been consistently associated with rRNA gene copy number. We found that the number of rRNA genes correlates with the rate at which phylogenetically diverse bacteria respond to resource availability. Soil bacteria that formed colonies rapidly upon exposure to a nutritionally complex medium contained an average of 5.5 copies of the small subunit rRNA gene, whereas bacteria that responded slowly contained an average of 1.4 copies. In soil microcosms pulsed with the herbicide 2,4-dichlorophenoxyacetic acid (2,4-D), indigenous populations of 2,4-D-degrading bacteria with multiple rRNA genes ( = 5.4) became dominant, whereas populations with fewer rRNA genes ( = 2.7) were favored in unamended controls. These findings demonstrate phenotypic effects associated with rRNA gene copy number that are indicative of ecological strategies influencing the structure of natural microbial communities.     

Microbial rRNA gene expression and co‐occurrence profiles associate with biokinetics and elemental composition in full‐scale anaerobic digesters

This study examined whether the abundance and expression of microbial 16S rRNA genes were associated with elemental concentrations and substrate conversion biokinetics in 20 full‐scale anaerobic digesters, including seven municipal sewage sludge (SS) digesters and 13 industrial codigesters. SS digester contents had higher methane production rates from acetate, propionate and phenyl acetate compared to industrial codigesters. SS digesters and industrial codigesters were distinctly clustered based on their elemental concentrations, with higher concentrations of NH₃‐N, Cl, K and Na observed in codigesters. Amplicon sequencing of 16S rRNA genes and reverse‐transcribed 16S rRNA revealed divergent grouping of microbial communities between mesophilic SS digesters, mesophilic codigesters and thermophilic digesters. Higher intradigester distances between Archaea 16S rRNA and rRNA gene profiles were observed in mesophilic codigesters, which also had the lowest acetate utilization biokinetics. Constrained ordination showed that microbial rRNA and rRNA gene profiles were significantly associated with maximum methane production rates from acetate, propionate, oleate and phenyl acetate, as well as concentrations of NH₃‐N, Fe, S, Mo and Ni. A co‐occurrence network of rRNA gene expression confirmed the three main clusters of anaerobic digester communities based on active populations. Syntrophic and methanogenic taxa were highly represented within the subnetworks, indicating that obligate energy‐sharing partnerships play critical roles in stabilizing the digester microbiome. Overall, these results provide new evidence showing that different feed substrates associate with different micronutrient compositions in anaerobic digesters, which in turn may influence microbial abundance, activity and function.     

Microbial Community Structure in Three Deep-Sea Carbonate Crusts

Carbonate crusts in marine environments can act as sinks for carbon dioxide. Therefore, understanding carbonate crust formation could be important for understanding global warming. In the present study, the microbial communities of three carbonate crust samples from deep-sea mud volcanoes in the eastern Mediterranean were characterized by sequencing 16S ribosomal RNA (rRNA) genes amplified from DNA directly retrieved from the samples. In combination with the mineralogical composition of the crusts and lipid analyses, sequence data were used to assess the possible role of prokaryotes in crust formation. Collectively, the obtained data showed the presence of highly diverse communities, which were distinct in each of the carbonate crusts studied. Bacterial 16S rRNA gene sequences were found in all crusts and the majority was classified as α-, γ-, and δ- Proteobacteria. Interestingly, sequences of Proteobacteria related to Halomonas and Halovibrio sp., which can play an active role in carbonate mineral formation, were present in all crusts. Archaeal 16S rRNA gene sequences were retrieved from two of the crusts studied. Several of those were closely related to archaeal sequences of organisms that have previously been linked to the anaerobic oxidation of methane (AOM). However, the majority of archaeal sequences were not related to sequences of organisms known to be involved in AOM. In combination with the strongly negative δ ¹³C values of archaeal lipids, these results open the possibility that organisms with a role in AOM may be more diverse within the Archaea than previously suggested. Different communities found in the crusts could carry out similar processes that might play a role in carbonate crust formation.     

Quantifying bacterial population dynamics in compost using 16S rRNA gene probes

Composting provides a dynamic setting for studying ecological topics such as succession, competition, and community stability in a relatively short period of time. This study used hierarchical small sub-unit-based rRNA gene probes to quantify the change in the relative abundance of phylogenetic groups common to compost in laboratory scale reactors. Bacterial 16S rRNA gene targets accounted for only 37% of all small subunit (SSU) rRNA genes initially, but increased to a maximum of 83% of the total at 84 h. The sum of rRNA genes detected using probes specific to Pseudomonas and low-G+C Gram-positive rRNA genes represented between 16% and 87% of the total. The lack of hybridization to the taxon-specific probes was most pronounced between 36 h and 60 h, when the pH was between 4.6 and 4.8. During this period the relative abundance of taxon-specific gene targets accounted for only 17–33% of the total bacterial rRNA gene targets. Pseudomonas-type 16S rRNA genes were the most abundant of the groups measured until 72 h. Those genes had their highest relative abundance at 12 h (78% of bacterial rRNA genes; 30% of all rRNA genes), after which time their relative abundance began to decline as the temperature increased. Prior to 72 h, 16S rRNA genes from low-G+C Gram-positive bacteria (LGC-GPB) represented less than 7% of the bacterial rRNA genes. However, by 84 h the relative abundance of LGC-GPB and Bacillus rRNA genes had increased to 60% and 18% of the bacterial rRNA gene targets, respectively (50% and 15% of all rRNA genes, respectively).     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

RiboFR-Seq: a novel approach to linking 16S rRNA amplicon profiles to metagenomes

16S rRNA amplicon analysis and shotgun metagenome sequencing are two main culture-independent strategies to explore the genetic landscape of various microbial communities. Recently, numerous studies have employed these two approaches together, but downstream data analyses were performed separately, which always generated incongruent or conflict signals on both taxonomic and functional classifications. Here we propose a novel approach, RiboFR-Seq (Ribosomal RNA gene flanking region sequencing), for capturing both ribosomal RNA variable regions and their flanking protein-coding genes simultaneously. Through extensive testing on clonal bacterial strain, salivary microbiome and bacterial epibionts of marine kelp, we demonstrated that RiboFR-Seq could detect the vast majority of bacteria not only in well-studied microbiomes but also in novel communities with limited reference genomes. Combined with classical amplicon sequencing and shotgun metagenome sequencing, RiboFR-Seq can link the annotations of 16S rRNA and metagenomic contigs to make a consensus classification. By recognizing almost all 16S rRNA copies, the RiboFR-seq approach can effectively reduce the taxonomic abundance bias resulted from 16S rRNA copy number variation. We believe that RiboFR-Seq, which provides an integrated view of 16S rRNA profiles and metagenomes, will help us better understand diverse microbial communities.     

Springtime dynamics,productivity and activity of prokaryotes in two Arctic fjords

In the Kongsfjorden–Krossfjorden system (Spitsbergen), increasing temperatures enhance glacier melting and concomitant intrusion of freshwater. These altered conditions affect the timing, intensity, and composition of the phytoplankton spring bloom in Kongsfjorden; yet, the effects on prokaryotes (bacteria and archaea) are not well understood. The aim of this study was to examine springtime prokaryote communities in both fjords as a function of hydrographic and phytoplankton variability. Prokaryote community composition was studied in two consecutive years by molecular fingerprinting of the 16S rRNA gene. In addition, we measured bacterial abundance, productivity (³H-Leucine uptake), and single-cell activity using catalyzed reporter deposition fluorescence in situ hybridization combined with microautoradiography. Differences in bacterial and archaeal communities were found between Kongsfjorden and Krossfjorden. Furthermore, an increase in productivity, abundance, and proportion of active bacterial cells was observed during the course of spring. Bacteroidetes were the most abundant bacterial group among the assessed taxa in both Kongsfjorden and Krossfjorden. Multivariate analysis of the microbial community fingerprints revealed a strong temporal shaping of both the bacterial and archaeal communities in addition to a spatial separation between the two fjords. A significant part of the observed bacterial variation could be explained by cyanobacterial biomass, as deduced from pigment analysis, and by phosphate concentration. Archaea were mainly controlled by abiotic factors. We speculate that the bacterial response to hydrographic changes and glacier meltwater is mediated through shifts in phytoplankton abundance and composition, whereas archaea are directly influenced by abiotic environmental variables.     

Coral reef invertebrate microbiomes correlate with the presence of photosymbionts

Coral reefs provide habitat for an array of marine invertebrates that host symbiotic microbiomes. Photosynthetic symbionts including Symbiodinium dinoflagellates and diatoms potentially influence the diversity of their host-associated microbiomes by releasing carbon-containing photosynthates and other organic compounds that fuel microbial metabolism. Here we used 16S ribosomal RNA (rRNA) gene amplicon pyrosequencing to characterise the microbiomes of 11 common Great Barrier Reef marine invertebrate species that host photosynthetic symbionts and five taxa in which they are absent. The presence of photosynthetic symbionts influenced the composition but not the species richness, evenness and phylogenetic diversity of invertebrate-associated microbiomes. Invertebrates without photosynthetic symbionts were dominated by Alphaproteobacteria, whereas those hosting photosynthetic symbionts were dominated by Gammaproteobacteria. Interestingly, many microbial species from photosymbiont-bearing invertebrates, including Oceanospirillales spp., Alteromonas spp., Pseudomonas spp., Halomonas spp., are implicated in the metabolism of dimethylsulfoniopropionate (DMSP). DMSP is produced in high concentrations by photosynthetic dinoflagellates and is involved in climate regulation by facilitating cloud formation. Microbiomes correlated with host taxa and replicate individuals from most sampled species grouped in distance-based redundancy analysis of retrieved 16S rRNA gene sequences. This study highlights the complex nature of invertebrate holobionts and confirms the importance of photosynthetic symbionts in structuring marine invertebrate bacterial communities.     

New PCR primers targeting hydrazine synthase and cytochrome c biogenesis proteins in anammox bacteria

PCR primers targeting genes encoding the two proteins of anammox bacteria, hydrazine synthase and cytochrome c biogenesis protein, were designed and tested in this study. Three different ecotypes of samples, namely ocean sediments, coastal wetland sediments, and wastewater treatment plant (WWTP) samples, were used to assess the primer efficiency and the community structures of anammox bacteria retrieved by 16S ribosomal RNA (rRNA) and the functional genes. Abundances of hzsB gene of anammox bacteria in South China Sea (SCS) samples were significantly correlated with 16S rRNA gene by qPCR method. And hzsB and hzsC gene primer pair hzsB364f-hzsB640r and hzsC745f-hzsC862r in combination with anammox bacterial 16S rRNA gene primers were recommended for quantifying anammox bacteria. Congruent with 16S rRNA gene-based community study, functional gene hzsB could also delineate the coastal-ocean distributing pattern, and seawater depth was positively associated with the diversity and abundance of anammox bacteria from shallow- to deep-sea. Both hzsC and ccsA genes could differentiate marine samples between deep and shallow groups of the Scalindua sp. clades. As for WWTP samples, non-Scalindua anammox bacteria reflected by hzsB, hzsC, ccsA, and ccsB gene-based libraries showed a similar distribution pattern with that by 16S rRNA gene. NH₄ ⁺ and NH₄ ⁺/Σ(NO₃ ⁻ + NO₂ ⁻) positively correlated with anammox bacteria gene diversity, but organic matter contents correlated negatively with anammox bacteria gene diversity in SCS. Salinity was positively associated with diversity indices of hzsC and ccsB gene-harboring anammox bacteria communities and could potentially differentiate the distribution patterns between shallow- and deep-sea sediment samples. SCS surface sediments harbored considerably diverse community of Scalindua. A new Mai Po clade representing coastal estuary wetland anammox bacteria group based on 16S rRNA gene phylogeny is proposed. Existence of anammox bacteria within wider coverage of genera in Mai Po wetland indicates this unique niche is very complex, and species of anammox bacteria are niche-specific with different physiological properties towards substrates competing and chemical tolerance capability.

     

Unusually high numbers of ribosomal RNA genes in copepods (Arthropoda: Crustacea) and their relationship to genome size.

We report on copy numbers of 18S ribosomal RNA genes in three species of copepods (Crustacea: Copepoda), two of which possess an unusual arrangement in which 5S genes are included within the 18S-5.8S-28S repeat unit. Slot blots of genomic and standard DNA were hybridized with an 18S rRNA gene probe constructed from one of the marine species and hybridization was quantified using chemiluminescence. Diploid 18S rRNA gene copy numbers are estimated as ca. 15 300 and 33 500 in the marine species Calanus finmarchicus (13.0 pg DNA in 2C adult nuclei) and C. glacialis (24.2 pg DNA), respectively, and ca. 840 and 730 in two freshwater populations of Mesocyclops edax (both ca. 3 pg DNA) from Virginia and Nova Scotia, respectively. The roughly proportional relationship between 2C somatic nuclear DNA contents and rRNA gene copy number in the sibling species C. finmarchicus and C. glacialis may reflect polytenic replication of entire genomes during abrupt speciation events. Copy numbers may also reflect differential losses during embryonic chromatin diminution.