Phylogeny and Rates of Molecular Evolution of Planktonic Foraminifera: SSU rDNA Sequences Compared to the Fossil Record

Planktonic foraminifera are marine protists, whose calcareous shells form oceanic sediments and are widely used for stratigraphic and paleoenvironmental analyses. The fossil record of planktonic foraminifera is compared here to their molecular phylogeny inferred from ribosomal DNA sequences. Eighteen partial SSU rDNA sequences from species representing all modern planktonic families (Globigerinidae, Hastigerinidae, Globorotaliidae, Candeinidae) were obtained and compared to seven sequences representing the major groups of benthic foraminifera. The phylogenetic analyses indicate a polyphyletic origin for the planktonic foraminifera. The Candeinidae, the Globorotaliidae, and the clade Globigerinidae + Hastigerinidae seem to have originated independently, at different epochs in the evolution of foraminifera. Inference of their relationships, however, is limited by substitution rates of heterogeneity. Rates of SSU rDNA evolution vary from 4.0 × 10⁻⁹ substitutions/site/year in the Globigerinidae to less than 1.0 × 10⁻⁹ substitutions/site/year in the Globorotaliidae. These variations may be related to different levels of adaptation to the planktonic mode of life. A clock-like evolution is observed among the Globigerinidae, for which molecular and paleontological data are congruent. Phylogeny of the Globorotaliidae is clearly biased by rapid rates of substitution in two species (G. truncatulinoides and G. menardii). Our study reveals differences in absolute rates of evolution at all taxonomic levels in planktonic foraminifera and demonstrates their effect on phylogenetic reconstructions. Received: 21 January 1997 / Accepted: 17 April 1997     

Molecular evidence for an independent origin of modern triserial planktonic foraminifera from benthic ancestors

Gallitellia vivans is the only Recent representative of the triserial planktonic foraminiferal family Guembelitriidae. The origin and evolution of this interesting albeit poorly known family are enigmatic. To elucidate the phylogenetic relationships between G. vivans and other planktonic foraminifera, we sequenced the small subunit ribosomal DNA (SSU rDNA) for comparison to our extensive database of planktonic and benthic species. Our analyses suggest that G. vivans represents a separate lineage of planktonic foraminifera, which branches close to the benthic rotaliids Stainforthia and Virgulinella. Both genera resemble Gallitellia in general morphological appearance, having elongate triserial tests at least in their early ontogenic stages. The divergence time of G. vivans is estimated at ca. 18 Ma (early Miocene), suggesting an origin independent from the Cretaceous and Paleogene triserial planktonic foraminifera. Our study thus indicates that modern triserial planktonic foraminifera are not related to the Cretaceous–Paleogene triserial species, and that the sporadic occurrences in the fossil record are not the result of poor preservation, but reflect multiple transitions from benthic to planktonic mode of life.     

Molecular evidence for an independent origin of modern triserial planktonic foraminifera from benthic ancestors

Gallitellia vivans is the only Recent representative of the triserial planktonic foraminiferal family Guembelitriidae. The origin and evolution of this interesting albeit poorly known family are enigmatic. To elucidate the phylogenetic relationships between G. vivans and other planktonic foraminifera, we sequenced the small subunit ribosomal DNA (SSU rDNA) for comparison to our extensive database of planktonic and benthic species. Our analyses suggest that G. vivans represents a separate lineage of planktonic foraminifera, which branches close to the benthic rotaliids Stainforthia and Virgulinella. Both genera resemble Gallitellia in general morphological appearance, having elongate triserial tests at least in their early ontogenic stages. The divergence time of G. vivans is estimated at ca. 18 Ma (early Miocene), suggesting an origin independent from the Cretaceous and Paleogene triserial planktonic foraminifera. Our study thus indicates that modern triserial planktonic foraminifera are not related to the Cretaceous–Paleogene triserial species, and that the sporadic occurrences in the fossil record are not the result of poor preservation, but reflect multiple transitions from benthic to planktonic mode of life.     

Phylogenetic analysis of the Metastrongyloidea (Nematoda: Strongylida) inferred from ribosomal RNA gene sequences

Phylogenetic relationships among nematodes of the strongylid superfamily Metastrongyloidea were analyzed using partial sequences from the large-subunit ribosomal RNA (LSU rRNA) and small-subunit ribosomal RNA (SSU rRNA) genes. Regions of nuclear ribosomal DNA (rDNA) were amplified by polymerase chain reaction, directly sequenced, aligned, and phylogenies inferred using maximum parsimony. Phylogenetic hypotheses inferred from the SSU rRNA gene supported the monophyly of representative taxa from each of the 7 currently accepted metastrongyloid families. Metastrongyloid taxa formed the sister group to representative trichostrongyloid sequences based on SSU data. Sequences from either the SSU or LSU RNA regions alone provided poor resolution for relationships within the Metastrongyloidea. However, a combined analysis using sequences from all rDNA regions yielded 3 equally parsimonious trees that represented the abursate Filaroididae as polyphyletic, Parafilaroides decorus as the sister species to the monophyletic Pseudaliidae, and a sister group relationship between Oslerus osleri and Metastrongylus salmi. Relationships among 3 members of the Crenosomatidae, and 1 representative of the Skrjabingylidae (Skrjabingylus chitwoodorum) were not resolved by these combined data. However, members of both these groups were consistently resolved as the sister group to the other metastrongyloid families. These relationships are inconsistent with traditional classifications of the Metastrongyloidea and existing hypotheses for their evolution.     

Early origin of foraminifera suggested by SSU rRNA gene sequences

Foraminifera are one of the largest groups of unicellular eukaryotes with probably the best known fossil record. However, the origin of foraminifera and their phylogenetic relationships with other eukaryotes are not well established. In particular, two recent reports, based on ribosomal RNA gene sequences, have reached strikingly different conclusions about foraminifera's evolutionary position within eukaryotes. Here, we present the complete small subunit (SSU) rRNA gene sequences of three species of foraminifera. Phylogenetic analysis of these sequences indicates that they branch very deeply in the eukaryotic evolutionary tree: later than those of the amitochondrial Archezoa, but earlier than those of the Euglenozoa and other mitochondria-bearing phyla. Foraminifera are clearly among the earliest eukaryotes with mitochondria, but because of the peculiar nature of their SSU genes we cannot be certain that they diverged first, as our data suggest.      

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

Higher-level phylogeny of Foraminifera inferred from the RNA polymerase II (RPB1) gene

Macroevolutionary relations among main lineages of Foraminifera have traditionally been inferred from the small subunit ribosomal genes (SSU rDNA). However, important discrepancies in the rates of SSU rDNA evolution between major lineages led to difficulties in accurate interpretation of SSU-based phylogenetic reconstructions. Recently, actin and beta-tubulin sequences have been used as alternative markers of foraminiferal phylogeny and their analyses globally confirm results obtained with SSU rDNA. In order to test new protein markers, we sequenced a fragment of the largest subunit of the RNA polymerase II (RPB1), a nuclear encoded single copy gene, for 8 foraminiferal species representing major orders of Foraminifera. Analyses of our data robustly confirm previous SSU rDNA and actin phylogenies and show (i) the paraphyly and ancestral position of monothalamid Foraminifera; (ii) the independent origin of miliolids; (iii) the monophyly of rotaliids, including buliminids and globigerinids; and (iv) the polyphyly of planktonic families Globigerinidae and Candeinidae. Additionally, the RPB1 phylogeny suggests Allogromiidae as the most ancestral foraminiferal lineage. In the light of our study, RPB1 appears as a valuable phylogenetic marker, particularly useful for groups of protists showing extreme variations of evolutionary rates in ribosomal genes.     

INTRAGENOMIC NUCLEOTIDE POLYMORPHISM AMONG SMALL SUBUNIT (18S) RDNA PARALOGS IN THE DIATOM GENUS SKELETONEMA (BACILLARIOPHYTA)1

Morphological features of the siliceous cell wall traditionally have been used to diagnose and classify species of diatoms, though an increasing number of studies distinguish new species, in part, by phylogenetic analysis of rDNA sequences. Intragenomic sequence variation is common among the hundreds to thousands of rDNA cistrons present within a genome, and this variation has strong potential to obscure species boundaries based on rDNA sequences. We screened six Skeletonema culture strains for intragenomic nucleotide polymorphisms in the small subunit (SSU) rDNA gene and found that all strains had polymorphic sites, with proportions ranging from 0.57% to 1.81%. In all cases, transitions accounted for more than 70% of nucleotide differences at polymorphic sites. Polymorphic sites were split nearly evenly in the SSU rRNA molecule between the base‐paired regions of helices (52%) and the unpaired regions of loops and bulges (48%). Phylogenetic analysis showed that SSU rDNA genotypes were monophyletic for two of the six culture strains examined. Genotypes from the other four culture strains either showed little or no phylogenetic structure compared with genotypes of other conspecific culture strains or had phylogenetic structure that was incongruent with existing species boundaries. Moderate to strong support for monophyly was recovered for four of the seven species included in the analysis. Phylogenetic results combined with the low sequence divergence of SSU rDNA genotypes within species suggest that concerted evolution has not proceeded to completion in these species and/or that the rate at which variation is being generated exceeds the rate at which concerted evolution is expunging variation.     

Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences

The Mollusca represent one of the most morphologically diverse animal phyla, prompting a variety of hypotheses on relationships between the major lineages within the phylum based upon morphological, developmental, and paleontological data. Analyses of small-ribosomal RNA (SSU rRNA) gene sequence have provided limited resolution of higher-level relationships within the Mollusca. Recent analyses suggest large-subunit (LSU) rRNA gene sequences are useful in resolving deep-level metazoan relationships, particularly when combined with SSU sequence. To this end, LSU (approximately 3.5 kb in length) and SSU (approximately 2 kb) sequences were collected for 33 taxa representing the major lineages within the Mollusca to improve resolution of intraphyletic relationships. Although the LSU and combined LSU+SSU datasets appear to hold potential for resolving branching order within the recognized molluscan classes, low bootstrap support was found for relationships between the major lineages within the Mollusca. LSU+SSU sequences also showed significant levels of rate heterogeneity between molluscan lineages. The Polyplacophora, Gastropoda, and Cephalopoda were each recovered as monophyletic clades with the LSU+SSU dataset. While the Bivalvia were not recovered as monophyletic clade in analyses of the SSU, LSU, or LSU+SSU, the Shimodaira-Hasegawa test showed that likelihood scores for these results did not differ significantly from topologies where the Bivalvia were monophyletic. Analyses of LSU sequences strongly contradict the widely accepted Diasoma hypotheses that bivalves and scaphopods are closely related to one another. The data are consistent with recent morphological and SSU analyses suggesting scaphopods are more closely related to gastropods and cephalopods than to bivalves. The dataset also presents the first published DNA sequences from a neomeniomorph aplacophoran, a group considered critical to our understanding of the origin and early radiation of the Mollusca.     

Phylogenetic analyses indicate that the 19'Hexanoyloxy-fucoxanthin-containing dinoflagellates have tertiary plastids of haptophyte origin

The three anomalously pigmented dinoflagellates Gymnodinium galatheanum, Gyrodinium aureolum, and Gymnodinium breve have plastids possessing 19'-hexanoyloxy-fucoxanthin as the major carotenoid rather than peridinin, which is characteristic of the majority of the dinoflagellates. Analyses of SSU rDNA from the plastid and the nuclear genome of these dinoflagellate species indicate that they have acquired their plastids via endosymbiosis of a haptophyte. The dinoflagellate plastid sequences appear to have undergone rapid sequence evolution, and there is considerable divergence between the three species. However, distance, parsimony, and maximum-likelihood phylogenetic analyses of plastid SSU rRNA gene sequences place the three species within the haptophyte clade. Pavlova gyrans is the most basal branching haptophyte and is the outgroup to a clade comprising the dinoflagellate sequences and those of other haptophytes. The haptophytes themselves are thought to have plastids of a secondary origin; hence, these dinoflagellates appear to have tertiary plastids. Both molecular and morphological data divide the plastids into two groups, where G. aureolum and G. breve have similar plastid morphology and G. galatheanum has plastids with distinctive features.     

PFR2: a curated database of planktonic foraminifera 18S ribosomal DNA as a resource for studies of plankton ecology,biogeography and evolution

Planktonic foraminifera (Rhizaria) are ubiquitous marine pelagic protists producing calcareous shells with conspicuous morphology. They play an important role in the marine carbon cycle, and their exceptional fossil record serves as the basis for biochronostratigraphy and past climate reconstructions. A major worldwide sampling effort over the last two decades has resulted in the establishment of multiple large collections of cryopreserved individual planktonic foraminifera samples. Thousands of 18S rDNA partial sequences have been generated, representing all major known morphological taxa across their worldwide oceanic range. This comprehensive data coverage provides an opportunity to assess patterns of molecular ecology and evolution in a holistic way for an entire group of planktonic protists. We combined all available published and unpublished genetic data to build PFR², the Planktonic foraminifera Ribosomal Reference database. The first version of the database includes 3322 reference 18S rDNA sequences belonging to 32 of the 47 known morphospecies of extant planktonic foraminifera, collected from 460 oceanic stations. All sequences have been rigorously taxonomically curated using a six‐rank annotation system fully resolved to the morphological species level and linked to a series of metadata. The PFR² website, available at http://pfr2.sb-roscoff.fr , allows downloading the entire database or specific sections, as well as the identification of new planktonic foraminiferal sequences. Its novel, fully documented curation process integrates advances in morphological and molecular taxonomy. It allows for an increase in its taxonomic resolution and assures that integrity is maintained by including a complete contingency tracking of annotations and assuring that the annotations remain internally consistent.     

Evolution of a Planktonic Foraminifer during Environmental Changes in the Tropical Oceans

Ecological adaptation to environmental changes is a strong driver of evolution, enabling speciation of pelagic plankton in the open ocean without the presence of effective physical barriers to gene flow. The tropical ocean environment, which plays an important role in shaping marine biodiversity, has drastically and frequently changed since the Pliocene. Nevertheless, the evolutionary history of tropical pelagic plankton has been poorly understood, as phylogeographic investigations are still in the developing state and paleontological approaches are insufficient to obtain a sequential record from the deep-sea sediments. The planktonic foraminifer Pulleniatina obliquiloculata is widely distributed in the tropical area throughout the world’s oceans, and its phylogeography is well established. It is thus one of the best candidates to examine how past environmental changes may have shifted the spatial distribution and affected the diversification of tropical pelagic plankton. Such an examination requires the divergence history of the planktonic foraminifers, yet the gene marker (partial small subunit (SSU) rDNA) previously used for phylogeographic studies was not powerful enough to achieve a high accuracy in estimating the divergence times. The present study focuses on improving the precision of divergence time estimates for the splits between sibling species (genetic types) of planktonic foraminifers by increasing the number of genes as well as the number of nucleotide bases used for molecular clock estimates. We have amplified the entire coding regions of two ribosomal RNA genes (SSU rDNA and large subunit (LSU) rDNA) of three genetic types of P. obliquiloculata and two closely related species for the first time and applied them to the Bayesian relaxed clock method. The comparison of the credible intervals of the four datasets consisting either of sequences of the partial SSU rDNA, the complete SSU rDNA, LSU rDNA, or a combination of both genes (SSU+LSU) clearly demonstrated that the two-gene dataset improved the accuracy of divergence time estimates. The P. obliquiloculata lineage diverged twice, first at the end of the Pliocene (3.1 Ma) and again in the middle Pleistocene (1.4 Ma). Both timings coincided with the environmental changes, which indirectly involved geographic separation of populations. The habitat of P. obliquiloculata was expanded toward the higher latitudinal zones during the stable warm periods and subsequently placed on the steep environmental gradients following the global cooling. Different environmental conditions in the stable warm tropics and unstable higher latitudes may have triggered ecological divergence among the populations, leading to adaptive differentiation and eventually speciation. A comprehensive analysis of divergence time estimates combined with phylogeography enabled us to reveal the evolutionary history of the pelagic plankton and to find the potential paleoenvironmental events, which could have changed their biogeography and ecology.     

Molecular versus Taxonomic Rates of Evolution in Planktonic Foraminifera

Neogene planktonic foraminifera are among the most widely used microfossils in the study of tempo and mode of evolution. Comparisons of taxonomic rates between the two major clades in this group have shown that the nonspinose globorotaliids have undergone a significantly more rapid evolutionary turnover than the spinose globigerinids (S. M. Stanleyet al., 1988,Paleobiology14, 235–249). In order to test if similar fluctuations are observed in molecular data, we have used different methods to calculate absolute and relative rates of substitutions based on 16 partial SSU rDNA sequences from representatives of both groups. According to our data, rates of substitution are relatively constant within the globigerinids with a mean value of 4.3 subst./site/10⁹years, but vary in the globorotaliid clade with three species having a rate of about 1 subst./site/10⁹years and two species evolving much faster with rates of more than 7 subst./site/10⁹years. Assuming that the fast rates result from recent accelerations, the globorotaliids have basically much slower molecular evolutionary rates than the globigerinids, in opposition to the fossil data.     

Planktic foraminiferal molecular evolution and their polyphyletic origins from benthic taxa

Phylogenetic analyses based on partial sequences of the small subunit (SSU) ribosomal (r) RNA gene have shown that the planktic and benthic foraminifera form a distinct monophyletic group within the eukaryotes. In order to determine the evolutionary relationships between benthic and planktic foraminifers, representatives of spinose and non-spinose planktic genera have been placed within a molecular SSU rDNA phylogeny containing sequences of the benthic suborders available to date. Our phylogenetic analysis shows that the planktic foraminifers are polyphyletic in origin, not evolving solely from a single ‘globigerinid-like’ lineage in the Mid-Jurassic, but derived from at least two ancestral benthic lines. The benthic ancestor of Neogloboquadrina dutertrei may have entered the plankton later than the Mid-Jurassic, and further investigation of related extant species should provide an indication of the timing of this event. The evolutionary origin of the non-spinose species Globorotalia menardii remains unclear. The divergences of the planktic spinose species generally support recent phylogenies based on the fossil record, which infer a radiation from a globigerinid common ancestor in the Mid- to Late Oligocene. The branching pattern indicates that there are possibly four distinct groups within the main spinose clade, with large evolutionary distances being observed between them. Globigerinoides conglobatus clusters strongly with Globigerinoides ruber and are divergent from Globigerinella siphonifera, Orbulina universa and Globigerinoides sacculifer.Conserved regions of the SSU rRNA gene show sufficient variation to discriminate foraminifers at the species level. Large genetic differences have been observed between the pink and white forms of Gs. ruber and between Ge. siphonifera Type I and II. The two types of Ge. siphonifera cannot be discriminated by traditional palaeontological methods, which has considerable implications for tracing fossil lineages and for the estimation of molecular evolutionary rates based upon the fossil record. The conserved regions show a high degree of sequence identity within a species, providing signature sequences for species identification. The variable regions of the gene may prove informative for population level studies in some species although complete sequence identity was observed in G. sacculifer and O. universa between specimens collected from the Caribbean and Western Pacific.     

Recovery and phylogenetic analysis of novel archaeal rRNA sequences from a deep-sea deposit feeder.

In 1992, two independent reports based on small-subunit rRNA gene (SSU rDNA) cloning revealed the presence of novel Archaea among marine bacterioplankton. Here, we report the presence of further novel Archaea SSU rDNA sequences recovered from the midgut contents of a deep-sea marine holothurian. Phylogenetic analyses show that these abyssal Archaea are a paraphyletic component of a highly divergent clade that also includes some planktonic sequences. Our data confirm that this clade is a deep-branching lineage in the tree of life.     

Evolution of the Diatoms (Bacillariophyta): IV. A Reconstruction of Their Age from Small Subunit rRNA Coding Regions and the Fossil Record

Small subunit ribosomal RNA (ssu rRNA) coding regions from 30 diatoms, 3 oomycetes, and 6 pelagophytes were used to construct linearized trees, maximum-likelihood trees, and neighbor-joining trees inferred from both unweighted and weighted distances. Stochastic accumulation of sequence substitutions among the diatoms was assessed with relative rate tests. Pennate diatoms evolved relatively slowly but within the limits set by a stochastic model; centric diatoms exceeded those limits. A rate distribution test was devised to identify those taxa showing an aberrant distribution of base substitutions within the ssu rRNA coding region. First appearance dates of diatom taxa from the fossil record were regressed against their corresponding branch lengths to infer the average and earliest possible age for the origin of the diatoms, the pennate diatoms, and the centric diatom order Thalassiosirales. Our most lenient age estimate (based on the median-evolving diatom taxon in the maximum-likelihood tree or on the average branch length in a linearized tree) suggests that their average age is approximately 164–166 Ma, which is close to their earliest fossil record. Both calculations suggest that it is unlikely that diatoms existed prior to 238–266 Ma. Rate variation among the diatoms' ssu rRNA coding regions and uncertainties associated with the origin of extant taxa in the fossil record contribute significantly to the variation in age estimates obtained. Different evolutionary models and the exclusion of fast or slow evolving taxa did not significantly affect age estimates; however, the inclusion of aberrantly fast evolving taxa did. Our molecular clock calibrations indicate that the rRNA coding regions in the diatoms are evolving at approximately 1% per 18 to 26 Ma, which is the fastest substitution rate reported in any pro- or eukaryotic group of organisms to date.     

Molecular Phylogeny and Morphology of Leannia veloxifera n. gen. et sp. Unveils a New Lineage of Monothalamous Foraminifera

Monothalamous (single‐chambered) foraminifera have long been considered as the “poor cousins” of multichambered species, which calcareous and agglutinated tests dominate in the fossil record. This view is currently changing with environmental DNA surveys showing that the monothalamids may be as diverse as hard‐shelled foraminifera. Yet, the majority of numerous molecular lineages revealed by eDNA studies remain anonymous. Here, we describe a new monothalamous species and genus isolated from the sample of sea grass collected in Gulf of Eilat (Red Sea). This new species, named Leannia veloxifera, is characterized by a tiny ovoid theca (about 50–100 μm) composed of thin organic wall, with two opposite apertures. The examined individuals are multinucleated and show very active reticulopodial movement. Phylogenetic analyses of SSU rDNA, actin, and beta‐tubulin (ß‐tubulin) show that the species represents a novel lineage branching separately from other monothalamous foraminifera. Interestingly, the SSU rDNA sequence of the new species is very similar to an environmental foraminiferal sequence from Bahamas, suggesting that the novel lineage may represent a group of shallow‐water tropical allogromiids, poorly studied until now.     

Can Abundance of Protists Be Inferred from Sequence Data: A Case Study of Foraminifera

Protists are key players in microbial communities, yet our understanding of their role in ecosystem functioning is seriously impeded by difficulties in identification of protistan species and their quantification. Current microscopy-based methods used for determining the abundance of protists are tedious and often show a low taxonomic resolution. Recent development of next-generation sequencing technologies offered a very powerful tool for studying the richness of protistan communities. Still, the relationship between abundance of species and number of sequences remains subjected to various technical and biological biases. Here, we test the impact of some of these biological biases on sequence abundance of SSU rRNA gene in foraminifera. First, we quantified the rDNA copy number and rRNA expression level of three species of foraminifera by qPCR. Then, we prepared five mock communities with these species, two in equal proportions and three with one species ten times more abundant. The libraries of rDNA and cDNA of the mock communities were constructed, Sanger sequenced and the sequence abundance was calculated. The initial species proportions were compared to the raw sequence proportions as well as to the sequence abundance normalized by rDNA copy number and rRNA expression level per species. Our results showed that without normalization, all sequence data differed significantly from the initial proportions. After normalization, the congruence between the number of sequences and number of specimens was much better. We conclude that without normalization, species abundance determination based on sequence data was not possible because of the effect of biological biases. Nevertheless, by taking into account the variation of rDNA copy number and rRNA expression level we were able to infer species abundance, suggesting that our approach can be successful in controlled conditions.