Genome skimming and exploration of DNA barcodes for Taiwan endemic cypresses

Cypresses are characterized by their longevity and valuable timber. In Taiwan, two endemic cypress species, Chamaecyparis formosensis and C. obtusa var. formosana, are threatened by prevalent illegal logging. A DNA barcode system is urgently needed for reforestation and conservation of these two cypresses. In this study, both plastomes and 35S rDNAs from 16, 10, and 6 individuals of C. formosensis, C. obtusa var. formosana, and C. obtusa var. obtusa were sequenced, respectively. We show that the loss of plastid trnT-GGU readily distinguishes C. formosensis from its congeneric species. We demonstrate that entire sequences of plastomes or 35S rDNAs are capable of correctly identifying cypress species and varieties, suggesting that they are effective super-barcodes. We also discover three short hypervariable loci (i.e., 3′ETS, ITS1, and trnH-psbA) that are promising barcodes for identifying cypress species and varieties. Moreover, nine species-specific indels of > 100 bp were detected in the cypress plastomes. These indels, together with the three aforementioned short barcodes, constitute an alternative and powerful barcode system crucial for identifying specimens that are fragmentary or contain degraded/poor DNA. Our sequenced data and barcode systems not only enrich the genetic reference for cypresses, but also contribute to future reforestation, conservation, and forensic investigations.

Plastome sequences help to improve the systematic position of trinerved Lindera species in the family Lauraceae

Lindera is a genus (c. 100 spp.) of trees belonging to the “core Laureae” group in the family Lauraceae. It is often confused with Litsea, and the systematics of the genus is unclear. Here, total 10 complete plastomes from nine trinerved Lindera species and another species Lindera obtusiloba (sect. Palminerviae Meissn.) were sequenced. Nine highly variable regions, trnH-GUG/psbA, psbM/trnD-GUC, petA/psbL, ndhF, trnL-UAG/ndhD, and ycf1, were identified among the 10 Lindera species. In addition, a total of 1,836 mutation events including six micro-inversions, 156 indels, and 1,674 substitutions, were also summarized. Comparing our sequences with other available plastomes in the “core Laureae,” we put forward that six hypervariable loci, trnH-GUG/psbA, ndhF, ndhF/rpl32, trl32/trnL-UAG, ndhD, and ycf1, could potentially be used as plastid barcode candidates for species identification. Further phylogenetic analyses were conducted using 49 complete Lauraceae plastomes. The results supported a close relationship among trinerved Lindera species and suggested an improved trinerved group comprising species of trinerved Lindera species and Iteadaphne caudate.

Genetic Polymorphism Study on Aedes albopictus of Different Geographical Regions Based on DNA Barcoding

Aedes albopictus is a very important vector for pathogens of many infectious diseases including dengue fever. In this study, we explored the genetic polymorphism of Aedes albopictus strains in different geographical regions using DNA barcoding of mitochondrial COI (MT-COI) gene. We collected MT-COI sequence of 106 Aedes albopictus mosquitos from 6 provinces in China including Fujian, Guangdong, Hainan, Yunnan, and Taiwan. The length of the sequences is 709bp with the content of A+T (67.7%) greater than that of G+C (32.3%). We identified mutations in 90 (13.68%) loci, of which 57 (63.33%) are transitions, 28 (31.11%) are transversions, and 5 (5.56%) are hypervariable loci. In addition, we obtained 42 haplotypes, 4 (9.52%) of which are shared among different populations. The haplotype diversity of Aedes albopictus is 0.882 and nucleotide diversity is 0.01017. Moreover, the pedigree network diagram shows that most haplotypes are under parallel evolution, suggesting a local expansion of Aedes albopictus in history. Finally, the Neighbor-Joining tree of MT-COI haplotypes reveals a certain correlation between haplotype clusters and geographical distribution, and there are differences among Aedes albopictus in different geographical regions. In conclusion, DNA barcoding of MT-COI gene is an effective method to study the genetic structure of Aedes albopictus.

Prediction of Local Transmission of Mycobacterium tuberculosis Isolates of a Predominantly Beijing Lineage by Use of a Variable-Number Tandem-Repeat Typing Method Incorporating a Consensus Set of Hypervariable Loci

ABSTRACT Strain genotyping based on the variable-number tandem repeat (VNTR) is widely applied for identifying the transmission of Mycobacterium tuberculosis. A consensus set of four hypervariable loci (1982, 3232, 3820, and 4120) has been proposed to improve the discrimination of Beijing lineage strains. Herein, we evaluated the utility of these four hypervariable loci for tracing local tuberculosis transmission in 981 cases over a 14-month period in Japan (2010 to 2011). We used six different VNTR systems, with or without the four hypervariable loci. Patient ages and weighted standard distances (a measure of the dispersion of genotype-clustered cases) were used as proxies for estimating local tuberculosis transmission. The highest levels of isolate discrimination were achieved with VNTR systems that incorporated the four hypervariable loci (i.e., the Japan Anti-Tuberculosis Association [JATA]18-VNTR, mycobacterial interspersed repetitive unit [MIRU]28-VNTR, and 24Beijing-VNTR). The clustering rates by JATA12-VNTR, MIRU15-VNTR, JATA15-VNTR, JATA18-VNTR, MIRU28-VNTR, and 24Beijing-VNTR systems were 52.2%, 51.0%, 39.0%, 24.1%, 23.1%, and 22.0%, respectively. As the discriminative power increased, the median weighted standard distances of the clusters tended to decrease (from 311 to 80 km, P < 0.001, Jonckheere-Terpstra trend test). Concurrently, the median ages of patients in the clusters tended to decrease (from 68 to 60 years, P < 0.001, Jonckheere-Terpstra trend test). These findings suggest that strain typing using the four hypervariable loci improves the prediction of active local tuberculosis transmission. The four-locus set can therefore contribute to the targeted control of tuberculosis in settings with high prevalence of Beijing lineage strains.

Evolutionary and functional implications of hypervariable loci within the skin virome

Localized genomic variability is crucial for the ongoing conflicts between infectious microbes and their hosts. An understanding of evolutionary and adaptive patterns associated with genomic variability will help guide development of vaccines and antimicrobial agents. While most analyses of the human microbiome have focused on taxonomic classification and gene annotation, we investigated genomic variation of skin-associated viral communities. We evaluated patterns of viral genomic variation across 16 healthy human volunteers. Human papillomavirus (HPV) and Staphylococcus phages contained 106 and 465 regions of diversification, or hypervariable loci, respectively. Propionibacterium phage genomes were minimally divergent and contained no hypervariable loci. Genes containing hypervariable loci were involved in functions including host tropism and immune evasion. HPV and Staphylococcus phage hypervariable loci were associated with purifying selection. Amino acid substitution patterns were virus dependent, as were predictions of their phenotypic effects. We identified diversity generating retroelements as one likely mechanism driving hypervariability. We validated these findings in an independently collected skin metagenomic sequence dataset, suggesting that these features of skin virome genomic variability are widespread. Our results highlight the genomic variation landscape of the skin virome and provide a foundation for better understanding community viral evolution and the functional implications of genomic diversification of skin viruses.

Evolutionary and Functional Implications of Hypervariable Loci Within the Skin Virome

ABSTRACTLocalized genomic variability is crucial for the ongoing conflicts between infectious microbes and their hosts. An understanding of evolutionary and adaptive patterns associated with genomic variability will help guide development of vaccines and anti-microbial agents. While most analyses of the human microbiome have focused on taxonomic classification and gene annotation, we investigated genomic variation of skin-associated viral communities. We evaluated patterns of viral genomic variation across 16 healthy human volunteers. HPV and Staphylococcus phages contained 106 and 465 regions of diversification, or hypervariable loci, respectively. Propionibacterium phage genomes were minimally divergent and contained no hypervariable loci. Genes containing hypervariable loci were involved in functions including host tropism and immune evasion. HPV and Staphylococcus phage hypervariable loci were associated with purifying selection. Amino acid substitution patterns were virus dependent, as were predictions of their phenotypic effects. We identified diversity generating retroelements as one likely mechanism driving hypervariability. We validated these findings in an independently collected skin metagenomic sequence dataset, suggesting that these features of skin virome genomic variability are widespread. Our results highlight the genomic variation landscape of the skin virome and provide a foundation for better understanding community viral evolution and the functional implications of genomic diversification of skin viruses.