Prediction of Protein-Destabilizing Polymorphisms by Manual Curation with Protein Structure

The relationship between sequence polymorphisms and human disease has been studied mostly in terms of effects of single nucleotide polymorphisms (SNPs) leading to single amino acid substitutions that change protein structure and function. However, less attention has been paid to more drastic sequence polymorphisms which cause premature termination of a protein’s sequence or large changes, insertions, or deletions in the sequence. We have analyzed a large set (n = 512) of insertions and deletions (indels) and single nucleotide polymorphisms causing premature termination of translation in disease-related genes. Prediction of protein-destabilization effects was performed by graphical presentation of the locations of polymorphisms in the protein structure, using the Genomes TO Protein (GTOP) database, and manual annotation with a set of specific criteria. Protein-destabilization was predicted for 44.4% of the nonsense SNPs, 32.4% of the frameshifting indels, and 9.1% of the non-frameshifting indels. A prediction of nonsense-mediated decay allowed to infer which truncated proteins would actually be translated as defective proteins. These cases included the proteins linked to diseases inherited dominantly, suggesting a relation between these diseases and toxic aggregation. Our approach would be useful in identifying potentially aggregation-inducing polymorphisms that may have pathological effects.     

Human Papillomavirus Type 6 and 11 Genetic Variants Found in 71 Oral and Anogenital Epithelial Samples from Australia

Genetic variation of 49 human papillomavirus (HPV) 6 and 22 HPV11 isolates from recurrent respiratory papillomatosis (RRP) (n = 17), genital warts (n = 43), anal cancer (n = 6) and cervical neoplasia cells (n = 5), was determined by sequencing the long control region (LCR) and the E6 and E7 genes. Comparative analysis of genetic variability was examined to determine whether different disease states resulting from HPV6 or HPV11 infection cluster into distinct variant groups. Sequence variation analysis of HPV6 revealed that isolates cluster into variants within previously described HPV6 lineages, with the majority (65%) clustering to HPV6 sublineage B1 across the three genomic regions examined. Overall 72 HPV6 and 25 HPV11 single nucleotide variations, insertions and deletions were observed within samples examined. In addition, missense alterations were observed in the E6/E7 genes for 6 HPV6 and 5 HPV11 variants. No nucleotide variations were identified in any isolates at the four E2 binding sites for HPV6 or HPV11, nor were any isolates found to be identical to the HPV6 lineage A or HPV11 sublineage A1 reference genomes. Overall, a high degree of sequence conservation was observed between isolates across each of the regions investigated for both HPV6 and HPV11. Genetic variants identified a slight association with HPV6 and anogenital lesions (p = 0.04). This study provides important information on the genetic diversity of circulating HPV 6 and HPV11 variants within the Australian population and supports the observation that the majority of HPV6 isolates cluster to the HPV6 sublineage B1 with anogenital lesions demonstrating an association with this sublineage (p = 0.02). Comparative analysis of Australian isolates for both HPV6 and HPV11 to those from other geographical regions based on the LCR revealed a high degree of sequence similarity throughout the world, confirming previous observations that there are no geographically specific variants for these HPV types.     

Digital Surveillance: A Novel Approach to Monitoring the Illegal Wildlife Trade

A dearth of information obscures the true scale of the global illegal trade in wildlife. Herein, we introduce an automated web crawling surveillance system developed to monitor reports on illegally traded wildlife. A resource for enforcement officials as well as the general public, the freely available website, http://www.healthmap.org/wildlifetrade, provides a customizable visualization of worldwide reports on interceptions of illegally traded wildlife and wildlife products. From August 1, 2010 to July 31, 2011, publicly available English language illegal wildlife trade reports from official and unofficial sources were collected and categorized by location and species involved. During this interval, 858 illegal wildlife trade reports were collected from 89 countries. Countries with the highest number of reports included India (n = 146, 15.6%), the United States (n = 143, 15.3%), South Africa (n = 75, 8.0%), China (n = 41, 4.4%), and Vietnam (n = 37, 4.0%). Species reported as traded or poached included elephants (n = 107, 12.5%), rhinoceros (n = 103, 12.0%), tigers (n = 68, 7.9%), leopards (n = 54, 6.3%), and pangolins (n = 45, 5.2%). The use of unofficial data sources, such as online news sites and social networks, to collect information on international wildlife trade augments traditional approaches drawing on official reporting and presents a novel source of intelligence with which to monitor and collect news in support of enforcement against this threat to wildlife conservation worldwide.     

Spatial Variation in Genetic Diversity and Natural Selection on the Thrombospondin-Related Adhesive Protein Locus of Plasmodium vivax (PvTRAP)

Thrombospondin-related adhesive protein (TRAP) of malaria parasites is essential for sporozoite motility and invasions into mosquito’s salivary gland and vertebrate’s hepatocyte; thereby, it is a promising target for pre-erythrocytic vaccine. TRAP of Plasmodium vivax (PvTRAP) exhibits sequence heterogeneity among isolates, an issue relevant to vaccine development. To gain insights into variation in the complete PvTRAP sequences of parasites in Thailand, 114 vivax malaria patients were recruited in 2006–2007 from 4 major endemic provinces bordering Myanmar (Tak in the northwest, n = 30 and Prachuap Khirikhan in the southwest, n = 25), Cambodia (Chanthaburi in the east, n = 29) and Malaysia (Yala and Narathiwat in the south, n = 30). In total, 26 amino acid substitutions were detected and 9 of which were novel, resulting in 44 distinct haplotypes. Haplotype and nucleotide diversities were lowest in southern P. vivax population while higher levels of diversities were observed in other populations. Evidences of positive selection on PvTRAP were demonstrated in domains II and IV and purifying selection in domains I, II and VI. Genetic differentiation was significant between each population except that between populations bordering Myanmar where transmigration was common. Regression analysis of pairwise linearized Fst and geographic distance suggests that P. vivax populations in Thailand have been isolated by distance. Sequence diversity of PvTRAP seems to be temporally stable over one decade in Tak province based on comparison of isolates collected in 1996 (n = 36) and 2006–2007. Besides natural selection, evidences of intragenic recombination have been supported in this study that could maintain and further generate diversity in this locus. It remains to be investigated whether amino acid substitutions in PvTRAP could influence host immune responses although several predicted variant T cell epitopes drastically altered the epitope scores. Knowledge on geographic diversity in PvTRAP constitutes an important basis for vaccine design provided that vaccination largely confers variant-specific immunity.     

Non-Immunogenicity of Overlapping Gag Peptides Pulsed on Autologous Cells after Vaccination of HIV Infected Individuals

Background

HIV Gag-specific CD4+ and CD8+ T-cell responses are important for HIV immune control. Pulsing overlapping Gag peptides on autologous lymphocytes (OPAL) has proven immunogenic and effective in reducing viral loads in multiple pigtail macaque studies, warranting clinical evaluation.

Methodology

We performed a phase I, single centre, placebo-controlled, double-blinded and dose-escalating study to evaluate the safety and preliminary immunogenicity of a novel therapeutic vaccine approach ‘OPAL-HIV-Gag(c)’. This vaccine is comprised of 120 15mer peptides, overlapping by 11 amino acids, spanning the HIV Gag C clade sequence proteome, pulsed on white blood cells enriched from whole blood using a closed system, followed by intravenous reinfusion. Patients with undetectable HIV viral loads (<50 copies/ml plasma) on HAART received four administrations at week 0, 4, 8 and 12, and were followed up for 12 weeks post-treatment. Twenty-three people were enrolled in four groups: 12 mg (n = 6), 24 mg (n = 7), 48 mg (n = 2) or matching placebo (n = 8) with 18 immunologically evaluable. T-cell immunogenicity was assessed by IFNγ ELIspot and intracellular cytokine staining (ICS).

Results

The OPAL-HIV-Gag(c) peptides were antigenic in vitro in 17/17 subjects. After vaccination with OPAL-HIV-Gag(c), 1/6 subjects at 12 mg and 1/6 subjects at 24 mg dose groups had a 2- and 3-fold increase in ELIspot magnitudes from baseline, respectively, of Gag-specific CD8+ T-cells at week 14, compared to 0/6 subjects in the placebo group. No Gag-specific CD4+ T-cell responses or overall change in Rev, Nef, Tat and CMV specific responses were detected. Marked, transient and self-limiting lymphopenia was observed immediately post-vaccination (4 hours) in OPAL-HIV-Gag(c) but not in placebo recipients, with median fall from 1.72 to 0.67 million lymphocytes/mL for active groups (P<0.001), compared to post-placebo from 1.70 to 1.56 lymphocytes/ml (P = 0.16).

Conclusion/Significance

Despite strong immunogenicity observed in several Macaca nemestrina studies using this approach, OPAL-HIV-Gag(c) was not significantly immunogenic in humans and improved methods of generating high-frequency Gag-specific T-cell responses are required.

Name of Registry

ClinicalTrials.gov, Registry number: {"type":"clinical-trial","attrs":{"text":"NCT01123915","term_id":"NCT01123915"}}NCT01123915, URL trial registry database: http://www.clinicaltrials.gov/ct2/results?term=OPAL-HIV-1001&Search=Search     

SIFT Indel: Predictions for the Functional Effects of Amino Acid Insertions/Deletions in Proteins

Indels in the coding regions of a gene can either cause frameshifts or amino acid insertions/deletions. Frameshifting indels are indels that have a length that is not divisible by 3 and subsequently cause frameshifts. Indels that have a length divisible by 3 cause amino acid insertions/deletions or block substitutions; we call these 3n indels. The new amino acid changes resulting from 3n indels could potentially affect protein function. Therefore, we construct a SIFT Indel prediction algorithm for 3n indels which achieves 82% accuracy, 81% sensitivity, 82% specificity, 82% precision, 0.63 MCC, and 0.87 AUC by 10-fold cross-validation. We have previously published a prediction algorithm for frameshifting indels. The rules for the prediction of 3n indels are different from the rules for the prediction of frameshifting indels and reflect the biological differences of these two different types of variations. SIFT Indel was applied to human 3n indels from the 1000 Genomes Project and the Exome Sequencing Project. We found that common variants are less likely to be deleterious than rare variants. The SIFT indel prediction algorithm for 3n indels is available at http://sift-dna.org/     

Molecular Characterization of Pneumococcal Isolates from Pets and Laboratory Animals

Background

Between 1986 and 2008 Streptococcus pneumoniae was isolated from 41 pets/zoo animals (guinea pigs (n = 17), cats (n = 12), horses (n = 4), dogs (n = 3), dolphins (n = 2), rat (n = 2), gorilla (n = 1)) treated in medical veterinary laboratories and zoos, and 44 laboratory animals (mastomys (multimammate mice; n = 32), mice (n = 6), rats (n = 4), guinea pigs (n = 2)) during routine health monitoring in an animal facility. S. pneumoniae was isolated from nose, lung and respiratory tract, eye, ear and other sites.

Methodology/Principal Findings

Carriage of the same isolate of S. pneumoniae over a period of up to 22 weeks was shown for four mastomys. Forty-one animals showed disease symptoms. Pneumococcal isolates were characterized by optochin sensitivity, bile solubility, DNA hybridization, pneumolysin PCR, serotyping and multilocus sequence typing. Eighteen of the 32 mastomys isolates (56%) were optochin resistant, all other isolates were optochin susceptible. All mastomys isolates were serotype 14, all guinea pig isolates serotype 19F, all horse isolates serotype 3. Rats had serotypes 14 or 19A, mice 33A or 33F. Dolphins had serotype 23F, the gorilla serotype 14. Cats and dogs had many different serotypes. Four isolates were resistant to macrolides, three isolates also to clindamycin and tetracyclin. Mastomys isolates were sequence type (ST) 15 (serotype 14), an ST/serotype combination commonly found in human isolates. Cats, dogs, pet rats, gorilla and dolphins showed various human ST/serotype combinations. Lab rats and lab mice showed single locus variants (SLV) of human STs, in human ST/serotype combinations. All guinea pig isolates showed the same completely new combination of known alleles. The horse isolates showed an unknown allele combination and three new alleles.

Conclusions/Significance

The isolates found in mastomys, mice, rats, cats, dogs, gorilla and dolphins are most likely identical to human pneumococcal isolates. Isolates from guinea pigs and horses appear to be specialized clones for these animals. Our data redraw attention to the fact that pneumococci are not strictly human pathogens. Pet animals that live in close contact to humans, especially children, can be infected by human isolates and also carriage of even resistant isolates is a realistic possibility.     

GINDEL: Accurate Genotype Calling of Insertions and Deletions from Low Coverage Population Sequence Reads

Insertions and deletions (indels) are important types of structural variations. Obtaining accurate genotypes of indels may facilitate further genetic study. There are a few existing methods for calling indel genotypes from sequence reads. However, none of these tools can accurately call indel genotypes for indels of all lengths, especially for low coverage sequence data. In this paper, we present GINDEL, an approach for calling genotypes of both insertions and deletions from sequence reads. GINDEL uses a machine learning approach which combines multiple features extracted from next generation sequencing data. We test our approach on both simulated and real data and compare with existing tools, including Genome STRiP, Pindel and Clever-sv. Results show that GINDEL works well for deletions larger than 50 bp on both high and low coverage data. Also, GINDEL performs well for insertion genotyping on both simulated and real data. For comparison, Genome STRiP performs less well for shorter deletions (50–200 bp) on both simulated and real sequence data from the 1000 Genomes Project. Clever-sv performs well for intermediate deletions (200–1500 bp) but is less accurate when coverage is low. Pindel only works well for high coverage data, but does not perform well at low coverage. To summarize, we show that GINDEL not only can call genotypes of insertions and deletions (both short and long) for high and low coverage population sequence data, but also is more accurate and efficient than other approaches. The program GINDEL can be downloaded at: http://sourceforge.net/p/gindel     

Angiotensin Converting Enzyme Inhibitor and HMG-CoA Reductase Inhibitor as Adjunct Treatment for Persons with HIV Infection: A Feasibility Randomized Trial

Background

Treatments that reduce inflammation and cardiovascular disease (CVD) risk among individuals with HIV infection receiving effective antiretroviral therapy (ART) are needed.

Design and Methods

We conducted a 2×2 factorial feasibility study of lisinopril (L) (10 mg daily) vs L-placebo in combination with pravastatin (P) (20 mg daily) vs P-placebo among participants receiving ART with undetectable HIV RNA levels, a Framingham 10 year risk score (FRS) ≥3%, and no indication for ACE-I or statin therapy. Tolerability and adherence were evaluated. Longitudinal mixed models assessed changes in blood pressure (BP), blood lipids, and inflammatory biomarkers from baseline through months 1 and 4.

Results

Thirty-seven participants were randomized and 34 [lisinopril/pravastatin (n = 9), lisinopril/P-placebo (n = 8), L-placebo/pravastatin (n = 9), L-placebo/P-placebo (n = 8)] attended at least one follow-up visit. Participants were 97% male, 41% white, 67% were current smokers, and 65% were taking a protease inhibitor. Median age was 48 years, CD4 count 483 cells/mm³, FRS 7.79%, total cholesterol 184 mg/dL, and LDL-C 95 mg/dL. There was no treatment difference for pravastatin vs P-placebo in total cholesterol, LDL-C, or any of the inflammatory biomarkers. Participants randomized to lisinopril vs. L-placebo had significant declines in diastolic BP (−3.3 mmHg, p = 0.05), hsCRP (−0.61 µg/mL, p = 0.02) and TNF-α (−0.17 pg/mL, p = 0.04). Participants taking lisinopril vs L-placebo were more likely to report missed doses (88 vs 35%; p = 0.001) and have adherence <90% by pill count (42 vs. 0%; p = 0.02). Few participants from either group reported side effects (n = 3 vs. n = 1).

Conclusions

The modest BP changes and decreased adherence with lisinopril and absence of lipid differences with pravastatin suggest future studies of these drug classes should consider a run-in period to assess adherence and use a different statin. Our results also indicate that ACE-I therapy may have anti-inflammatory benefits for ART-treated persons with HIV infection and this should be further evaluated.

Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT00982189","term_id":"NCT00982189"}}NCT00982189     

ACE Inhibitors Potently Reduce Vascular Inflammation,Results of an Open Proof-Of-Concept Study in the Abdominal Aortic Aneurysm

Background

Independent of their blood pressure lowering effect, ACE inhibitors are thought to reduce vascular inflammation. The clinical relevance of this effect is unclear with the current knowledge. Abdominal aortic aneurysms (AAA) are characterized by a broad, non-specific inflammatory response, and thus provide a clinical platform to evaluate the anti-inflammatory potential of ACE inhibitors.

Methods and Results

Eleven patients scheduled for open AAA repair received ramipril (5 mg/day) during 2–4 weeks preceding surgery. Aortic wall samples were collected during surgery, and compared to matched samples obtained from a biobank. An anti-inflammatory potential was evaluated in a comprehensive analysis that included immunohistochemistry, mRNA and protein analysis. A putative effect of ACE inhibitors on AAA growth was tested separately by comparing 18-month growth rate of patients on ACE inhibitors (n = 82) and those not taking ACE inhibitors (n = 204). Ramipril reduces mRNA expression of multiple pro-inflammatory cytokines such as IL-1β, IL-6, IL-8, TNF -α, Interferon-, and MCP-1, as well as aortic wall IL-8 and MCP-1 (P = 0.017 and 0.008, respectively) protein content. The is followed by clear effects on cell activation that included a shift towards anti-inflammatory macrophage (M2) subtype. Evaluation of data from the PHAST cohort did not indicate an effect of ACE inhibitors on 18-month aneurysm progression (mean difference at 18 months: −0.24 mm (95% CI: −0.90–0.45, P = NS).

Conclusions

ACE inhibition quenches multiple aspects of vascular inflammation in AAA. However, this does not translate into reduced aneurysm growth.

Trial Registration

Nederlands Trial Register 1345.     

Plasma Uric Acid Levels Correlate with Inflammation and Disease Severity in Malian Children with Plasmodium falciparum Malaria

Background

Plasmodium falciparum elicits host inflammatory responses that cause the symptoms and severe manifestations of malaria. One proposed mechanism involves formation of immunostimulatory uric acid (UA) precipitates, which are released from sequestered schizonts into microvessels. Another involves hypoxanthine and xanthine, which accumulate in parasitized red blood cells (RBCs) and may be converted by plasma xanthine oxidase to UA at schizont rupture. These two forms of ‘parasite-derived’ UA stimulate immune cells to produce inflammatory cytokines in vitro.

Methods and Findings

We measured plasma levels of soluble UA and inflammatory cytokines and chemokines (IL-6, IL-10, sTNFRII, MCP-1, IL-8, TNFα, IP-10, IFNγ, GM-CSF, IL-1β) in 470 Malian children presenting with uncomplicated malaria (UM), non-cerebral severe malaria (NCSM) or cerebral malaria (CM). UA levels were elevated in children with NCSM (median 5.74 mg/dl, 1.21-fold increase, 95% CI 1.09–1.35, n = 23, p = 0.0007) and CM (median 5.69 mg/dl, 1.19-fold increase, 95% CI 0.97–1.41, n = 9, p = 0.0890) compared to those with UM (median 4.60 mg/dl, n = 438). In children with UM, parasite density and plasma creatinine levels correlated with UA levels. These UA levels correlated with the levels of seven cytokines [IL-6 (r = 0.259, p<0.00001), IL-10 (r = 0.242, p<0.00001), sTNFRII (r = 0.221, p<0.00001), MCP-1 (r = 0.220, p<0.00001), IL-8 (r = 0.147, p = 0.002), TNFα (r = 0.132, p = 0.006) and IP-10 (r = 0.120, p = 0.012)]. In 39 children, UA levels were 1.49-fold (95% CI 1.34–1.65; p<0.0001) higher during their malaria episode [geometric mean titer (GMT) 4.67 mg/dl] than when they were previously healthy and aparasitemic (GMT 3.14 mg/dl).

Conclusions

Elevated UA levels may contribute to the pathogenesis of P. falciparum malaria by activating immune cells to produce inflammatory cytokines. While this study cannot identify the cause of elevated UA levels, their association with parasite density and creatinine levels suggest that parasite-derived UA and renal function may be involved. Defining pathogenic roles for parasite-derived UA precipitates, which we have not directly studied here, requires further investigation.

Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT00669084","term_id":"NCT00669084"}}NCT00669084     

Population-Specific Haplotype Association of the Postsynaptic Density Gene DLG4 with Schizophrenia,in Family-Based Association Studies

The post-synaptic density (PSD) of glutamatergic synapses harbors a multitude of proteins critical for maintaining synaptic dynamics. Alteration of protein expression levels in this matrix is a marked phenomenon of neuropsychiatric disorders including schizophrenia, where cognitive functions are impaired. To investigate the genetic relationship of genes expressed in the PSD with schizophrenia, a family-based association analysis of genetic variants in PSD genes such as DLG4, DLG1, PICK1 and MDM2, was performed, using Japanese samples (124 pedigrees, n = 376 subjects). Results showed a significant association of the rs17203281 variant from the DLG4 gene, with preferential transmission of the C allele (p = 0.02), although significance disappeared after correction for multiple testing. Replication analysis of this variant, found no association in a Chinese schizophrenia cohort (293 pedigrees, n = 1163 subjects) or in a Japanese case-control sample (n = 4182 subjects). The DLG4 expression levels between postmortem brain samples from schizophrenia patients showed no significant changes from controls. Interestingly, a five marker haplotype in DLG4, involving rs2242449, rs17203281, rs390200, rs222853 and rs222837, was enriched in a population specific manner, where the sequences A-C-C-C-A and G-C-C-C-A accumulated in Japanese (p = 0.0009) and Chinese (p = 0.0007) schizophrenia pedigree samples, respectively. However, this could not be replicated in case-control samples. None of the variants in other examined candidate genes showed any significant association in these samples. The current study highlights a putative role for DLG4 in schizophrenia pathogenesis, evidenced by haplotype association, and warrants further dense screening for variants within these haplotypes.     

Fitness Effects of Single Amino Acid Insertions and Deletions in TEM-1 β-Lactamase

Short insertions and deletions (InDels) are a common type of mutation found in nature and a useful source of variation in protein engineering. InDel events have important consequences in protein evolution, often opening new pathways for adaptation. However, much less is known about the effects of InDels compared to point mutations and amino acid substitutions. In particular, deep mutagenesis studies on the distribution of fitness effects of mutations have focused almost exclusively on amino acid substitutions. Here, we present a near-comprehensive analysis of the fitness effects of single amino acid InDels in TEM-1 β-lactamase. While we found InDels to be largely deleterious, partially overlapping deletion-tolerant and insertion-tolerant regions were observed throughout the protein, especially in unstructured regions and at the end of helices. The signal sequence of TEM-1 tolerated InDels more than the mature protein. Most regions of the protein tolerated insertions more than deletions, but a few regions tolerated deletions more than insertions. We examined the relationship between InDel tolerance and a variety of measures to help understand its origin. These measures included evolutionary variation in β-lactamases, secondary structure identity, tolerance to amino acid substitutions, solvent accessibility, and side-chain weighted contact number. We found secondary structure, weighted contact number, and evolutionary variation in class A beta-lactamases to be the somewhat predictive of InDel fitness effects.     

Effect of High- versus Low-Intensity Supervised Aerobic and Resistance Training on Modifiable Cardiovascular Risk Factors in Type 2 Diabetes; The Italian Diabetes and Exercise Study (IDES)

Background

While current recommendations on exercise type and volume have strong experimental bases, there is no clear evidence from large-sized studies indicating whether increasing training intensity provides additional benefits to subjects with type 2 diabetes.

Objective

To compare the effects of moderate-to-high intensity (HI) versus low-to-moderate intensity (LI) training of equal energy cost, i.e. exercise volume, on modifiable cardiovascular risk factors.

Design

Pre-specified sub-analysis of the Italian Diabetes and Exercise Study (IDES), a randomized multicenter prospective trial comparing a supervised exercise intervention with standard care for 12 months (2005–2006).

Setting

Twenty-two outpatient diabetes clinics across Italy.

Patients

Sedentary patients with type 2 diabetes assigned to twice-a-week supervised progressive aerobic and resistance training plus exercise counseling (n = 303).

Interventions

Subjects were randomized by center to LI (n = 142, 136 completed) or HI (n = 161, 152 completed) progressive aerobic and resistance training, i.e. at 55% or 70% of predicted maximal oxygen consumption and at 60% or 80% of predicted 1-Repetition Maximum, respectively, of equal volume.

Main Outcome Measure(s)

Hemoglobin (Hb) A_1c and other cardiovascular risk factors; 10-year coronary heart disease (CHD) risk scores.

Results

Volume of physical activity, both supervised and non-supervised, was similar in LI and HI participants. Compared with LI training, HI training produced only clinically marginal, though statistically significant, improvements in HbA_1c (mean difference −0.17% [95% confidence interval −0.44,0.10], P = 0.03), triglycerides (−0.12 mmol/l [−0.34,0.10], P = 0.02) and total cholesterol (−0.24 mmol/l [−0.46, −0.01], P = 0.04), but not in other risk factors and CHD risk scores. However, intensity was not an independent predictor of reduction of any of these parameters. Adverse event rate was similar in HI and LI subjects.

Conclusions

Data from the large IDES cohort indicate that, in low-fitness individuals such as sedentary subjects with type 2 diabetes, increasing exercise intensity is not harmful, but does not provide additional benefits on cardiovascular risk factors.

Trial Registration

www.ISRCTN.org ISRCTN-04252749.     

Gene Expression Profiling in the Type 1 Diabetes Rat Diaphragm

Background

Respiratory muscle contractile performance is impaired by diabetes, mechanisms of which included altered carbohydrate and lipid metabolism, oxidative stress and changes in membrane electrophysiology. The present study examined to what extent these cellular perturbations involve changes in gene expression.

Methodology/Principal Findings

Diaphragm muscle from streptozotocin-diabetic rats was analyzed with Affymetrix gene expression arrays. Diaphragm from diabetic rats had 105 genes with at least ±2-fold significantly changed expression (55 increased, 50 decreased), and these were assigned to gene ontology groups based on over-representation analysis using DAVID software. There was increased expression of genes involved in palmitoyl-CoA hydrolase activity (a component of lipid metabolism) (P = 0.037, n = 2 genes, fold change 4.2 to 27.5) and reduced expression of genes related to carbohydrate metabolism (P = 0.000061, n = 8 genes, fold change −2.0 to −8.5). Other gene ontology groups among upregulated genes were protein ubiquitination (P = 0.0053, n = 4, fold change 2.2 to 3.4), oxidoreductase activity (P = 0.024, n = 8, fold change 2.1 to 6.0), and morphogenesis (P = 0.012, n = 10, fold change 2.1 to 4.3). Other downregulated gene groups were extracellular region (including extracellular matrix and collagen) (P = 0.00032, n = 13, fold change −2.2 to −3.7) and organogenesis (P = 0.032, n = 7, fold change −2.1 to −3.7). Real-time PCR confirmed the directionality of changes in gene expression for 30 of 31 genes tested.

Conclusions/Significance

These data indicate that in diaphragm muscle type 1 diabetes increases expression of genes involved in lipid energetics, oxidative stress and protein ubiquitination, decreases expression of genes involved in carbohydrate metabolism, and has little effect on expression of ion channel genes. Reciprocal changes in expression of genes involved in carbohydrate and lipid metabolism may change the availability of energetic substrates and thereby directly modulate fatigue resistance, an important issue for a muscle like the diaphragm which needs to contract without rest for the entire lifetime of the organism.     

A Variant of the SLC10A2 Gene Encoding the Apical Sodium-Dependent Bile Acid Transporter Is a Risk Factor for Gallstone Disease

Background

Cholelithiasis is a multifactorial process and several mechanisms of gallstone formation have been postulated. As one of these mechanisms, a decreased expression of the ileal apical sodium-dependent bile acid transporter gene SLC10A2 in gallstone carriers was described previously. In this study the SLC10A2 gene was investigated to identify novel genetic variants and their association with gallstone formation.

Methodology/Principal Findings

Study subjects were selected with the presence or absence of gallstones confirmed by ultrasound and medical history. Genomic DNA was obtained from blood leukocytes. Sequence analysis was performed of all six exonic and flanking regions as well as of 2,400 base pairs of the SLC10A2 promoter in a cohort of gallstone carriers and control subjects from Stuttgart, Germany. Genotype frequencies of newly identified genetic variants (n = 6) and known single nucleotide polymorphisms (n = 24) were established using MALDI-TOF mass spectrometry. Six new genetic variants were found within the SLC10A2 gene. Although none of the variants was linked to gallstone disease in the Stuttgart cohort overall, the minor allele of SNP rs9514089 was more prevalent in male non-obese gallstone carriers (p = 0.06680, OR = 11.00). In a separate population from Aachen, Germany, the occurrence of rs9514089 was two-fold higher in gallstone patients (22%) than in corresponding controls (11%) (p = 0.00995, OR = 2.19). In the pooled Aachen/Stuttgart cohort rs9514089 was highly significantly linked to cholelithiasis (p = 0.00767, OR = 2.04). A more frequent occurrence was observed for male gallstone carriers (22%) compared to controls (9%) (p = 0.01017, OR = 2.99), for the total normal weight group (p = 0.00754, OR = 2.90), and for male non-obese gallstone patients (p = 0.01410, OR = 6.85). Moreover, for the minor allele of rs9514089 an association with low plasma cholesterol levels was found especially in gallstone carriers (p = 0.05).

Conclusions/Significance

We have identified SLC10A2 as a novel susceptibility gene for cholelithiasis in humans. Comprehensive statistical analysis provides strong evidence that rs9514089 is a genetic determinant especially in male non-obese gallstone carriers. The minor allele of rs9514089 is related to differences in plasma cholesterol levels among the subjects.     

The Unique hmuY Gene Sequence as a Specific Marker of Porphyromonas gingivalis

Porphyromonas gingivalis, a major etiological agent of chronic periodontitis, acquires heme from host hemoproteins using the HmuY hemophore. The aim of this study was to develop a specific P. gingivalis marker based on a hmuY gene sequence. Subgingival samples were collected from 66 patients with chronic periodontitis and 40 healthy subjects and the entire hmuY gene was analyzed in positive samples. Phylogenetic analyses demonstrated that both the amino acid sequence of the HmuY protein and the nucleotide sequence of the hmuY gene are unique among P. gingivalis strains/isolates and show low identity to sequences found in other species (below 50 and 56%, respectively). In agreement with these findings, a set of hmuY gene-based primers and standard/real-time PCR with SYBR Green chemistry allowed us to specifically detect P. gingivalis in patients with chronic periodontitis (77.3%) and healthy subjects (20%), the latter possessing lower number of P. gingivalis cells and total bacterial cells. Isolates from healthy subjects possess the hmuY gene-based nucleotide sequence pattern occurring in W83/W50/A7436 (n = 4), 381/ATCC 33277 (n = 3) or TDC60 (n = 1) strains, whereas those from patients typically have TDC60 (n = 21), W83/W50/A7436 (n = 17) and 381/ATCC 33277 (n = 13) strains. We observed a significant correlation between periodontal index of risk of infectiousness (PIRI) and the presence/absence of P. gingivalis (regardless of the hmuY gene-based sequence pattern of the isolate identified [r = 0.43; P = 0.0002] and considering particular isolate pattern [r = 0.38; P = 0.0012]). In conclusion, we demonstrated that the hmuY gene sequence or its fragments may be used as one of the molecular markers of P. gingivalis.     

IgG Antibodies to Cyclic Citrullinated Peptides Exhibit Profiles Specific in Terms of IgG Subclasses,Fc-Glycans and a Fab-Peptide Sequence

The Fc-glycan profile of IgG₁ anti-citrullinated peptide antibodies (ACPA) in rheumatoid arthritis (RA) patients has recently been reported to be different from non-ACPA IgG₁, a phenomenon which likely plays a role in RA pathogenesis. Herein we investigate the Fc-glycosylation pattern of all ACPA-IgG isotypes and simultaneously investigate in detail the IgG protein-chain sequence repertoire. IgG from serum or plasma (S/P, n = 14) and synovial fluid (SF, n = 4) from 18 ACPA-positive RA-patients was enriched using Protein G columns followed by ACPA-purification on cyclic citrullinated peptide-2 (CCP2)-coupled columns. Paired ACPA (anti-CCP2 eluted IgG) and IgG flow through (FT) fractions were analyzed by LC-MS/MS-proteomics. IgG peptides, isotypes and corresponding Fc-glycopeptides were quantified and interrogated using uni- and multivariate statistics. The Fc-glycans from the IgG₄ peptide EEQFNSTYR was validated using protein A column purification. Relative to FT-IgG₄, the ACPA-IgG₄ Fc-glycan-profile contained lower amounts (p = 0.002) of the agalacto and asialylated core-fucosylated biantennary form (FA2) and higher content (p = 0.001) of sialylated glycans. Novel differences in the Fc-glycan-profile of ACPA-IgG₁ compared to FT-IgG₁ were observed in the distribution of bisected forms (n = 5, p = 0.0001, decrease) and mono-antennnary forms (n = 3, p = 0.02, increase). Our study also confirmed higher abundance of FA2 (p = 0.002) and lower abundance of afucosylated forms (n = 4, p = 0.001) in ACPA-IgG₁ relative to FT-IgG₁ as well as lower content of IgG₂ (p = 0.0000001) and elevated content of IgG₄ (p = 0.004) in ACPA compared to FT. One λ-variable peptide sequence was significantly increased in ACPA (p = 0.0001). In conclusion, the Fc-glycan profile of both ACPA-IgG₁ and ACPA-IgG₄ are distinct. Given that IgG₁ and IgG₄ have different Fc-receptor and complement binding affinities, this phenomenon likely affects ACPA effector- and immune-regulatory functions in an IgG isotype-specific manner. These findings further highlight the importance of antibody characterization in relation to functional in vivo and in vitro studies.     

Human Papillomavirus 16 Non-European Variants Are Preferentially Associated with High-Grade Cervical Lesions

HPV16 accounts for 50–70% of cervical cancer cases worldwide. Characterization of HPV16 variants previously indicated that they differ in risks for viral persistence, progression to cervical precancer and malignant cancer. The aim of this study was to examine the association of severity of disease with HPV16 variants identified in specimens (n = 281) obtained from a Cervical Pathology and Colposcopy outpatient clinic in the University Hospital of Espírito Santo State, Southeastern Brazil, from April 2010 to November 2011. All cytologic and histologic diagnoses were determined prior to definitive treatment. The DNA was isolated using QIAamp DNA Mini Kit and HPV was detected by amplification with PGMY09/11 primers and positive samples were genotyped by RFLP analyses and reverse line blot. The genomes of the HPV16 positive samples were sequenced, from which variant lineages were determined. Chi² statistics was performed to test the association of HPV16 variants between case and control groups. The prevalence of HR-HPV types in <CIN1, CIN2 and CIN3+ were 33.7%, 84.4% and 91.6%, respectively. Thirty-eight of 49 (78%) HPV16 positive samples yielded HPV16 sequence information; of which, 32 complete genomes were sequenced and an additional 6 samples were partially sequenced. Phylogenetic analysis and patterns of variations identified 65.8% (n = 25) as HPV16 European (E) and 34.2% (n = 13) as non-European (NE) variants. Classification of disease into CIN3+ vs. <CIN3 indicated that NE types were associated with high-grade disease with an OR = 4.6 (1.07–20.2, p = 0.05). The association of HPV16 NE variants with an increased risk of CIN3+ is consistent with an HPV16 genetically determined enhanced oncogenicity. The prevalence of genetic variants of HPV16 is distributed across different geographical areas and with recent population admixture, only empiric data will provide information on the highest risk HPV16 variants within a given population.