Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.     

肺囊康定产生菌Glarea lozoyensis线粒体基因组序列的再分析

[目的] Glarea lozoyensis是抗真菌药物卡泊芬净的产生菌,其突变菌株ATCC 74030的线粒体基因组已被报道。我们此前的研究发现诱变剂能引起该菌某些细胞核基因的突变,但诱变剂是否也能引起线粒体DNA序列的改变并不清楚。[方法] 组装野生型菌株ATCC 20868的线粒体基因组,并与发表的突变型菌株ATCC 74030的线粒体基因组进行比较。通过PCR验证野生和突变菌株线粒体基因组间表现差异之处,并利用正确的线粒体基因组序列进行新的分析。[结果] 我们成功组装出野生型菌株ATCC 20868的线粒体基因组,通过比较其与发表的ATCC 74030的线粒体基因组序列,发现存在6处单核苷酸变异位点和2处具有长度差异的区域。然而,随后的PCR验证和序列比较并没有发现2个菌株间存在这些差异。最初观察到的碱基差异是因为发表的ATCC 74030线粒体基因组存在序列错误。有趣的是,在Glarea lozoyensis的线粒体基因组中,我们发现存在3个具有内含子的tRNA基因和1个rnpB基因。同时,该菌线粒体基因组中存在多种重复序列,在其线粒体和细胞核基因组间也存在明显的DNA片段重复事件。[结论] 诱变剂没有引起G. lozoyensis线粒体DNA的任何改变;发表的ATCC 74030的线粒体基因组存在序列错误。我们报道G. lozoyensis正确的线粒体基因组序列,并且发现该菌线粒体和细胞核基因组间频繁的基因交流。     

食线虫真菌洛斯里被毛孢线粒体基因组的再分析

【目的】鉴定洛斯里被毛孢OWVT-1菌株的线粒体基因组,验证公布的USA-87-5菌株线粒体基因组中的错误,对洛斯里被毛孢正确的线粒体基因组序列进行注释并开展不同被毛孢物种间的比较线粒体基因组学分析。【方法】借助DNA高通量测序数据并通过必要的Sanger测序组装OWVT-1的线粒体基因组。通过PCR验证OWVT-1与公布的USA-87-5线粒体基因组序列差异的真实性。利用多种生物信息方法分析和注释洛斯里被毛孢的线粒体基因组。【结果】公布的洛斯里被毛孢USA-87-5菌株的线粒体基因组存在几处序列错误,包括3处长片段的插入缺失和多处短片段的插入缺失。实际上,洛斯里被毛孢USA-87-5与OWVT-1菌株的线粒体基因组序列完全相同。该菌的线粒体基因组全长62949 bp,在7个基因中共插入13个内含子,部分内含子和基因间区显现出序列退化的特征。洛斯里被毛孢、明尼苏达被毛孢、线虫被毛孢的线粒体基因组具有较强的共线性关系。除一些独立的ORF外,核心蛋白编码基因、rRNA基因和tRNA基因的排列顺序非常保守。基因间区的长短是影响3种被毛孢线粒体基因组大小最主要的因素。【结论】公布的洛斯里被毛孢USA-87-5菌株线粒体基因组中存在序列错误。本文新报道了OWVT-1菌株的线粒体基因组,并进行注释和比较线粒体基因组学分析。     

Finding stories in noise: Mitochondrial portraits from RAD data

Mitochondrial DNA (mtDNA) has formed the backbone of phylogeographic research for many years; however, recent trends focus on genome‐wide analyses. One method proposed for calibrating inferences from noisy next‐generation data, such as RAD sequencing, is to compare these results with analyses of mitochondrial sequences. Most researchers using this approach appear to be unaware that many single nucleotide polymorphisms (SNPs) identified from genome‐wide sequence data are themselves mitochondrial, or assume that these are too few to bias analyses. Here, we demonstrate two methods for mining mitochondrial markers using RAD sequence data from three South African species of yellowfish, Labeobarbus. First, we use a rigorous SNP discovery pipeline using the program stacks , to identify variant sites in mtDNA, which we then combine into haplotypes. Second, we directly map sequence reads against a mitochondrial genome reference. This method allowed us to reconstruct up to 98% of the Labeobarbus mitogenome. We validated these mitogenome reconstructions through blast database searches and by comparison with cytochrome b gene sequences obtained through Sanger sequencing. Finally, we investigate the organismal consequences of these data including ancient genetic exchange and a recent translocation among populations of L. natalensis, as well as interspecific hybridization between L. aeneus and L. kimberleyensis.     

The complete mitochondrial genome sequence of the world's largest fish,the whale shark (Rhincodon typus), and its comparison with those of related shark species

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a “vulnerable” species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species.     

Mitogenome rearrangement in the cold-water scleractinian coral Lophelia pertusa (Cnidaria, Anthozoa) involves a long-term evolving group I intron

Group I introns are genetic insertion elements that invade host genomes in a wide range of organisms. In metazoans, however, group I introns are extremely rare, so far only identified within mitogenomes of hexacorals and some sponges. We sequenced the complete mitogenome of the cold-water scleractinian coral Lophelia pertusa, the dominating deep sea reef-building coral species in the North Atlantic Ocean. The mitogenome (16,150 bp) has the same gene content but organized in a unique gene order compared to that of other known scleractinian corals. A complex group I intron (6460 bp) inserted in the ND5 gene (position 717) was found to host seven essential mitochondrial protein genes and one ribosomal RNA gene. Phylogenetic analysis supports a vertical inheritance pattern of the ND5-717 intron among hexacoral mitogenomes with no examples of intron loss. Structural assessments of the Lophelia intron revealed an unusual organization that lacks the universally conserved ωG at the 3′ end, as well as a highly compact RNA core structure with overlapping ribozyme and protein coding capacities. Based on phylogenetic and structural analyses we reconstructed the evolutionary history of ND5-717, from its ancestral protist origin, through intron loss in some early metazoan lineages, and into a compulsory feature with functional implications in hexacorals.     

Variation of nuclear and mitochondrial DNAs in Korean and Chinese isolates of Clonorchis sinensis

We compared the DNA sequence difference of isolates of Clonorchis sinensis from one Korean (Kimhae) and two Chinese areas (Guangxi and Shenyang). The sequences of nuclear rDNA (18S, internal transcribed spacer 1 and 2: ITS1 and ITS2) and mitochondrial DNA (cytochrome c oxidase subunit 1: cox1) were compared. A very few intraspecific nucleotide substitution of the 18S, ITS1, ITS2 and cox1 was found among three isolates of C. sinensis and a few nucleotide insertion and deletion of ITS1 were detected. The 18S, ITS1, ITS2 and cox1 sequences were highly conserved among three isolates. These findings indicated that the Korean and two Chinese isolates are similar at the DNA sequence level.     

The De Novo Assembly of Mitochondrial Genomes of the Extinct Passenger Pigeon (Ectopistes migratorius) with Next Generation Sequencing

The information from ancient DNA (aDNA) provides an unparalleled opportunity to infer phylogenetic relationships and population history of extinct species and to investigate genetic evolution directly. However, the degraded and fragmented nature of aDNA has posed technical challenges for studies based on conventional PCR amplification. In this study, we present an approach based on next generation sequencing to efficiently sequence the complete mitochondrial genome (mitogenome) of two extinct passenger pigeons (Ectopistes migratorius) using de novo assembly of massive short (90 bp), paired-end or single-end reads. Although varying levels of human contamination and low levels of postmortem nucleotide lesion were observed, they did not impact sequencing accuracy. Our results demonstrated that the de novo assembly of shotgun sequence reads could be a potent approach to sequence mitogenomes, and offered an efficient way to infer evolutionary history of extinct species.     

The complete mitochondrial genome of the wild silkworm moth, Actias selene

The complete mitochondrial genome (mitogenome) of Actias selene (Lepidoptera: Saturniidae) was determined to be 15,236 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The arrangement of 13 PCGs was similar to that of other sequenced lepidopterans. The AT skew of the mitogenome of A. selene was slightly negative, indicating a higher number of T compared to A nucleotides. The nucleotide composition of the mitogenome of A. selene was also biased toward A+T nucleotides (78.91%). All PCGs were initiated by ATN codons, except for the gene encoding cytochrome c oxidase subunit 1 (cox1), which may be initiated by the TTAG, as observed in other lepidopterans. Three genes, including cox1, cox2, and nad5, had incomplete stop codons consisting of just a T. With an exception for trnS1(AGN), all the other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome of A. selene was 339 bp in length, and contains several features common to the Lepidopteras, including non-repetitive sequences, a conserved structure combining the motif ATAGA and an 18-bp poly-T stretch and a poly-A element upstream of trnM gene. Phylogenetic analysis showed that A. selene was close to Saturniidae.     

Navigating the tip of the genomic iceberg: Next-generation sequencing for plant systematics

? Premise of the study: Just as Sanger sequencing did more than 20 years ago, next-generation sequencing (NGS) is poised to revolutionize plant systematics. By combining multiplexing approaches with NGS throughput, systematists may no longer need to choose between more taxa or more characters. Here we describe a genome skimming (shallow sequencing) approach for plant systematics. ? Methods: Through simulations, we evaluated optimal sequencing depth and performance of single-end and paired-end short read sequences for assembly of nuclear ribosomal DNA (rDNA) and plastomes and addressed the effect of divergence on reference-guided plastome assembly. We also used simulations to identify potential phylogenetic markers from low-copy nuclear loci at different sequencing depths. We demonstrated the utility of genome skimming through phylogenetic analysis of the Sonoran Desert clade (SDC) of Asclepias (Apocynaceae). ? Key results: Paired-end reads performed better than single-end reads. Minimum sequencing depths for high quality rDNA and plastome assemblies were 40× and 30×, respectively. Divergence from the reference significantly affected plastome assembly, but relatively similar references are available for most seed plants. Deeper rDNA sequencing is necessary to characterize intragenomic polymorphism. The low-copy fraction of the nuclear genome was readily surveyed, even at low sequencing depths. Nearly 160000 bp of sequence from three organelles provided evidence of phylogenetic incongruence in the SDC. ? Conclusions: Adoption of NGS will facilitate progress in plant systematics, as whole plastome and rDNA cistrons, partial mitochondrial genomes, and low-copy nuclear markers can now be efficiently obtained for molecular phylogenetics studies.