Investigating the Relationship Between Rare Genetic Variants and Fibrosis in Pediatric Nonalcoholic Fatty Liver Disease

Background and Aims

Nonalcoholic Fatty Liver Disease (NAFLD) is a complex human disease. Common genetic variation in the patatin-like phospholipase domain containing 3 (PNPLA3) and transmembrane 6 superfamily member 2 (TM6SF2) genes have been associated with an increased risk of developing NAFLD, nonalcoholic steatohepatitis (NASH), and fibrosis in adults. The role of rare genetic variants in the development and progression of NAFLD in children is not well known. We aimed to explore the role of rare genetic variants in pediatric patients with advanced fibrosis.

Methods

Whole exome sequencing data was generated for 229 pediatric patients diagnosed with NAFLD recruited from the NASH Clinical Research Network (NASH CRN). Case-control single variant and gene-based collapsing analyses were used to test for rare variants that were enriched or depleted within the pediatric NAFLD cohort specifically for advanced fibrosis (cases) versus those without fibrosis (controls) or six other histologic characteristics. Exome data from non-NAFLD population controls were also used for additional analyses. All results were adjusted for multiple testing using a Bonferroni correction.

Results

No genome-wide significant associations were found between rare variation and presence of advanced fibrosis or NASH, nor the severity of steatosis, inflammation, or hepatocellular ballooning. Significantly, no enrichment of rare variants in PNPLA3 or TM6SF2 was observed across phenotypes.

Conclusion

In a cohort of children with histologically proven NAFLD, no genome-wide significant associations were found between rare genetic variation and advanced fibrosis or six other histologic features. Of particular interest was the lack of association with genes of interest in adults: PNPLA3 and TM6SF2, though limitations in sample size may reduce the ability to detect associations, particularly with rare variation.

The best of both worlds: Blending cutting-edge research with clinical processes for a productive exome clinic

Genomic medicine has been transformed by next-generation sequencing (NGS), inclusive of exome sequencing (ES) and genome sequencing (GS). Currently, ES is offered widely in clinical settings, with a less prevalent alternative model consisting of hybrid programs that incorporate research ES along with clinical patient workflows. We were among the earliest to implement a hybrid ES clinic, have provided diagnoses to 45% of probands, and have identified several novel candidate genes. Our program is enabled by a cost-effective investment by the health system and is unique in encompassing all the processes that have been variably included in other hybrid/clinical programs. These include careful patient selection, utilization of a phenotype-agnostic bioinformatics pipeline followed by manual curation of variants and phenotype integration by clinicians, close collaborations between the clinicians and the bioinformatician, pursuit of interesting variants, communication of results to patients in categories that are predicated upon the certainty of a diagnosis, and tracking changes in results over time and the underlying mechanisms for such changes. Due to its effectiveness, scalability to GS and its resource efficiency, specific elements of our paradigm can be incorporated into existing clinical settings, or the entire hybrid model can be implemented within health systems that have genomic medicine programs, to provide NGS in a scientifically rigorous, yet pragmatic setting.

The diagnostic yield of exome sequencing in liver diseases from a curated gene panel

Exome sequencing (ES) has been used in a variety of clinical settings but there are limited data on its utility for diagnosis and/or prediction of monogenic liver diseases. We developed a curated list of 502 genes for monogenic disorders associated with liver phenotypes and analyzed ES data for these genes in 758 patients with chronic liver diseases (CLD). For comparison, we examined ES data in 7856 self-declared healthy controls (HC), and 2187 patients with chronic kidney disease (CKD). Candidate pathogenic (P) or likely pathogenic (LP) variants were initially identified in 19.9% of participants, most of which were attributable to previously reported pathogenic variants with implausibly high allele frequencies. After variant annotation and filtering based on population minor allele frequency (MAF ≤ 10-4 for dominant disorders and MAF ≤ 10-3 for recessive disorders), we detected a significant enrichment of P/LP variants in the CLD cohort compared to the HC cohort (X2 test OR 5.00, 95% CI 3.06-8.18, p value = 4.5e-12). A second-level manual annotation was necessary to capture true pathogenic variants that were removed by stringent allele frequency and quality filters. After these sequential steps, the diagnostic rate of monogenic disorders was 5.7% in the CLD cohort, attributable to P/LP variants in 25 genes. We also identified concordant liver disease phenotypes for 15/22 kidney disease patients with P/LP variants in liver genes, mostly associated with cystic liver disease phenotypes. Sequencing results had many implications for clinical management, including familial testing for early diagnosis and management, preventative screening for associated comorbidities, and in some cases for therapy. Exome sequencing provided a 5.7% diagnostic rate in CLD patients and required multiple rounds of review to reduce both false positive and false negative findings. The identification of concordant phenotypes in many patients with P/LP variants and no known liver disease also indicates a potential for predictive testing for selected monogenic liver disorders.

A genome-wide screen for variants influencing certolizumab pegol response in a moderate to severe rheumatoid arthritis population

Certolizumab pegol (CZP) is a PEGylated Fc-free tumor necrosis factor (TNF) inhibitor antibody approved for use in the treatment of rheumatoid arthritis (RA), Crohn’s disease, psoriatic arthritis, axial spondyloarthritis and psoriasis. In a clinical trial of patients with severe RA, CZP improved disease symptoms in approximately half of patients. However, variability in CZP efficacy remains a problem for clinicians, thus, the aim of this study was to identify genetic variants predictive of CZP response. We performed a genome-wide association study (GWAS) of 302 RA patients treated with CZP in the REALISTIC trial to identify common single nucleotide polymorphisms (SNPs) associated with treatment response. Whole-exome sequencing was also performed for 74 CZP extreme responders and non-responders within the same population, as well as 1546 population controls. No common SNPs or rare functional variants were significantly associated with CZP response, though a non-significant enrichment in the RA-implicated KCNK5 gene was observed. Two SNPs near spondin-1 and semaphorin-4G approached genome-wide significance. The results of the current study did not provide an unambiguous predictor of CZP response.

Association of CYP2B6 Single-Nucleotide Polymorphisms Altering Efavirenz Metabolism With Hepatitis C Virus (HCV) Treatment Relapse Among Human Immunodeficiency Virus/HCV-Coinfected African Americans Receiving Ledipasvir/Sofosbuvir in the ION-4 Trial

In the ION-4 trial, hepatitis C virus relapse was rare, occurring only in African American patients, 80% receiving efavirenz for human immunodeficiency virus infection. We observed no indication that CYP2B6 polymorphisms associated with increased plasma efavirenz exposure explained the relapses.

Genome-wide association study (GWAS) of human host factors influencing viral severity of herpes simplex virus type 2 (HSV-2)

Herpes simplex virus type 2 (HSV-2) is an incurable viral infection with severity ranging from asymptomatic to frequent recurrences. The viral shedding rate has been shown as a reproducible HSV-2 severity end point that correlates with lesion rates. We used a genome-wide association study (GWAS) to investigate the role of common human genetic variation in HSV-2 severity. We performed a GWAS on 223 HSV-2-positive participants of European ancestry. Severity was measured by viral shedding rate, as defined by the percent of days PCR+ for HSV-2 DNA over at least 30 days. Analyses were performed under linear regression models, adjusted for age, sex, and ancestry. There were no genome-wide significant (p < 5E-08) associations with HSV-2 viral shedding rate. The top nonsignificant SNP (rs75932292, p = 6.77E-08) associated with HSV-2 viral shedding was intergenic, with the nearest known biologically interesting gene (ABCA1) ~130 kbp downstream. Several other SNPs approaching significance were in or near genes with viral or neurological associations, including four SNPs in KIF1B. The current study is the first comprehensive genome-wide investigation of human genetic variation in virologic severity of established HSV-2 infection. However, no significant associations were observed with HSV-2 virologic severity, leaving the exact role of human variation in HSV-2 severity unclear.

Exome sequencing of an adolescent with nonalcoholic fatty liver disease identifies a clinically actionable case of Wilson disease

Diagnostic whole-exome sequencing has proven highly successful in a range of rare diseases, particularly early-onset genetic conditions. In more common conditions, however, exome sequencing for diagnostic purposes remains the exception. Here we describe a patient initially diagnosed with a common, complex liver disease, nonalcoholic fatty liver disease (NAFLD), who was determined to have Wilson disease (WD) upon research-related exome sequencing. The patient presented as a 14.5-yr-old adolescent with chronically elevated aminotransferases, normal ceruloplasmin, and histologic examination consistent with NAFLD with advanced fibrosis. He was enrolled in a large longitudinal study of patients with NAFLD and was found to have WD by exome sequencing performed 4 yr later. This new diagnosis, confirmed clinically by 24 h urine copper quantification, led to a change in the therapy from lifestyle counseling to directed treatment with d-penicillamine, a copper chelating agent. In this case, the likelihood of making the correct diagnosis and thereby choosing the appropriate treatment was increased by exome sequencing and careful interpretation. This example illustrates the utility of exome sequencing diagnostically in more common conditions not currently considered as targets for genome-wide evaluation and adds to a growing body of evidence that patients diagnosed with more common conditions often in fact have rarer genetically determined syndromes that have escaped clinical detection.

Killer cell immunoglobulin-like receptor 3DL1 variation modifies HLA-B*57 protection against HIV-1

HLA-B*57 control of HIV involves enhanced CD8+ T cell responses against infected cells, but extensive heterogeneity exists in the level of HIV control among B*57+ individuals. Using whole-genome sequencing of untreated B*57+ HIV-1-infected controllers and noncontrollers, we identified a single variant (rs643347A/G) encoding an isoleucine-to-valine substitution at position 47 (I47V) of the inhibitory killer cell immunoglobulin-like receptor KIR3DL1 as the only significant modifier of B*57 protection. The association was replicated in an independent cohort and across multiple outcomes. The modifying effect of I47V was confined to B*57:01 and was not observed for the closely related B*57:03. Positions 2, 47, and 54 tracked one another nearly perfectly, and 2 KIR3DL1 allotypes differing only at these 3 positions showed significant differences in binding B*57:01 tetramers, whereas the protective allotype showed lower binding. Thus, variation in an immune NK cell receptor that binds B*57:01 modifies its protection. These data highlight the exquisite specificity of KIR-HLA interactions in human health and disease.

Evaluating the Impact of Functional Genetic Variation on HIV-1 Control

Background

Previous genetic association studies of human immunodeficiency virus-1 (HIV-1) progression have focused on common human genetic variation ascertained through genome-wide genotyping.

Methods

We sought to systematically assess the full spectrum of functional variation in protein coding gene regions on HIV-1 progression through exome sequencing of 1327 individuals. Genetic variants were tested individually and in aggregate across genes and gene sets for an influence on HIV-1 viral load.

Results

Multiple single variants within the major histocompatibility complex (MHC) region were observed to be strongly associated with HIV-1 outcome, consistent with the known impact of classical HLA alleles. However, no single variant or gene located outside of the MHC region was significantly associated with HIV progression. Set-based association testing focusing on genes identified as being essential for HIV replication in genome-wide small interfering RNA (siRNA) and clustered regularly interspaced short palindromic repeats (CRISPR) studies did not reveal any novel associations.

Conclusions

These results suggest that exonic variants with large effect sizes are unlikely to have a major contribution to host control of HIV infection.

No Evidence for Association of β-Defensin Genomic Copy Number with HIV Susceptibility, HIV Load during Clinical Latency, or Progression to AIDS

Common single-nucleotide variation in the host accounts for 25% of the variability in the plasma levels of HIV during the clinical latency stage (viral load set point). However, the role of rare variants and copy number variants remains relatively unexplored. Previous work has suggested copy number variation of a cluster of β-defensin genes affects HIV load in treatment-naïve sub-Saharan Africans and rate of response to antiretroviral treatment. Here we analyse a total of 1827 individuals from two cohorts of HIV-infected individuals from Europe and sub-Saharan Africa to investigate the role of β-defensin copy number variation on HIV load at set point. We find no evidence for association of copy number with viral load. We also compare distribution of β-defensin copy number between European cases and controls and find no differences, arguing against a role of β-defensin copy number in HIV acquisition. Taken together, our data argue against an effect of copy number variation of the β-defensin region in the spontaneous control of HIV infection.

Immune perturbations in HIV-1-infected individuals who make broadly neutralizing antibodies

Induction of broadly neutralizing antibodies (bnAbs) is a goal of HIV-1 vaccine development. bnAbs occur in some HIV-1-infected individuals and frequently have characteristics of autoantibodies. We have studied cohorts of HIV-1-infected individuals who made bnAbs and compared them with those who did not do so, and determined immune traits associated with the ability to produce bnAbs. HIV-1-infected individuals with bnAbs had a higher frequency of blood autoantibodies, a lower frequency of regulatory CD4⁺ T cells, a higher frequency of circulating memory T follicular helper CD4⁺ cells, and a higher T regulatory cell level of programmed cell death-1 expression compared with HIV-1-infected individuals without bnAbs. Thus, induction of HIV-1 bnAbs may require vaccination regimens that transiently mimic immunologic perturbations in HIV-1-infected individuals.

Testing for risk and protective trends in genetic analyses of HIV acquisition

Host genetics studies of HIV-1 acquisition are critically important for the identification of new targets for drug and vaccine development. Analyses of such studies typically focus on pairwise comparisons of three different groups: HIV-1 positive individuals, HIV-1 high-risk seronegative individuals, and population controls. Because there is a clear expectation of how gene frequencies of risk or protective alleles would be ordered in the three groups, we are able to construct a statistical framework that offers a consistent increase in power over a wide-range of the magnitude of risk/protective effects. In this paper, we develop tests that constrain the alternative hypothesis to appropriately reflect risk or protective trends jointly across the three groups and show that they lead to a substantial increase in power over the naive pairwise approach. We develop both likelihood-ratio and score statistics that test for genetic effects across the three groups while constraining the alternative hypothesis to reflect biologically motivated trends of risk or protection. The asymptotic distribution of both statistics (likelihood ratio and score) is derived. We investigate the performance of our approach via extensive simulation studies using a biologically motivated model of HIV-1 acquisition, and find that our proposed approach leads to an increase in power of roughly 10-28%. We illustrate our approach with an analysis of the effect of the CCR5Δ32 mutation on HIV acquisition.

Association study of common genetic variants and HIV-1 acquisition in 6,300 infected cases and 7,200 controls

Multiple genome-wide association studies (GWAS) have been performed in HIV-1 infected individuals, identifying common genetic influences on viral control and disease course. Similarly, common genetic correlates of acquisition of HIV-1 after exposure have been interrogated using GWAS, although in generally small samples. Under the auspices of the International Collaboration for the Genomics of HIV, we have combined the genome-wide single nucleotide polymorphism (SNP) data collected by 25 cohorts, studies, or institutions on HIV-1 infected individuals and compared them to carefully matched population-level data sets (a list of all collaborators appears in Note S1 in Text S1). After imputation using the 1,000 Genomes Project reference panel, we tested approximately 8 million common DNA variants (SNPs and indels) for association with HIV-1 acquisition in 6,334 infected patients and 7,247 population samples of European ancestry. Initial association testing identified the SNP rs4418214, the C allele of which is known to tag the HLA-B*57:01 and B*27:05 alleles, as genome-wide significant (p = 3.6×10⁻¹¹). However, restricting analysis to individuals with a known date of seroconversion suggested that this association was due to the frailty bias in studies of lethal diseases. Further analyses including testing recessive genetic models, testing for bulk effects of non-genome-wide significant variants, stratifying by sexual or parenteral transmission risk and testing previously reported associations showed no evidence for genetic influence on HIV-1 acquisition (with the exception of CCR5Δ32 homozygosity). Thus, these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Comparing the frequency differences between common DNA variants in disease-affected cases and in unaffected controls has been successful in uncovering the genetic component of multiple diseases. This approach is most effective when large samples of cases and controls are available. Here we combine information from multiple studies of HIV infected patients, including more than 6,300 HIV+ individuals, with data from 7,200 general population samples of European ancestry to test nearly 8 million common DNA variants for an impact on HIV acquisition. With this large sample we did not observe any single common genetic variant that significantly associated with HIV acquisition. We further tested 22 variants previously identified by smaller studies as influencing HIV acquisition. With the exception of a deletion polymorphism in the CCR5 gene (CCR5Δ32) we found no convincing evidence to support these previous associations. Taken together these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Host genetics of HIV acquisition and viral control

Since the discovery of HIV as the cause of AIDS, numerous insights have been gained from studies of its natural history and epidemiology. It has become clear that there are substantial interindividual differences in the risk of HIV acquisition and course of disease. Meanwhile, the field of human genetics has undergone a series of rapid transitions that have fundamentally altered the approach to studying HIV host genetics. We aim to describe the field as it has transitioned from the era of candidate-gene studies and the era of genome-wide association studies (GWAS) to its current state in the infancy of comprehensive sequencing. In some ways the field has come full circle, having evolved from being driven almost exclusively by our knowledge of immunology, to a bias-free GWAS approach, to a point where our ability to catalogue human variation far outstrips our ability to biologically interpret it.

Group A Streptococcus emm gene types in pharyngeal isolates, Ontario, Canada, 2002-2010

Determination of emm variations may help improve vaccine design.

Group A Streptococcus (GAS) is a human-adapted pathogen that causes a variety of diseases, including pharyngitis and invasive infections. GAS strains are categorized by variation in the nucleotide sequence of the gene (emm) that encodes the M protein. To identify the emm types of GAS strains causing pharyngitis in Ontario, Canada, we sequenced the hypervariable region of the emm gene in 4,635 pharyngeal GAS isolates collected during 2002–2010. The most prevalent emm types varied little from year to year. In contrast, fine-scale geographic analysis identified inter-site variability in the most common emm types. Additionally, we observed fluctuations in yearly frequency of emm3 strains from pharyngitis patients that coincided with peaks of emm3 invasive infections. We also discovered a striking increase in frequency of emm89 strains among isolates from patients with pharyngitis and invasive disease. These findings about the epidemiology of GAS are potentially useful for vaccine research.

Polymorphisms in regulator of protease B (RopB) alter disease phenotype and strain virulence of serotype M3 group A Streptococcus

Whole-genome sequencing of serotype M3 group A streptococci (GAS) from oropharyngeal and invasive infections in Ontario recently showed that the gene encoding regulator of protease B (RopB) is highly polymorphic in this population. To test the hypothesis that ropB is under diversifying selective pressure among all serotype M3 GAS strains, we sequenced this gene in 1178 strains collected from different infection types, geographic regions, and time periods. The results confirmed our hypothesis and discovered a significant association between mutant ropB alleles, decreased activity of its major regulatory target SpeB, and pharyngitis. Additionally, isoallelic strains with ropB polymorphisms were significantly less virulent in a mouse model of necrotizing fasciitis. These studies provide a model strategy for applying whole-genome sequencing followed by deep single-gene sequencing to generate new insight to the rapid evolution and virulence regulation of human pathogens.

Naturally occurring single amino acid replacements in a regulatory protein alter streptococcal gene expression and virulence in mice

Infection with different strains of the same species of bacteria often results in vastly different clinical outcomes. Despite extensive investigation, the genetic basis of microbial strain-specific virulence remains poorly understood. Recent whole-genome sequencing has revealed that SNPs are the most prevalent form of genetic diversity among different strains of the same species of bacteria. For invasive serotype M3 group A streptococci (GAS) strains, the gene encoding regulator of proteinase B (RopB) has the highest frequency of SNPs. Here, we have determined that ropB polymorphisms alter RopB function and modulate GAS host-pathogen interactions. Sequencing of ropB in 171 invasive serotype M3 GAS strains identified 19 distinct ropB alleles. Inactivation of the ropB gene in strains producing distinct RopB variants had dramatically divergent effects on GAS global gene expression. Additionally, generation of isoallelic GAS strains differing only by a single amino acid in RopB confirmed that variant proteins affected transcript levels of the gene encoding streptococcal proteinase B, a major RopB-regulated virulence factor. Comparison of parental, RopB-inactivated, and RopB isoallelic strains in mouse infection models demonstrated that ropB polymorphisms influence GAS virulence and disease manifestations. These data detail a paradigm in which unbiased, whole-genome sequence analysis of populations of clinical bacterial isolates creates new avenues of productive investigation into the pathogenesis of common human infections.

RNASEL and RNASEL-inhibitor variation and prostate cancer risk in Afro-Caribbeans

BACKGROUND

Afro-Caribbeans from Tobago are at high risk of developing prostate cancer. This elevated risk of prostate cancer is shared by populations of African ancestry living in diverse environments in the Western hemisphere. Variation in the ribonuclease L (RNASEL) gene has recently been reported to be associated with an increased risk of prostate cancer. However, whether RNASEL variation contributes to the increased risk of prostate cancer observed in populations of African ancestry remains unclear.

METHODS

We resequenced the positional candidate gene RNASEL in 48 prostate cancer cases and genotyped the previously reported R462Q and D541E polymorphisms in 230 prostate cancer cases and 458 controls. We also examined the inhibitor of RNASEL (ABCE1) for variation associated with prostate cancer risk.

RESULTS

We found no evidence of association between R462Q and D541E polymorphisms and prostate cancer risk in our case/control analysis. A novel variant (K294E) was identified in a single heterozygous individual with prostate cancer. We also observed a 20 bp insertion/deletion polymorphism 1,109 bp upstream of the initiation codon, but this variant was not associated with prostate cancer. We identified 16 single nucleotide polymorphisms in the ABCE1 gene, only 3 of which had a minor allele frequency >5%. A common A/G transition -1,071 bp from the transcriptional start site was genotyped and showed no evidence of association with prostate cancer.

CONCLUSIONS

Our results suggest that common variation in the putative prostate cancer susceptibility gene, RNASEL, or its inhibitor does not contribute significantly to prostate cancer risk in this Afro-Caribbean population.

Androgen receptor CAG repeat polymorphism is associated with fat-free mass in men

The human androgen receptor (AR) gene contains a CAG (glutamine) repeat polymorphism in exon 1 that is inversely associated with transcriptional activity of the AR. We studied the association of AR CAG repeat length, fat-free mass (FFM), and testosterone in two independent cohorts: 294 Caucasian men, aged 55-93 yr, from the Study of Osteoporotic Risk in Men (STORM), and 202 Caucasian volunteers (112 men and 90 women), aged 19-90 yr, from the Baltimore Longitudinal Study of Aging (BLSA). Subjects were genotyped to determine the number of AR CAG repeats and grouped as carrying either < 22 or > or =22 repeats. Whole body soft tissue composition was measured by dual-energy X-ray absorptiometry. Men with greater CAG repeat number exhibited significantly greater total FFM than those with fewer CAG repeats in both cohorts (STORM: 59.2 +/- 0.3 vs. 58.0 +/- 0.4 kg, P = 0.02; BLSA: 57.2 +/- 1.1 vs. 53.8 +/- 1.1 kg, P = 0.04). Similar results were observed for total FFM normalized to height. No differences were seen in women in the BLSA cohort. In the BLSA cohort, serum testosterone levels were higher in subjects with greater repeat number (P = 0.003). This same pattern approached significance in the STORM cohort (P = 0.07). In conclusion, the androgen receptor CAG repeat polymorphism is associated with FFM in men in two independent cohorts. Additional studies are needed to confirm this observation and to clarify the mechanisms involved.

Vitamin D Receptor Genotype Is Associated With Fat-Free Mass and Sarcopenia in Elderly Men

We investigated the association of vitamin D receptor (VDR) genotype with fat-free mass (FFM) in a cohort of 302 older (aged 58-93 years) Caucasian men who underwent body composition analysis by dual-energy X-ray absorptiometry, and completed questionnaires addressing comorbidities, physical activity, and dietary intake. All participants were genotyped for a VDR translation start site (FokI) polymorphism [FF (37.7%), Ff (48.4%), and ff (13.9%)] and the previously studied BsmI polymorphism [BB (24.9%), Bb (37.7%), and bb (37.4%)]. The BsmI polymorphism was not associated with FFM in any analysis; however, the FokI polymorphism was significantly associated with total FFM, appendicular FFM, and relative (kg/m(2)) appendicular FFM (all p <.05), with the FF group demonstrating significantly lower FFM than the Ff and ff groups (e.g., total FFM: FF = 57.6 +/- 0.4, Ff = 59.4 +/- 0.4, ff = 59.4 +/- 0.7 kg; p <.02). Age-adjusted logistic regression revealed a 2.17-fold higher risk for sarcopenia (defined previously as appendicular FFM <7.26 kg/m(2)) in FF homozygotes (95% CI [confidence interval] = 1.19-3.85; p =.03) compared to men with one or more f alleles. The VDR translation start site (FokI) polymorphism is significantly associated with FFM and sarcopenia in this cohort of older Caucasian men.

ELAC2 and prostate cancer risk in Afro-Caribbeans of Tobago

To test the hypothesis that variation in the putative prostate cancer susceptibility gene ELAC2 contributes to the elevated risk of prostate cancer in Afro-Caribbean males from Tobago, we genotyped the S217L and A514T polymorphisms, previously reported to be associated with prostate cancer risk in a large sample of cases and controls. The frequency of the high-risk Leu allele at the S217L site was the same in cases and controls. Both cases and controls were homozygous for the low-risk Ala allele at the A514T site. In addition, we sequenced the exons and 3'- and 5'-flanking regions of ELAC2 in 24 individuals with histologically confirmed prostate cancer. We identified 17 new single nucleotide polymorphisms. An A(-1196)T polymorphism, which alters a predicted TATA box consensus sequence, was tested in cases and controls, and no significant difference in allele or genotype frequencies was observed. The absence of ELAC2 mutations and lack of association between polymorphisms in ELAC2 and prostate cancer in cases and controls leads us to conclude that ELAC2 does not contribute significantly to the elevated prevalence of prostate cancer in Afro-Caribbean males of Tobago.