Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Solving ambiguities in contig assembly of Idiomarina loihiensis L2TR chromosome by in silico analyses.

Nucleotide composition analyses of bacterial genomes such as cumulative GC skew highlight the atypical, strongly asymmetric architecture of the recently published chromosome of Idiomarina loihiensis L2TR, suggesting that an inversion of a 600-kb chromosomal segment occurred. The presence of 3.4-kb inverted repeated sequences at the borders of the putative rearrangement supports this hypothesis. Reverting in silico this segment restores (1) a symmetric chromosome architecture; (2) the co-orientation of transcription of all rRNA operons with DNA replication; and (3) a better conservation of gene order between this chromosome and other gamma-proteobacterial ones. Finally, long-range PCRs encompassing the ends of the 600-kb segment reveal the existence of the reverted configuration but not of the published one. This demonstrates how cumulative nucleotide-skew analyses can validate genome assemblies.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Predicting stroke through genetic risk functions: the CHARGE Risk Score Project.

BACKGROUND AND PURPOSE: Beyond the Framingham Stroke Risk Score, prediction of future stroke may improve with a genetic risk score (GRS) based on single-nucleotide polymorphisms associated with stroke and its risk factors. METHODS: The study includes 4 population-based cohorts with 2047 first incident strokes from 22,720 initially stroke-free European origin participants aged ≥55 years, who were followed for up to 20 years. GRSs were constructed with 324 single-nucleotide polymorphisms implicated in stroke and 9 risk factors. The association of the GRS to first incident stroke was tested using Cox regression; the GRS predictive properties were assessed with area under the curve statistics comparing the GRS with age and sex, Framingham Stroke Risk Score models, and reclassification statistics. These analyses were performed per cohort and in a meta-analysis of pooled data. Replication was sought in a case-control study of ischemic stroke. RESULTS: In the meta-analysis, adding the GRS to the Framingham Stroke Risk Score, age and sex model resulted in a significant improvement in discrimination (all stroke: Δjoint area under the curve=0.016, P=2.3×10(-6); ischemic stroke: Δjoint area under the curve=0.021, P=3.7×10(-7)), although the overall area under the curve remained low. In all the studies, there was a highly significantly improved net reclassification index (P<10(-4)). CONCLUSIONS: The single-nucleotide polymorphisms associated with stroke and its risk factors result only in a small improvement in prediction of future stroke compared with the classical epidemiological risk factors for stroke.

Integrating detection probabilities in species distribution models of amphibians breeding in Mediterranean temporary ponds

Aim  The imperfect detection of species may lead to erroneous conclusions about species-environment relationships. Accuracy in species detection usually requires temporal replication at sampling sites, a time-consuming and costly monitoring scheme. Here, we applied a lower-cost alternative based on a double-sampling approach to incorporate the reliability of species detection into regression-based species distribution modelling.Location  Doñana National Park (south-western Spain).Methods  Using species-specific monthly detection probabilities, we estimated the detection reliability as the probability of having detected the species given the species-specific survey time. Such reliability estimates were used to account explicitly for data uncertainty by weighting each absence. We illustrated how this novel framework can be used to evaluate four competing hypotheses as to what constitutes primary environmental control of amphibian distribution: breeding habitat, aestivating habitat, spatial distribution of surrounding habitats and/or major ecosystems zonation. The study was conducted on six pond-breeding amphibian species during a 4-year period.Results  Non-detections should not be considered equivalent to real absences, as their reliability varied considerably. The occurrence of Hyla meridionalis and Triturus pygmaeus was related to a particular major ecosystem of the study area, where suitable habitat for these species seemed to be widely available. Characteristics of the breeding habitat (area and hydroperiod) were of high importance for the occurrence of Pelobates cultripes and Pleurodeles waltl. Terrestrial characteristics were the most important predictors of the occurrence of Discoglossus galganoi and Lissotriton boscai, along with spatial distribution of breeding habitats for the last species.Main conclusions  We did not find a single best supported hypothesis valid for all species, which stresses the importance of multiscale and multifactor approaches. More importantly, this study shows that estimating the reliability of non-detection records, an exercise that had been previously seen as a naïve goal in species distribution modelling, is feasible and could be promoted in future studies, at least in comparable systems.

Association of the hemochromatosis gene with pazopanib-induced transaminase elevation in renal cell carcinoma.

BACKGROUND & AIMS: Pazopanib has demonstrated clinical benefit in patients with advanced renal cell carcinoma (RCC) and is generally well tolerated. However, transaminase elevations have commonly been observed. This 2-stage study sought to identify genetic determinants of alanine transaminase (ALT) elevations in pazopanib-treated white patients with RCC.¦METHODS: Data from two separate clinical studies were used to examine the association of genetic polymorphisms with maximum on-treatment ALT levels.¦RESULTS: Of 6852 polymorphisms in 282 candidate genes examined in an exploratory dataset of 115 patients, 92 polymorphisms in 40 genes were significantly associated with ALT elevation (p<0.01). Two markers (rs2858996 and rs707889) in the HFE gene, which are not yet known to be associated with hemochromatosis, showed evidence for replication. Because of multiple comparisons, there was a 12% likelihood the replication occurred by chance. These two markers demonstrated strong linkage disequilibrium (r(2)=0.99). In the combined dataset, median (25-75th percentile) maximum ALT values were 1.2 (0.7-1.9), 1.1 (0.8-2.5), and 5.4 (1.9-7.6)×ULN for rs2858996 GG (n=148), GT (n=82), and TT (n=1 2) genotypes, respectively. All 12 TT patients had a maximum ALT>ULN, and 8 (67%) had ALT≥3×ULN. The odds ratio (95% CI) for ALT≥3×ULN for TT genotype was 39.7 (2.2-703.7) compared with other genotypes. As a predictor of ALT≥3×ULN, the TT genotype had a negative predictive value of 0.83 and positive predictive value of 0.67. No TT patients developed liver failure.¦CONCLUSIONS: The rs2858996/rs707889 polymorphisms in the HFE gene may be associated with reversible ALT elevation in pazo-panib-treated patients with RCC.

Large TCR diversity of virus-specific CD8 T cells provides the mechanistic basis for massive TCR renewal after antigen exposure.

Ex vivo analysis of virus-specific CD8 T cell populations by anchored PCR has shown that the CD8 TCR repertoire was less oligoclonal (seven to nine clonotypes per individual epitope) than previously thought. In the current study, TCR diversity was investigated by assessing both the overall TCR β-chain variable regions usage as well as the CDR3 regions in ex vivo-isolated CMV- and EBV-specific CD8 T cells from 27 healthy donors. The average number of clonotypes specific to most single viral epitopes comprised between 14 and 77. Changes in the CD8 TCR repertoire were also longitudinally assessed under conditions of HIV-1 chronic infection (i.e., in patients with suppressed virus replication and after treatment interruption and Ag re-exposure). The results showed that a large renewal (≤80%) of the TRB repertoire occurred after Ag re-exposure and was eventually associated with an increased T cell recognition functional avidity. These results demonstrate that the global CD8 TCR repertoire is much more diverse (≤9-fold) than previously estimated and provide the mechanistic basis for supporting massive repertoire renewal during chronic virus infection and Ag re-exposure.

Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis

Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10(-8)) loci, some including known iron-related genes (HFE, SLC40A1, TF, TFR2, TFRC, TMPRSS6) and others novel (ABO, ARNTL, FADS2, NAT2, TEX14). SNPs at ARNTL, TF, and TFR2 affect iron markers in HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease.

Genome-wide association study to identify common variants associated with brachial circumference: a meta-analysis of 14 cohorts.

Brachial circumference (BC), also known as upper arm or mid arm circumference, can be used as an indicator of muscle mass and fat tissue, which are distributed differently in men and women. Analysis of anthropometric measures of peripheral fat distribution such as BC could help in understanding the complex pathophysiology behind overweight and obesity. The purpose of this study is to identify genetic variants associated with BC through a large-scale genome-wide association scan (GWAS) meta-analysis. We used fixed-effects meta-analysis to synthesise summary results across 14 GWAS discovery and 4 replication cohorts comprising overall 22,376 individuals (12,031 women and 10,345 men) of European ancestry. Individual analyses were carried out for men, women, and combined across sexes using linear regression and an additive genetic model: adjusted for age and adjusted for age and BMI. We prioritised signals for follow-up in two-stages. We did not detect any signals reaching genome-wide significance. The FTO rs9939609 SNP showed nominal evidence for association (p<0.05) in the age-adjusted strata for men and across both sexes. In this first GWAS meta-analysis for BC to date, we have not identified any genome-wide significant signals and do not observe robust association of previously established obesity loci with BC. Large-scale collaborations will be necessary to achieve higher power to detect loci underlying BC.

Type I interferons protect T cells against NK cell attack mediated by the activating receptor NCR1.

Direct type I interferon (IFN) signaling on T cells is necessary for the proper expansion, differentiation, and survival of responding T cells following infection with viruses prominently inducing type I IFN. The reasons for the abortive response of T cells lacking the type I IFN receptor (Ifnar1(-/-)) remain unclear. We report here that Ifnar1(-/-) T cells were highly susceptible to natural killer (NK) cell-mediated killing in a perforin-dependent manner. Depletion of NK cells prior to lymphocytic choriomeningitis virus (LCMV) infection completely restored the early expansion of Ifnar1(-/-) T cells. Ifnar1(-/-) T cells had elevated expression of natural cytotoxicity triggering receptor 1 (NCR1) ligands upon infection, rendering them targets for NCR1 mediated NK cell attack. Thus, direct sensing of type I IFNs by T cells protects them from NK cell killing by regulating the expression of NCR1 ligands, thereby revealing a mechanism by which T cells can evade the potent cytotoxic activity of NK cells.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

ITMIG consensus statement on the use of the WHO histological classification of thymoma and thymic carcinoma: refined definitions, histological criteria, and reporting.

INTRODUCTION: The 2004 version of the World Health Organization classification subdivides thymic epithelial tumors into A, AB, B1, B2, and B3 (and rare other) thymomas and thymic carcinomas (TC). Due to a morphological continuum between some thymoma subtypes and some morphological overlap between thymomas and TC, a variable proportion of cases may pose problems in classification, contributing to the poor interobserver reproducibility in some studies. METHODS: To overcome this problem, hematoxylin-eosin-stained and immunohistochemically processed sections of prototypic, "borderland," and "combined" thymomas and TC (n = 72) were studied by 18 pathologists at an international consensus slide workshop supported by the International Thymic Malignancy Interest Group. RESULTS: Consensus was achieved on refined criteria for decision making at the A/AB borderland, the distinction between B1, B2, and B3 thymomas and the separation of B3 thymomas from TCs. "Atypical type A thymoma" is tentatively proposed as a new type A thymoma variant. New reporting strategies for tumors with more than one histological pattern are proposed. CONCLUSION: These guidelines can set the stage for reproducibility studies and the design of a clinically meaningful grading system for thymic epithelial tumors.

Multiplex liquid chromatography-tandem mass spectrometry assay for simultaneous therapeutic drug monitoring of ribavirin, boceprevir, and telaprevir.

New directly acting antivirals (DAAs) that inhibit hepatitis C virus (HCV) replication are increasingly used for the treatment of chronic hepatitis C. A marked pharmacokinetic variability and a high potential for drug-drug interactions between DAAs and numerous drug classes have been identified. In addition, ribavirin (RBV), commonly associated with hemolytic anemia, often requires dose adjustment, advocating for therapeutic drug monitoring (TDM) in patients under combined antiviral therapy. However, an assay for the simultaneous analysis of RBV and DAAs constitutes an analytical challenge because of the large differences in polarity among these drugs, ranging from hydrophilic (RBV) to highly lipophilic (telaprevir [TVR]). Moreover, TVR is characterized by erratic behavior on standard octadecyl-based reversed-phase column chromatography and must be separated from VRT-127394, its inactive C-21 epimer metabolite. We have developed a convenient assay employing simple plasma protein precipitation, followed by high-performance liquid chromatography coupled to tandem mass spectrometry (HPLC-MS/MS) for the simultaneous determination of levels of RBV, boceprevir, and TVR, as well as its metabolite VRT-127394, in plasma. This new, simple, rapid, and robust HPLC-MS/MS assay offers an efficient method of real-time TDM aimed at maximizing efficacy while minimizing the toxicity of antiviral therapy.

Genetic loci for retinal arteriolar microcirculation.

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10(-8). This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10(-12) in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.

Backseat drivers: the hidden influence of microbial viruses on disease.

Because viral replication depends on the vigour of its host, many viruses have evolved incentives of fitness to pay their keep. When the viral host is a human pathogen, these fitness factors can surface as virulence: creating a Russian doll of pathogenesis where pathogens within pathogens complicate the disease process. Microbial viruses can even be independently immunogenic, as we recently reported for leishmania-virus. Thus, the incidence of this 'hyperpathogenism' is becoming an important clinical consideration and by appreciating the microbial-virus as a backseat driver of human disease, we could exploit its presence as a diagnostic biomarker and molecular target for therapeutic intervention. Here we discuss the prevalence of clinically relevant hyperpathogenism as well as the environmental sanctuaries that breed it.