Distinguishing low frequency mutations from RT-PCR and sequence errors in viral deep sequencing data

Background

RNA viruses have high mutation rates and exist within their hosts as large, complex and heterogeneous populations, comprising a spectrum of related but non-identical genome sequences. Next generation sequencing is revolutionising the study of viral populations by enabling the ultra deep sequencing of their genomes, and the subsequent identification of the full spectrum of variants within the population. Identification of low frequency variants is important for our understanding of mutational dynamics, disease progression, immune pressure, and for the detection of drug resistant or pathogenic mutations. However, the current challenge is to accurately model the errors in the sequence data and distinguish real viral variants, particularly those that exist at low frequency, from errors introduced during sequencing and sample processing, which can both be substantial.

Results

We have created a novel set of laboratory control samples that are derived from a plasmid containing a full-length viral genome with extremely limited diversity in the starting population. One sample was sequenced without PCR amplification whilst the other samples were subjected to increasing amounts of RT and PCR amplification prior to ultra-deep sequencing. This enabled the level of error introduced by the RT and PCR processes to be assessed and minimum frequency thresholds to be set for true viral variant identification. We developed a genome-scale computational model of the sample processing and NGS calling process to gain a detailed understanding of the errors at each step, which predicted that RT and PCR errors are more likely to occur at some genomic sites than others. The model can also be used to investigate whether the number of observed mutations at a given site of interest is greater than would be expected from processing errors alone in any NGS data set. After providing basic sample processing information and the site’s coverage and quality scores, the model utilises the fitted RT-PCR error distributions to simulate the number of mutations that would be observed from processing errors alone.

Conclusions

These data sets and models provide an effective means of separating true viral mutations from those erroneously introduced during sample processing and sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1456-x) contains supplementary material, which is available to authorized users.     

A specific protein disorder catalyzer of HIV-1 Nef

The HIV-1 auxiliary protein Nef is required for the onset and progression of AIDS in HIV-1-infected persons. Here, we have deciphered the mode of action of a second-generation inhibitor of Nef, DLC27-14, presenting a competitive IC(50) of ～16 μM measured by MALDI-TOF experiments. Thermal protein denaturation experiments revealed a negative effect on stability of Nef in the presence of a saturating concentration of the inhibitor. The destabilizing action of DLC27-14 was confirmed by a HIV protease-based experiment, in which the protease sensitivity of DLC27-14-bound Nef was three times as high as that of apo Nef. The only compatible docking modes of action for DLC27-14 suggest that DLC27-14 promotes an opening of two α-helices that would destabilize the Nef core domain. DLC27-14 thus acts as a specific protein disorder catalyzer that destabilizes the folded conformation of the protein. Our results open novel avenues toward the development of next-generation Nef inhibitors.     

Mendelian breeding units versus standard sampling strategies: Mitochondrial DNA variation in southwest Sardinia

We report a sampling strategy based on Mendelian Breeding Units (MBUs), representing an interbreeding group of individuals sharing a common gene pool. The identification of MBUs is crucial for case-control experimental design in association studies. The aim of this work was to evaluate the possible existence of bias in terms of genetic variability and haplogroup frequencies in the MBU sample, due to severe sample selection. In order to reach this goal, the MBU sampling strategy was compared to a standard selection of individuals according to their surname and place of birth. We analysed mitochondrial DNA variation (first hypervariable segment and coding region) in unrelated healthy subjects from two different areas of Sardinia: the area around the town of Cabras and the western Campidano area. No statistically significant differences were observed when the two sampling methods were compared, indicating that the stringent sample selection needed to establish a MBU does not alter original genetic variability and haplogroup distribution. Therefore, the MBU sampling strategy can be considered a useful tool in association studies of complex traits.     

Copper exposure effects on yeast mitochondrial proteome

Mitochondria play an important role on the entire cellular copper homeostatic mechanisms. Alteration of cellular copper levels may thus influence mitochondrial proteome and its investigation represents an important contribution to the general understanding of copper-related cellular effects. In these study we have performed an organelle targeted proteomic investigation focusing our attention on the effect of non-lethal 1mM copper concentration on Saccharomyces cerevisiae mitochondrial proteome. Functional copper effects on yeast mitochondrial proteome were evaluated by using both 2D electrophoresis (2-DE) and liquid chromatography coupled with tandem mass spectrometry. Proteomic data have been then analyzed by different unsupervised meta-analysis approaches that highlight the impairment of mitochondrial functions and the activation of oxidative stress response. Interestingly, our data have shown that stress response generated by 1mM copper treatment determines the activation of S. cerevisiae survival pathway. To investigate these findings we have treated yeast cells responsiveness to copper with hydrogen peroxide and observed a protective role of this metal. These results are suggestive of a copper role in the protection from oxidative stress possibly due to the activation of mechanisms involved in cellular survival and growth.     

Principles and patterns in the interaction between mono-heme cytochrome c and its partners in electron transfer processes

Cytochromes c are very widespread proteins that play key roles in the electron transfer events associated to a wide variety of physiological redox processes. The function of cytochromes c is, at the broad level, to interact with different partners in order to allow electrons to flow from one protein to another. Here, we focused our attention on the protein-protein interactions that involve mono-heme cytochrome c domains in order to identify possible general vs. specific patterns of intermolecular interactions at the structural level. We observed that a number of physico-chemical properties are statistically different in transient vs. permanent and fused complexes. These include the extent of the protein interface area, the amino acid composition and the packing density at the interface. The understanding of the features of transient redox complexes is of particular importance because of the difficulty of obtaining co-crystals that preserve the physiologically relevant configuration. In addition, we identified three different structural modes of interaction that cover all the structurally characterized cytochrome c interactions except one. The mode of interaction does not correlate with the nature of the complex (transient, permanent, fused). Regardless of the mode of interaction, the distance between the heme iron and the partner metal center or organic cofactor center of mass is typically around 19-20 ? for complexes permitting direct electron transfer between the two sites.     

In-cell NMR in E. coli to monitor maturation steps of hSOD1

In-cell NMR allows characterizing the folding state of a protein as well as posttranslational events at molecular level, in the cellular context. Here, the initial maturation steps of human copper, zinc superoxide dismutase 1 are characterized in the E. coli cytoplasm by in-cell NMR: from the apo protein, which is partially unfolded, to the zinc binding which causes its final quaternary structure. The protein selectively binds only one zinc ion, whereas in vitro also the copper site binds a non-physiological zinc ion. However, no intramolecular disulfide bridge formation occurs, nor copper uptake, suggesting the need of a specific chaperone for those purposes.     

Pyrrolizidine alkaloids from Anchusa strigosa and their antifeedant activity

The pyrrolizidine alkaloid (PA) content of flowers, leaves, and roots of Anchusa strigosa (Boraginaceae) was analysed by ESI-LC-MS. Six PAs, including two new natural compounds, were detected, characterized by NMR spectroscopy, and quantified in each plant organ. The results indicated that the highest total concentration of PAs was in the leaves (23.63 mg/g of dried part), followed by the flowers (19.77 mg/g), and finally by the roots (1.80 mg/g). All PAs isolated were subjected to Spodoptera exigua and Pieris brassicae larvae. Feeding activity by both herbivore species using a bioassay was inhibited up to circa 75% depending on PA and applied concentration.     

The effect of deficit irrigation on seasonal variations of plant water use in Olea europaea L.

A field experiment on olive trees (Olea europaea L.) was designed with the objective to search for an optimum irrigation scheduling by analyzing the possible effects of deficit irrigation. Treatments were: a non-irrigated control (rainfed) and three treatments that received seasonal water amount equivalent to 33 and 66% of crop evapotranspiration (ET_C) in the period August–September (respectively 33II and 66II), and 66% of (ET_C) from late May to early October (66I-II). Atmospheric evaporative demand and soil moisture conditions were regularly monitored. Irrigation effects on plant water relations were characterized throughout a growing season. Whole-plant water use, in deficit irrigated (66I-II) and rainfed olive trees, was determined using a xylem sap flow method (compensation heat-pulse technique). The magnitude of variations in water use and the seasonal dynamic of water relations varied among treatments, suggesting that olive trees were strongly responsive to both irrigation amount and time. Physiological parameters responded to variations in tree water status, soil moisture conditions and atmospheric evaporative demand. All measurements of tree water status were highly correlated with one another. There was a considerable degree of agreement between daily transpiration deduced from heat-pulse velocity and that determined by calibration using the water balance technique. Deficit irrigation during the whole summer (66I-II) resulted in improved plant water relations with respect to other watering regimes; while, severe regulated deficit irrigation differentiated only slightly 33II treatment from rainfed plants. Nevertheless, regulated deficit irrigation of olive trees after pit hardening (66II) could be recommended, at least in soil, cultivar and environmental conditions of this study.     

Asymptomatic Helicobacter pylori infection increases asymmetric dimethylarginine levels in healthy subjects

BACKGROUND: Chronic infections have been demonstrated to be early factors of atherosclerosis and cardiovascular diseases, and their relevance increases when they are caused by agents with extremely broad spectrum of disease outcome such as Helicobacter pylori. The consequent endothelial impairment leads to a reduced bioavailability of nitric oxide. Increasing evidences have pointed out that the endogenous inhibitor of nitric oxide synthase, asymmetric dimethylarginine, defined as a risk factor for cardiovascular disease, may increase in infections and plays an important role impairing the vascular functions of the endothelium. Starting from these findings, we aim to investigate whether H. pylori may affect asymmetric dimethylarginine levels. MATERIALS AND METHODS: The study was carried out on a group of 186 subjects (age 46.2 +/- 14.9 years). We evaluated asymmetric dimethylarginine, symmetric dimethylarginine, L-arginine, presence of H. pylori by 13C-urea breath test, and the main parameters of glyco and lipo metabolic balance. RESULTS: Increased levels of asymmetric dimethylarginine were found in H. pylori-positive subjects with respect to H. pylori-negative subjects (0.46 x/ / 1.13 versus 0.42 x/ / 1.23 mol/l, p < .001, respectively). No differences were detected in L-arginine levels between the two groups. Multiple regression analysis performed in H. pylori-positive subjects and H. pylori-negative subjects showed profound differences in the variables related to asymmetric dimethylarginine (R2 = 66.9%, p < .01 versus 34.3%, p < .01, respectively) and symmetric dimethylarginine (R2 = 39.2%, p < .01 versus 20.6%, p = .09, respectively) levels. CONCLUSIONS: Our data clearly demonstrate that H. pylori infection increases asymmetric dimethylarginine levels. Moreover, this infection causes a profound metabolic modification that alters the role of the known determinants of asymmetric dimethylarginine levels. We conclude that H. pylori infection must be taken into account as a cause of increased asymmetric dimethylarginine levels and that the eradication of H. pylori may therefore lead to a decrease in asymmetric dimethylarginine levels, which is a further reason for the reduction of the risk for cardiovascular disease in this large portion of population.     

Comparative analysis of the ADAM and ADAMTS families

The "A Disintegrin And Metalloproteinase" (ADAM) protein family and the "A Disintegrin-like And Metalloproteinase with ThromboSpondin motifs" (ADAMTS) protein family are two related families of human proteins. The similarities and differences between these two families have been investigated using phylogenetic trees and homology modeling. The phylogenetic analysis indicates that the two families are well differentiated, even when only the common metalloprotease domain is taken into account. Within the ADAM family, several proteins are lacking the binding motif for the catalytic zinc in the active site and thus presumably lack any catalytic activity. These proteins tend to cluster within the ADAM phylogenetic tree and are expressed in specific tissues, suggesting a functional differentiation. The present analysis allows us to propose the following: (i) ADAMTS proteins have a conserved role in the human organism as proteases, with some differentiation in terms of substrate specificity; (ii) ADAM proteins can act as proteases and/or mediators of intermolecular interactions; (iii) proteolytically active ADAMs tend to be more ubiquitously expressed than the inactive ones.