Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

An optimized protocol for retina single-cell RNA sequencing.

Purpose: Single-cell RNA sequencing (scRNA-seq) is a powerful technique used to explore gene expression at the single cell level. However, appropriate preparation of samples is essential to obtain the most information out of this transformative technology. Generating high-quality single-cell suspensions from the retina is critical to preserve the native expression profile that will ensure meaningful transcriptome data analysis. Methods: We modified the conditions for rapid and optimal dissociation of retina sample preparation. We also included additional filtering steps in data analysis for retinal scRNA-seq. Results: We report a gentle method for dissociation of the mouse retina that minimizes cell death and preserves cell morphology. This protocol also results in detection of higher transcriptional complexity. In addition, the modified computational pipeline leads to better-quality single-cell RNA-sequencing data in retina samples. We also demonstrate the advantages and limitations of using fresh versus frozen retinas to prepare cell or nuclei suspensions for scRNA-seq. Conclusions: We provide a simple yet robust and reproducible protocol for retinal scRNA-seq analysis, especially for comparative studies.

Association of Aflatoxin and Gallbladder Cancer.

BACKGROUND & AIMS: Aflatoxin, which causes hepatocellular carcinoma, may also cause gallbladder cancer. We investigated whether patients with gallbladder cancer have higher exposure to aflatoxin than patients with gallstones. METHODS: We measured aflatoxin B1 (AFB1)-lysine adducts in plasma samples from the Shanghai Biliary Tract Cancer case-control study, conducted from 1997 through 2001. We calculated age- and sex-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) and the population-attributable fraction for 209 patients with gallbladder cancer and gallstones vs 250 patients with gallstones without cancer (controls). In 54 patients with gallbladder cancer, tumor tissue was examined for the R249S mutation in TP53, associated with aflatoxin exposure, through targeted sequencing. RESULTS: The AFB1-lysine adduct was detected in 67 (32%) of 209 patients with gallbladder cancer and 37 (15%) of the 250 controls (χ2 P < .0001), almost threefold more patients with gallbladder cancer than controls (OR, 2.71; 95% CI, 1.70-4.33). Among participants with detectable levels of AFB1-lysine, the median level of AFB1-lysine was 5.4 pg/mg in those with gallbladder cancer, compared with 1.2 pg/mg in controls. For patients in the fourth quartile of AFB1-lysine level vs the first quartile, the OR for gallbladder cancer was 7.61 (95% CI, 2.01-28.84). None of the 54 gallbladder tumors sequenced were found to have the R249S mutation in TP53. The population-attributable fraction for cancer related to aflatoxin was 20% (95% CI, 15%-25%). CONCLUSIONS: In a case-control study of patients with gallbladder cancer and gallstones vs patients with gallstones without cancer, we associated exposure to aflatoxin (based on plasma level of AFB1-lysine) with gallbladder cancer. Gallbladder cancer does not appear associate with the R249S mutation in TP53. If aflatoxin is a cause of gallbladder cancer, it may have accounted for up to 20% of the gallbladder cancers in Shanghai, China, during the study period, and could account for an even higher proportion in high-risk areas. If our findings are verified, reducing aflatoxin exposure might reduce the incidence of gallbladder cancer.

Multiple rare variants in high-risk pancreatic cancer-related genes may increase risk for pancreatic cancer in a subset of patients with and without germline CDKN2A mutations.

The risk of pancreatic cancer (PC) is increased in melanoma-prone families but the causal relationship between germline CDKN2A mutations and PC risk is uncertain, suggesting the existence of non-CDKN2A factors. One genetic possibility involves patients having mutations in multiple high-risk PC-related genes; however, no systematic examination has yet been conducted. We used next-generation sequencing data to examine 24 putative PC-related genes in 43 PC patients with and 23 PC patients without germline CDKN2A mutations and 1001 controls. For each gene and the four pathways in which they occurred, we tested whether PC patients (overall or CDKN2A+ and CDKN2A- cases separately) had an increased number of rare nonsynonymous variants. Overall, we identified 35 missense variants in PC patients, 14 in CDKN2A+ and 21 in CDKN2A- PC cases. We found nominally significant associations for mismatch repair genes (MLH1, MSH2, MSH6, PMS2) in all PC patients and for ATM, CPA1, and PMS2 in CDKN2A- PC patients. Further, nine CDKN2A+ and four CDKN2A- PC patients had rare potentially deleterious variants in multiple PC-related genes. Loss-of-function variants were only observed in CDKN2A- PC patients, with ATM having the most pathogenic variants. Also, ATM variants (n = 5) were only observed in CDKN2A- PC patients with a family history that included digestive system tumors. Our results suggest that a subset of PC patients may have increased risk because of germline mutations in multiple PC-related genes.

Whole exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We used whole exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mu-tant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

The principal genetic determinants for nasopharyngeal carcinoma in China involve the HLA class I antigen recognition groove.

Nasopharyngeal carcinoma (NPC) is an epithelial malignancy facilitated by Epstein-Barr Virus infection. Here we resolve the major genetic influences for NPC incidence using a genome-wide association study (GWAS), independent cohort replication, and high-resolution molecular HLA class I gene typing including 4,055 study participants from the Guangxi Zhuang Autonomous Region and Guangdong province of southern China. We detect and replicate strong association signals involving SNPs, HLA alleles, and amino acid (aa) variants across the major histocompatibility complex-HLA-A, HLA -B, and HLA -C class I genes (P(HLA-A-aa-site-62) = 7.4 × 10(-29); P (HLA-B-aa-site-116) = 6.5 × 10(-19); P (HLA-C-aa-site-156) = 6.8 × 10(-8) respectively). Over 250 NPC-HLA associated variants within HLA were analyzed in concert to resolve separate and largely independent HLA-A, -B, and -C gene influences. Multivariate logistical regression analysis collapsed significant associations in adjacent genes spanning 500 kb (OR2H1, GABBR1, HLA-F, and HCG9) as proxies for peptide binding motifs carried by HLA- A*11:01. A similar analysis resolved an independent association signal driven by HLA-B*13:01, B*38:02, and B*55:02 alleles together. NPC resistance alleles carrying the strongly associated amino acid variants implicate specific class I peptide recognition motifs in HLA-A and -B peptide binding groove as conferring strong genetic influence on the development of NPC in China.

Genome-wide association study implicates PARD3B-based AIDS restriction.

BACKGROUND: Host genetic variation influences human immunodeficiency virus (HIV) infection and progression to AIDS. Here we used clinically well-characterized subjects from 5 pretreatment HIV/AIDS cohorts for a genome-wide association study to identify gene associations with rate of AIDS progression. METHODS: European American HIV seroconverters (n = 755) were interrogated for single-nucleotide polymorphisms (SNPs) (n = 700,022) associated with progression to AIDS 1987 (Cox proportional hazards regression analysis, co-dominant model). RESULTS: Association with slower progression was observed for SNPs in the gene PARD3B. One of these, rs11884476, reached genome-wide significance (relative hazard = 0.3; P =3. 370 × 10(-9)) after statistical correction for 700,022 SNPs and contributes 4.52% of the overall variance in AIDS progression in this study. Nine of the top-ranked SNPs define a PARD3B haplotype that also displays significant association with progression to AIDS (hazard ratio, 0.3; P = 3.220 × 10(-8)). One of these SNPs, rs10185378, is a predicted exonic splicing enhancer; significant alteration in the expression profile of PARD3B splicing transcripts was observed in B cell lines with alternate rs10185378 genotypes. This SNP was typed in European cohorts of rapid progressors and was found to be protective for AIDS 1993 definition (odds ratio, 0.43, P = .025). CONCLUSIONS: These observations suggest a potential unsuspected pathway of host genetic influence on the dynamics of AIDS progression.

Genetic associations of variants in genes encoding HIV-dependency factors required for HIV-1 infection.

BACKGROUND: High-throughput genome-wide techniques have facilitated the identification of previously unknown host proteins involved in cellular human immunodeficiency virus (HIV) infection. Recently, 3 independent studies have used small interfering RNA technology to silence each gene in the human genome to determine the importance of each in HIV infection. Genes conferring a significant effect were termed HIV-dependency factors (HDFs). METHODS: We assembled high-density panels of 6380 single-nucleotide polymorphisms (SNPs) in 278 HDF genes and tested for genotype associations with HIV infection and AIDS progression in 1633 individuals from clinical AIDS cohorts. RESULTS: After statistical correction for multiple tests, significant associations with HIV acquisition were found for SNPs in 2 genes, NCOR2 and IDH1. Weaker associations with AIDS progression were revealed for SNPs within the TM9SF2 and EGFR genes. CONCLUSIONS: This study independently verifies the influence of NCOR2 and IDH1 on HIV transmission, and its findings suggest that variation in these genes affects susceptibility to HIV infection in exposed individuals.

Genetic variants in nuclear-encoded mitochondrial genes influence AIDS progression.

BACKGROUND: The human mitochondrial genome includes only 13 coding genes while nuclear-encoded genes account for 99% of proteins responsible for mitochondrial morphology, redox regulation, and energetics. Mitochondrial pathogenesis occurs in HIV patients and genetically, mitochondrial DNA haplogroups with presumed functional differences have been associated with differential AIDS progression. METHODOLOGY/PRINCIPAL FINDINGS: Here we explore whether single nucleotide polymorphisms (SNPs) within 904 of the estimated 1,500 genes that specify nuclear-encoded mitochondrial proteins (NEMPs) influence AIDS progression among HIV-1 infected patients. We examined NEMPs for association with the rate of AIDS progression using genotypes generated by an Affymetrix 6.0 genotyping array of 1,455 European American patients from five US AIDS cohorts. Successfully genotyped SNPs gave 50% or better haplotype coverage for 679 of known NEMP genes. With a Bonferroni adjustment for the number of genes and tests examined, multiple SNPs within two NEMP genes showed significant association with AIDS progression: acyl-CoA synthetase medium-chain family member 4 (ACSM4) on chromosome 12 and peroxisomal D3,D2-enoyl-CoA isomerase (PECI) on chromosome 6. CONCLUSIONS: Our previous studies on mitochondrial DNA showed that European haplogroups with presumed functional differences were associated with AIDS progression and HAART mediated adverse events. The modest influences of nuclear-encoded mitochondrial genes found in the current study add support to the idea that mitochondrial function plays a role in AIDS pathogenesis.

A high-density admixture map for disease gene discovery in african americans.

Admixture mapping (also known as "mapping by admixture linkage disequilibrium," or MALD) provides a way of localizing genes that cause disease, in admixed ethnic groups such as African Americans, with approximately 100 times fewer markers than are required for whole-genome haplotype scans. However, it has not been possible to perform powerful scans with admixture mapping because the method requires a dense map of validated markers known to have large frequency differences between Europeans and Africans. To create such a map, we screened through databases containing approximately 450000 single-nucleotide polymorphisms (SNPs) for which frequencies had been estimated in African and European population samples. We experimentally confirmed the frequencies of the most promising SNPs in a multiethnic panel of unrelated samples and identified 3011 as a MALD map (1.2 cM average spacing). We estimate that this map is approximately 70% informative in differentiating African versus European origins of chromosomal segments. This map provides a practical and powerful tool, which is freely available without restriction, for screening for disease genes in African American patient cohorts. The map is especially appropriate for those diseases that differ in incidence between the parental African and European populations.

High-density map of short tandem repeats across the human major histocompatibility complex.

The human genome contains one short tandem repeat (STR) roughly every 2,000 base pairs. They are particularly useful markers for gene mapping and disease association studies due to their high degree of polymorphism and ubiquitous frequency throughout the genome. The major histocompatibility complex (MHC) has been the focus of many disease association studies, and the recent availability of the entire sequence of the complex has logarithmically expanded the density of potential markers for fine mapping disease loci. Here we present a complete assessment of the available STRs within a 3.8-Mb genomic segment encompassing the MHC. Of 443 potential STRs identified by computer analysis and tested for variation in a single sample containing pooled DNA from 36 individuals, 249 polymorphic STRs located throughout the complex were identified. The class of repeat (di-, tri-, etc.), precise nucleotide position, position relative to known genes, PCR conditions, and D6S numbers for the 249 polymorphic STRs are provided as a resource for selecting appropriate markers to use in future studies of MHC molecular genetics and disease association.