Comprehensive Analysis to Improve the Validation Rate for Single Nucleotide Variants Detected by Next-Generation Sequencing

Next-generation sequencing (NGS) has enabled the high-throughput discovery of germline and somatic mutations. However, NGS-based variant detection is still prone to errors, resulting in inaccurate variant calls. Here, we categorized the variants detected by NGS according to total read depth (TD) and SNP quality (SNPQ), and performed Sanger sequencing with 348 selected non-synonymous single nucleotide variants (SNVs) for validation. Using the SAMtools and GATK algorithms, the validation rate was positively correlated with SNPQ but showed no correlation with TD. In addition, common variants called by both programs had a higher validation rate than caller-specific variants. We further examined several parameters to improve the validation rate, and found that strand bias (SB) was a key parameter. SB in NGS data showed a strong difference between the variants passing validation and those that failed validation, showing a validation rate of more than 92% (filtering cutoff value: alternate allele forward [AF]≥20 and AF<80 in SAMtools, SB<–10 in GATK). Moreover, the validation rate increased significantly (up to 97–99%) when the variant was filtered together with the suggested values of mapping quality (MQ), SNPQ and SB. This detailed and systematic study provides comprehensive recommendations for improving validation rates, saving time and lowering cost in NGS analyses.     

Trends in First-Line Antiretroviral Therapy in Asia: Results from the TREAT Asia HIV Observational Database

Background

Antiretroviral therapy (ART) has evolved rapidly since its beginnings. This analysis describes trends in first-line ART use in Asia and their impact on treatment outcomes.

Methods

Patients in the TREAT Asia HIV Observational Database receiving first-line ART for ≥6 months were included. Predictors of treatment failure and treatment modification were assessed.

Results

Data from 4662 eligible patients was analysed. Patients started ART in 2003–2006 (n = 1419), 2007–2010 (n = 2690) and 2011–2013 (n = 553). During the observation period, tenofovir, zidovudine and abacavir use largely replaced stavudine. Stavudine was prescribed to 5.8% of ART starters in 2012/13. Efavirenz use increased at the expense of nevirapine, although both continue to be used extensively (47.5% and 34.5% of patients in 2012/13, respectively). Protease inhibitor use dropped after 2004. The rate of treatment failure or modification declined over time (22.1 [95%CI 20.7–23.5] events per 100 patient/years in 2003–2006, 15.8 [14.9–16.8] in 2007–2010, and 11.6 [9.4–14.2] in 2011–2013). Adjustment for ART regimen had little impact on the temporal decline in treatment failure rates but substantially attenuated the temporal decline in rates of modification due to adverse event. In the final multivariate model, treatment modification due to adverse event was significantly predicted by earlier period of ART initiation (hazard ratio 0.52 [95%CI 0.33–0.81], p = 0.004 for 2011–2013 versus 2003–2006), older age (1.56 [1.19–2.04], p = 0.001 for ≥50 years versus <30years), female sex (1.29 [1.11–1.50], p = 0.001 versus male), positive hepatitis C status (1.33 [1.06–1.66], p = 0.013 versus negative), and ART regimen (11.36 [6.28–20.54], p<0.001 for stavudine-based regimens versus tenofovir-based).

Conclusions

The observed trends in first-line ART use in Asia reflect changes in drug availability, global treatment recommendations and prescriber preferences over the past decade. These changes have contributed to a declining rate of treatment modification due to adverse event, but not to reductions in treatment failure.     

The species range‐size patterns for vascular plants of Seorak Mountain (Korea): Relationship between group of life forms and phytogeography affinity along the elevational gradient

Research on species richness patterns and the advanced elevational Rapoport rule (ERR) has been widespread in recent years; however, there is a lack of such research for the temperate mountainous regions in northeast Asia. Here, we collected plant species from the Seorak Mountain in northeast Asia through field surveys. The species were divided into 11 groups according to the life‐form types and phytogeography affinities of each species. The ERR was evaluated using Steven''s method and by examining the species richness patterns of each group. The species richness patterns revealed a positive multimodal pattern along the elevation gradient, but phytogeography affinities (increasing trend) and life‐form analysis (unimodal) exhibited different patterns. The elevation gradients (1,350 m for the mean elevation–range relationships), which are affected by the boundary effect and different life forms, did not consistently support the ERR. However, herbs as well as rare, endemic, and red list species showed consistent support for the ERR, which could be attributed to the influence by phytogeography affinities. Therefore, the results from Seorak Mountain showed that the ERR was not consistent for different plant life forms in the same area; however, phytogeography affinities could support and explain ERR.     

Genome-wide copy number analysis uncovers a new HSCR gene: NRG3

Hirschsprung disease (HSCR) is a congenital disorder characterized by aganglionosis of the distal intestine. To assess the contribution of copy number variants (CNVs) to HSCR, we analysed the data generated from our previous genome-wide association study on HSCR patients, whereby we identified NRG1 as a new HSCR susceptibility locus. Analysis of 129 Chinese patients and 331 ethnically matched controls showed that HSCR patients have a greater burden of rare CNVs (p = 1.50 × 10(-5)), particularly for those encompassing genes (p = 5.00 × 10(-6)). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64 × 10(-3)). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36 × 10(-5)) and 2 de novo duplications among patients and two deletions among controls. Importantly, NRG3 is a paralog of NRG1. Stratification of patients by presence/absence of HSCR-associated syndromes showed that while syndromic-HSCR patients carried significantly longer CNVs than the non-syndromic or controls (p = 1.50 × 10(-5)), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00 × 10(-6)) or the syndromic counterpart. Our results suggest a role for NRG3 in HSCR etiology and provide insights into the relative contribution of structural variants in both syndromic and non-syndromic HSCR. This would be the first genome-wide catalog of copy number variants identified in HSCR.     

Development of a monitoring vector for Leuconostoc mesenteroides using the green fluorescent protein gene

The vector pCW5 with plasmid pC7, originally isolated in Lactobacillus paraplantarum C7 derived from kimchi, was constructed using a p32 strong promoter, the pC7 replicon, and green fluorescent protein (GFP) as the reporter. The constructed vector was transformed into E. coli and Leuconostoc mesenteroides, and GFP expression detected using a Western blot analysis. GFP fluorescence was recognized in E. coli and Leuconostoc mesenteroides using a confocal microscope. In addition, GFP fluorescence was also clearly detected in several industrially important lactic acid bacteria (LAB), including Lactobacillus bulgaricus, Lactobacillus paraplantarum, and Lactobacillus plantarum. Thus, pCW5 was shown to be effective for Leuconostoc mesenteroides when using GFP as the reporter, and it can also be used as a broad-host-range vector for other lactic acid bacteria.     

Role of NKG2D in Obesity-Induced Adipose Tissue Inflammation and Insulin Resistance

The early events that initiate inflammation in the adipose tissue during obesity are not well defined. It is unclear whether the recruitment of CD8 T cells to the adipose tissue during onset of obesity occurs through antigen-dependent or -independent processes. We have previously shown that interaction between NKG2D (natural-killer group 2, member D) and its ligand Rae-1ε is sufficient to recruit cytotoxic T lymphocytes to the pancreas and induce insulitis. Here, we tested whether NKG2D–NKG2D ligand interaction is also involved in obesity-induced adipose tissue inflammation and insulin resistance. We observed a significant induction of NKG2D ligand expression in the adipose tissue of obese mice, especially during the early stages of obesity. However, mice lacking NKG2D developed similar levels of insulin resistance and adipose tissue inflammation compared to control mice when placed on a high-fat diet. Moreover, overexpression of Rae-1ε in the adipose tissue did not increase immune cell infiltration to the adipose tissue either in the setting of a normal or high-fat diet. These results indicate that, unlike in the pancreas, NKG2D–NKG2D ligand interaction does not play a critical role in obesity-induced inflammation in the adipose tissue.     

The Postnatal Development of the Cerebellum— A fMRI and Silver Study

Summary The aim of this study was to evaluate the postnatal development of the cerebella of the pig and to compare this with the activation of the fMRI. The cells in the cerebella were studied by silver technique and the activation of the fMRI in the cerebella was initiated by flexion and extension of the hind paw. Our results showed an increase of the branching of the cells of the cerebellar cortex postnatally, coordinated with registration of fMRI active sites in the cerebella at 6-month postnatal. We concluded that the full maturation of the cerebella was around 6-month postnatal in the pig.