EGFR Gene Variants Are Associated with Specific Somatic Aberrations in Glioma

A number of gene variants have been associated with an increased risk of developing glioma. We hypothesized that the reported risk variants may be associated with tumor genomic instability. To explore potential correlations between germline risk variants and somatic genetic events, we analyzed matched tumor and blood samples from 95 glioma patients by means of SNP genotyping. The generated genotype data was used to calculate genome-wide allele-specific copy number profiles of the tumor samples. We compared the copy number profiles across samples and found two EGFR gene variants (rs17172430 and rs11979158) that were associated with homozygous deletion at the CDKN2A/B locus. One of the EGFR variants (rs17172430) was also associated with loss of heterozygosity at the EGFR locus. Our findings were confirmed in a separate dataset consisting of matched blood and tumor samples from 300 glioblastoma patients, compiled from publically available TCGA data. These results imply there is a functional effect of germline EGFR variants on tumor progression.     

Variants in the ASB10 Gene Are Associated with Primary Open Angle Glaucoma

Background

Recently nonsynonymous coding variants in the ankyrin repeats and suppressor of cytokine signaling box-containing protein 10 (ASB10) gene were found to be associated with primary open angle glaucoma (POAG) in cohorts from Oregon and Germany, but this finding was not confirmed in an independent cohort from Iowa. The aim of the current study was to assess the role of ASB10 gene variants in Pakistani glaucoma patients.

Methods

Sanger sequencing of the coding exons and splice junctions of the ASB10 gene was performed in 30 probands of multiplex POAG families, 208 sporadic POAG patients and 151 healthy controls from Pakistan. Genotypic associations of individual variants with POAG were analyzed with the Fisher’s exact or Chi-square test.

Results

In total 24 variants were identified in POAG probands and sporadic patients, including 11 novel variants and 13 known variants. 13 of the variants were nonsynonymous, 6 were synonymous, and 5 were intronic. Three nonsynonymous variants (p.Arg49Cys, p.Arg237Gly, p.Arg453Cys) identified in the probands were not segregating in the respective families. This is not surprising since glaucoma is a multifactorial disease, and multiple factors are likely to be involved in the disease manifestation in these families. However a nonsynonymous variant, p.Arg453Cys (rs3800791), was found in 6 sporadic POAG patients but not in controls, suggesting that it infers increased risk for the disease. In addition, one synonymous variant was found to be associated with sporadic POAG: p.Ala290Ala and the association of the variant with POAG remained significant after correction for multiple testing (uncorrected p-value 0.002, corrected p-value 0.047). The cumulative burden of rare, nonsynonymous variants was significantly higher in sporadic POAG patients compared to control individuals (p-value 0.000006).

Conclusions

Variants in ASB10 were found to be significantly associated with sporadic POAG in the Pakistani population. This supports previous findings that sequence variants in the ASB10 gene may act as a risk factor for glaucoma.     

Variants in Neuropeptide Y Receptor 1 and 5 Are Associated with Nutrient-Specific Food Intake and Are Under Recent Selection in Europeans

There is a large variation in caloric intake and macronutrient preference between individuals and between ethnic groups, and these food intake patterns show a strong heritability. The transition to new food sources during the agriculture revolution around 11,000 years ago probably created selective pressure and shaped the genome of modern humans. One major player in energy homeostasis is the appetite-stimulating hormone neuropeptide Y, in which the stimulatory capacity may be mediated by the neuropeptide Y receptors 1, 2 and 5 (NPY1R, NPY2R and NPY5R). We assess association between variants in the NPY1R, NPY2R and NPY5R genes and nutrient intake in a cross-sectional, single-center study of 400 men aged 40 to 80 years, and we examine whether genomic regions containing these genes show signatures of recent selection in 270 HapMap individuals (90 Africans, 90 Asians, and 90 Caucasians) and in 846 Dutch bloodbank controls. Our results show that derived alleles in NPY1R and NPY5R are associated with lower carbohydrate intake, mainly because of a lower consumption of mono- and disaccharides. We also show that carriers of these derived alleles, on average, consume meals with a lower glycemic index and glycemic load and have higher alcohol consumption. One of these variants shows the hallmark of recent selection in Europe. Our data suggest that lower carbohydrate intake, consuming meals with a low glycemic index and glycemic load, and/or higher alcohol consumption, gave a survival advantage in Europeans since the agricultural revolution. This advantage could lie in overall health benefits, because lower carbohydrate intake, consuming meals with a low GI and GL, and/or higher alcohol consumption, are known to be associated with a lower risk of chronic diseases.     

Common Genetic Variants in FOXP2 Are Not Associated with Individual Differences in Language Development

Much of our current knowledge regarding the association of FOXP2 with speech and language development comes from singleton and small family studies where a small number of rare variants have been identified. However, neither genome-wide nor gene-specific studies have provided evidence that common polymorphisms in the gene contribute to individual differences in language development in the general population. One explanation for this inconsistency is that previous studies have been limited to relatively small samples of individuals with low language abilities, using low density gene coverage. The current study examined the association between common variants in FOXP2 and a quantitative measure of language ability in a population-based cohort of European decent (n = 812). No significant associations were found for a panel of 13 SNPs that covered the coding region of FOXP2 and extended into the promoter region. Power analyses indicated we should have been able to detect a QTL variance of 0.02 for an associated allele with MAF of 0.2 or greater with 80% power. This suggests that, if a common variant associated with language ability in this gene does exist, it is likely of small effect. Our findings lead us to conclude that while genetic variants in FOXP2 may be significant for rare forms of language impairment, they do not contribute appreciably to individual variation in the normal range as found in the general population.     

Variants in the Dopamine-4-Receptor Gene Promoter Are Not Associated with Sensation Seeking in Skiers

Sensation seeking is a personality trait that has been associated with disinhibited behaviours including substance use and gambling, but also with high-risk sport practices including skydiving, paragliding, and downhill skiing. Twin studies have shown that sensation seeking is moderately heritable, and candidate genes encoding components involved in dopaminergic transmission have been investigated as contributing to this type of behaviour. To determine whether variants in the regulatory regions of the dopamine-4-receptor gene (DRD4) influenced sport-specific sensation seeking, we analyzed five polymorphisms (−1106T/C, −906T/C, −809G/A, −291C/T, 120-bp duplication) in the promoter region of the gene in a cohort of skiers and snowboarders (n = 599) that represented a broad range of sensation seeking behaviours. We grouped subjects by genotype at each of the five loci and compared impulsive sensation seeking and domain-specific (skiing) sensation seeking between groups. There were no significant associations between genotype(s) and general or domain-specific sensation seeking in the skiers and snowboarders, suggesting that while DRD4 has previously been implicated in sensation seeking, the promoter variants investigated in this study do not contribute to sensation seeking in this athlete population.     

Common Genetic Variants of Surfactant Protein-D (SP-D) Are Associated with Type 2 Diabetes

Context

Surfactant protein-D (SP-D) is a primordial component of the innate immune system intrinsically linked to metabolic pathways. We aimed to study the association of single nucleotide polymorphisms (SNPs) affecting SP-D with insulin resistance and type 2 diabetes (T2D).

Research Design and Methods

We evaluated a common genetic variant located in the SP-D coding region (rs721917, Met³¹Thr) in a sample of T2D patients and non-diabetic controls (n = 2,711). In a subset of subjects (n = 1,062), this SNP was analyzed in association with circulating SP-D concentrations, insulin resistance, and T2D. This SNP and others were also screened in the publicly available Genome Wide Association (GWA) database of the Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC).

Results

We found the significant association of rs721917 with circulating SP-D, parameters of insulin resistance and T2D. Indeed, G carriers showed decreased circulating SP-D (p = 0.004), decreased fasting glucose (p = 0.0002), glycated hemoglobin (p = 0.0005), and 33% (p = 0.002) lower prevalence of T2D, estimated under a dominant model, especially among women. Interestingly, these differences remained significant after controlling for origin, age, gender, and circulating SP-D. Moreover, this SNP and others within the SP-D genomic region (i.e. rs10887344) were significantly associated with quantitative measures of glucose homeostasis, insulin sensitivity, and T2D, according to GWAS datasets from MAGIC.

Conclusions

SP-D gene polymorphisms are associated with insulin resistance and T2D. These associations are independent of circulating SP-D concentrations.     

Common Variants in Left/Right Asymmetry Genes and Pathways Are Associated with Relative Hand Skill

Humans display structural and functional asymmetries in brain organization, strikingly with respect to language and handedness. The molecular basis of these asymmetries is unknown. We report a genome-wide association study meta-analysis for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)] (n = 728). The most strongly associated variant, rs7182874 (P = 8.68×10⁻⁹), is located in PCSK6, further supporting an association we previously reported. We also confirmed the specificity of this association in individuals with RD; the same locus was not associated with relative hand skill in a general population cohort (n = 2,666). As PCSK6 is known to regulate NODAL in the development of left/right (LR) asymmetry in mice, we developed a novel approach to GWAS pathway analysis, using gene-set enrichment to test for an over-representation of highly associated variants within the orthologs of genes whose disruption in mice yields LR asymmetry phenotypes. Four out of 15 LR asymmetry phenotypes showed an over-representation (FDR≤5%). We replicated three of these phenotypes; situs inversus, heterotaxia, and double outlet right ventricle, in the general population cohort (FDR≤5%). Our findings lead us to propose that handedness is a polygenic trait controlled in part by the molecular mechanisms that establish LR body asymmetry early in development.     

Genetic Variants in MUC4 Gene Are Associated with Lung Cancer Risk in a Chinese Population

Mucin MUC4, which is encoded by the MUC4 gene, plays an important role in epithelial cell proliferation and differentiation. Aberrant MUC4 overexpression is associated with invasive tumor proliferation and poor outcome in epithelial cancers. Collectively, the existing evidence suggests that MUC4 has tumor-promoter functions. In this study, we performed a case-control study of 1,048 incident lung cancer cases and 1,048 age- and sex frequency-matched cancer-free controls in a Chinese population to investigate the role of MUC4 gene polymorphism in lung cancer etiology. We identified nine SNPs that were significantly associated with increased lung cancer risk (P = 0.0425 for rs863582, 0.0333 for rs842226, 0.0294 for rs842225, 0.0010 for rs2550236, 0.0149 for rs2688515, 0.0191 for rs 2641773, 0.0058 for rs3096337, 0.0077 for rs859769, and 0.0059 for rs842461 in an additive model). Consistent with these single-locus analysis results, the haplotype analyses revealed an adverse effect of the haplotype “GGC” of rs3096337, rs859769, and rs842461 on lung cancer. Both the haplotype and diplotype “CTGAGC” of rs863582, rs842226, rs2550236, rs842225, and rs2688515 had an adverse effect on lung cancer, which is also consistent with the single-locus analysis. Moreover, we observed statistically significant interactions for rs863582 and rs842461 in heavy smokers. Our results suggest that MUC4 gene polymorphisms and their interaction with smoking may contribute to lung cancer etiology.     

Two Variants of the C-Reactive Protein Gene Are Associated with Risk of Pre-Eclampsia in an American Indian Population

Background

The etiology of pre-eclampsia (PE) is unknown; but it is accepted that normal pregnancy represents a distinctive challenge to the maternal immune system. C-reactive protein is a prominent component of the innate immune system; and we previously reported an association between PE and the CRP polymorphism, rs1205. Our aim was to explore the effects of additional CRP variants. The IBC (Cardiochip) genotyping microarray focuses on candidate genes and pathways related to the pathophysiology of cardiovascular disease.

Methods

This study recruited 140 cases of PE and 270 matched controls, of which 95 cases met criteria as severe PE, from an American Indian community. IBC array genotypes from 10 suitable CRP SNPs were analyzed. A replication sample of 178 cases and 427 controls of European ancestry was also genotyped.

Results

A nominally significant difference (p value <0.05) was seen in the distribution of discordant matched pairs for rs3093068; and Bonferroni corrected differences (P<0.005) were seen for rs876538, rs2794521, and rs3091244. Univariate conditional logistic regression odds ratios (OR) were nominally significant for rs3093068 and rs876538 models only. Multivariate logistic models with adjustment for mother''s age, nulliparity and BMI attenuated the effect (OR 1.58, P = 0.066, 95% CI 0.97–2.58) for rs876538 and (OR 2.59, P = 0.050, 95% CI 1.00–6.68) for rs3093068. An additive risk score of the above two risk genotypes shows a multivariate adjusted OR of 2.04 (P = 0.013, 95% CI 1.16–3.56). The replication sample also demonstrated significant association between PE and the rs876538 allele (OR = 1.55, P = 0.01, 95% CI 2.16–1.10). We also show putative functionality for the rs876538 and rs3093068 CRP variants.

Conclusion

The CRP variants, rs876538 and rs3093068, previously associated with other cardiovascular disease phenotypes, show suggestive association with PE in this American Indian population, further supporting a possible role for CRP in PE.     

Common Variants of KCNJ10 Are Associated with Susceptibility and Anti-Epileptic Drug Resistance in Chinese Genetic Generalized Epilepsies

To explore genetic mechanism of genetic generalized epilepsies (GGEs) is challenging because of their complex heritance pattern and genetic heterogeneity. KCNJ10 gene encodes Kir4.1 channels and plays a major role in modulating resting membrane potentials in excitable cells. It may cause GGEs if mutated. The purpose of this study was to investigate the possible association between KCNJ10 common variants and the susceptibility and drug resistance of GGEs in Chinese population. The allele-specific MALDI–TOF mass spectrometry method was used to assess 8 single nucleotide polymorphisms (SNPs) of KCNJ10 in 284 healthy controls and 483 Chinese GGEs patients including 279 anti-epileptic drug responsive patients and 204 drug resistant patients. We found the rs6690889 TC+TT genotypes were lower frequency in the GGEs group than that in the healthy controls (6.7% vs 9.5%, p = 0.01, OR = 0.50[0.29–0.86]). The frequency of rs1053074 G allele was lower in the childhood absence epilepsy (CAE) group than that in the healthy controls (28.4% vs 36.2%, p = 0.01, OR = 0.70[0.53–0.93]). The frequency of rs12729701 G allele and AG+GG genotypes was lower in the CAE group than that in the healthy controls (21.2% vs 28.4%, p = 0.01, OR = 0.74[0.59–0.94] and 36.3% vs 48.1%, p = 0.01, OR = 0.83[0.72–0.96], respectively). The frequency of rs12402969 C allele and the CC+CT genotypes were higher in the GGEs drug responsive patients than that in the drug resistant patients (9.3% vs 5.6%, OR = 1.73[1.06–2.85], p = 0.026 and 36.3% vs 48.1%, p = 0.01, OR = 0.83[0.72–0.96], respectively). This study identifies potential SNPs of KCNJ10 gene that may contribute to seizure susceptibility and anti-epileptic drug resistance.     

Common Variants of FTO Are Associated with Childhood Obesity in a Cross-Sectional Study of 3,126 Urban Indian Children

Background

FTO variants are robustly associated with obesity and related traits in many population and shown to have variable impact during life course. Although studies have shown association of FTO variants with adiposity in adult Indian, its association in Indian children is yet to be confirmed.

Methods

Here we examined association of FTO variants (rs9939609 and rs8050136) with obesity and related anthropometric and biochemical traits in 3,126 Indian children (aged 11–17 years) including 2,230 normal-weight and 896 over-weight/obese children. We also compared effects observed in the present study with that observed in previous studies on South Asian adults and children of other ethnic groups.

Results

The variant rs9939609 showed significant association with risk of obesity [OR = 1.21, P = 2.5×10⁻³] and its measures BMI, weight, waist circumference and hip circumference [β range = 0.11 to 0.14 Z-score units; P range = 1.3×10⁻⁴ to 1.6×10⁻⁷] in children. The observed effect sizes in Indian children were similar to those reported for European children. Variant rs9939609 explained 0.88% of BMI variance in Indian children. The effect sizes of rs9939609 on BMI and WC were ∼2 fold higher in children than adults. Interestingly rs9939609 was also associated with serum levels of thyroid stimulating hormone (TSH) [β = 0.10 Z-score, P = 5.8×10⁻³]. The other variant rs8050136 was in strong linkage disequilibrium with rs9939609 (r² = 0.97) and provided similar association results.

Conclusion

The study provides first report of association of FTO variants with obesity and related anthropometric traits in Indian children with higher impact in children compared to adults. We also demonstrated association of FTO variant with serum levels of TSH, indicating putative influence of FTO in hypothalamic-pituitary-thyroid axis.     

Massive Withdrawal Symptoms and Affective Vulnerability Are Associated with Variants of the CHRNA4 Gene in a Subgroup of Smokers

Heterogeneous phenotypes of complex disorders pose a great challenge for genetic association studies and for the development of personalized treatment strategies. Cluster analysis of phenotypic data has been recently proposed as a reliable auxiliary method for such studies. A cohort of 236 treatment-seeking smokers was investigated after overnight nicotine abstinence. Alpha4 nicotinic acetylcholine receptor (nAChR) subunit-related phenotypes were assessed by the Fagerström Test for Nicotine Dependence (FTND), exhaled carbon monoxide (CO) measurements, the Minnesota Nicotine Withdrawal Scale (MNWS) and the Zung Self-Rating Depression Scale (ZSDS). Seven tag SNPs (single-nucleotide polymorphisms) across CHRNA4 (the gene encoding alpha4 subunit of the nicotinic acetylcholine receptor) were genotyped and two-step cluster analysis was used for phenotypic cluster characterization. Haplotype estimation was determined by HapStat module of R 2.0 software. Three different phenotypic clusters were identified and the C3 cluster was characterized by the highest ZSDS and MNWS scores compared to others. Furthermore, lifetime prevalence of major depression was significantly higher in the C3 cluster (p = 0.019). In genetic association tests, this cluster was also significantly associated with rs3787138 genotypes (p = 0.004) while haplotype analyses of three SNPs (rs3787138, rs1044396, rs3787140) revealed that the risk for C3 phenotype was almost three times higher in GCC haplotype carriers compared to others (p_perm = 0.013). This is the first report on a significant association between CHRNA4 variants and a subgroup of smokers characterized by massive withdrawal symptoms and affective vulnerability. Identification of such a phenotypic cluster can be a pivotal step for further pharmacogenetic studies on ligands of the alpha4 nAChR subunit. Our results suggest that performing cluster analysis in genetic association studies can be proposed for complex disorders.     

Coding Variants of the FMO3 Gene Are Associated with the Risk of Chronic Kidney Disease: A Case-Control Study

Background:Chronic kidney disease (CKD) is a global health concern involving roughly one-tenth of developed countries'' populations. The flavin-containing dimethylaniline monooxygenase 3 (FMO3) gene encodes an enzyme that catalyzes trimethylamine N-oxide (TMAO), a toxin in CKD sufferers. This preliminary study aims to evaluate the association between coding region variations of FMO3, rs2266782G/A (E158K), rs2266780A/G (E308G), and rs1736557G/A (V257M), and the susceptibility to CKD.Methods:A total of 356 participants were enrolled, including 157 patients diagnosed with CKD and 199 age-matched healthy individuals. Genotyping of FMO3 gene variations was performed via PCR-RFLP and ARMS-PCR methods.Results:Our findings revealed a significant association between rs2266780A/G and rs1736557G/A and CKD under different genetic models. Compared to the GGG haplotype of rs2266782/rs1736557/rs2266780, the GAG, GAA, AAG, and AAA haplotype combinations conferred an increased risk of CKD in our population. Interaction analysis revealed that some genotype combinations, including GA/AA/AA, AA/AA/AA, GA/AA/GA, and GG/AG/AA, dramatically increased CKD risk in the Iranian population. No correlation was found between FMO3 polymorphisms and CKD stages.DiscussionThese observations highlight the potential impact of coding variants of the FMO3 gene on the onset of CKD. Further investigations into expanded populations and diverse races are needed to confirm our findings.Key Words: Chronic kidney disease, FMO3, Genetic variant, Single-nucleotide polymorphism, Trimethylamine N-oxide     

Common Variants Related to Serum Uric Acid Concentrations Are Associated with Glucose Metabolism and Insulin Secretion in a Chinese Population

Background

Elevated serum uric acid concentration is an independent risk factor and predictor of type 2 diabetes (T2D). Whether the uric acid-associated genes have an impact on T2D remains unclear. We aimed to investigate the effects of the uric acid-associated genes on the risk of T2D as well as glucose metabolism and insulin secretion.

Method

We recruited 2,199 normal glucose tolerance subjects from the Shanghai Diabetes Study I and II and 2,999 T2D patients from the inpatient database of Shanghai Diabetes Institute. Fifteen single nucleotide polymorphisms (SNPs) mapped in or near 11 loci (PDZK1, GCKR, LRP2, SLC2A9, ABCG2, LRRC16A, SLC17A1, SLC17A3, SLC22A11, SLC22A12 and SF1) were genotyped and serum biochemical parameters related to uric acid and T2D were determined.

Results

SF1 rs606458 showed strong association to T2D in both males and females (p = 0.034 and 0.0008). In the males, LRRC16A was associated with 2-h insulin and insulin secretion (p = 0.009 and 0.009). SLC22A11 was correlated with HOMA-B and insulin secretion (p = 0.048 and 0.029). SLC2A9 rs3775948 was associated with 2-h glucose (p = 0.043). In the females, LRP2 rs2544390 and rs1333049 showed correlations with fasting insulin, HOMA-IR and insulin secretion (p = 0.028, 0.033 and 0.052 and p = 0.034, 0.047 and 0.038, respectively). SLC2A9 rs11722228 was correlated with 2-h glucose, 2-h insulin and insulin secretion (p = 0.024, 0.049 and 0.049, respectively).

Conclusions

Our results indicated that the uric acid-associated genes have an impact on the risk of T2D, glucose metabolism and insulin secretion in a Chinese population.     

New Gene Variants Associated with the Risk of Chronic HBV Infection

Some patients with chronic hepatitis B virus(HBV) infection failed to clear HBV, even persistently continue to produce antibodies to HBV. Here we performed a two stage genome wide association study in a cohort of Chinese patients designed to discover single nucleotide variants that associate with HBV infection and clearance of HBV. The first stage involved genome wide exome sequencing of 101 cases(HBsAg plus anti-HBs positive) compared with 102 control patients(antiHBs positive, HBsAg negative). Over 80% of individual sequences displayed 209 sequence coverage. Adapters,uncertain bases [10% or low-quality base calls([50%) were filtered and compared to the human reference genome hg19. In the second stage, 579 chronic HBV infected cases and 439 HBV clearance controls were sequenced with selected genes from the first stage. Although there were no significant associated gene variants in the first stage, two significant gene associations were discovered when the two stages were assessed in a combined analysis. One association showed rs506121-‘‘T' allele [within the dedicator of cytokinesis 8(DOCK8) gene] was higher in chronic HBV infection group than that in clearance group(P = 0.002, OR = 0.77, 95% CI [0.65, 0.91]). The second association involved rs2071676—A allele within the Carbonic anhydrase(CA9) gene that was significantly elevated in chronic HBV infection group compared to the clearance group(P = 0.0003, OR = 1.35, 95% CI [1.15, 1.58]). Upon replication these gene associations would suggest the influence of DOCK8 and CA9 as potential risk genetic factors in the persistence of HBV infection.     

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.