Evaluating the efficacy of the new Ion PGM Hi-Q Sequencing Kit applied to bacterial genomes

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Comparative polyphasic characterization of Weissella strains isolated from beaked whale and rainbow trout (Oncorhynchus mykiss): confirmation of Weissella ceti sp. nov. and description of the novel Weissella tructae sp. nov. isolated from farmed rainbow trout

The weissellosis agent bacterium (WS08T = CBMAI 2730) was isolated from diseased rainbow trout (Oncorhynchus mykiss) in Brazil. The whole genome sequence of this strain was compared with the Mexican W-1 strain, also isolated from diseased rainbow trout, and with the Weissella ceti type strain CECT 7719 T (= 1119-1A-09 T = CCUG 59653 T), recovered from the beaked whale. Digital DNA-DNA hybridization pairwise analyses scored 98.7% between the Mexican W-1 and Brazilian WS08T but just 24.4% for both fish isolates compared to the W. ceti type strain CECT 7719 T. The 16S rRNA gene sequence comparisons with isolates of W. ceti, available at GenBank, were conducted. All rainbow trout-pathogenic isolates grouped close (97% bootstrap confirmation), but when this group was compared to the W. ceti type strain CECT 7719 T the similarity varied from 78.9 to 79.1%. Phenotypic assays were also conducted, and the W. ceti type strain diverged from WS08T and W-1 in the hydrolysis of aesculin, D-mannose, and potassium gluconate and in the hydrolysis of hippurate. Moreover, WS08T and W-1 showed weak growth at 5 °C whereas no growth was observed for W. ceti CECT 7719 T. The major fatty acids (> 10% total fatty acids) presented by WS08T and W-1 were summed feature 8 (C18:1 ω7c/C18:1 ω6c), summed feature 3 (C16:1 ω6c/C16:1ω7c), and C16:0. The results of phylogenetic and phenotypic analyses clearly differentiated the W. ceti CECT 7719 T type strain from the assessed pathogenic strains obtained from rainbow trout. Therefore, Weissella strains isolated from rainbow trout, here represented by strain WS08T (= CBMAI 2730), should be known as members of a novel species for which the name Weissella tructae sp. nov. is proposed.

A Hybrid Microfluidic Electronic Sensing Platform for Life Science Applications

This paper presents a novel hybrid microfluidic electronic sensing platform, featuring an electronic sensor incorporated with a microfluidic structure for life science applications. This sensor with a large sensing area of 0.7 mm² is implemented through a foundry process called Open-Gate Junction FET (OG-JFET). The proposed OG-JFET sensor with a back gate enables the charge by directly introducing the biological and chemical samples on the top of the device. This paper puts forward the design and implementation of a PDMS microfluidic structure integrated with an OG-JFET chip to direct the samples toward the sensing site. At the same time, the sensor’s gain is controlled with a back gate electrical voltage. Herein, we demonstrate and discuss the functionality and applicability of the proposed sensing platform using a chemical solution with different pH values. Additionally, we introduce a mathematical model to describe the charge sensitivity of the OG-JFET sensor. Based on the results, the maximum value of transconductance gain of the sensor is ~1 mA/V at Vgs = 0, which is decreased to ~0.42 mA/V at Vgs = 1, all in Vds = 5. Furthermore, the variation of the back-gate voltage from 1.0 V to 0.0 V increases the sensitivity from ~40 mV/pH to ~55 mV/pH. As per the experimental and simulation results and discussions in this paper, the proposed hybrid microfluidic OG-JFET sensor is a reliable and high-precision measurement platform for various life science and industrial applications.

Prokaryotic Diversity and Metabolically Active Communities in Brines from Two Perennially Ice-Covered Antarctic Lakes

The genomic diversity of bacteria and archaea in brines (BC1, BC2, and BC3) from two adjacent and perennially frozen Antarctic lakes (L16 and L-2) in the Boulder Clay (BC) area was investigated together with the metabolically active fraction of both communities, by analyzing the bulk rRNA as a general marker of metabolic activity. Although similar bacterial and archaeal assemblages were observed at phylum level, differences were encountered when considering the distribution in species. Overall, the total bacterial communities were dominated by Bacteroidetes. A massive occurrence of flavobacterial sequences was observed within the metabolically active bacterial communities of the BC1 brine, whereas the active fractions in BC2 and BC3 strongly differed from the bulk communities being dominated by Betaproteobacteria (mainly Hydrogenophaga members). The BC lakes also hosted sequences of the most thermally tolerant archaea, also related to well-known hyperthermophiles. Interestingly, RNA sequences of the hyperthermophilic genus Ferroglobus were retrieved in all brine samples. Finally, a high abundance of the strictly anaerobic methanogens (such as Methanosarcina members) within the active community suggests that anoxic conditions might occur in the lake brines. Our findings indicate perennially ice-covered Antarctic lakes as plausible terrestrial candidates for the study of the potential for extant life on different bodies of our solar system.

High Prevalence of Genetically Related Clostridium Difficile Strains at a Single Hemato-Oncology Ward Over 10 Years

Aims: Clostridium difficile (C. difficile) infection (CDI) is the main cause of healthcare-associated infectious diarrhea. We used whole-genome sequencing (WGS) to measure the prevalence and genetic variability of C. difficile at a single hemato-oncology ward over a 10 year period. Methods: Between 2008 and 2018, 2077 stool samples were obtained from diarrheal patients hospitalized at the Department of Lymphoma; of these, 618 were positive for toxin A/B. 140 isolates were then subjected to WGS on Ion Torrent PGM sequencer. Results: 36 and 104 isolates were recovered from 36 to 46 patients with single and multiple CDIs, respectively. Of these, 131 strains were toxigenic. Toxin gene profiles tcdA(+);tcdB(+);cdtA/cdtB(+) and tcdA(+);tcdB(+);cdtA/cdtB(-) were identified in 122 and nine strains, respectively. No isolates showed reduced susceptibility to metronidazole and vancomycin. All tested strains were resistant to ciprofloxacin, and 72.9, 42.9, and 72.9% of strains were resistant to erythromycin, clindamycin, or moxifloxacin, respectively. Multi-locus sequence typing (MLST) identified 23 distinct sequence types (STs) and two unidentified strains. Strains ST1 and ST42 represented 31 and 30.1% of all strains tested, respectively. However, while ST1 was detected across nearly all years studied, ST42 was detected only from 2009 to 2011. Conclusion: The high proportion of infected patients in 2008-2011 may be explained by the predominance of more transmissible and virulent C. difficile strains. Although this retrospective study was not designed to define outbreaks of C. difficile, the finding that most isolates exhibited high levels of genetic relatedness suggests nosocomial acquisition.

Genetic diversity of Listeria monocytogenes strains in ruminant abortion and rhombencephalitis cases in comparison with the natural environment

BACKGROUND

Listeria monocytogenes is the causative agent of listeriosis, a serious disease affecting both animals and humans. Here, multilocus sequence typing (MLST) was used to characterize the genetic diversity of Listeria monocytogenes strains isolated from the natural environment and animal clinical cases in Europe. The prevalence of clonal complexes (CCs) obtained was compared according to (i) the origin of isolation - clinical cases vs. natural environment - and (ii) the clinical form of animal listeriosis - rhombencephalitis vs. abortion. To this aim, two datasets were constructed. The clinical dataset consisted of 350 animal clinical isolates originating from France and Slovenia and supplemented with isolates from Switzerland and Great Britain. The natural environment dataset consisted of 253 isolates from the natural environment originating from Slovenia and supplemented with isolates from nine other European countries.

RESULTS

For the clinical cases, CC1, CC4-CC217 and CC412 were the most prevalent in rhombencephalitis and CC1, CC37 and CC4-CC217 in abortion. The hypervirulent CC1 and CC4-CC217 prevailed in both datasets. These results indicated that livestock is constantly exposed to hypervirulent CCs. CC1 was significantly associated with a clinical origin, whereas CC9, CC29 and CC14 were associated with the natural environment. CC1 was predominant among rhombencephalitis cases both in cattle and small ruminants, and its prevalence did not differ significantly between these two groups. A novel association of CC37 and CC6 with abortion cases was revealed.

CONCLUSIONS

Here, we show that CC1 and CC4-CC217 are prevalent in isolates of environmental and animal clinical origin, suggesting that ruminants are frequently exposed to hypervirulent CCs. The presence of CC4 in two mastitis cases calls for further attention due to direct threat to the consumer. We showed several associations between CCs and the origin of isolation or clinical form of listeriosis, e.g. CC37 and CC6 with abortion. This study improves our understanding of the population structure of L. monocytogenes isolates from the natural environment and animal clinical cases. Moreover, it provides a basis for future studies aiming to determine the underlying mechanisms of phenotypic traits of interest.

Microbial Assemblages in Pressurized Antarctic Brine Pockets (Tarn Flat, Northern Victoria Land): A Hotspot of Biodiversity and Activity

Two distinct pressurized hypersaline brine pockets (named TF4 and TF5), separated by a thin ice layer, were detected below an ice-sealed Antarctic lake. Prokaryotic (bacterial and archaeal) diversity, abundances (including virus-like particles) and metabolic profiles were investigated by an integrated approach, including traditional and new-generation methods. Although similar diversity indices were computed for both Bacteria and Archaea, distinct bacterial and archaeal assemblages were observed. Bacteroidetes and Gammaproteobacteria were more abundant in the shallowest brine pocket, TF4, and Deltaproteobacteria, mainly represented by versatile sulphate-reducing bacteria, dominated in the deepest, TF5. The detection of sulphate-reducing bacteria and methanogenic Archaea likely reflects the presence of a distinct synthrophic consortium in TF5. Surprisingly, members assigned to hyperthermophilic Crenarchaeota and Euryarchaeota were common to both brines, indicating that these cold habitats host the most thermally tolerant Archaea. The patterns of microbial communities were different, coherently with the observed microbiological diversity between TF4 and TF5 brines. Both the influence exerted by upward movement of saline brines from a sub-surface anoxic system and the possible occurrence of an ancient ice remnant from the Ross Ice Shelf were the likely main factors shaping the microbial communities.

High-Throughput Targeted Repeat Element Bisulfite Sequencing (HT-TREBS)

High-throughput targeted repeat element bisulfite sequencing (HT-TREBS) is designed to assay the methylation level of individual retrotransposon loci of a targeted family, in a locus-specific manner, and on a genome-wide scale. Briefly, genomic DNA is sheared and ligated to Ion Torrent A adaptors (with methylated cytosines), followed by bisulfite-conversion, and amplification with primers designed to bind the targeted retrotransposon. Since the primers carry the Ion Torrent P1 adaptor as a 5'-extension, the amplified library is ready to be size-selected and sequenced on a next-generation sequencing platform. Once sequenced, each retrotransposon is mapped to a particular genomic locus, which is achieved through ensuring at least a 10-bp overlap with flanking unique sequence, followed by the calculation of methylation levels of the mapped retrotransposon using a BiQ Analyzer HT. A complete protocol for library construction as well as the bioinformatics for HT-TREBS is described in this chapter.

Large-scale genomic analyses reveal the population structure and evolutionary trends of Streptococcus agalactiae strains in Brazilian fish farms

Streptococcus agalactiae is a major pathogen and a hindrance on tilapia farming worldwide. The aims of this work were to analyze the genomic evolution of Brazilian strains of S. agalactiae and to establish spatial and temporal relations between strains isolated from different outbreaks of streptococcosis. A total of 39 strains were obtained from outbreaks and their whole genomes were sequenced and annotated for comparative analysis of multilocus sequence typing, genomic similarity and whole genome multilocus sequence typing (wgMLST). The Brazilian strains presented two sequence types, including a newly described ST, and a non-typeable lineage. The use of wgMLST could differentiate each strain in a single clone and was used to establish temporal and geographical correlations among strains. Bayesian phylogenomic analysis suggests that the studied Brazilian population was co-introduced in the country with their host, approximately 60 years ago. Brazilian strains of S. agalactiae were shown to be heterogeneous in their genome sequences and were distributed in different regions of the country according to their genotype, which allowed the use of wgMLST analysis to track each outbreak event individually.

Genomic GC-Content Affects the Accuracy of 16S rRNA Gene Sequencing Based Microbial Profiling due to PCR Bias

Profiling of microbial community composition is frequently performed by partial 16S rRNA gene sequencing on benchtop platforms following PCR amplification of specific hypervariable regions within this gene. Accuracy and reproducibility of this strategy are two key parameters to consider, which may be influenced during all processes from sample collection and storage, through DNA extraction and PCR based library preparation to the final sequencing. In order to evaluate both the reproducibility and accuracy of 16S rRNA gene based microbial profiling using the Ion Torrent PGM platform, we prepared libraries and performed sequencing of a well-defined and validated 20-member bacterial DNA mock community on five separate occasions and compared results with the expected even distribution. In general the applied method had a median coefficient of variance of 11.8% (range 5.5–73.7%) for all 20 included strains in the mock community across five separate sequencing runs, with underrepresented strains generally showing the largest degree of variation. In terms of accuracy, mock community species belonging to Proteobacteria were underestimated, whereas those belonging to Firmicutes were mostly overestimated. This could be explained partly by premature read truncation, but to larger degree their genomic GC-content, which correlated negatively with the observed relative abundances, suggesting a PCR bias against GC-rich species during library preparation. Increasing the initial denaturation time during the PCR amplification from 30 to 120 s resulted in an increased average relative abundance of the three mock community members with the highest genomic GC%, but did not significantly change the overall evenness of the community distribution. Therefore, efforts should be made to optimize the PCR conditions prior to sequencing in order to maximize accuracy.