Genetic Analyses Reveal a Role for Vitamin D Insufficiency in HCV-Associated Hepatocellular Carcinoma Development.

BACKGROUND: Vitamin D insufficiency has been associated with the occurrence of various types of cancer, but causal relationships remain elusive. We therefore aimed to determine the relationship between genetic determinants of vitamin D serum levels and the risk of developing hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC). METHODOLOGYPRINCIPAL FINDINGS: Associations between CYP2R1, GC, and DHCR7 genotypes that are determinants of reduced 25-hydroxyvitamin D (25[OH]D3) serum levels and the risk of HCV-related HCC development were investigated for 1279 chronic hepatitis C patients with HCC and 4325 without HCC, respectively. The well-known associations between CYP2R1 (rs1993116, rs10741657), GC (rs2282679), and DHCR7 (rs7944926, rs12785878) genotypes and 25(OH)D3 serum levels were also apparent in patients with chronic hepatitis C. The same genotypes of these single nucleotide polymorphisms (SNPs) that are associated with reduced 25(OH)D3 serum levels were found to be associated with HCV-related HCC (P = 0.07 [OR = 1.13, 95% CI = 0.99-1.28] for CYP2R1, P = 0.007 [OR = 1.56, 95% CI = 1.12-2.15] for GC, P = 0.003 [OR = 1.42, 95% CI = 1.13-1.78] for DHCR7; ORs for risk genotypes). In contrast, no association between these genetic variations and liver fibrosis progression rate (P>0.2 for each SNP) or outcome of standard therapy with pegylated interferon-α and ribavirin (P>0.2 for each SNP) was observed, suggesting a specific influence of the genetic determinants of 25(OH)D3 serum levels on hepatocarcinogenesis. CONCLUSIONSSIGNIFICANCE: Our data suggest a relatively weak but functionally relevant role for vitamin D in the prevention of HCV-related hepatocarcinogenesis.

Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Discovery of amino acid variants in the human glucose-dependent insulinotropic polypeptide (GIP) receptor: the impact on the pancreatic beta cell responses and functional expression studies in Chinese hamster fibroblast cells.

The two incretins, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1), are insulinotropic factors released from the small intestine to the blood stream in response to oral glucose ingestion. The insulinotropic effect of GLP-1 is maintained in patients with Type II (non-insulin-dependent) diabetes mellitus, whereas, for unknown reasons, the effect of GIP is diminished or lacking. We defined the exon-intron boundaries of the human GIP receptor, made a mutational analysis of the gene and identified two amino acid substitutions, A207 V and E354Q. In an association study of 227 Caucasian Type II diabetic patients and 224 matched glucose tolerant control subjects, the allelic frequency of the A207 V polymorphism was 1.1% in Type II diabetic patients and 0.7% in control subjects (p = 0.48), whereas the allelic frequency of the codon 354 polymorphism was 24.9% in Type II diabetic patients versus 23.2% in control subjects. Interestingly, the glucose tolerant subjects (6% of the population) who were homozygous for the codon 354 variant had on average a 14% decrease in fasting serum C-peptide concentration (p = 0.01) and an 11% decrease in the same variable 30 min after an oral glucose load (p = 0.03) compared with subjects with the wild-type receptor. Investigation of the function of the two GIP receptor variants in Chinese hamster fibroblasts showed, however, that the GIP-induced cAMP formation and the binding of GIP to cells expressing the variant receptors were not different from the findings in cells expressing the wildtype GIP receptor. In conclusion, amino acid variants in the GIP receptor are not associated with random Type II diabetes in patients of Danish Caucasian origin or with altered GIP binding and GIP-induced cAMP production when stably transfected in Chinese hamster fibroblasts. The finding of an association between homozygosity for the codon 354 variant and reduced fasting and post oral glucose tolerance test (OGTT) serum C-peptide concentrations, however, calls for further investigations and could suggest that GIP even in the fasting state regulates the beta-cell secretory response.

BlastR--fast and accurate database searches for non-coding RNAs.

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.

Current genetic data do not improve the prediction of type 2 diabetes mellitus: the CoLaus study.

CONTEXT: Several genetic risk scores to identify asymptomatic subjects at high risk of developing type 2 diabetes mellitus (T2DM) have been proposed, but it is unclear whether they add extra information to risk scores based on clinical and biological data. OBJECTIVE: The objective of the study was to assess the extra clinical value of genetic risk scores in predicting the occurrence of T2DM. DESIGN: This was a prospective study, with a mean follow-up time of 5 yr. SETTING AND SUBJECTS: The study included 2824 nondiabetic participants (1548 women, 52 ± 10 yr). MAIN OUTCOME MEASURE: Six genetic risk scores for T2DM were tested. Four were derived from the literature and two were created combining all (n = 24) or shared (n = 9) single-nucleotide polymorphisms of the previous scores. A previously validated clinic + biological risk score for T2DM was used as reference. RESULTS: Two hundred seven participants (7.3%) developed T2DM during follow-up. On bivariate analysis, no differences were found for all but one genetic score between nondiabetic and diabetic participants. After adjusting for the validated clinic + biological risk score, none of the genetic scores improved discrimination, as assessed by changes in the area under the receiver-operating characteristic curve (range -0.4 to -0.1%), sensitivity (-2.9 to -1.0%), specificity (0.0-0.1%), and positive (-6.6 to +0.7%) and negative (-0.2 to 0.0%) predictive values. Similarly, no improvement in T2DM risk prediction was found: net reclassification index ranging from -5.3 to -1.6% and nonsignificant (P ≥ 0.49) integrated discrimination improvement. CONCLUSIONS: In this study, adding genetic information to a previously validated clinic + biological score does not seem to improve the prediction of T2DM.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

The Growing Importance of CNVs: New Insights for Detection and Clinical Interpretation.

Differences between genomes can be due to single nucleotide variants, translocations, inversions, and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 500 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease. Hence there is a need for better-tailored and more robust tools for the detection and genome-wide analyses of CNVs. While a link between a given CNV and a disease may have often been established, the relative CNV contribution to disease progression and impact on drug response is not necessarily understood. In this review we discuss the progress, challenges, and limitations that occur at different stages of CNV analysis from the detection (using DNA microarrays and next-generation sequencing) and identification of recurrent CNVs to the association with phenotypes. We emphasize the importance of germline CNVs and propose strategies to aid clinicians to better interpret structural variations and assess their clinical implications.

Decreased PM10 exposure attenuates age-related lung function decline: genetic variants in p53, p21, and CCND1 modify this effect.

BACKGROUND: Decreasing exposure to airborne particulates was previously associated with reduced age-related decline in lung function. However, whether the benefit from improved air quality depends on genetic background is not known. Recent evidence points to the involvement of the genes p53 and p21 and of the cell cycle control gene cyclin D1 (CCND1) in the response of bronchial cells to air pollution. OBJECTIVE: We determined in 4,326 participants of the Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults (SAPALDIA) whether four single-nucleotide polymorphisms in three genes [CCND1 (rs9344 [P242P], rs667515), p53 (rs1042522 [R72P]), and p21 (rs1801270 [S31R])] modified the previously observed attenuation of the decline in the forced expiratory flow between 25% and 75% of the forced vital capacity (FEF(25-75)) associated with improved air quality. METHODS: Subjects of the prospective population-based SAPALDIA cohort were assessed in 1991 and 2002 by spirometry, questionnaires, and biological sample collection for genotyping. We assigned spatially resolved concentrations of particulate matter with aerodynamic diameter < or = 10 microm (PM(10)) to each participant's residential history 12 months before the baseline and follow-up assessments. RESULTS: The effect of diminishing PM(10) exposure on FEF(25-75) decline appeared to be modified by p53 R72P, CCND1 P242P, and CCND1 rs667515. For example, a 10-microg/m(3) decline in average PM(10) exposure over an 11-year period attenuated the average annual decline in FEF(25-75) by 21.33 mL/year (95% confidence interval, 10.57-32.08) among participants homozygous for the CCND1 (P242P) GG genotype, by 13.72 mL/year (5.38-22.06) among GA genotypes, and by 6.00 mL/year (-4.54 to 16.54) among AA genotypes. CONCLUSIONS: Our results suggest that cell cycle control genes may modify the degree to which improved air quality may benefit respiratory function in adults.

A large study of androgen receptor germline variants and their relation to sex hormone levels and prostate cancer risk. Results from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium

Background: Androgens are key regulators of prostate gland maintenance and prostate cancer growth, and androgen deprivation therapy has been the mainstay of treatment for advanced prostate cancer for many years. A long-standing hypothesis has been that inherited variation in the androgen receptor (AR) gene plays a role in prostate cancer initiation. However, studies to date have been inconclusive and often suffered from small sample sizes. Objective and Methods: We investigated the association of AR sequence variants with circulating sex hormone levels and prostate cancer risk in 6058 prostate cancer cases and 6725 controls of Caucasian origin within the Breast and Prostate Cancer Cohort Consortium. We genotyped a highly polymorphic CAG microsatellite in exon 1 and six haplotype tagging single nucleotide polymorphisms and tested each genetic variant for association with prostate cancer risk and with sex steroid levels. Results: We observed no association between AR genetic variants and prostate cancer risk. However, there was a strong association between longer CAG repeats and higher levels of testosterone (P = 4.73 × 10−5) and estradiol (P = 0.0002), although the amount of variance explained was small (0.4 and 0.7%, respectively). Conclusions: This study is the largest to date investigating AR sequence variants, sex steroid levels, and prostate cancer risk. Although we observed no association between AR sequence variants and prostate cancer risk, our results support earlier findings of a relation between the number of CAG repeats and circulating levels of testosterone and estradiol.

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels.

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9)

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95% confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

High prevalence and low E6 genetic variability of human papillomavirus 58 in women with cervical cancer and precursor lesions in Southeast Mexico

Infection with some genotypes of human papillomavirus (HPV) is the most important risk factor associated with cervical cancer (CC). Throughout the world, HPV type 58 prevalence varies from one region to another; it is higher in women from certain countries in Asia and Latin America, such as China and Mexico. Although intratypic variants have been reported on a few occasions, our knowledge about HPV 58 genetic variation remains limited. Therefore, this work aims to (i) determine the prevalence of HPV type 58 amongst Mexican women with invasive CC or precursor lesions and (ii) identify HPV 58 sequence variants. One hundred and forty five colposcopy clinic patients were studied. Genotyping of HPV 16, 18 and 58 was determined by specific nested PCR and HPV 58 variants were detected by direct sequencing. The general prevalence of HPV was 51.7% (75/145). HPV 16 was found in 30.6% (23/75) and HPV 58 in 24% (18/75) of the patients. HPV 18 was not identified in patients with cervical intraepithelial neoplasia (CIN) grade I; it was only found in those with CIN II, with a prevalence of 6.8% (3/44). In patients with CC, the prevalence of HPV 16 and 58 was 78.9%. Regarding HPV 58 variants, 94.4% of the HPV 58 sequences were identical to the prototype strain, whereas one sample showed changes at a single nucleotide. This study demonstrates a high prevalence of HPV 58 and a low genetic variability of E6 sequences amongst Mexican colposcopy patients.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Allelic variation in the VMD2 gene in best disease and age-related macular degeneration.

PURPOSE: To assess the allelic variation of the VMD2 gene in patients with Best disease and age-related macular degeneration (AMD). METHODS: Three hundred twenty-one AMD patients, 192 ethnically similar control subjects, 39 unrelated probands with familial Best disease, and 57 unrelated probands with the ophthalmoscopic findings of Best disease but no family history were screened for sequence variations in the VMD2 gene by single-strand conformation polymorphism (SSCP) analysis. Amplimers showing a bandshift were reamplified and sequenced bidirectionally. In addition, the coding regions of the VMD2 gene were completely sequenced in six probands with familial Best disease who showed no SSCP shift. RESULTS: Forty different probable or possible disease-causing mutations were found in one or more Best disease or AMD patients. Twenty-nine of these variations are novel. Of the 39 probands with familial Best disease, mutations were detected in all 39 (33 by SSCP and 6 by DNA sequencing). SSCP screening of the 57 probands with a clinical diagnosis of Best disease but no family history revealed 16 with mutations. Mutations were found in 5 of 321 AMD patients (1.5%), a fraction that was not significantly greater than in control individuals (0/192, 0%). CONCLUSIONS: Patients with the clinical diagnosis of Best disease are significantly more likely to have a mutation in the VMD2 gene if they also have a positive family history. These findings suggest that a small fraction of patients with the clinical diagnosis of AMD may actually have a late-onset variant of Best disease, whereas at the same time, a considerable fraction of isolated patients with the ophthalmoscopic features of Best disease are probably affected with some other macular disease.

Change in mutation patterns of Plasmodium vivax dihydrofolate reductase (Pvdhfr) and dihydropteroate synthase (Pvdhps) in P. vivax isolates from malaria endemic areas of Thailand

Malaria is the most important public health problem in several countries. In Thailand, co-infections of Plasmodium vivax and Plasmodium falciparum are common. We examined the prevalence and patterns of mutations in P. vivax dihydrofolate reductase (Pvdhfr) and P. vivax dihydropteroate synthase (Pvdhps) in 103 blood samples collected from patients with P. vivax infection who had attended the malaria clinic in Mae Sot, Tak Province during 2009 and 2010. Using nested polymerase chain reaction-restriction fragment length polymorfism, we examined single nucleotide polymorphisms-haplotypes at amino acid positions 13, 33, 57, 58, 61, 117 and 173 of Pvdhfr and 383 and 553 of Pvdhps. All parasite isolates carried mutant Pvdhfr alleles, of which the most common alleles were triple mutants (99%). Eight different types of Pvdhfr and combination alleles were found, as follows: 57I/58R/117T, 57I/58R/117T, 57I/58R/117T/N, 57L/58R/117T, 57L/58R/117T, 58R/61M/117N, 58R/61M/117N and 13L/57L/58R/117T. The most common Pvdhfr alleles were 57I/58R/117T (77.7%), 57I/58R/117T/N (1%), 57L/58R/117T (5.8%) and 58R/61M/117N (14.5%). The most common Pvdhfr alleles were 57I/58R/117T (77.7%), 57I/58R/117T/N (1%), 57L/58R/117T (5.8%) and 58R/61M/117N (14.5%). Additionally, we recovered one isolate of a carrying a quadruple mutant allele, 13L/57L/58R/117T. The most prevalent Pvdhps allele was a single mutation in amino acid 383 (82.5%), followed by the wild-type A383/A553 (17.5%) allele. Results suggest that all P. vivax isolates in Thailand carry some combination of mutations in Pvdhfr and Pvdhps. Our findings demonstrate that development of new antifolate drugs effective against sulfadoxine-pyrimethamine-resistant P. vivax is required.