A comprehensive analysis of common genetic variation in prolactin (PRL) and PRL receptor (PRLR) genes in relation to plasma prolactin levels and breast cancer risk: the Multiethnic Cohort

<p>Abstract</p> <p>Background</p> <p>Studies in animals and humans clearly indicate a role for prolactin (PRL) in breast epithelial proliferation, differentiation, and tumorigenesis. Prospective epidemiological studies have also shown that women with higher circulating PRL levels have an increase in risk of breast cancer, suggesting that variability in PRL may also be important in determining a woman's risk.</p> <p>Methods</p> <p>We evaluated genetic variation in the PRL and PRL receptor (PRLR) genes as predictors of plasma PRL levels and breast cancer risk among African-American, Native Hawaiian, Japanese-American, Latina, and White women in the Multiethnic Cohort Study (MEC). We selected single nucleotide polymorphisms (SNPs) from both the public (dbSNP) and private (Celera) databases to construct high density SNP maps that included up to 20 kilobases (kb) upstream of the transcription initiation site and 10 kb downstream of the last exon of each gene, for a total coverage of 59 kb in PRL and 210 kb in PRLR. We genotyped 80 SNPs in PRL and 173 SNPs in PRLR in a multiethnic panel of 349 unaffected subjects to characterize linkage disequilibrium (LD) and haplotype patterns. We sequenced the coding regions of PRL and PRLR in 95 advanced breast cancer cases (19 of each racial/ethnic group) to uncover putative functional variation. A total of 33 and 60 haplotype "tag" SNPs (tagSNPs) that allowed for high predictability (R_h²≥ 0.70) of the common haplotypes in PRL and PRLR, respectively, were then genotyped in a multiethnic breast cancer case-control study of 1,615 invasive breast cancer cases and 1,962 controls in the MEC. We also assessed the association of common genetic variation with circulating PRL levels in 362 postmenopausal controls without a history of hormone therapy use at blood draw. Because of the large number of comparisons being performed we used a relatively stringent type I error criteria (p < 0.0005) for evaluating the significance of any single association to correct for performing approximately 100 independent tests, close to the number of tagSNPs genotyped for both genes.</p> <p>Results</p> <p>We observed no significant associations between PRL and PRLR haplotypes or individual SNPs in relation to breast cancer risk. A nominally significant association was noted between prolactin levels and a tagSNP (tagSNP 44, rs2244502) in intron 1 of PRL. This SNP showed approximately a 50% increase in levels between minor allele homozygotes vs. major allele homozygotes. However, this association was not significant (p = 0.002) using our type I error criteria to correct for multiple testing, nor was this SNP associated with breast cancer risk (p = 0.58).</p> <p>Conclusion</p> <p>In this comprehensive analysis covering 59 kb of the PRL locus and 210 kb of the PRLR locus, we found no significant association between common variation in these candidate genes and breast cancer risk or plasma PRL levels. The LD characterization of PRL and PRLR in this multiethnic population provide a framework for studying these genes in relation to other disease outcomes that have been associated with PRL, as well as for larger studies of plasma PRL levels.</p

What has GWAS done for HLA and disease associations?

The major histocompatibility complex (MHC) is located in chromosome 6p21 and contains crucial regulators of immune response, including human leucocyte antigen (HLA) genes, alongside other genes with nonimmunological roles. More recently, a repertoire of noncoding RNA genes, including expressed pseudogenes, has also been identified. The MHC is the most gene dense and most polymorphic part of the human genome. The region exhibits haplotype-specific linkage disequilibrium patterns, contains the strongest cis- and trans-eQTLs/meQTLs in the genome and is known as a hot spot for disease associations. Another layer of complexity is provided to the region by the extreme structural variation and copy number variations. While the HLA-B gene has the highest number of alleles, the HLA-DR/DQ subregion is structurally most variable and shows the highest number of disease associations. Reliance on a single reference sequence has complicated the design, execution and analysis of GWAS for the MHC region and not infrequently, the MHC region has even been excluded from the analysis of GWAS data. Here, we contrast features of the MHC region with the rest of the genome and highlight its complexities, including its functional polymorphisms beyond those determined by single nucleotide polymorphisms or single amino acid residues. One of the several issues with customary GWAS analysis is that it does not address this additional layer of polymorphisms unique to the MHC region. We highlight alternative approaches that may assist with the analysis of GWAS data from the MHC region and unravel associations with all functional polymorphisms beyond single SNPs. We suggest that despite already showing the highest number of disease associations, the true extent of the involvement of the MHC region in disease genetics may not have been uncovered