COPE: an accurate k-mer-based pair-end reads connection tool to facilitate genome assembly

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the diversity in RNA size and chemical modification makes it difficult to capture all RNAs in a cell. We developed a method that combines quasi-random priming with template switching to construct sequencing libraries from RNA molecules of any length and with any type of 3' modifications, allowing for the sequencing of virtually all RNA species. Our ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to identify and quantify all classes of coding and non-coding RNAs. With LIDAR, we comprehensively characterized the transcriptomes of mouse embryonic stem cells, neural progenitor cells, mouse tissues, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared with traditional ligation-dependent sequencing methods and uncovered tDRs with blocked 3' ends that had previously escaped detection. Therefore, LIDAR can capture all RNAs in a sample and uncover RNA species with potential regulatory functions.

CAREx: context-aware read extension of paired-end sequencing data

BACKGROUND

Commonly used next generation sequencing machines typically produce large amounts of short reads of a few hundred base-pairs in length. However, many downstream applications would generally benefit from longer reads.

RESULTS

We present CAREx-an algorithm for the generation of pseudo-long reads from paired-end short-read Illumina data based on the concept of repeatedly computing multiple-sequence-alignments to extend a read until its partner is found. Our performance evaluation on both simulated data and real data shows that CAREx is able to connect significantly more read pairs (up to 99 % for simulated data) and to produce more error-free pseudo-long reads than previous approaches. When used prior to assembly it can achieve superior de novo assembly results. Furthermore, the GPU-accelerated version of CAREx exhibits the fastest execution times among all tested tools.

CONCLUSION

CAREx is a new MSA-based algorithm and software for producing pseudo-long reads from paired-end short read data. It outperforms other state-of-the-art programs in terms of (i) percentage of connected read pairs, (ii) reduction of error rates of filled gaps, (iii) runtime, and (iv) downstream analysis using de novo assembly. CAREx is open-source software written in C++ (CPU version) and in CUDA/C++ (GPU version). It is licensed under GPLv3 and can be downloaded at ( https://github.com/fkallen/CAREx ).

Determinants of genetic diversity in Neotropical salamanders (Plethodontidae: Bolitoglossini)

Genetic diversity is the raw material of evolution, yet the reasons why it varies among species remain poorly understood. While studies at deeper phylogenetic scales point to the influence of life history traits on genetic diversity, it appears to be more affected by population size but less predictable at shallower scales. We used proxies for population size, mutation rate, direct selection, and linked selection to test factors affecting genetic diversity within a diverse assemblage of Neotropical salamanders, which vary widely for these traits. We estimated genetic diversity of noncoding loci using ddRADseq and coding loci using RNAseq for an assemblage of Neotropical salamanders distributed from northern Mexico to Costa Rica. Using ddRADseq loci, we found no significant association with genetic diversity, while for RNAseq data we found that environmental heterogeneity and proxies of population size predict a substantial portion of the variance in genetic diversity across species. Our results indicate that diversity of coding loci may be more predictable than that of noncoding loci, which appears to be mostly unpredictable at shallower phylogenetic scales. Our results suggest that coding loci may be more appropriate for genetic diversity estimates used in conservation planning because of the lack of any association between the variables we used and genetic diversity of noncoding loci.

The co-presence of polystyrene nanoplastics and ofloxacin demonstrates combined effects on the structure, assembly, and metabolic activities of marine microbial community

Nanoplastic is increasing in environments and can address toxic effects on various organisms. Particle size, concentration, and surface functionalization most influence nanoplastic toxicity. Besides, nanoplastic can adsorb other contaminants (e.g., antibiotics) to aggravate its adverse effects. The combined effects of nanoplastics and antibiotics on planktonic/benthic microbial communities, however, are still largely unknown. In this study, the combined effects of polystyrene nanoplastic and ofloxacin on the structure, assembly, and metabolic activities of marine microbial communities were investigated based on amplicon sequencing data. The results mainly demonstrate that: (1) nanoplastic and ofloxacin have greater impacts on prokaryotic communities than eukaryotic ones; (2) niche breadths of planktonic prokaryotes and benthic eukaryotes were shrank with both high nanoplastic and ofloxacin concentrations; (3) increased ofloxacin mainly reduces nodes/edges of co-occurrence networks, while nanoplastic centralizes network modularity; (4) increased nanoplastic under high ofloxacin concentration induces more differential prokaryotic pathways in planktonic communities, while benthic communities are less influenced. The present work indicates that co-presence of nanoplastics and ofloxacin has synergistic combined effects on community structure shifts, niche breadth shrinking, network simplifying, and differential prokaryotic pathways inducing in marine microbial communities, suggesting nanoplastics and its combined impacts with other pollutions should be paid with more concerns.

Hybridization and Transgressive Evolution Generate Diversity in an Adaptive Radiation of Anolis Lizards

Interspecific hybridization may act as a major force contributing to the evolution of biodiversity. Although generally thought to reduce or constrain divergence between 2 species, hybridization can, paradoxically, promote divergence by increasing genetic variation or providing novel combinations of alleles that selection can act upon to move lineages toward new adaptive peaks. Hybridization may, then, play a key role in adaptive radiation by allowing lineages to diversify into new ecological space. Here, we test for signatures of historical hybridization in the Anolis lizards of Puerto Rico and evaluate 2 hypotheses for the role of hybridization in facilitating adaptive radiation-the hybrid swarm origins hypothesis and the syngameon hypothesis. Using whole genome sequences from all 10 species of Puerto Rican anoles, we calculated D and f-statistics (from ABBA-BABA tests) to test for introgression across the radiation and employed multispecies network coalescent methods to reconstruct phylogenetic networks that allow for hybridization. We then analyzed morphological data for these species to test for patterns consistent with transgressive evolution, a phenomenon in which the trait of a hybrid lineage is found outside of the range of its 2 parents. Our analyses uncovered strong evidence for introgression at multiple stages of the radiation, including support for an ancient hybrid origin of a clade comprising half of the extant Puerto Rican anole species. Moreover, we detected significant signals of transgressive evolution for 2 ecologically important traits, head length and toepad width, the latter of which has been described as a key innovation in Anolis. [Adaptive radiation; introgression; multispecies network coalescent; phenotypic evolution; phylogenetic network; reticulation; syngameon; transgressive segregation.].

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

B Lymphocyte Development in the Bursa of Fabricius of Young Broilers is Influenced by the Gut Microbiota

Chickens have been used as a valuable and traditional model for studies on basic immunology. B lymphocytes were first identified in the bursa of Fabricius (BF) of broilers. The microbiota is important for immune system development and function. However, the effect of the microbiota on mediating B cell development and its regulatory mechanism is poorly elucidated. Here, we show that the gut microbiota is associated with the development of bursal B cells in young chickens. Changing patterns of both the alpha diversity and the expression of the B cell marker Bu-1α in the gut microbiota were related to the ages of chickens at different growth phases. Further correlation analysis revealed the marked correlation between the relative abundances of Intestinimonas, Bilophila, Parasutterella, Bacteroides, Helicobacter, Campylobacter, and Mucispirillum and the expression of Bu-1α. In antibiotic-treated chickens, BF and B cell development had aberrations as the relative abundance of the microbiota in early life decreased. These findings were consistent with Spearman's correlation results. Single-cell transcriptome analysis indicated that the heterogeneity in the cellular composition and developmental trajectory of bursal B cells from antibiotic-treated chickens was large. We found a novel subpopulation of unnamed B cells and identified Taf1 as a new pivotal regulator of B cell lineage differentiation. Therefore, we provide novel insights into the regulatory role of the gut microbiota in B cell development in early life and the maturation of host humoral immunity. IMPORTANCE In this study, we used young broilers to investigate the relationship between their gut microbiota and bursal B cell development. We characterized the important variables, microbes, B cells, and immunoglobulins during the posthatch development of birds. We also identified several candidate taxa in the cecal contents associated with B cells. Our study provides a rich resource and cell-cell cross talk model supporting B cell differentiation from the bursa in vitro at single-cell resolution. Furthermore, we determined a new pivotal regulator (Taf1) of B cell differentiation. We believe that our study makes a significant contribution to the literature because our findings may elucidate the role of the gut microbiota in B cell differentiation. This study also serves as a basis for developing new strategies that modulate B cell differentiation to prevent diseases.

Application of T-cell receptor repertoire as a novel monitor in dynamic tracking and assessment: A cohort-study based on RA patients

T-cell receptor repertoire (TCRR) sequencing has been widely applied in many fields as a novel tool. This study explored characteristics of TCRR in detail with a cohort of 598 rheumatoid arthritis (RA) patients before and after anti-rheumatic treatments. We highlighted the abnormal TCRR distribution in RA characterized by decreased diversity and increased proportion of hyperexpanded clones (HECs), which was potentially attributed to skewed usage of global V/J segments but not a few certain ones. Enriched motifs analysis in RA community demonstrated the huge heterogeneity of CDR3 sequences, so that individual factors are strongly recommended to be taken into consideration when it comes to clinical application of TCRR. Disease-modifying antirheumatic drugs (DMARDs) can regulate immune system through recovery of TCRR richness to relieve symptoms. Remarkably, sensitive gene profile and advantageous gene profile were identified in this study as new biomarkers for different DMARDs regimens.

CMash: fast, multi-resolution estimation of k-mer-based Jaccard and containment indices

Motivation

K-mer-based methods are used ubiquitously in the field of computational biology. However, determining the optimal value of k for a specific application often remains heuristic. Simply reconstructing a new k-mer set with another k-mer size is computationally expensive, especially in metagenomic analysis where datasets are large. Here, we introduce a hashing-based technique that leverages a kind of bottom-m sketch as well as a k-mer ternary search tree (KTST) to obtain k-mer-based similarity estimates for a range of k values. By truncating k-mers stored in a pre-built KTST with a large k=kmax value, we can simultaneously obtain k-mer-based estimates for all k values up to k_max. This truncation approach circumvents the reconstruction of new k-mer sets when changing k values, making analysis more time and space-efficient.

Results

We derived the theoretical expression of the bias factor due to truncation. And we showed that the biases are negligible in practice: when using a KTST to estimate the containment index between a RefSeq-based microbial reference database and simulated metagenome data for 10 values of k, the running time was close to 10× faster compared to a classic MinHash approach while using less than one-fifth the space to store the data structure.

Availability and implementation

A python implementation of this method, CMash, is available at https://github.com/dkoslicki/CMash. The reproduction of all experiments presented herein can be accessed via https://github.com/KoslickiLab/CMASH-reproducibles.

Supplementary information

Supplementary data are available at Bioinformatics online.

Multi-omics analysis reveals gut microbiota-induced intramuscular fat deposition via regulating expression of lipogenesis-associated genes

The gut microbiome has great effects on the digestion, absorption, and metabolism of lipids. However, the microbiota composition that can alter the fat deposition and the meat quality of pigs remains unclear. Here, we used Laiwu (LW) pigs (a native Chinese breed with higher intramuscular fat) compared with commercial crossbreed Duroc × (Landrace × Yorkshire) (DLY) pigs to investigate the effects of microbiota on meat quality, especially in intramuscular fat content. A total of 32 DLY piglets were randomly allotted to 4 groups and transplanted with fecal microbiota from healthy LW pigs. The results indicated that the high dose of fecal microbiota transplantation (HFMT) selectively enhanced fat deposition in longissimus dorsi (P < 0.05) but decreased backfat thickness (P < 0.05) compared with control group. HFMT significantly altered meat color and increased feed conversation ratio (P < 0.05). Furthermore, the multi-omics analysis revealed that Bacteroides uniformis, Sphaerochaeta globosa, Hydrogenoanaerobacterium saccharovorans, and Pyramidobacter piscolens are the core species which can regulate lipid deposition. A total of 140 male SPF C57BL/6j mice were randomly allotted into 7 groups and administrated with these 4 microbes alone or consortium to validate the relationships between microbiota and lipid deposition. Inoculating the bacterial consortium into mice increased intramuscular fat content (P < 0.05) compared with control mice. Increased expressions of lipogenesis-associated genes including cluster of differentiation 36 (Cd36), diacylglycerol O-acyltransferase 2 (Dgat2), and fatty acid synthase (FASN) were observed in skeletal muscle in the mice with mixed bacteria compared with control mice. Together, our results suggest that the gut microbiota may play an important role in regulating the lipid deposition in the muscle of pigs and mice.

Effects of Multispecies Probiotic on Intestinal Microbiota and Mucosal Barrier Function of Neonatal Calves Infected With E. coli K99

Altered gut microbiota are implicated in inflammatory neonatal calf diarrhea caused by E. coli K99. Beneficial probiotics are used to modulate gut microbiota. However, factors that mediate host-microbe interactions remain unclear. We evaluated the effects of a combination of multispecies probiotics (MSP) on growth, intestinal epithelial development, intestinal immune function and microbiota of neonatal calves infected with E. coli K99. Twelve newborn calves were randomly assigned as follows: C (control, without MSP); D (E. coli O78:K99 + gentamycin); and P (E. coli O78:K99 + supplemental MSP). All groups were studied for 21 d. MSP supplementation significantly (i) changed fungal Chao1 and Shannon indices of the intestine compared with group D; (ii) reduced the relative abundance of Bacteroides and Actinobacteria, while increasing Bifidobacteria, Ascomycetes, and Saccharomyces, compared with groups C and D; (iii) improved duodenal and jejunal mucosal SIgA and total Short Chain Fatty Acids (SCFA) concentrations compared with group D; (iv) increased relative ZO-1 and occludin mRNA expression in jejunal mucosa compared with group D; and (v) enhanced intestinal energy metabolism and defense mechanisms of calves by reducing HSP90 expression in E. coli K99, thereby alleviating the inflammatory response and promoting recovery of mucosal function. Our research may provide direct theoretical support for future applications of MSP in ruminant production.

BigFiRSt: A Software Program Using Big Data Technique for Mining Simple Sequence Repeats From Large-Scale Sequencing Data

Background

Simple Sequence Repeats (SSRs) are short tandem repeats of nucleotide sequences. It has been shown that SSRs are associated with human diseases and are of medical relevance. Accordingly, a variety of computational methods have been proposed to mine SSRs from genomes. Conventional methods rely on a high-quality complete genome to identify SSRs. However, the sequenced genome often misses several highly repetitive regions. Moreover, many non-model species have no entire genomes. With the recent advances of next-generation sequencing (NGS) techniques, large-scale sequence reads for any species can be rapidly generated using NGS. In this context, a number of methods have been proposed to identify thousands of SSR loci within large amounts of reads for non-model species. While the most commonly used NGS platforms (e.g., Illumina platform) on the market generally provide short paired-end reads, merging overlapping paired-end reads has become a common way prior to the identification of SSR loci. This has posed a big data analysis challenge for traditional stand-alone tools to merge short read pairs and identify SSRs from large-scale data.

Results

In this study, we present a new Hadoop-based software program, termed BigFiRSt, to address this problem using cutting-edge big data technology. BigFiRSt consists of two major modules, BigFLASH and BigPERF, implemented based on two state-of-the-art stand-alone tools, FLASH and PERF, respectively. BigFLASH and BigPERF address the problem of merging short read pairs and mining SSRs in the big data manner, respectively. Comprehensive benchmarking experiments show that BigFiRSt can dramatically reduce the execution times of fast read pairs merging and SSRs mining from very large-scale DNA sequence data.

Conclusions

The excellent performance of BigFiRSt mainly resorts to the Big Data Hadoop technology to merge read pairs and mine SSRs in parallel and distributed computing on clusters. We anticipate BigFiRSt will be a valuable tool in the coming biological Big Data era.

Benefits of merging paired-end reads before pre-processing environmental metagenomics data

BACKGROUND

High throughput sequencing of environmental DNA has applications in biodiversity monitoring, taxa abundance estimation, understanding the dynamics of community ecology, and marine species studies and conservation. Environmental DNA, especially, marine eDNA, has a fast degradation rate. Aside from the good quality reads, the data could have a significant number of reads that fall slightly below the default PHRED quality threshold of 30 on sequencing. For quality control, trimming methods are employed, which generally precede the merging of the read pairs. However, in the case of eDNA, a significant percentage of reads within the acceptable quality score range are also dropped.

METHODS

To infer the ideal merge tool that is sensitive to eDNA, two Hiseq paired-end eDNA datasets were utilized to study the merging by the tools - FLASH (Fast Length Adjustment of SHort reads), PANDAseq, COPE, BBMerge, and VSEARCH without preprocessing. We assessed these tools on the following parameters: Time taken to process, the quality, and the number of merged reads. Trimmomatic, a widely-used preprocessing tool, was also assessed by preprocessing the datasets at different parameters for the two approaches of preprocessing: Sliding Window and Maximum Information. The preprocessed read pairs were then merged using the ideal merge tool identified earlier.

RESULTS

FLASH is the most efficient merge tool balancing data conservation, quality of reads, and processing time. We compared Trimmomatic's two quality trimming options with increasing strictness with FLASH's direct merge. The raw reads processed with Trimmomatic then merged, yielded a significant drop in reads compared to the direct merge. An average of 29% of reads was dropped when directly merged with FLASH. Maximum Information option resulted in 30.7% to 68.05% read loss with lowest and highest stringency parameters, respectively. The Sliding Window approach conserves approximately 10% more reads at a PHRED score of 25 set as the threshold for a window of size 4. The lowered PHRED cut off conserves about 50% of the reads that could potentially be informative. We noted no significant reduction of data while optimizing the number of reads read in a window with the ideal quality (Q) score.

CONCLUSIONS

Losing reads can negatively impact the downstream processing of the environmental data, especially for sequence alignment studies. The quality trim-first-merge-later approach can significantly decrease the number of reads conserved. However, direct merging of pair-end reads using FLASH conserved more than 60% of the reads. Therefore, direct merging of the paired-end reads can prevent potential removal of informative reads that do not comply by the trimming tool's strict checks. FLASH to be an efficient tool in conserving reads while carrying out quality trimming in moderation. Overall, our results show that merging paired-end reads of eDNA data before trimming can conserve more reads.

Transition of Serotype 35B Pneumococci From Commensal to Prevalent Virulent Strain in Children

In our community-based prospective cohort study in young children, we observed a significant increase in pneumococcal serotype 35B nasopharyngeal (NP) commensal colonization during the 2011–2014 timeframe, but these strains were not associated with disease. Beginning in 2015 and continuing through to the present, the serotype 35B virulence changed, and it became the dominant bacteria isolated and associated with pneumococcal acute otitis-media (AOM) in our cohort. We performed comparative analyses of 250 35B isolates obtained from 140 children collected between 2006 and 2019. Changes in prevalence, clonal-complex composition, and antibiotic resistance were analyzed. Seventy-two (29%) of 35B isolates underwent whole-genome sequencing to investigate genomic changes associated with the shift in virulence that resulted in increased rates of 35B-associated AOM disease. 35B strains that were commensals and AOM disease-causing were mainly associated with sequence type (ST) 558. Antibiotic concentrations of β-lactams and ofloxacin necessary to inhibit growth of 35B strains rose significantly (2006–2019) (p<0.005). However, only isolates from the 35B/ST558 showed significant increases in MIC₅₀ of penicillin and ofloxacin between the years 2006–2014 and 2015–2019 (p=0.007 and p<0.0001). One hundred thirty-eight SNPs located in 34 different genes were significantly associated with post-2015 strains. SNPs were found in nrdG (metal binding, 10%); metP and metN (ABC transporter, 9%); corA (Mg²⁺ transporter, 6%); priA (DNA replication, 5%); and on the enzymic gene ldcB (LD-carboxypeptidase, 3%). Pneumococcal serotype 35B strains was a common NP commensal during 2010–2014. In 2015, a shift in increasing number of AOM cases occurred in young children caused by 35B, that was associated with changes in genetic composition and antibiotic susceptibility.

Phylogeographic analysis delimits three evolutionary significant units of least chipmunks in North America and identifies unique genetic diversity within the imperiled Peñasco population

Although least chipmunks (Neotamias minimus) are a widely distributed North American species of least concern, the southernmost population, N. m. atristriatus (Peñasco least chipmunk), is imperiled and a candidate for federal listing as a subspecies. We conducted a phylogeographic analysis across the N. minimus range to assess genomic differentiation and distinctiveness of the N. m. atristriatus population. Additionally, we leveraged the historical component of sampling to conduct a temporal analysis of N. minimus genetic diversity and also considered climate change effects on range persistence probability by projecting a species distribution model into the IPCC5 RCP 2.6 and 8.5 scenarios. We identified three geographically structured groups (West, North, and South) that were supported by both mitochondrial and nuclear data. N. m. atristriatus grouped within a unique South subclade but were not reciprocally monophyletic from N. m. operarius, and nuclear genome analyses did not separate N. m. atristriatus, N. m. caryi, and N. m. operarius. Thus, while least chipmunks in the Southwest represent an evolutionary significant unit, subspecies distinctions were not supported and listing of the Peñasco population as a Distinct Population Segment of N. m. operarius may be warranted. Our results also support consideration of populations with North and West mitogenomes as two additional evolutionary significant units. We found that N. minimus genetic diversity declined by ~87% over the last century, and our models predicted substantial future habitat contraction, including the loss of the full contemporary ranges of N. m. atristriatus, N. m. arizonensis, and N. m. chuskaensis.

Shifts in the protist community associated with an anticyclonic gyre in the Alboran Sea (Mediterranean Sea)

The diversity of protists was researched in the Alboran Sea (SW Mediterranean Sea) by means of high-throughput sequencing technologies based on the amplification of the V9 region of 18S rRNA. Samples were collected at different depths in seven stations following an environmental gradient from a coastal upwelling zone to the core of an oligotrophic anticyclonic gyre (AG). Sampling was performed during summer, when the water column was stratified. The superphyla Alveolata, Stramenopila and Rhizaria accounted for 84% of the total operational taxonomic units (OTUs). The most diverse groups were Dinophyceae (21% of OTUs), Marine Alveolates-II (MALV-II; 20%), Ciliophora (9%) and MALV-I (6%). In terms of read abundance, the predominant groups were Dinophyceae (29%), Bacillariophyta (14%), MALV-II (11%) and Ciliophora (11%). Samples were clustered into three groups according to the sampling depth and position. The shallow community in coastal stations presented distinguishable patterns of diatoms and ciliates compared with AG stations. These results indicate that there was a strong horizontal coupling between phytoplankton and ciliate communities. Abundance of Radiolaria and Syndiniales increased with depth. Our analyses demonstrate that the stratification disruption produced by the AG caused shifts in the trophic ecology of the plankton assemblages inducing a transition from bottom-up to top-down control.

OTUs and ASVs Produce Comparable Taxonomic and Diversity from Shrimp Microbiota 16S Profiles Using Tailored Abundance Filters

The interplay between shrimp immune system, its environment, and microbiota contributes to the organism’s homeostasis and optimal production. The metagenomic composition is typically studied using 16S rDNA profiling by clustering amplicon sequences into operational taxonomic units (OTUs) and, more recently, amplicon sequence variants (ASVs). Establish the compatibility of the taxonomy, α, and β diversity described by both methods is necessary to compare past and future shrimp microbiota studies. Here, we used identical sequences to survey the V3 16S hypervariable-region using 97% and 99% OTUs and ASVs to assess the hepatopancreas and intestine microbiota of L. vannamei from two ponds under standardized rearing conditions. We found that applying filters to retain clusters >0.1% of the total abundance per sample enabled a consistent taxonomy comparison while preserving >94% of the total reads. The three sets turned comparable at the family level, whereas the 97% identity OTU set produced divergent genus and species profiles. Interestingly, the detection of organ and pond variations was robust to the clustering method’s choice, producing comparable α and β-diversity profiles. For comparisons on shrimp microbiota between past and future studies, we strongly recommend that ASVs be compared at the family level to 97% identity OTUs or use 99% identity OTUs, both using tailored frequency filters.

Effects of microbial organic fertilizers on Astragalus membranaceus growth and rhizosphere microbial community

Purpose

The application of excessive chemical fertilizers during the cultivation of Astragalus membranaceus leads to a decline in the quality of this medicinal plant as well as the soil’s sustainable productivity. In this study, we developed a special microbial organic fertilizer for A. membranaceus and investigated its effects on plant growth and rhizosphere microbial communities.

Methods

The root biomass and main active components of A. membranaceus in different growth stages were measured to assess the impacts of microbial organic manure on plant growth. Meanwhile, 16S rRNA and ITS1 amplicons were amplified and high-throughput sequencing was performed to detect the dynamic impacts of microbial organic manure on rhizosphere microbial communities.

Result

The results demonstrated that microbial organic manure significantly increased wet and dry weights of A. membranaceus seedlings and the accumulation of two effective components (flavonoids and saponin) in bacterial fertilizer treatment groups are significant higher than the control group. Research on rhizosphere microbial flora shows that the number and polymorphism of bacteria and fungi were decreased after the application of special fertilizer during the rapid growth period of plant and then gradually increased with seedling growth. The community structure of bacteria was regulated after the application of special fertilizer, and the beneficial bacteria for plant growth are enriched. Functional profiles prediction showed that significant shifts in metabolic functions impacting KEGG pathways of the microbial fertilizer treatment groups are related to metabolism and biosynthesis.

Conclusion

The results indicate that the microbial organic manure can improve A. membranaceus growth by providing appropriate nutrients and regulating the rhizosphere microbial community which has good potential in ecological cultivation of A. membranaceus.

Demography, genetic diversity and expansion load in the colonizing species Leontodon longirostris (Asteraceae) throughout its native range

Unravelling the evolutionary processes underlying range expansions is fundamental to understand the distribution of organisms, as well as to predict their future responses to environmental change. Predictions for range expansions include a loss of genetic diversity and an accumulation of deleterious alleles along the expansion axis, which can decrease fitness at the range-front (expansion load). In plants, empirical studies supporting expansion load are scarce, and its effects remain to be tested outside a few model species. Leontodon longirostris is a colonizing Asteraceae with a widespread distribution in the Western Mediterranean, providing a particularly interesting system to gain insight into the factors that can enhance or mitigate expansion load. In this study, we produced a first genome draft for the species, covering 418 Mbp (~53% of the genome). Although incomplete, this draft was suitable to design a targeted sequencing of ~1.5 Mbp in 238 L. longirostris plants from 21 populations distributed along putative colonization routes in the Iberian Peninsula. Inferred demographic history supports a range expansion from southern Iberia around 40,000 years ago, reaching northern Iberia around 25,000 years ago. The expansion was accompanied by a loss of genetic diversity and a significant increase in the proportion of putatively deleterious mutations. However, levels of expansion load in L. longirostris were smaller than those found in other plant species, which can be explained, at least partially, by its high dispersal ability, the self-incompatible mating system, and the fact that the expansion occurred along a strong environmental cline.

Fecal bacterial community of finishing beef steers fed ruminally protected and non-protected active dried yeast

Our previous study suggested that supplementation of high-grain diets with ruminally protected and non-protected active dried yeast (ADY) may potentially reduce manure pathogen excretion by feedlot cattle. We hypothesized that feeding ruminally protected ADY might change the fecal bacterial community of finishing cattle. The objective of this study was to investigate the effects of feeding ruminally protected and non-protected ADY to finishing beef steers on their fecal bacterial community. Fresh fecal samples were collected on day 56 from 50 steers fed one of five treatments: 1) control (no monensin, tylosin, or ADY), 2) antibiotics (ANT, 330 mg monensin + 110 mg tylosin·steer-1d-1), 3) ADY (1.5 g·steer-1d-1), 4) encapsulated ADY (EDY; 3 g·steer-1d-1), and 5) a mixture of ADY and EDY (MDY; 1.5 g ADY + 3 g EDY·steer-1d-1). Bacterial DNA was extracted from fecal samples and sequenced using a MiSeq high-throughput sequencing platform. A total number of 2,128,772 high-quality V4 16S rRNA sequences from 50 fecal samples were analyzed, and 1,424 operational taxonomic units (OTU) were detected based on 97% nucleotide sequence identity among reads, with 769 OTU shared across the five treatments. Alpha diversity indices, including species observed, Chao estimate, abundance-based coverage estimator, Shannon, Simpson, and coverage, did not differ among treatments, and principal coordinate analysis revealed a high similarity among treatments without independent distribution. Bacteroidetes and Firmicutes were dominant phyla in the fecal bacterial community for all treatments, with a tendency (P < 0.10) for greater relative abundance of Bacteroidetes but lesser Firmicutes with ANT, EDY, and MDY compared with control steers. Prevotella was the dominant genus in all treatments and steers supplemented with ANT, EDY, and MDY had greater (P < 0.05) relative abundance of Prevotella than control steers, but lesser (P < 0.03) relative abundance of Oscillospira. No differences between ADY and control were observed for the aforementioned variables. Fecal starch contents were not different among treatments, but the relative abundance of Bacteroidetes, as well as Prevotella at genera level, tended (P < 0.06) to be positively correlated to fecal starch content. We conclude that supplementing ruminally protected or non-protected ADY or ANT had no effect on diversity and richness of fecal bacteria of finishing beef cattle, whereas feeding protected ADY or ANT to finishing beef steers altered the dominant fecal bacteria at phylum and genus levels. Therefore, supplementation of ruminally protected ADY may potentially improve intestinal health by stimulating the relative abundance of Prevotella.

High-throughput sequencing reveals the diversity of TCR β chain CDR3 repertoire in patients with severe acne

Acne is a common chronic inflammatory skin disease, and the inflammation immune response runs through all stages of acne lesions. In this study, we use a combination of multiplex-PCR and high-throughput sequencing technologies to analyze T cell receptor β chain CDR3 (complementarity-determining region 3) in peripheral blood isolated from severe acne patients. Once compared with healthy controls, we propose to identify acne-relevant CDR3 peptides. Our results reveal that the diversity of T cell receptor β chain (TRB) CDR3 sequences in the peripheral blood of the severe acne vulgaris (SA) group differed from that of the control group. In addition, we find 10 TRB CDR3 sequences, amino acid sequences and V-J combinations with significantly different expressions between the SA group and the non-acne (NA) group (P < 0.0001). These findings may contribute to a better understanding of the role of immunity in the pathogenesis of acne and may serve as biomarkers for evaluating risk or prognosis of severe acne disease in future.

Doing More with Less: A Comparison of 16S Hypervariable Regions in Search of Defining the Shrimp Microbiota

The shrimp has become the most valuable traded marine product in the world, and its microbiota plays an essential role in its development and overall health status. Massive high-throughput sequencing techniques using several hypervariable regions of the 16S rRNA gene are broadly applied in shrimp microbiota studies. However, it is essential to consider that the use of different hypervariable regions can influence the obtained data and the interpretation of the results. The present study compares the shrimp microbiota structure and composition obtained by three types of amplicons: one spanning both the V3 and V4 hypervariable regions (V3V4), one for the V3 region only (V3), and one for the V4 region only (V4) using the same experimental and bioinformatics protocols. Twenty-four samples from hepatopancreas and intestine were sequenced and evaluated using the GreenGenes and silva reference databases for clustering and taxonomic classification. In general, the V3V4 regions resulted in higher richness and diversity, followed by V3 and V4. All three regions establish an apparent clustering effect that discriminates between the two analyzed organs and describe a higher richness for the intestine and a higher diversity for the hepatopancreas samples. Proteobacteria was the most abundant phyla overall, and Cyanobacteria was more common in the intestine, whereas Firmicutes and Actinobacteria were more prevalent in hepatopancreas samples. Also, the genus Vibrio was significantly abundant in the intestine, as well as Acinetobacter and Pseudomonas in the hepatopancreas suggesting these taxa as markers for their respective organs independently of the sequenced region. The use of a single hypervariable region such as V3 may be a low-cost alternative that enables an adequate description of the shrimp microbiota, allowing for the development of strategies to continually monitor the microbial communities and detect changes that could indicate susceptibility to pathogens under real aquaculture conditions while the use of the full V3V4 regions can contribute to a more in-depth characterization of the microbial composition.

Updates on Legume Genome Sequencing

Legumes (Fabaceae) are agronomically and economically one of the most important crops. Because legumes serve as a source of food, feed, and industrial materials, many studies in the field of legume genomics, including genome sequencing, have been conducted over the last decade. Here, we update the progress in genome sequencing of legume crops, including soybean (Glycine max [L.] Merr.), mung bean (V. radiata var. radiata), adzuki bean (V. angularis var. angularis), common bean (Phaseolus vulgaris L.), pigeon pea (Cajanus cajan), chickpea (Cicer arietinum), and peanut (Arachis hypogaea). Since the publication of the first reference genome sequence of each species, many accessions have been resequenced to study genetic diversity, speciation, and polyploidization in the legume lineage.

Exploring functional core bacteria in fermentation of a traditional Chinese food, Aspergillus-type douchi

Douchi is a type of traditional Chinese flavoring food that has been used for thousands of years and is produced by multispecies solid-state fermentation. However, the correlation between the flavor, the microbiota, and the functional core microbiota in Aspergillus-type douchi fermentation remains unclear. In this study, Illumina MiSeq sequencing and chromatography were used to investigate the bacterial community and flavor components in Aspergillus-type douchi fermentation. The dominant phyla were Firmicutes, Proteobacteria, and Actinobacteria, and the dominant genera were Weissella, Bacillus, Anaerosalibacter, Lactobacillus, Staphylococcus, and Enterococcus. A total of 58 flavor components were detected during fermentation, including two alcohols, 14 esters, five pyrazines, three alkanes, four aldehydes, three phenols, six acids, and five other compounds. Bidirectional orthogonal partial least square modeling showed that Corynebacterium_1, Lactococcus, Atopostipes, Peptostreptococcus, norank_o__AKYG1722, Truepera, Gulosibacter, norank_f__Actinomycetaceae, and unclassified_f__Rhodobacteraceae are the functional core microbiota responsible for the formation of the flavor components during douchi fermentation. This is the first study to investigate the functional core microbiota in douchi fermentation using Illumina MiSeq sequencing and chromatographic techniques. Our findings extend our understanding of the relationships between flavor, the microbiota, and the functional core microbiota during Aspergillus-type douchi fermentation.

Metatranscriptomic and metabolite profiling reveals vertical heterogeneity within a Zygnema green algal mat from Svalbard (High Arctic)

Within streptophyte green algae Zygnematophyceae are the sister group to the land plants that inherited several traits conferring stress protection. Zygnema sp., a mat-forming alga thriving in extreme habitats, was collected from a field site in Svalbard, where the bottom layers are protected by the top layers. The two layers were investigated by a metatranscriptomic approach and GC-MS-based metabolite profiling. In the top layer, 6569 genes were significantly upregulated and 149 were downregulated. Upregulated genes coded for components of the photosynthetic apparatus, chlorophyll synthesis, early light-inducible proteins, cell wall and carbohydrate metabolism, including starch-degrading enzymes. An increase in maltose in the top layer and degraded starch grains at the ultrastructural levels corroborated these findings. Genes involved in amino acid, redox metabolism and DNA repair were upregulated. A total of 29 differentially accumulated metabolites (out of 173 identified ones) confirmed higher metabolic turnover in the top layer. For several of these metabolites, differential accumulation matched the transcriptional changes of enzymes involved in associated pathways. In summary, the findings support the hypothesis that in a Zygnema mat the top layer shields the bottom layers from abiotic stress factors such as excessive irradiation.

Cold Acclimation Improves the Desiccation Stress Resilience of Polar Strains of Klebsormidium (Streptophyta)

Biological soil crusts (BSCs) are complex communities of autotrophic, heterotrophic, and saprotrophic (micro)organisms. In the polar regions, these biocrust communities have essential ecological functions such as primary production, nitrogen fixation, and ecosystem engineering while coping with extreme environmental conditions (temperature, desiccation, and irradiation). The microalga Klebsormidium is commonly found in BSCs all across the globe. The ecophysiological resilience of various Klebsormidium species to desiccation and other stresses has been studied intensively. Here we present the results of transcriptomic analyses of two different Klebsormidium species, K. dissectum and K. flaccidum, isolated from Antarctic and Arctic BSCs. We performed desiccation stress experiments at two different temperatures mimicking fluctuations associated with global change. Cultures grown on agar plates were desiccated on membrane filters at 10% relative air humidity until the photosynthetic activity as reflected in the effective quantum yield of photosystem II [Y(II)] ceased. For both species, the response to dehydration was much faster at the higher temperature. At the transcriptome level both species responded more strongly to the desiccation stress at the higher temperature suggesting that adaptation to cold conditions enhanced the resilience of both algae to desiccation stress. Interestingly, the two different species responded differently to the applied desiccation stress with respect to the number as well as function of genes showing differential gene expression. The portion of differentially expressed genes shared between both taxa was surprisingly low indicating that both Klebsormidium species adapted independently to the harsh conditions of Antarctica and the Arctic, respectively. Overall, our results indicate that environmental acclimation has a great impact on gene expression and the response to desiccation stress in Klebsormidium.

A Microbiota-Derived Bacteriocin Targets the Host to Confer Diarrhea Resistance in Early-Weaned Piglets

Alternatives to antibiotics for preventing diarrhea in early-weaned farm animals are sorely needed. CM piglets (a native Chinese breed) are more resistant to early-weaning stress-induced diarrhea than the commercial crossbred LY piglets. Transferring fecal microbiota, but not saline, from healthy CM into LY piglets by oral administration prior to early weaning conferred diarrhea resistance. By comparing the relative abundance of intestinal microbiota in saline and microbiota transferred LY piglets, we identified and validated Lactobacillus gasseri LA39 and Lactobacillus frumenti as two bacterial species that mediate diarrhea resistance. Diarrhea resistance depended on the bacterial secretory circular peptide gassericin A, a bacteriocin. The binding of gassericin A to Keratin 19 (KRT19) on the plasma membrane of intestinal epithelial cells was essential for enhancement of fluid absorption and decreased secretion. These findings suggest the use of L. gasseri LA39 and L. frumenti as antibiotic alternatives for preventing diarrhea in mammals.

Viability and Composition Validation of Commercial Probiotic Products by Selective Culturing Combined with Next-Generation Sequencing

The consumption of dietary supplements to treat health complications or to improve overall health conditions has become a globally increasing trend that leads to the development of a large number of health-related novel products and expands the associated manufacturing industries around the world. In the current study, we applied selective culturing combined with next-generation sequencing to examine the microbial viability in terms of its culturability on culture medium, composition, and possible contamination in the selected 17 commercial probiotic products sold in the mainland China market. Additionally, the relative abundance of each individual bacterial content was also evaluated by using the generated sequencing reads. The tested probiotic product samples were subjected to Illumina HiSeq-2000 sequencing platform and thoroughly analyzed by the in-house developed bioinformatics pipeline. The comprehensive culturing and sequencing analysis revealed both viability and composition inaccuracy among the several tested probiotic products, however, no contaminant was identified during the analysis. Among the total, five probiotic products (29.41%) were found with an inaccurate or lower colony-forming unit (CFU) counts on culture media while four probiotic products (23.52%) have inaccurately labeled classification. This study provides an ideal qualitative and quantitative assessment approach, which can be used as a diagnostic tool for the accurate assessment of commercial probiotic supplements.

Effects of a probiotic drink containing Lactobacillus casei strain Shirota on dental plaque microbiota

Objective

To investigate the effects of a probiotic drink containing Lactobacillus casei strain Shirota on dental plaque microbiota.

Methods

Caries-free young adults were administered a probiotic drink containing L. casei Shirota (Yakult) three times on day 1, then once daily for 28 days. Dental plaque samples were collected and analysed by 16S rRNA sequencing before (day 1), during (day 2), and one day following intervention (day 30).

Results

Out of samples from 10 included participants, 256 814 sequences passed through quality control, clustered into 170–234 different ‘species-level’ phylotypes. Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria and Proteobacteria were the most abundant phyla. Dental plaque microbiota composition significantly changed at different taxonomic levels following probiotic consumption. At genus level, the relative abundance of Veillonella and Kingella increased significantly following intervention, while that of Leptotrichia reduced significantly during intervention, but recovered to baseline level at day 30. Alpha diversity and overall structure of the dental plaque microbiota was not significantly impacted by the probiotic.

Conclusion

Yakult intake changed the abundance of some bacteria related to caries, suggesting that the change of composition may be beneficial to oral health, while the overall microbiota structure remained unaffected.

Genomic sequencing of Troides aeacus nucleopolyhedrovirus (TraeNPV) from golden birdwing larvae (Troides aeacus formosanus) to reveal defective Autographa californica NPV genomic features

Background

The golden birdwing butterfly (Troides aeacus formosanus) is a rarely observed species in Taiwan. Recently, a typical symptom of nuclear polyhedrosis was found in reared T. aeacus larvae. From the previous Kimura-2 parameter (K-2-P) analysis based on the nucleotide sequence of three genes in this isolate, polh, lef-8 and lef-9, the underlying virus did not belong to any known nucleopolyhedrovirus (NPV) species. Therefore, this NPV was provisionally named “TraeNPV”. To understand this NPV, the nucleotide sequence of the whole TraeNPV genome was determined using next-generation sequencing (NGS) technology.

Results

The genome of TraeNPV is 125,477 bp in length with 144 putative open reading frames (ORFs) and its GC content is 40.45%. A phylogenetic analysis based on the 37 baculoviral core genes suggested that TraeNPV is a Group I NPV that is closely related to Autographa californica nucleopolyhedrovirus (AcMNPV). A genome-wide analysis showed that TraeNPV has some different features in its genome compared with other NPVs. Two novel ORFs (Ta75 and Ta139), three truncated ORFs (pcna, he65 and bro) and one duplicated ORF (38.7 K) were found in the TraeNPV genome; moreover, there are fewer homologous regions (hrs) than there are in AcMNPV, which shares eight hrs within the TraeNPV genome. TraeNPV shares similar genomic features with AcMNPV, including the gene content, gene arrangement and gene/genome identity, but TraeNPV lacks 15 homologous ORFs from AcMNPV in its genome, such as ctx, host cell-specific factor 1 (hcf-1), PNK/PNL, vp15, and apsup, which are involved in the auxiliary functions of alphabaculoviruses.

Conclusions

Based on these data, TraeNPV would be clarified as a new NPV species with defective AcMNPV genomic features. The precise relationship between TraeNPV and other closely related NPV species were further investigated. This report could provide comprehensive information on TraeNPV for evolutionary insights into butterfly-infected NPV.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5713-2) contains supplementary material, which is available to authorized users.

Impact of Introduction of Arbuscular Mycorrhizal Fungi on the Root Microbial Community in Agricultural Fields

Arbuscular mycorrhizal (AM) fungi are important members of the root microbiome and may be used as biofertilizers for sustainable agriculture. To elucidate the impact of AM fungal inoculation on indigenous root microbial communities, we used high-throughput sequencing and an analytical pipeline providing fixed operational taxonomic units (OTUs) as an output to investigate the bacterial and fungal communities of roots treated with a commercial AM fungal inoculum in six agricultural fields. AM fungal inoculation significantly influenced the root microbial community structure in all fields. Inoculation changed the abundance of indigenous AM fungi and other fungal members in a field-dependent manner. Inoculation consistently enriched several bacterial OTUs by changing the abundance of indigenous bacteria and introducing new bacteria. Some inoculum-associated bacteria closely interacted with the introduced AM fungi, some of which belonged to the genera Burkholderia, Cellulomonas, Microbacterium, Sphingomonas, and Streptomyces and may be candidate mycorrhizospheric bacteria that contribute to the establishment and/or function of the introduced AM fungi. Inoculated AM fungi also co-occurred with several indigenous bacteria with putative beneficial traits, suggesting that inoculated AM fungi may recruit specific taxa to confer better plant performance. The bacterial families Methylobacteriaceae, Acetobacteraceae, Armatimonadaceae, and Alicyclobacillaceae were consistently reduced by the inoculation, possibly due to changes in the host plant status caused by the inoculum. To the best of our knowledge, this is the first large-scale study to investigate interactions between AM fungal inoculation and indigenous root microbial communities in agricultural fields.

Species-level bacterial community profiling of the healthy sinonasal microbiome using Pacific Biosciences sequencing of full-length 16S rRNA genes

BACKGROUND

Pan-bacterial 16S rRNA microbiome surveys performed with massively parallel DNA sequencing technologies have transformed community microbiological studies. Current 16S profiling methods, however, fail to provide sufficient taxonomic resolution and accuracy to adequately perform species-level associative studies for specific conditions. This is due to the amplification and sequencing of only short 16S rRNA gene regions, typically providing for only family- or genus-level taxonomy. Moreover, sequencing errors often inflate the number of taxa present. Pacific Biosciences' (PacBio's) long-read technology in particular suffers from high error rates per base. Herein, we present a microbiome analysis pipeline that takes advantage of PacBio circular consensus sequencing (CCS) technology to sequence and error correct full-length bacterial 16S rRNA genes, which provides high-fidelity species-level microbiome data.

RESULTS

Analysis of a mock community with 20 bacterial species demonstrated 100% specificity and sensitivity with regard to taxonomic classification. Examination of a 250-plus species mock community demonstrated correct species-level classification of > 90% of taxa, and relative abundances were accurately captured. The majority of the remaining taxa were demonstrated to be multiply, incorrectly, or incompletely classified. Using this methodology, we examined the microgeographic variation present among the microbiomes of six sinonasal sites, by both swab and biopsy, from the anterior nasal cavity to the sphenoid sinus from 12 subjects undergoing trans-sphenoidal hypophysectomy. We found greater variation among subjects than among sites within a subject, although significant within-individual differences were also observed. Propiniobacterium acnes (recently renamed Cutibacterium acnes) was the predominant species throughout, but was found at distinct relative abundances by site.

CONCLUSIONS

Our microbial composition analysis pipeline for single-molecule real-time 16S rRNA gene sequencing (MCSMRT, https://github.com/jpearl01/mcsmrt ) overcomes deficits of standard marker gene-based microbiome analyses by using CCS of entire 16S rRNA genes to provide increased taxonomic and phylogenetic resolution. Extensions of this approach to other marker genes could help refine taxonomic assignments of microbial species and improve reference databases, as well as strengthen the specificity of associations between microbial communities and dysbiotic states.

Genomic analysis of MHC-based mate choice in the monogamous California mouse

Variation at Major Histocompatibility Complex (MHC) genes is thought to be an important mechanism underlying mate choice in vertebrates, with individuals typically predicted to prefer MHC-dissimilar reproductive partners. However, analyses based on individual MHC loci have generated contradictory results regarding the role of these genes in mate-choice decisions. To provide a more comprehensive assessment of relationships between MHC variation and mating behavior, we used an exome capture strategy to characterize variability at 13 MHC loci, 312 innate immune system genes, and 1044 nonimmune genes in 25 obligate monogamous pairs of California mice (Peromyscus californicus) from 2 free-living populations of this species in Monterey County, California. Pairwise genotypic comparisons and analyses of SNP-based allelic differences failed to detect disassortative mating based on MHC variability; reproductive partners were not more dissimilar than randomly generated male-female pairs at MHC, innate or nonimmune loci. Within populations, individuals tended to be more closely related at MHC genes than at innate or nonimmune genes. Consistent with the functional role of immunogenes, the 2 study populations were highly differentiated at MHC and innate genes but not at nonimmune loci. Collectively, our results suggest that MHC genetic variation in California mice reflects local differences in pathogen exposure rather than disassortative mating based on variability at MHC Class I and II genes.

Vertical Distribution of Microbial Eukaryotes From Surface to the Hadal Zone of the Mariana Trench

Marine microbial eukaryotes are ubiquitous, comprised of phylogenetically diverse groups and play key roles in microbial food webs and global biogeochemical cycling. However, their vertical distribution in the deep sea has received little attention. In this study, we investigated the composition and diversity of the eukaryotes of both 0.2–3 μm and >3 μm size fractions from the surface to the hadal zone (8727 m) of the Mariana Trench using Illumina MiSeq sequencing for the 18S rDNA. The microbial eukaryotic community structure differed substantially across size fractions and depths. Operational taxonomic unit (OTU) richness in the >3 μm fraction was higher than that in the 0.2–3 μm fraction at the same depth. For the 0.2–3 μm fraction, sequences of Retaria (Rhizaria) were most abundant in the surface water (53.5%). Chrysophyceae (Stramenopiles) sequences dominated mostly in the samples from water depths below 1795 m. For the >3 μm fraction, sequences of Dinophyceae (Alveolata) were most abundant in surface waters (49.3%) and remained a significant proportion of total sequences at greater depths (9.8%, on average). Retaria sequences were abundant in samples of depths ≥1000 m. Amoebozoa and Apusozoa sequences were enriched in the hadal sample, comprising 38 and 20.4% of total sequences, respectively. Fungi (Opisthokonta) sequences were most abundant at 1759 m in both size fractions. Strong positive associations were found between Syndiniales (mainly MALV-I and MALV-II) and Retaria while negative associations were shown between MALV-II and Fungi in a co-occurrence analysis. This study compared the community structure of microbial eukaryotes in different zones in the deep sea and identified a distinct hadal community in the larger size fraction, suggesting the uniqueness of the eukaryotes in the biosphere in the Mariana Trench.

Yeast Probiotics Shape the Gut Microbiome and Improve the Health of Early-Weaned Piglets

Weaning is one of the most stressful challenges in the pig's life, which contributes to dysfunctions of intestinal and immune system, disrupts the gut microbial ecosystem, and therefore compromises the growth performance and health of piglets. To mitigate the negative impact of the stress on early-weaned piglets, effective measures are needed to promote gut health. Toward this end, we tamed a Saccharomyces cerevisiae strain and developed a probiotic Duan-Nai-An, which is a yeast culture of the tamed S. cerevisiae on egg white. In this study, we tested the performance of Duan-Nai-An on growth and health of early-weaned piglets and analyzed its impact on fecal microbiota. The results showed that Duan-Nai-An significantly improved weight gain and feed intake, and reduced diarrhea and death of early-weaned piglets. Analysis of the gut microbiota showed that the bacterial community was shaped by Duan-Nai-An and maintained as a relatively stable structure, represented by a higher core OTU number and lower unweighted UniFrac distances across the early weaned period. However, fungal community was not significantly shaped by the yeast probiotics. Notably, 13 bacterial genera were found to be associated with Duan-Nai-An feeding, including Enterococcus, Succinivibrio, Ruminococcus, Sharpea, Desulfovibrio, RFN20, Sphaerochaeta, Peptococcus, Anaeroplasma, and four other undefined genera. These findings suggest that Duan-Nai-An has the potential to be used as a feed supplement in swine production.

PIPEBAR and OverlapPER: tools for a fast and accurate DNA barcoding analysis and paired-end assembly

Background

Taxonomic identification of plants and insects is a hard process that demands expert taxonomists and time, and it’s often difficult to distinguish on morphology only. DNA barcodes allow a rapid species discovery and identification and have been widely used for taxonomic identification by targeting known gene regions that permit to discriminate these species. DNA barcode sequence analysis is usually carried out with processes and tools that still demand a high interaction with the user or researcher. To reduce at most such interaction, we proposed PIPEBAR, a pipeline for DNA chromatograms analysis of Sanger platform sequencing, ensuring high quality consensus sequences along with efficient running time. We also proposed a paired-end reads assembly tool, OverlapPER, which is used in sequence or independently of PIPEBAR.

Results

PIPEBAR is a command line tool to automatize the processing of large number of trace files. It is accurate as the proprietary Geneious tool and faster than most popular software for barcoding analysis. It is 7 times faster than Geneious and 14 times faster than SeqTrace for processing hundreds of barcoding sequences. OverlapPER is a novel tool for overlapping paired-end reads accurately that accepts both substitution and indel errors and returns both overlapped and non-overlapped regions between a pair of reads. OverlapPER obtained the best results compared to currently used tools when merging 1,000,000 simulated paired-end reads.

Conclusions

PIPEBAR and OverlapPER run on most operating systems and are freely available, along with supporting code and documentation, at https://sourceforge.net/projects/PIPEBAR/ and https://sourceforge.net/projects/overlapper-reads/.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2307-y) contains supplementary material, which is available to authorized users.

The roles of the TLR/NF‑κB signaling pathway in the mutual interactions between the lung and the large intestine

The 'exterior-interior relationship between the lung and the large intestine' is a classical basic theory in Traditional Chinese Medicine. The present study aimed to investigate the roles of the toll like receptor/nuclear factor‑κB (TLR/NF‑κB) signaling pathway in the mutual interactions between the lung and the large intestine. A rat model of allergic asthma complicated with intestinal flora disorder was established by oral administration of Candida albicans and intraperitoneal injection with ovalbumin. The number of inflammatory cells and expression levels immunoglobulin (Ig)E, secretory IgA, interleukin (IL)‑4 and interferon‑γ in serum and bronchoalveolar lavage fluid were subsequently measured. Bacterial colonies and expression of 16S ribosomal DNA were studied in feces samples and pathological alterations of lung tissues were identified. Furthermore, the expression levels of genes associated with the TLR/NF‑κB signaling pathway in the lung and intestinal tissues were determined by reverse transcription‑quantitative polymerase chain reaction. The results of the present study indicated that, in the rat model of allergic asthma complicated with intestinal flora disorder, the expression levels of IL‑4 and IgE, and the numbers of inflammatory cells and C. albicans increased, and marked inflammatory cell infiltration was observed in lung tissues, suggesting that the animal model was successfully established. Furthermore, the present results revealed the mRNA expression levels of genes associated with the TLR/NF‑κB signaling (including myeloid differentiation primary response 88, TNF receptor associated factor 6 and β‑arrestin) were upregulated in both of the lung and intestinal tissues of the model group rats. Collectively, the results demonstrated that the TLR/NF‑κB signaling may serve roles in the mutual interactions between the lung and the large intestine, and TLR and NF‑κB may be potential targets for the treatment of lung diseases complicated with intestinal disorders.

Recent advances in sequence assembly: principles and applications

The application of advanced sequencing technologies and the rapid growth of various sequence data have led to increasing interest in DNA sequence assembly. However, repeats and polymorphism occur frequently in genomes, and each of these has different impacts on assembly. Further, many new applications for sequencing, such as metagenomics regarding multiple species, have emerged in recent years. These not only give rise to higher complexity but also prevent short-read assembly in an efficient way. This article reviews the theoretical foundations that underlie current mapping-based assembly and de novo-based assembly, and highlights the key issues and feasible solutions that need to be considered. It focuses on how individual processes, such as optimal k-mer determination and error correction in assembly, rely on intelligent strategies or high-performance computation. We also survey primary algorithms/software and offer a discussion on the emerging challenges in assembly.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

MERIT reveals the impact of genomic context on sequencing error rate in ultra-deep applications

BACKGROUND

Rapid progress in high-throughput sequencing (HTS) and the development of novel library preparation methods have improved the sensitivity of detecting mutations in heterogeneous samples, specifically in high-depth (> 500×) clinical applications. However, HTS methods are bounded by their technical and theoretical limitations and sequencing errors cannot be completely eliminated. Comprehensive quantification of the background noise can highlight both the efficiency and the limitations of any HTS methodology, and help differentiate true mutations at low abundance from artifacts.

RESULTS

We introduce MERIT (Mutation Error Rate Inference Toolkit), designed for in-depth quantification of erroneous substitutions and small insertions and deletions. MERIT incorporates an all-inclusive variant caller and considers genomic context, including the nucleotides immediately at 5 'and 3 ', thereby establishing error rates for 96 possible substitutions as well as four single-base and 16 double-base indels. We applied MERIT to ultra-deep sequencing data (1,300,000 ×) obtained from the amplification of multiple clinically relevant loci, and showed a significant relationship between error rates and genomic contexts. In addition to observing significant difference between transversion and transition rates, we identified variations of more than 100-fold within each error type at high sequencing depths. For instance, T >G transversions in trinucleotide GTCs occurred 133.5 ± 65.9 more often than those in ATAs. Similarly, C >T transitions in GCGs were observed at 73.8 ± 10.5 higher rate than those in TCTs. We also devised an in silico approach to determine the optimal sequencing depth, where errors occur at rates similar to those of expected true mutations. Our analyses showed that increasing sequencing depth might improve sensitivity for detecting some mutations based on their genomic context. For example, T >G rate of error in GTCs did not change when sequenced beyond 10,000 ×; in contrast, T >G rate in TTAs consistently improved even at above 500,000 ×.

CONCLUSIONS

Our results demonstrate significant variation in nucleotide misincorporation rates, and suggest that genomic context should be considered for comprehensive profiling of specimen-specific and sequencing artifacts in high-depth assays. This data provide strong evidence against assigning a single allele frequency threshold to call mutations, for it can result in substantial false positive as well as false negative variants, with important clinical consequences.

Dissection of niche competition between introduced and indigenous arbuscular mycorrhizal fungi with respect to soybean yield responses

Arbuscular mycorrhizal (AM) fungi associate with most land plants and deliver phosphorus to the host. Identification of biotic/abiotic factors that determine crop responses to AM fungal inoculation is an essential step for successful application of the fungi in sustainable agriculture. We conducted three field trials on soybean with a commercial inoculum and developed a new molecular tool to dissect interactions between the inoculum and indigenous fungi on the MiSeq sequencing platform. Regression analysis indicated that sequence read abundance of the inoculum fungus was the most significant factor that determined soybean yield responses to the inoculation, suggesting that dominance of the inoculum fungus is a necessary condition for positive yield responses. Agricultural practices (fallow/cropping in the previous year) greatly affected the colonization levels (i.e. read abundances) of the inoculum fungus via altering the propagule density of indigenous AM fungi. Analysis of niche competition revealed that the inoculum fungus competed mainly with the indigenous fungi that are commonly distributed in the trial sites, probably because their life-history strategy is the same as that of the inoculum fungus. In conclusion, we provide a new framework for evaluating the significance of environmental factors towards successful application of AM fungi in agriculture.

Detection of Tick-Borne Pathogens in Lambs Undergoing Prophylactic Treatment Against Ticks on Two Swedish Farms

Tick-borne pathogens (TBPs), especially Anaplasma phagocytophilum, cause disease in grazing livestock. Tick prophylaxis is, therefore, a routine practice in sheep flocks in Sweden, especially in central, southern, and coastal areas of the country where ixodid ticks (Ixodes ricinus and Haemaphysalis punctata) are present. In the present study, the status of infection by A. phagocytophilum and other TBPs in lambs treated with tick prophylaxis has been assessed serologically and with polymerase chain reaction (PCR). Blood samples (n = 78) from lambs (n = 20) subjected to regular tick prophylactic treatment (flumethrin, Bayticol^®) at two sites in different regions in Sweden (Östergötland, Gotland) were collected on four occasions from May until July 2013. The severity of clinical signs in Anaplasma-infected animals is known to differ between these two regions. In total, 20% of blood samples were PCR-positive for A. phagocytophilum. Serological analyses showed that 33% of all collected samples were positive for A. phagocytophilum, while 2.5% were positive for Borrelia burgdorferi s.l. and 13% for tick-borne encephalitis virus (TBEV). Percentages of lambs positive were 75 and 45% for A. phagocytophilum antibodies and DNA, respectively, while 10 and 45% were serologically positive for B. burgdorferi s.l. and TBEV, respectively. Sequencing of partial 16S rRNA genes from Anaplasma PCR positive samples revealed presence of A. phagocytophilum in all animals in Östergötland, while sequences consistent with A. phagocytophilum as well as A. capra and A. bovis were found on the island of Gotland. This is the first report of the occurrence of the latter two species in Sweden.

From Short Reads to Chromosome-Scale Genome Assemblies

A high-quality, annotated genome assembly is the foundation for many downstream studies. However, obtaining such an assembly is a complex, reiterative process that requires the assimilation of high-quality data and combines different approaches and data types. While some software packages incorporating multiple steps of genome assembly are commercially available, they may not be flexible enough to be routinely applied to all organisms, particularly to nonmodel species such as pathogenic oomycetes and fungi. If researchers understand and apply the most appropriate, currently available tools for each step, it is possible to customize parameters and optimize results for their organism of study. Based on our experience of de novo assembly and annotation of several oomycete species, this chapter provides a modular workflow from processing of raw reads, to initial assembly generation, through optimization, chromosome-scale scaffolding and annotation, outlining input and output data as well as examples and alternative software used for each step. The accompanying Notes provide background information for each step as well as alternative options. The final result of this workflow could be an annotated, high-quality, validated, chromosome-scale assembly or a draft assembly of sufficient quality to meet specific needs of a project.

Enhanced Desiccation Tolerance in Mature Cultures of the Streptophytic Green Alga Zygnema circumcarinatum Revealed by Transcriptomics

Desiccation tolerance is commonly regarded as one of the key features for the colonization of terrestrial habitats by green algae and the evolution of land plants. Extensive studies, focused mostly on physiology, have been carried out assessing the desiccation tolerance and resilience of the streptophytic genera Klebsormidium and Zygnema. Here we present transcriptomic analyses of Zygnema circumcarinatum exposed to desiccation stress. Cultures of Z. circumcarinatum grown in liquid medium or on agar plates were desiccated at ∼86% relative air humidity until the effective quantum yield of PSII [Y(II)] ceased. In general, the response to dehydration was much more pronounced in Z. circumcarinatum cultured in liquid medium for 1 month compared with filaments grown on agar plates for 7 and 12 months. Culture on solid medium enables the alga to acclimate to dehydration much better and an increase in desiccation tolerance was clearly correlated to increased culture age. Moreover, gene expression analysis revealed that photosynthesis was strongly repressed upon desiccation treatment in the liquid culture while only minor effects were detected in filaments cultured on agar plates for 7 months. Otherwise, both samples showed induction of stress protection mechanisms such as reactive oxygen species scavenging (early light-induced proteins, glutathione metabolism) and DNA repair as well as the expression of chaperones and aquaporins. Additionally, Z. circumcarinatum cultured in liquid medium upregulated sucrose-synthesizing enzymes and strongly induced membrane modifications in response to desiccation stress. These results corroborate the previously described hardening and associated desiccation tolerance in Zygnema in response to seasonal fluctuations in water availability.

Filling reference gaps via assembling DNA barcodes using high-throughput sequencing-moving toward barcoding the world

Over the past decade, biodiversity researchers have dedicated tremendous efforts to constructing DNA reference barcodes for rapid species registration and identification. Although analytical cost for standard DNA barcoding has been significantly reduced since early 2000, further dramatic reduction in barcoding costs is unlikely because Sanger sequencing is approaching its limits in throughput and chemistry cost. Constraints in barcoding cost not only led to unbalanced barcoding efforts around the globe, but also prevented high-throughput sequencing (HTS)-based taxonomic identification from applying binomial species names, which provide crucial linkages to biological knowledge. We developed an Illumina-based pipeline, HIFI-Barcode, to produce full-length Cytochrome c oxidase subunit I (COI) barcodes from pooled polymerase chain reaction amplicons generated by individual specimens. The new pipeline generated accurate barcode sequences that were comparable to Sanger standards, even for different haplotypes of the same species that were only a few nucleotides different from each other. Additionally, the new pipeline was much more sensitive in recovering amplicons at low quantity. The HIFI-Barcode pipeline successfully recovered barcodes from more than 78% of the polymerase chain reactions that didn't show clear bands on the electrophoresis gel. Moreover, sequencing results based on the single molecular sequencing platform Pacbio confirmed the accuracy of the HIFI-Barcode results. Altogether, the new pipeline can provide an improved solution to produce full-length reference barcodes at about one-tenth of the current cost, enabling construction of comprehensive barcode libraries for local fauna, leading to a feasible direction for DNA barcoding global biomes.

The performance and archaeal community shifts in a modified anaerobic baffled reactor treating sweet potato starch wastewater at ambient temperatures

A conventional anaerobic baffled reactors (ABRs) treating high strength sweet potato starch wastewater at ambient temperatures resulted in acidification and bad performances. After modification, the acidification was remitted and COD removal efficiencies reached 92.73% at high temperatures and were maintained at 71.19% at low temperatures. Moreover, as much as 1.014 ± 0.056 L CH₄/L/d were collected at Stage III. The q-PCR results revealed that the largest methanogen populations emerged at Stage III as well, which was 5.29 × 10⁸ mcrA copies per milliliter sludge. A comparable shift in the archaeal community structure at different stages and acetoclastic methanogens Methanosaeta predominated the archaeal community in every compartment in Stages I (63.73%) and II (48.63%). Finally, the net energy gains analysis at mesophilic, thermophilic, and ambient temperature revealed that modified ABR at ambient temperature was not only economical but also profitable and could generated 3.68 KJ energy per gram COD removed.