Unveiling the Biodiversity of Deep-Sea Nematodes through Metabarcoding: Are We Ready to Bypass the Classical Taxonomy?

Nematodes inhabiting benthic deep-sea ecosystems account for >90% of the total metazoan abundances and they have been hypothesised to be hyper-diverse, but their biodiversity is still largely unknown. Metabarcoding could facilitate the census of biodiversity, especially for those tiny metazoans for which morphological identification is difficult. We compared, for the first time, different DNA extraction procedures based on the use of two commercial kits and a previously published laboratory protocol and tested their suitability for sequencing analyses of 18S rDNA of marine nematodes. We also investigated the reliability of Roche 454 sequencing analyses for assessing the biodiversity of deep-sea nematode assemblages previously morphologically identified. Finally, intra-genomic variation in 18S rRNA gene repeats was investigated by Illumina MiSeq in different deep-sea nematode morphospecies to assess the influence of polymorphisms on nematode biodiversity estimates. Our results indicate that the two commercial kits should be preferred for the molecular analysis of biodiversity of deep-sea nematodes since they consistently provide amplifiable DNA suitable for sequencing. We report that the morphological identification of deep-sea nematodes matches the results obtained by metabarcoding analysis only at the order-family level and that a large portion of Operational Clustered Taxonomic Units (OCTUs) was not assigned. We also show that independently from the cut-off criteria and bioinformatic pipelines used, the number of OCTUs largely exceeds the number of individuals and that 18S rRNA gene of different morpho-species of nematodes displayed intra-genomic polymorphisms. Our results indicate that metabarcoding is an important tool to explore the diversity of deep-sea nematodes, but still fails in identifying most of the species due to limited number of sequences deposited in the public databases, and in providing quantitative data on the species encountered. These aspects should be carefully taken into account before using metabarcoding in quantitative ecological research and monitoring programmes of marine biodiversity.     

Evaluating high-throughput sequencing as a method for metagenomic analysis of nematode diversity

Nematodes play an important role in ecosystem processes, yet the relevance of nematode species diversity to ecology is unknown. Because nematode identification of all individuals at the species level using standard techniques is difficult and time-consuming, nematode communities are not resolved down to the species level, leaving ecological analysis ambiguous. We assessed the suitability of massively parallel sequencing for analysis of nematode diversity from metagenomic samples. We set up four artificial metagenomic samples involving 41 diverse reference nematodes in known abundances. Two samples came from pooling polymerase chain reaction products amplified from single nematode species. Two additional metagenomic samples consisted of amplified products of DNA extracted from pooled nematode species. Amplified products involved two rapidly evolving ~400-bp sections coding for the small and large subunit of rRNA. The total number of reads ranged from 4159 to 14771 per metagenomic sample. Of these, 82% were > 199 bp in length. Among the reads > 199 bp, 86% matched the referenced species with less than three nucleotide differences from a reference sequence. Although neither rDNA section recovered all nematode species, the use of both loci improved the detection level of nematode species from 90 to 97%. Overall, results support the suitability of massively parallel sequencing for identification of nematodes. In contrast, the frequency of reads representing individual species did not correlate with the number of individuals in the metagenomic samples, suggesting that further methodological work is necessary before it will be justified for inferring the relative abundances of species within a nematode community.     

Reproducibility of read numbers in high-throughput sequencing analysis of nematode community composition and structure

Although nematodes are the most abundant metazoan animals on Earth, their diversity is largely unknown. To overcome limitations of traditional approaches (labour, time, and cost) for assessing biodiversity of nematode species in environmental samples, we have previously examined the suitability of high-throughput sequencing for assessing species level diversity with a set of control experiments employing a mixture of nematodes of known number and with known sequences for target diagnostic loci. Those initial experiments clearly demonstrated the suitability of the approach for identification of nematode taxa but lacked the replicate experiments necessary to evaluate reproducibility of the approach. Here, we analyze reads generated from three different PCR amplifications and three different sequencing reactions to examine the differential PCR amplification, the possibility of emulsion PCR artefacts, and differences between sequencing reactions. Our results suggest that both qualitative and quantitative data are consistent and highly reproducible. Variation associated with in-house PCR amplification or emPCR and sequencing are present but the representation of each nematode is very consistent from experiment to experiment and supports the use of read counts to estimate relative abundance of taxa in a metagenetic sample.     

Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies

The advent of next generation sequencing has coincided with a growth in interest in using these approaches to better understand the role of the structure and function of the microbial communities in human, animal, and environmental health. Yet, use of next generation sequencing to perform 16S rRNA gene sequence surveys has resulted in considerable controversy surrounding the effects of sequencing errors on downstream analyses. We analyzed 2.7×10(6) reads distributed among 90 identical mock community samples, which were collections of genomic DNA from 21 different species with known 16S rRNA gene sequences; we observed an average error rate of 0.0060. To improve this error rate, we evaluated numerous methods of identifying bad sequence reads, identifying regions within reads of poor quality, and correcting base calls and were able to reduce the overall error rate to 0.0002. Implementation of the PyroNoise algorithm provided the best combination of error rate, sequence length, and number of sequences. Perhaps more problematic than sequencing errors was the presence of chimeras generated during PCR. Because we knew the true sequences within the mock community and the chimeras they could form, we identified 8% of the raw sequence reads as chimeric. After quality filtering the raw sequences and using the Uchime chimera detection program, the overall chimera rate decreased to 1%. The chimeras that could not be detected were largely responsible for the identification of spurious operational taxonomic units (OTUs) and genus-level phylotypes. The number of spurious OTUs and phylotypes increased with sequencing effort indicating that comparison of communities should be made using an equal number of sequences. Finally, we applied our improved quality-filtering pipeline to several benchmarking studies and observed that even with our stringent data curation pipeline, biases in the data generation pipeline and batch effects were observed that could potentially confound the interpretation of microbial community data.     

Phylogenetic analysis of 70 kD heat shock protein sequences suggests a chimeric origin for the eukaryotic cell nucleus

Background: The evolutionary relationships between archaebacteria, eubacteria and eukaryotic cells are of central importance in biology. The current view is that each of these three groups of organisms constitutes a monophyletic domain, and that eukaryotic cells have evolved from an archaebacterial ancestor. Recent studies on a number of highly conserved protein sequences do not, however, support this view and raise important questions concerning the evolutionary relationships between all extant organisms, particularly regarding the origin of eukaryotic cells.Results We have used sequences of 70 kD heat shock protein (hsp70) — the most conserved protein found to date in all species — to examine the evolutionary relationship between various species. We have obtained two new archaebacterial hsp70 sequences from the species, Thermoplasma acidophilum and Halobacterium cutirubrum. A global comparison of hsp70 sequences, including our two new sequences, shows that all known archaebacterial homologs share a number of sequence signatures with the Gram-positive group of bacteria that are not found in any other prokaryotic or eukaryotic species. In contrast, the eukaryotic homologs are shown to share a number of unique sequence features with the Gram-negative bacteria that are not present in any archaebacteria. Detailed phylogenetic analyses of hsp70 sequences strongly support a specific evolutionary relationship between archaebacteria and Gram-positive bacteria on the one hand, and Gram-negative bacteria and eukaryotes on the other. The phylogenetic analyses also indicate a polyphyletic branching of archaebacteria within the Gram-positive species. The possibility that the observed relationships are due to horizontal gene transfers can be excluded on the basis of sequence characteristics of different groups of homologs.Conclusion Our results do not support the view that archaebacteria constitute a monophyletic domain, but instead suggest a close evolutionary linkage between archaebacteria and Gram-positive bacteria. Furthermore, in contrast to the presently accepted view, eukaryotic hsp70s show a close and specific relationship to those from Gram-negative species. To explain the phylogenies based on different gene sequences, a chimeric model for the origin of the eukaryotic cell nucleus involving fusion between an archaebacterium and a Gram-negative eubacterium is proposed. Several predictions from the chimeric model are discussed.     

Quantitative and qualitative assessment of pollen DNA metabarcoding using constructed species mixtures

Pollen DNA metabarcoding—marker‐based genetic identification of potentially mixed‐species pollen samples—has applications across a variety of fields. While basic species‐level pollen identification using standard DNA barcode markers is established, the extent to which metabarcoding (a) correctly assigns species identities to mixes (qualitative matching) and (b) generates sequence reads proportionally to their relative abundance in a sample (quantitative matching) is unclear, as these have not been assessed relative to known standards. We tested the quantitative and qualitative robustness of metabarcoding in constructed pollen mixtures varying in species richness (1–9 species), taxonomic relatedness (within genera to across class) and rarity (5%–100% of grains), using Illumina MiSeq with the markers rbcL and ITS2. Qualitatively, species composition determinations were largely correct, but false positives and negatives occurred. False negatives were typically driven by lack of a barcode gap or rarity in a sample. Species richness and taxonomic relatedness, however, did not strongly impact correct determinations. False positives were likely driven by contamination, chimeric sequences and/or misidentification by the bioinformatics pipeline. Quantitatively, the proportion of reads for each species was only weakly correlated with its relative abundance, in contrast to suggestions from some other studies. Quantitative mismatches are not correctable by consistent scaling factors, but instead are context‐dependent on the other species present in a sample. Together, our results show that metabarcoding is largely robust for determining pollen presence/absence but that sequence reads should not be used to infer relative abundance of pollen grains.     

Sample richness and genetic diversity as drivers of chimera formation in nSSU metagenetic analyses

Eukaryotic diversity in environmental samples is often assessed via PCR-based amplification of nSSU genes. However, estimates of diversity derived from pyrosequencing environmental data sets are often inflated, mainly because of the formation of chimeric sequences during PCR amplification. Chimeras are hybrid products composed of distinct parental sequences that can lead to the misinterpretation of diversity estimates. We have analyzed the effect of sample richness, evenness and phylogenetic diversity on the formation of chimeras using a nSSU data set derived from 454 Roche pyrosequencing of replicated, large control pools of closely and distantly related nematode mock communities, of known intragenomic identity and richness. To further investigate how chimeric molecules are formed, the nSSU gene secondary structure was analyzed in several individuals. For the first time in eukaryotes, chimera formation proved to be higher in both richer and more genetically diverse samples, thus providing a novel perspective of chimera formation in pyrosequenced environmental data sets. Findings contribute to a better understanding of the nature and mechanisms involved in chimera formation during PCR amplification of environmentally derived DNA. Moreover, given the similarities between biodiversity analyses using amplicon sequencing and those used to assess genomic variation, our findings have potential broad application for identifying genetic variation in homologous loci or multigene families in general.     

Influence of bacterivorous nematodes on the decomposition of cordgrass

The influence of bacterivorous nematodes (Diplolaimelloides meyli, Diplolaimelloides oschei, Diplolaimella dievengatensis, Panagrolaimus paetzoldi) on the decomposition of a macrophyte (Spartina anglica) in an aquatic environment was investigated by using laboratory microcosm experiments. Several earlier studies have shown enhancement of the decomposition process in the presence of nematodes, but nematode species-specific effects were never tested. In this study four bacterivorous nematode species were applied separately to microcosms to investigate such species-specific influences.No stimulation of the decomposition process nor of the microbial community was observed in the presence of the nematodes, both were highest in the absence of nematodes. However, clear differences were found between nematode treatments. P. paetzoldi reached much higher numbers than the other species, causing a decrease in microbial activity, probably due to (over)grazing. Remarkably this low microbial activity did not result in a slow-down of the decomposition process compared to the other nematode treatments, raising the question whether P. paetzoldi might be able to directly assimilate detrital compounds. Other nematode species reached much lower densities, but nevertheless an influence on the decomposition process was observed. However, this experiment does not support the view that bacterivorous nematodes enhance decomposition rate.The experimental results show that in nematode communities the use of functional groups is inadequate for biodiversity studies. The four nematode species used in this study belong to the same functional group, but are clearly not functionally redundant since they all have a different influence on the cordgrass decomposition. This suggests that the relationship between nematode species diversity and ecosystem functioning may be idiosyncratic.     

Systematic Characteristic Exploration of the Chimeras Generated in Multiple Displacement Amplification through Next Generation Sequencing Data Reanalysis

Background

The chimeric sequences produced by phi29 DNA polymerase, which are named as chimeras, influence the performance of the multiple displacement amplification (MDA) and also increase the difficulty of sequence data process. Despite several articles have reported the existence of chimeric sequence, there was only one research focusing on the structure and generation mechanism of chimeras, and it was merely based on hundreds of chimeras found in the sequence data of E. coli genome.

Method

We finished data mining towards a series of Next Generation Sequencing (NGS) reads which were used for whole genome haplotype assembling in a primary study. We established a bioinformatics pipeline based on subsection alignment strategy to discover all the chimeras inside and achieve their structural visualization. Then, we artificially defined two statistical indexes (the chimeric distance and the overlap length), and their regular abundance distribution helped illustrate of the structural characteristics of the chimeras. Finally we analyzed the relationship between the chimera type and the average insertion size, so that illustrate a method to decrease the proportion of wasted data in the procedure of DNA library construction.

Results/Conclusion

131.4 Gb pair-end (PE) sequence data was reanalyzed for the chimeras. Totally, 40,259,438 read pairs (6.19%) with chimerism were discovered among 650,430,811 read pairs. The chimeric sequences are consisted of two or more parts which locate inconsecutively but adjacently on the chromosome. The chimeric distance between the locations of adjacent parts on the chromosome followed an approximate bimodal distribution ranging from 0 to over 5,000 nt, whose peak was at about 250 to 300 nt. The overlap length of adjacent parts followed an approximate Poisson distribution and revealed a peak at 6 nt. Moreover, unmapped chimeras, which were classified as the wasted data, could be reduced by properly increasing the length of the insertion segment size through a linear correlation analysis.

Significance

This study exhibited the profile of the phi29MDA chimeras by tens of millions of chimeric sequences, and helped understand the amplification mechanism of the phi29 DNA polymerase. Our work also illustrated the importance of NGS data reanalysis, not only for the improvement of data utilization efficiency, but also for more potential genomic information.     

Inventory of the benthic eukaryotic diversity in the oldest European lake

We have profound knowledge on biodiversity on Earth including plants and animals. In the recent decade, we have also increased our understanding on microorganisms in different hosts and the environment. However, biodiversity is not equally well studied among different biodiversity groups and Earth''s systems with eukaryotes in freshwater sediments being among the least known. In this study, we used high‐throughput sequencing of the 18S rRNA gene to investigate the entire diversity of benthic eukaryotes in three distinct habitats (littoral sediment and hard substrate, profundal sediment) of Lake Ohrid, the oldest European lake. Eukaryotic sequences were dominated by annelid and arthropod animals (54% of all eukaryotic reads) and protists (Ochrophyta and Ciliophora; together 40% of all reads). Eukaryotic diversity was 15% higher in the deep profundal than on either near‐surface hard substrates or littoral sediments. The three habitats differed in their taxonomic and functional community composition. Specifically, heterotrophic organisms accounted for 92% of the reads in the profundal, whereas phototrophs accounted for 43% on the littoral hard substrate. The profundal community was the most homogeneous, and its network was the most complex, suggesting its highest stability among the sampled habitats.     

Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

Deep-sea nematode biodiversity in the Mediterranean basin: testing for longitudinal, bathymetric and energetic gradients

The knowledge of the processes controlling the spatial distribution of species diversity is one of the main challenges of the present ecological research. Spatial patterns of benthic biodiversity in the deep sea are poorly known in comparison with other ecosystems and this limits our understanding of the mechanisms controlling the distribution and maintenance of high biodiversity in the largest ecosystems of our biosphere. Although the Mediterranean basin covers <1% of the world ocean surface, none the less it hosts >7.5% of the global biodiversity. The high biogeographic complexity and the presence of steep ecological gradients contribute in making the Mediterranean a region of very high diversity. Here we report the results of an investigation on the patterns of nematode biodiversity in the deep-Mediterranean Sea, in relation with bathymetric, longitudinal and energetic gradients. Our results indicate that benthic biodiversity in the deep-Mediterranean decreases significantly with increasing depth. Moreover, at equally deep sites, nematode diversity decreased from the western to the eastern basin and longitudinal gradients were evident when comparing sites at 4000-m depth, with 3000-m depth. The analysis of the available energy (measured as labile organic matter content of the sediments) suggests that biodiversity patterns are not controlled by the amounts of food resources, but instead bio-availability is the key factor. A more detailed analysis revealed an extremely high deep-sea beta-diversity (turnover diversity), both among sites at different depths as well as at similar depths of different longitude or within the same basin. This new finding has not only important implications on the estimates of the overall regional diversity (gamma diversity), but also suggests the presence of high biogeographic complexity in the deep benthic domain of the Mediterranean Sea.     

The fate of female donor blastodermal cells in male chimeric chickens.

In previous experiments in our laboratories, chickens that are chimeric in their gamete, melanocyte, and blood cell populations have been produced by injection of dispersed stage X blastodermal donor cells into the subgerminal cavity of stage X recipient embryos. In some experiments, donor cells were transfected with reporter gene constructs prior to injection as a preliminary step in the production of transgenic birds. Chimerism was assessed by test mating, observation of plumage, and DNA fingerprinting. Methods were sought that would provide a relatively rapid analysis of the spatial distribution of descendants of donor cells in chimeras to assess the efficacy of various methods of chimera construction. To date, the sex of donor and recipient embryos was not known and, therefore, numerous mixed sex chimeras must have been constructed by chance, since donor cells were usually collected from several embryos rather than from individual embryos. The presence of female-derived cells was determined by in situ hybridization using a W-chromosome-specific DNA probe, using smears of washed erythrocytes from 16 phenotypically male chimeric chickens ranging in age from 4 days to 42 months posthatching. The proportion of female cells detected in the erythrocyte samples was zero (eight samples) or very low (0.020-0.083%), although 1% of the erythrocytes from a phenotypically male chick that was killed 4 days after hatch were female-derived. The low proportions of female-derived cells were surprising, considering that most of these chimeras had been produced by the injection of cells pooled from several donor embryos and most recipients had been exposed to gamma irradiation prior to injection, thus dramatically enhancing the level of incorporation of donor cells into the resulting chimeras. By contrast, 0-100% of the erythrocytes were female-derived in blood samples taken at 10 days of incubation from the chorioallantois of seven phenotypically normal male embryos that resulted from the injection of blastodermal cells pooled from five embryos into irradiated recipient embryos. Approximately 70% of the erythrocytes in a blood sample from a phenotypically normal female chimeric embryo were female-derived, and 100% of the erythrocytes examined from an intersex embryo bearing a right testis and a left ovary were female-derived. These results indicate that female-derived cells can contribute to the formation of erythropoietic tissue during the early development of what will become a phenotypically male chimeric embryo. It would appear, therefore, that female-derived cells are blocked in development or destroyed, or certain male-female combinations of cells may be lethal prior to hatching.(ABSTRACT TRUNCATED AT 400 WORDS)     

The use of random chimeragenesis to study structure/function properties of rat and human P450c17

The microsomal 17alpha-hydroxylase/17,20-lyase cytochrome P450 (P450c17) catalyzes the 17alpha-hydroxylase reaction required to produce cortisol, the major glucocorticoid in many species and the 17,20-lyase activity required for the production of androgens in all species. Utilizing the technique of random chimeragenesis we have attempted to map regions of primary sequence that contribute to the species-specific biochemical differences between rat and human P450c17. We have previously reported significant differences between rat and human P450c17 in their activities, stability and substrate-dependent coupling efficiencies even though they share 68% amino acid identity. Identification of the regions of primary sequence that contribute to each of these properties would be helpful in understanding the structure/function relationships in this enzyme. A single plasmid containing the cDNAs encoding both enzymes in a tandem orientation was constructed. This plasmid was linearized at unique restriction sites and used to transform Escherichia coli. A three-step screening protocol identified five chimeras with a uniform distribution of 5' rat and 3' human sequence. All chimeric proteins yield the characteristic reduced-CO difference spectra, indicating proper folding. The chimeras exhibit a range of stability and activities that are not consistent with the degree of parental primary sequence. A chimera containing 301 N-terminal rat P450c17 amino acids and lacking the rat P450c17 phenylalanine 343, had the highest lyase activity. Generation of these functional rat/human chimeras suggests that the tertiary structures of rat and human P450c17 are sufficiently conserved to allow proper folding of chimeric enzymes. However, the properties of these chimeras did not permit identification of a region of primary sequence that contributes to a species-specific property of rat and human P450c17. Stability of these chimeras and insight into the presence of secondary structural elements is discussed.     

Nematode mitochondrial metagenomics: A new tool for biodiversity analysis

DNA barcoding approaches have greatly increased our understanding of biodiversity on the planet, and metabarcoding is widely used for classifying members of the phylum Nematoda. However, loci typically utilized in metabarcoding studies are often unable to resolve closely related species or are unable to recover all taxa present in a sample due to inadequate PCR primer binding. Mitochondrial metagenomics (mtMG) is an alternative approach utilizing shotgun sequencing of total DNA to recover the mitochondrial genomes of all species present in samples. However, this approach requires a comprehensive reference database for identification and currently available mitochondrial sequences for nematodes are highly dominated by sequences from the order Rhabditida, and excludes many clades entirely. Here, we analysed the efficacy of mtMG for the recovery of nematode taxa and the generation of mitochondrial genomes. We first developed a curated reference database of nematode mitochondrial sequences and expanded it with 40 newly sequenced taxa. We then tested the mito-metagenomics approach using a series of nematode mock communities consisting of morphologically identified nematode species representing various feeding traits, life stages, and phylogenetic relationships. We were able to identify all but two species through the de novo assembly of COX1 genes. We were also able to recover additional mitochondrial protein coding genes (PCGs) for 23 of the 24 detected species including a full array of 12 PCGs from five of the species. We conclude that mtMG offers a potential for the effective recovery of nematode biodiversity but remains limited by the breadth of the reference database.     

Metabarcoding reveals high diversity of benthic foraminifera linked to water masses circulation at coastal Svalbard

Arctic marine biodiversity is undergoing rapid changes due to global warming and modifications of oceanic water masses circulation. These changes have been demonstrated in the case of mega- and macrofauna, but much less is known about their impact on the biodiversity of smaller size organisms, such as foraminifera that represent a main component of meiofauna in the Arctic. Several studies analyzed the distribution and diversity of Arctic foraminifera. However, all these studies are based exclusively on the morphological identification of specimens sorted from sediment samples. Here, we present the first assessment of Arctic foraminifera diversity based on metabarcoding of sediment DNA samples collected in fjords and open sea areas in the Svalbard Archipelago. We obtained a total of 5,968,786 reads that represented 1384 amplicon sequence variants (ASVs). More than half of the ASVs (51.7%) could not be assigned to any group in the reference database suggesting a high genetic novelty of Svalbard foraminifera. The sieved and unsieved samples resolved comparable communities, sharing 1023 ASVs, comprising over 97% of reads. Our analyses show that the foraminiferal assemblage differs between the localities, with communities distinctly separated between fjord and open sea stations. Each locality was characterized by a specific assemblage, with only a small overlap in the case of open sea areas. Our study demonstrates a clear pattern of the influence of water masses on the structure of foraminiferal communities. The stations situated on the western coast of Svalbard that are strongly influenced by warm and salty Atlantic water (AW) are characterized by much higher diversity than stations in the northern and eastern part, where the impact of AW is less pronounced. This high diversity and specificity of Svalbard foraminifera associated with water mass distribution indicate that the foraminiferal metabarcoding data can be very useful for inferring present and past environmental conditions in the Arctic.     

澜沧江两条支流捕食线虫真菌水陆分布差异及关联性

为探究微生物在流域中的水、陆分布差异和相互关系,对澜沧江两条支流捕食线虫真菌多样性水陆分布进行了调查研究。在枯水期,以澜沧江一级支流沘江和黑惠江为研究区域,系统布设12个采样点(水流交汇点),每个采样点采集水、陆对称样品各5份,共采集土样120份。结合传统分离纯化、形态学及分子生物学方法筛选和鉴定菌株,按照《Nematode-Trapping Fungi》进行分类,共获得2属13种88株捕食线虫真菌;其中,陆地样品中共分离到2属11种45株,水体底泥中分离出1属10种43株,检出率分别为41.67%和53.33%。结果表明,澜沧江两条支流捕食线虫真菌在物种、属、群落结构3个层面上均存在水陆差异,也相互联系;陆地土壤可能是流域内捕食线虫真菌多样性的源,水流是其重要的传播因子。在流域生态系统内,水陆间的扩散限制和水流的连通性都是维持微生物物种多样性的重要机制。     

Barcoding marine nematodes: an improved set of nematode 18S rRNA primers to overcome eukaryotic co-interference

Nematodes form an important component of many benthic marine ecosystems and DNA barcoding approaches could provide an insight into nematode community composition from different environments globally. We have amplified nematode 18S rRNA sequences using standard nematode18S rRNA primers from environmental DNA extracted from intertidal sediment collected from New Jersey coast, USA to test whether the published marine nematode 18S rRNA sequences from GenBank and EMBL databases can effectively assign unknown nematode sequences into genus or species level. Most of the sequenced clones showed some degree of identities with published marine nematode 18S rRNA sequences. However, relatively very few of the sequences could be assigned even to genus level based on sequence assignment rule. In addition, other eukaryotic 18S rRNA sequences were found to be co-amplified with commonly used nematode 18S rRNA primers. We found that the majority of the current nematode 18S rRNA primers will co-amplify other eukaryotes if environmental DNA is the target template. We therefore designed a new set of nematode 18S rRNA primers and evaluated them using environmental DNA in intertidal sediment from the New Jersey coast. In total, 40 clones were screened and subsequently sequenced and all the sequences showed varying degree of identities with published nematode 18S rRNA sequences from GenBank and EMBL databases, and no obvious eukaryotic co-amplicons were detected with new primers. Only 13 out of 40 clones amplified with the new primer set showed 100% identity to published Daptonema and Metachromadora 18S rRNA sequences. The current molecular databases for nematodes are dominated by sequences from NW Europe and need to be more extensively populated with new full length 18S rRNA nematode sequences collected from different biogeographic locations. The new primers developed in this study, in combination with an updated nematode 18S rRNA sequence database, would help us to better investigate and understand the diversity and community composition of free-living marine nematodes based on DNA barcoding approaches during biodiversity or biomonitoring surveys on a global-scale.     

Comparing the effectiveness of metagenomics and metabarcoding for diet analysis of a leaf‐feeding monkey (Pygathrix nemaeus)

Faecal samples are of great value as a non‐invasive means to gather information on the genetics, distribution, demography, diet and parasite infestation of endangered species. Direct shotgun sequencing of faecal DNA could give information on these simultaneously, but this approach is largely untested. Here, we used two faecal samples to characterize the diet of two red‐shanked doucs langurs (Pygathrix nemaeus) that were fed known foliage, fruits, vegetables and cereals. Illumina HiSeq produced ~74 and 67 million paired reads for these samples, of which ~10 000 (0.014%) and ~44 000 (0.066%), respectively, were of chloroplast origin. Sequences were matched against a database of available chloroplast ‘barcodes’ for angiosperms. The results were compared with ‘metabarcoding’ using PCR amplification of the P6 loop of trnL. Metagenomics identified seven and nine of the likely 16 diet plants while six and five were identified by metabarcoding. Metabarcoding produced thousands of reads consistent with the known diet, but the barcodes were too short to identify several plant species to genus. Metagenomics utilized multiple, longer barcodes that combined had greater power of identification. However, rare diet items were not recovered. Read numbers for diet species in metagenomic and metabarcoding data were correlated, indicating that both are useful for determining relative sequence abundance. Metagenomic reads were uniformly distributed across the chloroplast genomes; thus, if chloroplast genomes were used as reference, the precision of identifications and species recovery would improve further. Metagenomics also recovered the host mitochondrial genome and numerous intestinal parasite sequences in addition to generating data useful for characterizing the microbiome.     

A tertiary plastid uses genes from two endosymbionts

The origin and subsequent spread of plastids by endosymbiosis had a major environmental impact and altered the course of a great proportion of eukaryotic biodiversity. The ancestor of dinoflagellates contained a secondary plastid that was acquired in an ancient endosymbiotic event, where a eukaryotic cell engulfed a red alga. This is known as secondary endosymbiosis and has happened several times in eukaryotic evolution. Certain dinoflagellates, however, are unique in having replaced this secondary plastid in an additional (tertiary) round of endosymbiosis. Most plastid proteins are encoded in the nucleus of the host and are targeted to the organelle. When secondary or tertiary endosymbiosis takes place, it is thought that these genes move from nucleus to nucleus, so the plastid retains the same proteome. We have conducted large-scale expressed sequence tag (EST) surveys from Karlodinium micrum, a dinoflagellate with a tertiary haptophyte-derived plastid, and two haptophytes, Isochrysis galbana and Pavlova lutheri. We have identified all plastid-targeted proteins, analysed the phylogenetic origin of each protein, and compared their plastid-targeting transit peptides. Many plastid-targeted genes in the Karlodinium nucleus are indeed of haptophyte origin, but some genes were also retained from the original plastid (showing the two plastids likely co-existed in the same cell), in other cases multiple isoforms of different origins exist. We analysed plastid-targeting sequences and found the transit peptides in K.micrum are different from those found in either dinoflagellates or haptophytes, pointing to a plastid with an evolutionarily chimeric proteome, and a massive remodelling of protein trafficking during plastid replacement.