The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

Ancestry adjustment improves genome-wide estimates of regional intolerance.

Genomic regions subject to purifying selection are more likely to carry disease-causing mutations than regions not under selection. Cross species conservation is often used to identify such regions but with limited resolution to detect selection on short evolutionary timescales such as that occurring in only one species. In contrast, genetic intolerance looks for depletion of variation relative to expectation within a species, allowing species-specific features to be identified. When estimating the intolerance of noncoding sequence, methods strongly leverage variant frequency distributions. As the expected distributions depend on ancestry, if not properly controlled for, ancestral population source may obfuscate signals of selection. We demonstrate that properly incorporating ancestry in intolerance estimation greatly improved variant classification. We provide a genome-wide intolerance map that is conditional on ancestry and likely to be particularly valuable for variant prioritization.

A Stat1 bound enhancer promotes Nampt expression and function within tumor associated macrophages.

Tumor associated macrophage responses are regulated by distinct metabolic states that affect their function. However, the ability of specific signals in the local tumor microenvironment to program macrophage metabolism remains under investigation. Here, we identify NAMPT, the rate limiting enzyme in NAD salvage synthesis, as a target of STAT1 during cellular activation by interferon gamma, an important driver of macrophage polarization and antitumor responses. We demonstrate that STAT1 occupies a conserved element within the first intron of Nampt, termed Nampt-Regulatory Element-1 (NRE1). Through disruption of NRE1 or pharmacological inhibition, a subset of M1 genes is sensitive to NAMPT activity through its impact on glycolytic processes. scRNAseq is used to profile in vivo responses by NRE1-deficient, tumor-associated leukocytes in melanoma tumors through the creation of a unique mouse strain. Reduced Nampt and inflammatory gene expression are present in specific myeloid and APC populations; moreover, targeted ablation of NRE1 in macrophage lineages results in greater tumor burden. Finally, elevated NAMPT expression correlates with IFNγ responses and melanoma patient survival. This study identifies IFN and STAT1-inducible Nampt as an important factor that shapes the metabolic program and function of tumor associated macrophages.

Genomic characterization of MICA gene using multiple next generation sequencing platforms: A validation study.

We have developed a protocol regarding the genomic characterization of the MICA gene by next generation sequencing (NGS). The amplicon includes the full length of the gene and is about 13 kb. A total of 156 samples were included in the study. Ninety-seven of these samples were previously characterized at MICA by legacy methods (Sanger or sequence specific oligonucleotide) and were used to evaluate the accuracy, precision, specificity, and sensitivity of the assay. An additional 59 DNA samples of unknown ethnicity volunteers from the United States were only genotyped by NGS. Samples were chosen to contain a diverse set of alleles. Our NGS approach included a first round of sequencing on the Illumina MiSeq platform and a second round of sequencing on the MinION platform by Oxford Nanopore Technology (ONT), on selected samples for the purpose of either characterizing new alleles or setting phase among multiple polymorphisms to resolve ambiguities or generate complete sequence for alleles that were only partially reported in the IMGT/HLA database. Complete consensus sequences were generated for every allele sequenced with ONT, extending from the 5' untranslated region (UTR) to the 3' UTR of the MICA gene. Thirty-two MICA sequences were submitted to the IMGT/HLA database including either new alleles or filling up the gaps (exonic, intronic and/or UTRs) of already reported alleles. Some of the challenges associated with the characterization of these samples are discussed.

Utilizing nanopore sequencing technology for the rapid and comprehensive characterization of eleven HLA loci; addressing the need for deceased donor expedited HLA typing.

The comprehensive characterization of human leukocyte antigen (HLA) genomic sequences remains a challenging problem. Despite the significant advantages of next-generation sequencing (NGS) in the field of Immunogenetics, there has yet to be a single solution for unambiguous, accurate, simple, cost-effective, and timely genotyping necessary for all clinical applications. This report demonstrates the benefits of nanopore sequencing introduced by Oxford Nanopore Technologies (ONT) for HLA genotyping. Samples (n = 120) previously characterized at high-resolution three-field (HR-3F) for 11 loci were assessed using ONT sequencing paired to a single-plex PCR protocol (Holotype) and to two multiplex protocols OmniType (Omixon) and NGSgo®-MX6-1 (GenDx). The results demonstrate the potential of nanopore sequencing for delivering accurate HR-3F typing with a simple, rapid, and cost-effective protocol. The protocol is applicable to time-sensitive applications, such as deceased donor typings, enabling better assessments of compatibility and epitope analysis. The technology also allows significantly shorter turnaround time for multiple samples at a lower cost. Overall, the nanopore technology appears to offer a significant advancement over current next-generation sequencing platforms as a single solution for all HLA genotyping needs.

Resolving MiSeq-Generated Ambiguities in HLA-DPB1 Typing by Using the Oxford Nanopore Technology.

The technical limitations of current next-generation sequencing technologies, combined with an ever-increasing number of human leukocyte antigen (HLA) alleles, form the basis for the additional ambiguities encountered at an increasing rate in clinical practice. HLA-DPB1 characterization, particularly, generates a significant percentage of ambiguities (25.5%), posing a challenge for accurate and unambiguous HLA-DPB1 genotyping. Phasing of exonic heterozygous positions between exon 2 and all other downstream exons has been the major cause of ambiguities. In this study, the Oxford Nanopore MinION, a third-generation sequencing technology, was used to resolve the phasing. The accurate MiSeq sequencing data, combined with the long reads obtained from the MinION platform, allow for the resolution of the tested ambiguities.

Next-generation HLA typing of 382 International Histocompatibility Working Group reference B-lymphoblastoid cell lines: Report from the 17th International HLA and Immunogenetics Workshop.

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m5C, Ψ, and m1A at single-base resolution transcriptome-wide. With RBS-Seq, m5C and m1A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1-2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg2+-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m5C, m1A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m5C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m5C sites in noncoding and coding RNAs and the absence of substantial m1A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

FoxA1 and FoxA2 drive gastric differentiation and suppress squamous identity in NKX2-1-negative lung cancer.

Changes in cancer cell identity can alter malignant potential and therapeutic response. Loss of the pulmonary lineage specifier NKX2-1 augments the growth of KRAS-driven lung adenocarcinoma and causes pulmonary to gastric transdifferentiation. Here, we show that the transcription factors FoxA1 and FoxA2 are required for initiation of mucinous NKX2-1-negative lung adenocarcinomas in the mouse and for activation of their gastric differentiation program. Foxa1/2 deletion severely impairs tumor initiation and causes a proximal shift in cellular identity, yielding tumors expressing markers of the squamocolumnar junction of the gastrointestinal tract. In contrast, we observe downregulation of FoxA1/2 expression in the squamous component of both murine and human lung adenosquamous carcinoma. Using sequential in vivo recombination, we find that FoxA1/2 loss in established KRAS-driven neoplasia originating from SPC-positive alveolar cells induces keratinizing squamous cell carcinomas. Thus, NKX2-1, FoxA1 and FoxA2 coordinately regulate the growth and identity of lung cancer in a context-specific manner.

The Lineage-Defining Transcription Factors SOX2 and NKX2-1 Determine Lung Cancer Cell Fate and Shape the Tumor Immune Microenvironment.

The major types of non-small-cell lung cancer (NSCLC)-squamous cell carcinoma and adenocarcinoma-have distinct immune microenvironments. We developed a genetic model of squamous NSCLC on the basis of overexpression of the transcription factor Sox2, which specifies lung basal cell fate, and loss of the tumor suppressor Lkb1 (SL mice). SL tumors recapitulated gene-expression and immune-infiltrate features of human squamous NSCLC; such features included enrichment of tumor-associated neutrophils (TANs) and decreased expression of NKX2-1, a transcriptional regulator that specifies alveolar cell fate. In Kras-driven adenocarcinomas, mis-expression of Sox2 or loss of Nkx2-1 led to TAN recruitment. TAN recruitment involved SOX2-mediated production of the chemokine CXCL5. Deletion of Nkx2-1 in SL mice (SNL) revealed that NKX2-1 suppresses SOX2-driven squamous tumorigenesis by repressing adeno-to-squamous transdifferentiation. Depletion of TANs in SNL mice reduced squamous tumors, suggesting that TANs foster squamous cell fate. Thus, lineage-defining transcription factors determine the tumor immune microenvironment, which in turn might impact the nature of the tumor.

Characterization of 108 Genomic DNA Reference Materials for 11 Human Leukocyte Antigen Loci: A GeT-RM Collaborative Project.

The highly polymorphic human leukocyte antigen (HLA) genes, located in the human major histocompatibility complex, encode the class I and II antigen-presenting molecules, which are centrally involved in the immune response. HLA typing is used for several clinical applications, such as transplantation, pharmacogenetics, and diagnosis of autoimmune disease. HLA typing is highly complex because of the homology of HLA genes and pseudogenes and the extensive polymorphism in the population. The Centers for Disease Control and Prevention established the Genetic Testing Reference Materials Coordination Program (GeT-RM) in partnership with the genetics community to improve the availability of genomic DNA reference materials necessary for quality assurance of genetic laboratory testing. The GeT-RM together with three clinical laboratories and the Coriell Cell Repositories have characterized genomic DNA obtained from a panel of 108 cell lines for all HLA classic polymorphic loci: HLA-A, B, C, DRB1, DRB3, DRB4, DRB5, DQA1, DQB1, DPA1, and DPB1. The goal was to develop a publicly available and renewable source of well-characterized genomic DNA reference materials to support molecular HLA typing assay development, validation, and verification, quality control, and proficiency testing. These genomic DNA samples are publicly available from the National Institutes of General Medical Science Repository at the Coriell Cell Repositories.

Reversible LSD1 inhibition with HCI-2509 induces the p53 gene expression signature and disrupts the MYCN signature in high-risk neuroblastoma cells.

Lysine-Specific Demethylase 1 (LSD1) over-expression correlates with poorly differentiated neuroblastoma and predicts poor outcome despite multimodal therapy. We have studied the efficacy of reversible and specific LSD1 inhibition with HCI-2509 in neuroblastoma cell lines and particularly the effect of HCI-2509 on the transcriptomic profile in MYCN amplified NGP cells. Cell survival assays show that HCI-2509 is cytotoxic to poorly differentiated neuroblastoma cell lines in low micromole or lower doses. Transcriptional profiling of NGP cells treated with HCI-2509 shows a significant effect on p53, cell cycle, MYCN and hypoxia pathway gene sets. HCI-2509 results in increased histone methyl marks and p53 levels along with cell cycle arrest in the G2/M phase and inhibition of colony formation of NGP cells. Our findings indicate that LSD1 inhibition with HCI-2509 has a multi-target effect in neuroblastoma cell lines, mediated in part via p53. MYCN-amplified neuroblastoma cells have a targeted benefit as HCI-2509 downregulates the MYCN upregulated gene set.

Active BRAF-V600E is the key player in generation of a sessile serrated polyp-specific DNA methylation profile.

BACKGROUND: Sessile serrated polyps (SSPs) have emerged as important precursors for a large number of sporadic colorectal cancers. They are difficult to detect during colonoscopy due to their flat shape and the excessive amounts of secreted mucin that cover the polyps. The underlying genetic and epigenetic basis for the emergence of SSPs is largely unknown with existing genetic studies confined to a limited number of oncogenes and tumor suppressors. A full characterization of the genetic and epigenetic landscape of SSPs would provide insight into their origin and potentially offer new biomarkers useful for detection of SSPs in stool samples. METHODS: We used a combination of genome-wide mutation detection, exome sequencing and DNA methylation profiling (via methyl-array and whole-genome bisulfite sequencing) to analyze multiple samples of sessile serrated polyps and compared these to familial adenomatous polyps. RESULTS: Our analysis revealed BRAF-V600E as the sole recurring somatic mutation in SSPs with no additional major genetic mutations detected. The occurrence of BRAF-V600E was coincident with a unique DNA methylation pattern revealing a set of DNA methylation markers showing significant (~3 to 30 fold) increase in their methylation levels, exclusively in SSP samples. These methylation patterns effectively distinguished sessile serrated polys from adenomatous polyps and did so more effectively than parallel gene expression profiles. CONCLUSIONS: This study provides an important example of a single oncogenic mutation leading to reproducible global DNA methylation changes. These methylated markers are specific to SSPs and could be of important clinical relevance for the early diagnosis of SSPs using non-invasive approaches such as fecal DNA testing.

miR-155 promotes FLT3-ITD-induced myeloproliferative disease through inhibition of the interferon response.

FLT3-ITD+ acute myeloid leukemia (AML) accounts for ∼25% of all AML cases and is a subtype that carries a poor prognosis. microRNA-155 (miR-155) is specifically overexpressed in FLT3-ITD+ AML compared with FLT3 wild-type (FLT3-WT) AML and is critical for the growth of FLT3-ITD+ AML cells in vitro. However, miR-155's role in regulating FLT3-ITD-mediated disease in vivo remains unclear. In this study, we used a genetic mouse model to determine whether miR-155 influences the development of FLT3-ITD-induced myeloproliferative disease. Results indicate that miR-155 promotes FLT3-ITD-induced myeloid expansion in the bone marrow, spleen, and peripheral blood. Mechanistically, miR-155 increases proliferation of the hematopoietic stem and progenitor cell compartments by reducing the growth-inhibitory effects of the interferon (IFN) response, and this involves targeting of Cebpb. Consistent with our observations in mice, primary FLT3-ITD+ AML clinical samples have significantly higher miR-155 levels and a lower IFN response compared with FLT3-WT AML samples. Further, inhibition of miR-155 in FLT3-ITD+ AML cell lines using CRISPR/Cas9, or primary FLT3-ITD+ AML samples using locked nucleic acid antisense inhibitors, results in an elevated IFN response and reduces colony formation. Altogether, our data reveal that miR-155 collaborates with FLT3-ITD to promote myeloid cell expansion in vivo and that this involves a multitarget mechanism that includes repression of IFN signaling.

MYC Drives Progression of Small Cell Lung Cancer to a Variant Neuroendocrine Subtype with Vulnerability to Aurora Kinase Inhibition.

Loss of the tumor suppressors RB1 and TP53 and MYC amplification are frequent oncogenic events in small cell lung cancer (SCLC). We show that Myc expression cooperates with Rb1 and Trp53 loss in the mouse lung to promote aggressive, highly metastatic tumors, that are initially sensitive to chemotherapy followed by relapse, similar to human SCLC. Importantly, MYC drives a neuroendocrine-low "variant" subset of SCLC with high NEUROD1 expression corresponding to transcriptional profiles of human SCLC. Targeted drug screening reveals that SCLC with high MYC expression is vulnerable to Aurora kinase inhibition, which, combined with chemotherapy, strongly suppresses tumor progression and increases survival. These data identify molecular features for patient stratification and uncover a potential targeted treatment approach for MYC-driven SCLC.

MEK Inhibitors Reverse Growth of Embryonal Brain Tumors Derived from Oligoneural Precursor Cells.

Malignant brain tumors are the leading cause of cancer-related deaths in children. Primitive neuroectodermal tumors of the CNS (CNS-PNETs) are particularly aggressive embryonal tumors of unknown cellular origin. Recent genomic studies have classified CNS-PNETs into molecularly distinct subgroups that promise to improve diagnosis and treatment; however, the lack of cell- or animal-based models for these subgroups prevents testing of rationally designed therapies. Here, we show that a subset of CNS-PNETs co-express oligoneural precursor cell (OPC) markers OLIG2 and SOX10 with coincident activation of the RAS/MAPK (mitogen-activated protein kinase) pathway. Modeling NRAS activation in embryonic OPCs generated malignant brain tumors in zebrafish that closely mimic the human oligoneural/NB-FOXR2 CNS-PNET subgroup by histology and comparative oncogenomics. The zebrafish CNS-PNET model was used to show that MEK inhibitors selectively eliminate Olig2+/Sox10+ CNS-PNET tumors in vivo without impacting normal brain development. Thus, MEK inhibitors represent a promising rationally designed therapy for children afflicted with oligoneural/NB-FOXR2 CNS-PNETs.

Genome-Wide CRISPR-Cas9 Screen Identifies MicroRNAs That Regulate Myeloid Leukemia Cell Growth.

Mammalian microRNA expression is dysregulated in human cancer. However, the functional relevance of many microRNAs in the context of tumor biology remains unclear. Using CRISPR-Cas9 technology, we performed a global loss-of-function screen to simultaneously test the functions of individual microRNAs and protein-coding genes during the growth of a myeloid leukemia cell line. This approach identified evolutionarily conserved human microRNAs that suppress or promote cell growth, revealing that microRNAs are extensively integrated into the molecular networks that control tumor cell physiology. miR-155 was identified as a top microRNA candidate promoting cellular fitness, which we confirmed with two distinct miR-155-targeting CRISPR-Cas9 lentiviral constructs. Further, we performed anti-correlation functional profiling to predict relevant microRNA-tumor suppressor gene or microRNA-oncogene interactions in these cells. This analysis identified miR-150 targeting of p53, a connection that was experimentally validated. Taken together, our study describes a powerful genetic approach by which the function of individual microRNAs can be assessed on a global level, and its use will rapidly advance our understanding of how microRNAs contribute to human disease.

Discordant Haplotype Sequencing Identifies Functional Variants at the 2q33 Breast Cancer Risk Locus.

The findings from genome-wide association studies hold enormous potential for novel insight into disease mechanisms. A major challenge in the field is to map these low-risk association signals to their underlying functional sequence variants (FSV). Simple sequence study designs are insufficient, as the vast numbers of statistically comparable variants and a limited knowledge of noncoding regulatory elements complicate prioritization. Furthermore, large sample sizes are typically required for adequate power to identify the initial association signals. One important question is whether similar sample sizes need to be sequenced to identify the FSVs. Here, we present a proof-of-principle example of an extreme discordant design to map FSVs within the 2q33 low-risk breast cancer locus. Our approach employed DNA sequencing of a small number of discordant haplotypes to efficiently identify candidate FSVs. Our results were consistent with those from a 2,000-fold larger, traditional imputation-based fine-mapping study. To prioritize further, we used expression-quantitative trait locus analysis of RNA sequencing from breast tissues, gene regulation annotations from the ENCODE consortium, and functional assays for differential enhancer activities. Notably, we implicate three regulatory variants at 2q33 that target CASP8 (rs3769823, rs3769821 in CASP8, and rs10197246 in ALS2CR12) as functionally relevant. We conclude that nested discordant haplotype sequencing is a promising approach to aid mapping of low-risk association loci. The ability to include more efficient sequencing designs into mapping efforts presents an opportunity for the field to capitalize on the potential of association loci and accelerate translation of association signals to their underlying FSVs. Cancer Res; 76(7); 1916-25. ©2016 AACR.