Droplet digital PCR as a tool for investigating dynamics of cryptic symbionts

Using digital PCR to predict ciliate abundance from ribosomal RNA gene copy numbers

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.

Heterogeneous distribution of Cardinium in whitefly populations is associated with host nuclear genetic background

Inherited bacterial symbionts are very common in arthropods, but infection frequency can vary widely among populations. Experiments and interpopulation comparisons suggest that host genetic background might be important in explaining this variation. Our extensive field investigation showed that the infection pattern of the facultative symbiont Cardinium was heterogeneous across geographical populations of the invasive whitefly Bemisia tabaci Mediterranean (MED) in China, with genetic nuclear differences evident in 2 of the populations: 1 with a low infection rate (SD line) and 1 with a high infection rate (HaN line). However, whether the heterogeneous frequency of Cardinium is associated with the host genetic background remains poorly understood. Here, we compared the fitness of the Cardinium-infected and uninfected sublines with similar nuclear genetic backgrounds from SD and HaN lines, respectively, and further determine whether host extranuclear or nuclear genotype influenced the Cardinium-host phenotype by performing 2 new introgression series of 6 generations between SD and HaN lines (i.e., Cardinium-infected females of SD were backcrossed with uninfected males of HaN, and vice versa). The results showed that Cardinium provides marginal fitness benefits in the SD line, whereas Cardinium provides strong fitness benefits in the HaN line. Further, both Cardinium and the Cardinium-host nuclear interaction influence the fecundity and pre-adult survival rate of B. tabaci, whereas the extranuclear genotype does not. In conclusion, our results provide evidence that Cardinium-mediated fitness effects were closely associated with the host genetic background, which provides a fundamental basis for understanding the mechanism underlying the heterogeneous distribution of Cardinium in B. tabaci MED populations across China.

Coinfection patterns of two marine apicomplexans are not associated with genetic diversity of their polychaete host

Coinfections of two or more parasites within one host are more of a rule than an exception in nature. Interactions between coinfecting parasites can greatly affect their abundance and prevalence. Characteristics of the host, such as genetic diversity, can also affect the infection dynamics of coinfecting parasites. Here, we investigate for the first time the association of coinfection patterns of two marine apicomplexans, Rhytidocystis sp. and Selenidium pygospionis, with the genetic diversity of their host, the polychaete Pygospio elegans, from natural populations. Host genetic diversity was determined with seven microsatellite loci and summarized as allelic richness, inbreeding coefficient, and individual heterozygosity. We detected nonsignificant correlations between infection loads and both individual host heterozygosity and population genetic diversity. Prevalence and infection load of Rhytidocystis sp. were higher than those of S. pygospionis, and both varied spatially. Coinfections were common, and almost all hosts infected by S. pygospionis were also infected by Rhytidocystis sp. Rhytidocystis sp. infection load was significantly higher in dual infections. Our results suggest that factors other than host genetic diversity might be more important in marine apicomplexan infection patterns and experimental approaches would be needed to further determine how interactions between the apicomplexans and their host influence infection.