Parallel tagged amplicon sequencing of relatively long PCR products using the Illumina HiSeq platform and transcriptome assembly

Changes in transcriptome regulations of a marine rotifer Brachionus plicatilis under methylmercury stress

Mercury (Hg), a heavy metal pollutant worldwide, can be transformed into methylmercury (MeHg) by various aquatic microorganisms in water, thus accumulating along the aquatic food chain and posing a particular challenge to human health. Zooplankton plays a crucial role in aquatic ecosystems and serves as a major component of the food chain. To evaluate the effects of MeHg on the rotifer Brachionus plicatilis and reveal the underlying mechanism of these effects, we exposed B. plicatilis to MeHg by either direct immersion or by feeding with MeHg-poisoned Chlorella pyrenoidesa, respectively, and conducted a transcriptomic analysis. The results showed that B. plicatilis directly exposed to MeHg by immersion showed significant enrichment of the glutathione metabolism pathway for detoxification of MeHg. In addition, the exposure to MeHg by feeding induced a significant enrichment of lysosome and notch signaling pathways of rotifers, supporting the hypothesis that MeHg can induce autophagy dysfunction in cells and disturb the nervous system of rotifers. In two different routes of MeHg exposure, the pathway of cytochrome P450 in rotifers showed significant enrichment for resisting MeHg toxicity. Our results suggest further studies on the potential mechanism and biological responses of MeHg toxicity in other links of the aquatic food chain.

Diversification of freshwater crabs on the sky islands in the Hengduan Mountains Region, China

The numerous naturally-fragmented sky islands (SIs) in the Hengduan Mountains Region (HMR) of southwestern China constitute discontinuous landscapes where montane habitats are isolated by dry-hot valleys which have fostered exceptional species diversification and endemicity. However, studies documenting the crucial role of SI on the speciation dynamics of native freshwater organisms are scarce. Here we used a novel set of comprehensive genetic markers (24 nuclear DNA sequences and complete mitogenomes), morphological characters, and biogeographical information to reveal the evolutionary history and speciation mechanisms of a group of small-bodied montane potamids in the genus Tenuipotamon. Our results provide a robustly supported phylogeny, and suggest that the vicariance events of these montane crabs correlate well with the emergence of SIs due to the uplift of the HMR during the Late Oligocene. Furthermore, ancestrally, mountain ridges provided corridors for the dispersal of these montane crabs that led to the colonization of moist montane-specific habitats, aided by past climatic conditions that were the crucial determinants of their evolutionary history. The present results illustrated that the mechanisms isolating SIs are reinforced by the harsh-dry isolating climatic features of dry-hot valleys separating SIs and continue to affect local diversification. This offers insights into the causes of the high biodiversity and endemism shown by the freshwater crabs of the HMR-SIs in southwestern China.

Lactobacillus acidophilus JYLA-126 Ameliorates Obesity-Associated Metabolic Disorders by Positively Regulating the AMPK Signaling Pathway Through the Gut-Liver Axis

Obesity is a chronic metabolic disease worldwide and is considered a major health problem in contemporary society. Lactobacillus acidophilus have demonstrated beneficial effects on obesity, but the specific mechanism of how it exerts beneficial effects has not been elucidated. Here, we found that L. acidophilus JYLA-126 had good biological properties for intestinal health, such as antioxidation, acid tolerance, bile salt tolerance, antimicrobial activity, and gut colonization. We further identified that supplementation of L. acidophilus JYLA-126 obese mice possessed a dose-dependent amelioration of body weight, intestinal imbalance, and metabolic disorders compared to HFD-induced mice. Mechanistically, the excellent slimming effect of L. acidophilus JYLA-126 was achieved mainly by reversing HFD-induced gut dysbiosis, inhibiting inflammatory factors and balancing the homeostasis of the gut-liver axis. Specifically, L. acidophilus JYLA-126 improved hepatic glycogen synthesis, lowered oxidative stress, and facilitated lipid metabolism by regulating AMPK signaling pathway-related protein expression to restore the overall metabolic level. Accordingly, L. acidophilus JYLA-126 promoted energy uptake efficiency in obese mice, resulting in significant expression of uncoupling protein 1 (UCP1) protein in brown adipose tissue (BAT), and markedly reduced the size of adipocytes. These findings indicate that the anti-obesity activity of L. acidophilus JYLA-126 correlates with activation of the AMPK signaling pathway through improved gut-liver interactions.

Combined transcriptomic and metabolomic analyses of high temperature stress response of quinoa seedlings

BACKGROUND

Quinoa (Chenopodium quinoa Willd.) originates in high altitude areas, such as the Andes, and has some inherent characteristics of cold, drought, and salinity tolerance, but is sensitive to high temperature.

RESULTS

To gain insight into the response mechanism of quinoa to high temperature stress, we conducted an extensive targeted metabolomic study of two cultivars, Dianli-3101 and Dianli-3051, along with a combined transcriptome analysis. A total of 794 metabolites and 54,200 genes were detected, in which the genes related to photosynthesis were found down-regulated at high temperatures, and two metabolites, lipids and flavonoids, showed the largest changes in differential accumulation. Further analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and transcription factors revealed that quinoa inhibits photosynthesis at high temperatures, and the possible strategies being used for high temperature stress management are regulation of heat stress transcription factors (HSFs) to obtain heat tolerance, and regulation of purine metabolism to enhance stress signals for rapid response to high temperature stress. The tolerant genotype could have an enhanced response through lower purine levels. The induction of the stress response could be mediated by HSF transcription factors. The results of this study may provide theoretical references for understanding the response mechanism of quinoa to high temperature stress, and for screening potential high temperature tolerant target genes and high temperature tolerant strains.

CONCLUSIONS

These findings reveal the regulation of the transcription factor family HSF and the purinergic pathway in response to high temperature stress to improve quinoa varieties with high temperature tolerance.

Assessing the Diet of a Predator Using a DNA Metabarcoding Approach

The diet of top predators is vital information needed to determine their ecological function and for their conservation management. However, the elusive habit and low population density of many predators constrains determination of their diets. While the morphological identification of scat contents is the traditional method, DNA metabarcoding has lately proven a more efficient and accurate method of identifying prey taxa. We applied DNA metabarcoding to analyzing the diet of the Eurasian otter (Lutra lutra), a top predator in freshwater ecosystems, using 12S and 16S rRNA mitochondrial primers target vertebrate prey. Diet did not vary among different data removal thresholds of 0.1, 1, 3, and 5%, comprising fishes (>90%), amphibians and birds (>2%), and occasionally mammals (<2%). Both 12S and 16S primers revealed similar otter diets, indicating that a single set of primers with a higher threshold is cost-effective for detecting the main prey taxa. Using 12S primers and a 5% threshold, we found no seasonal variation of otter diet in the Tangjiahe National Nature Reserve. A different prey community was found outside the reserve, which resulted in different prey composition for otters. However, prey taxon richness was not different between otters in- and outside the reserve. Otters preferred Schizothorax spp., the largest-sized fish species in the reserve, whereas they mainly preyed on Triplophysa bleekeri, a small-sized fish species, outside the reserve. Otters’ flexible feeding strategy reflect their high adaptability. However, greater human disturbance outside the reserve may present significant challenges to otters by altering prey communities and reducing prey profitability. Combining fecal DNA metabarcoding and local fish survey will provide opportunities for more detailed studies on the impact of different levels of human disturbances on prey communities and otters.

Transcriptome and Metabolome Analysis of the Synthesis Pathways of Allelochemicals in Eupatorium adenophorum

Eupatorium adenophorum (Crofton weed) is an invasive weed in more than 30 countries. It inhibits the growth of surrounding plants by releasing allelochemicals during its invasion. However, the synthetic pathways and molecular mechanisms of its allelochemicals have been rarely reported. In this study, the related genes and pathways of allelochemicals in E. adenophorum were analyzed. Transcriptome analysis showed that differentially expressed genes (DEGs) were mainly enriched in the phenylpropanoid biosynthetic pathway and flavonoid biosynthetic pathway. Thirty-three DEGs involved in the synthesis of allelochemicals were identified, and 30 DEGs showed significant differences in blades and stems. Six allelochemicals were identified from blades and stems by ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Correlation analysis of genes and metabolites showed a strong correlation between the five genes and allelochemicals. In addition, this study supplemented the biosynthetic pathway of Eupatorium adenophorum B (HHO). It was found that acyclic sesquiterpene synthase (NES), δ-cadinene synthase (TPS), and cytochrome P450 (P450) were involved in the synthesis of HHO. These findings provide a dynamic spectrum consisting of allelochemical metabolism and a coexpression network of allelochemical synthesis genes in E. adenophorum.

Transcriptomic Analysis Exploring the Molecular Mechanisms of Hanchuan Zupa Granules in Alleviating Asthma in Rat

Objective

To investigate the molecular mechanisms of HCZP treatment of asthma.

Materials and Methods

Thirty Sprague Dawley (SD) rats were divided into normal, asthma, and HCZP groups (n = 10). The asthma model was sensitized by 1 mg ovalbumin (OVA)/aluminum hydroxide Al(OH)₃mixture and then challenged with 1% aerosolized OVA for four weeks. Rats in the HCZP group received 10.08 g/kg/d HCZP for four weeks during OVA challenge. Then, lung tissues of rats in each group were collected for RNA sequencing. Moreover, the expression level of some core genes was detected by using western blotting and immunohistochemistry.

Results

Inflammatory cell infiltration and pathological damage of the lungs improved in the HCZP group. Compared with the asthma group (0.049 ± 0.002 mm²/mm; 0.036 ± 0.006 mm²/mm; and 0.014 ± 0.001 mm²/mm), total wall thickness (0.042 ± 0.001 mm²/mm), inner wall thickness (0.013 ± 0.001 mm²/mm), and smooth muscle layer thickness (0.012 ± 0.001 mm²/mm) significantly decreased in the HCZP group. Bioinformatics analysis showed that hub genes such as bradykinin receptor B2 (Bdkrb2) and CD4 molecule (Cd4) had different expression patterns between model and HCZP groups. Two transcription factors, forkhead box Q1 (Foxq1) and nuclear factor of activated T cells 2 (Nfatc2), served important regulatory roles in asthma. Compared with the model group, Bdkrb2 protein expression increased and Nfatc2 protein expression decreased in the HCZP group. Discussion and Conclusion. HCZP could alleviate asthma via regulating the expression of several hub genes, which might serve as therapeutic targets for asthma. However, the mechanism of these genes will be studied in the future.

Plant-herbivorous insect networks: who is eating what revealed by long barcodes using high-throughput sequencing and Trinity assembly

Interactions between plants and insects are among the most important life functions for all organism at a particular natural community. Usually a large number of samples are required to identify insect diets in food web studies. Previously, Sanger sequencing and next generation sequencing (NGS) with short DNA barcodes were used, resulting in low species-level identification; meanwhile the costs of Sanger sequencing are expensive for metabarcoding together with more samples. Here, we present a fast and effective sequencing strategy to identify larvae of Lepidoptera and their diets at the same time without increasing the cost on Illumina platform in a single HiSeq run, with long-multiplex-metabarcoding (COI for insects, rbcL, matK, ITS and trnL for plants) obtained by Trinity assembly (SHMMT). Meanwhile, Sanger sequencing (for single individuals) and NGS (for polyphagous) were used to verify the reliability of the SHMMT approach. Furthermore, we show that SHMMT approach is fast and reliable, with most high-quality sequences of five DNA barcodes of 63 larvae individuals (54 species) recovered (full length of 100% of the COI gene and 98.3% of plant DNA barcodes) using Trinity assembly (up-sized to 1015 bp). For larvae diets identification, 95% are reliable; the other 5% failed because their guts were empty. The diets identified by SHMMT approach are 100% consistent with the host plants that the larvae were feeding on during our collection. Our study demonstrates that SHMMT approach is reliable and cost-effective for insect-plants network studies. This will facilitate insect-host plant studies that generally contain a huge number of samples.

New insights into the phylogeny and evolution of lady beetles (Coleoptera: Coccinellidae) by extensive sampling of genes and species

Ladybirds (family Coccinellidae) are one of the most diverse groups of beetles and globally comprise over 6000 species. Despite their scientific and economic significance, the taxonomy of Coccinellidae remains unstable, and we still know little about their evolutionary history. By using a small number of genes, previous phylogenetic analyses have not reliably resolved the relationships among major ladybird lineages. In this study, we sequenced 94 nuclear protein-coding genes for 214 species of Coccinellidae and 14 outgroups, covering 90 genera and 35 tribes. We found that nucleotide compositional heterogeneity is present among ladybird tribes so that phylogenetic inference at the amino acid level is more reliable than at the DNA level. Based on the maximum likelihood analyses of the amino acid dataset, we recognize three subfamilies in Coccinellidae: Microweiseinae, Monocoryninae stat. nov., and Coccinellinae. The subfamily relationships are strongly supported as (Microweiseinae, (Monocoryninae stat. nov., Coccinellinae)). The tribes of ladybirds are mostly monophyletic, except Ortaliini, Sticholotidini, Scymnini, and Coccidulini. The phylogenetic relationships among tribes of Coccinellinae are still not well resolved, with many nodes weakly supported. Our divergence time analysis suggests that the crown group of extant lady beetles arose in the Early Cretaceous ~ 143 million years ago (Mya) and experienced a rapid diversification during the Late Cretaceous (120-70 Mya). We hypothesize that the boom of angiosperms in the Late Cretaceous promoted the diversification of herbivorous Sternorrhyncha insects, especially aphids, which in turn drove the rapid radiation of predatory lady beetles. In summary, our work provides a comprehensive time-calibrated phylogeny of Coccinellidae that provides a sound framework for revising their classification and understanding the origin of their biodiversity.

A large-scale systematic framework of Chinese snakes based on a unified multilocus marker system

Snakes are one of the most diverse groups of terrestrial vertebrates, with approximately 3500 extant species. A robust phylogeny and taxonomy of snakes is crucial for us to know, study and protect them. For a large group such as snakes, broad-scale phylogenetic reconstructions largely rely on data integration. Increasing the compatibility of the data from different researches is thus important, which can be facilitated by standardization of the loci used in systematic analyses. In this study, we proposed a unified multilocus marker system for snake systematics by conflating 5 mitochondrial markers, 19 vertebrate-universal nuclear protein coding (NPC) markers and 72 snake-specific noncoding intron markers. This marker system is an addition to the large squamate conserved locus set (SqCL) for studies preferring a medium-scale data set. We applied this marker system to over 440 snake samples and constructed the currently most comprehensive systematic framework of the snakes in China. Robust snake phylogenetic relationships were recovered at both deep and shallow evolutionary depths, demonstrating the usefulness of this multilocus marker system. Discordance was revealed by a parallel comparison between the snake tree based on the multilocus marker system and that based on only the mitochondrial loci, highlighting the necessity of using multiple types of markers to better understand the snake evolutionary histories. The divergence times of different snake groups were estimated with the nuclear data set. Our comprehensive snake tree not only confirms many important nodes inferred in previous studies but also contributes new insights into many snake phylogenetic relationships. Suggestions are made for the current Chinese snake taxonomy.

Sequence capture across large phylogenetic scales by using pooled PCR-generated baits: A case study of Lepidoptera

Sequence capture across large phylogenetic scales is not easy because hybridization capture is only effective when the genetic distance between the bait and target is small. Here, we propose a simple but effective strategy to tackle this issue: pooling DNA from a number of selected representative species of different clades to prepare PCR-generated baits to minimize the genetic distance between the bait and target. To demonstrate the utility of this strategy, we newly developed a set of universal nuclear markers (including 94 nuclear protein-coding genes) for Lepidoptera, a superdiverse insect group. We used a DNA pool from six lepidopteran species (representing six superfamilies) to prepare PCR baits for the 94 markers. These homemade PCR baits were used to capture sequence data from 43 species of 17 lepidopteran families, and 94% of the target loci were recovered. We constructed two data sets from the obtained data (one containing ~90 kb target coding sequences and the other containing ~120 kb target + flanking coding sequences). Both data sets yielded highly similar and well-resolved trees with 90% of nodes having >95% bootstrap support. Our capture experiment indicated that using DNA mixtures pooled from different clade-representative species of Lepidoptera to prepare PCR baits can reliably capture a large number of targeted nuclear markers across different Lepidoptera lineages. We hope that this newly developed nuclear marker set will serve as a new phylogenetic tool for Lepidoptera phylogenetics, and the PCR bait preparation strategy can facilitate the application of sequence capture techniques by researchers to accelerate data collection.

Evolutionary history of Coleoptera revealed by extensive sampling of genes and species

Beetles (Coleoptera) are the most diverse and species-rich group of insects, and a robust, time-calibrated phylogeny is fundamental to understanding macroevolutionary processes that underlie their diversity. Here we infer the phylogeny and divergence times of all major lineages of Coleoptera by analyzing 95 protein-coding genes in 373 beetle species, including ~67% of the currently recognized families. The subordinal relationships are strongly supported as Polyphaga (Adephaga (Archostemata, Myxophaga)). The series and superfamilies of Polyphaga are mostly monophyletic. The species-poor Nosodendridae is robustly recovered in a novel position sister to Staphyliniformia, Bostrichiformia, and Cucujiformia. Our divergence time analyses suggest that the crown group of extant beetles occurred ~297 million years ago (Mya) and that ~64% of families originated in the Cretaceous. Most of the herbivorous families experienced a significant increase in diversification rate during the Cretaceous, thus suggesting that the rise of angiosperms in the Cretaceous may have been an 'evolutionary impetus' driving the hyperdiversity of herbivorous beetles.

Characterization of Brachyspira communities from clinical cases of swine mucohaemorrhagic diarrhea through deep sequencing of the NADH oxidase (nox) gene

Swine dysentery is traditionally associated with Brachyspira hyodysenteriae, but the re-emergence of Brachyspira-associated disease in North America associated with a novel causative species, B. hampsonii, is now a concern for swine producers. The pathogenesis of Brachyspira-associated disease is not completely understood, and it is not known whether mixed infections of Brachyspira spp. are important in disease development. Deep sequencing of partial sequences of the nox gene amplified with genus-specific primers was used to detect Brachyspira spp. in 55 fecal samples from clinical cases of mucohaemorrhagic diarrhea in pigs from Western Canada that had been identified as positive for one or more Brachyspira species using established diagnostic tests. Synthetic mixtures of Brachyspira genomic DNA were included in the study to define detection limits for the technique and identify biases in detection of different species. Multiple species were detected in all clinical cases for which sufficient nox sequence data were generated (n = 47), indicating that mixed species Brachyspira infections are common, although in most cases, one species accounted for at least half of the sequences identified. In all cases, the species detected in the original diagnostic investigation of each case was also detected by nox sequencing. Results from synthetic communities indicated that the method was highly reproducible, but also indicated potential PCR bias against B. hampsonii genomovar I. Deep sequencing of the nox gene target is a suitable method for simultaneous detection of multiple Brachyspira species in clinical case material that may offer advantages over current, more targeted diagnostic approaches for investigating the significance of mixed infections in disease development.

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein-coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single-copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next-generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p-distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.

Phylogenomics reveals rapid, simultaneous diversification of three major clades of Gondwanan frogs at the Cretaceous-Paleogene boundary

Frogs (Anura) are one of the most diverse groups of vertebrates and comprise nearly 90% of living amphibian species. Their worldwide distribution and diverse biology make them well-suited for assessing fundamental questions in evolution, ecology, and conservation. However, despite their scientific importance, the evolutionary history and tempo of frog diversification remain poorly understood. By using a molecular dataset of unprecedented size, including 88-kb characters from 95 nuclear genes of 156 frog species, in conjunction with 20 fossil-based calibrations, our analyses result in the most strongly supported phylogeny of all major frog lineages and provide a timescale of frog evolution that suggests much younger divergence times than suggested by earlier studies. Unexpectedly, our divergence-time analyses show that three species-rich clades (Hyloidea, Microhylidae, and Natatanura), which together comprise ∼88% of extant anuran species, simultaneously underwent rapid diversification at the Cretaceous-Paleogene (K-Pg) boundary (KPB). Moreover, anuran families and subfamilies containing arboreal species originated near or after the KPB. These results suggest that the K-Pg mass extinction may have triggered explosive radiations of frogs by creating new ecological opportunities. This phylogeny also reveals relationships such as Microhylidae being sister to all other ranoid frogs and African continental lineages of Natatanura forming a clade that is sister to a clade of Eurasian, Indian, Melanesian, and Malagasy lineages. Biogeographical analyses suggest that the ancestral area of modern frogs was Africa, and their current distribution is largely associated with the breakup of Pangaea and subsequent Gondwanan fragmentation.

Replacing Sanger with Next Generation Sequencing to improve coverage and quality of reference DNA barcodes for plants

We estimate the global BOLD Systems database holds core DNA barcodes (rbcL + matK) for about 15% of land plant species and that comprehensive species coverage is still many decades away. Interim performance of the resource is compromised by variable sequence overlap and modest information content within each barcode. Our model predicts that the proportion of species-unique barcodes reduces as the database grows and that 'false' species-unique barcodes remain >5% until the database is almost complete. We conclude the current rbcL + matK barcode is unfit for purpose. Genome skimming and supplementary barcodes could improve diagnostic power but would slow new barcode acquisition. We therefore present two novel Next Generation Sequencing protocols (with freeware) capable of accurate, massively parallel de novo assembly of high quality DNA barcodes of >1400 bp. We explore how these capabilities could enhance species diagnosis in the coming decades.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.