Development and Validation of a Breast Cancer Polygenic Risk Score on the Basis of Genetic Ancestry Composition.

PURPOSE: Polygenic risk scores (PRSs) for breast cancer (BC) risk stratification have been developed primarily in women of European ancestry. Their application to women of non-European ancestry has lagged because of the lack of a formal approach to incorporate genetic ancestry and ancestry-dependent variant frequencies and effect sizes. Here, we propose a multiple-ancestry PRS (MA-PRS) that addresses these issues and may be useful in the development of equitable PRSs across other cancers and common diseases. MATERIALS AND METHODS: Women referred for hereditary cancer testing were divided into consecutive cohorts for development (n = 189,230) and for independent validation (n = 89,126). Individual genetic composition as fractions of three reference ancestries (African, East Asian, and European) was determined from ancestry-informative single-nucleotide polymorphisms. The MA-PRS is a combination of three ancestry-specific PRSs on the basis of genetic ancestral composition. Stratification of risk was evaluated by multivariable logistic regression models controlling for family cancer history. Goodness-of-fit analysis compared expected with observed relative risks by quantiles of the MA-PRS distribution. RESULTS: In independent validation, the MA-PRS was significantly associated with BC risk in the full cohort (odds ratio, 1.43; 95% CI, 1.40 to 1.46; P = 8.6 × 10-308) and within each major ancestry. The top decile of the MA-PRS consistently identified patients with two-fold increased risk of developing BC. Goodness-of-fit tests showed that the MA-PRS was well calibrated and predicted BC risk accurately in the tails of the distribution for both European and non-European women. CONCLUSION: The MA-PRS uses genetic ancestral composition to expand the utility of polygenic risk prediction to non-European women. Inclusion of genetic ancestry in polygenic risk prediction presents an opportunity for more personalized treatment decisions for women of varying and mixed ancestries.

Homologous recombination deficiency (HRD) status predicts response to standard neoadjuvant chemotherapy in patients with triple-negative or BRCA1/2 mutation-associated breast cancer.

PURPOSE: Defects in the homologous recombination (HR) DNA repair pathway sensitize tumors to therapeutics that target this pathway. A significant proportion of triple-negative breast cancers (TNBC) carry HR defects. The HRD assay is highly associated with sensitivity to neoadjuvant platinum-based chemotherapy in TNBC. Standard chemotherapy consists of some combination of an anthracycline, cyclophosphamide, and taxane. This study assesses the association of HR deficiency status with response to standard neoadjuvant chemotherapy in TNBC or BRCA1/2 mutation-associated breast cancer. METHODS: Tumor samples were retrospectively obtained from 45 TNBC patients and 2 BRCA1/2 mutant, hormone receptor-positive/HER2-negative breast cancer patients who received anthracycline- and/or taxane-based neoadjuvant chemotherapy at Stanford University or Cedars-Sinai Medical Centers. The HRD score and tumor BRCA1/2 mutation status were determined from baseline tumor biopsies. HR deficient tumors were those with a HRD score of ≥ 42 or a tumor BRCA1/2 mutation. Response was categorized by the residual cancer burden (RCB) index. RESULTS: HR deficient patients were more likely to achieve a pathologic complete response (pCR) compared with non-deficient patients (OR 13.06, CI 1.52-11.241, p = 0.0028). Among BRCA1/2 mutation wild-type patients, HR deficient patients were more likely to achieve a pCR (OR 16, 95% CI 1.65-160.41, p = 0.0041) compared with HR non-deficient patients. Further, HRD scores were highly concordant pre- and post-therapy (Spearman correlation > 99%). CONCLUSIONS: HR deficiency status is significantly associated with response to standard neoadjuvant chemotherapy in TNBC. This observation is consistent with the mechanisms of action of doxorubicin and cyclophosphamide as DNA damaging agents.

Homologous Recombination Deficiency (HRD) Score Predicts Response to Platinum-Containing Neoadjuvant Chemotherapy in Patients with Triple-Negative Breast Cancer.

PURPOSE: BRCA1/2-mutated and some sporadic triple-negative breast cancers (TNBC) have DNA repair defects and are sensitive to DNA-damaging therapeutics. Recently, three independent DNA-based measures of genomic instability were developed on the basis of loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale state transitions (LST). EXPERIMENTAL DESIGN: We assessed a combined homologous recombination deficiency (HRD) score, an unweighted sum of LOH, TAI, and LST scores, in three neoadjuvant TNBC trials of platinum-containing therapy. We then tested the association of HR deficiency, defined as HRD score ≥42 or BRCA1/2 mutation, with response to platinum-based therapy. RESULTS: In a trial of neoadjuvant platinum, gemcitabine, and iniparib, HR deficiency predicted residual cancer burden score of 0 or I (RCB 0/I) and pathologic complete response (pCR; OR = 4.96, P = 0.0036; OR = 6.52, P = 0.0058). HR deficiency remained a significant predictor of RCB 0/I when adjusted for clinical variables (OR = 5.86, P = 0.012). In two other trials of neoadjuvant cisplatin therapy, HR deficiency predicted RCB 0/I and pCR (OR = 10.18, P = 0.0011; OR = 17.00, P = 0.0066). In a multivariable model of RCB 0/I, HR deficiency retained significance when clinical variables were included (OR = 12.08, P = 0.0017). When restricted to BRCA1/2 nonmutated tumors, response was higher in patients with high HRD scores: RCB 0/I P = 0.062, pCR P = 0.063 in the neoadjuvant platinum, gemcitabine, and iniparib trial; RCB 0/I P = 0.0039, pCR P = 0.018 in the neoadjuvant cisplatin trials. CONCLUSIONS: HR deficiency identifies TNBC tumors, including BRCA1/2 nonmutated tumors more likely to respond to platinum-containing therapy. Clin Cancer Res; 22(15); 3764-73. ©2016 AACR.

Phase II neoadjuvant clinical trial of carboplatin and eribulin in women with triple negative early-stage breast cancer (NCT01372579).

The purpose of this study is to evaluate the efficacy and safety of neoadjuvant treatment with carboplatin and eribulin in patients with early-stage triple negative breast cancer (TNBC), and to explore biomarkers based on DNA and protein expression profiles as predictors of response. Patients with histologically confirmed early-stage TNBC received carboplatin AUC 6 iv every 21 days, and eribulin 1.4 mg/m(2) day 1 and day 8 every 21 days for four cycles. The primary endpoint of the study was pathologic complete response (pCR), with secondary endpoints including clinical response and safety of the combination. Exploratory studies assessed DNA-based biomarkers [homologous recombination deficiency (HRD) score, and HR deficiency status (HRD score + BRCA1/BRCA2 mutation status)], protein-based biomarkers (Ki67, TP53, androgen receptor, Cyclin E, CDK2, Cyclin D, CDK4, Pin1 and Smad3), and clinical pretreatment factors as predictors of pCR. 13/30 (43.3 %) patients enrolled in the study achieved pCR. 24 (80.0 %) had a clinical complete or partial response. The combination was safe with mostly grade 1 and 2 toxicities. HRD score (P = 0.0024) and HR deficiency status (P = 0.0012) significantly predicted pCR. Pretreatment cytoplasmic CDK2 was also associated with pCR (P = 0.021). Significant differences in pre- versus post-treatment expression levels of nuclear Cyclin D (P = 0.020), nuclear CDK4 (P = 0.0030), and nuclear Smad3 (P = 0.015) were detected. The combination of carboplatin and eribulin is safe and efficacious in the treatment of early-stage TNBC. HRD score, HR deficiency status, and cytoplasmic CDK2 predicted pCR in this patient population.

Phase II Study of Gemcitabine, Carboplatin, and Iniparib As Neoadjuvant Therapy for Triple-Negative and BRCA1/2 Mutation-Associated Breast Cancer With Assessment of a Tumor-Based Measure of Genomic Instability: PrECOG 0105.

PURPOSE: This study was designed to assess efficacy, safety, and predictors of response to iniparib in combination with gemcitabine and carboplatin in early-stage triple-negative and BRCA1/2 mutation-associated breast cancer. PATIENTS AND METHODS: This single-arm phase II study enrolled patients with stage I to IIIA (T ≥ 1 cm) estrogen receptor-negative (≤ 5%), progesterone receptor-negative (≤ 5%), and human epidermal growth factor receptor 2-negative or BRCA1/2 mutation-associated breast cancer. Neoadjuvant gemcitabine (1,000 mg/m(2) intravenously [IV] on days 1 and 8), carboplatin (area under curve of 2 IV on days 1 and 8), and iniparib (5.6 mg/kg IV on days 1, 4, 8, and 11) were administered every 21 days for four cycles, until the protocol was amended to six cycles. The primary end point was pathologic complete response (no invasive carcinoma in breast or axilla). All patients underwent comprehensive BRCA1/2 genotyping, and homologous recombination deficiency was assessed by loss of heterozygosity (HRD-LOH) in pretreatment core breast biopsies. RESULTS: Among 80 patients, median age was 48 years; 19 patients (24%) had germline BRCA1 or BRCA2 mutations; clinical stage was I (13%), IIA (36%), IIB (36%), and IIIA (15%). Overall pathologic complete response rate in the intent-to-treat population (n = 80) was 36% (90% CI, 27 to 46). Mean HRD-LOH scores were higher in responders compared with nonresponders (P = .02) and remained significant when BRCA1/2 germline mutations carriers were excluded (P = .021). CONCLUSION: Preoperative combination of gemcitabine, carboplatin, and iniparib is active in the treatment of early-stage triple-negative and BRCA1/2 mutation-associated breast cancer. The HRD-LOH assay was able to identify patients with sporadic triple-negative breast cancer lacking a BRCA1/2 mutation, but with an elevated HRD-LOH score, who achieved a favorable pathologic response. Confirmatory controlled trials are warranted.

Association of BRCA1/2 defects with genomic scores predictive of DNA damage repair deficiency among breast cancer subtypes.

INTRODUCTION: Homologous recombination (HR) DNA repair is of clinical relevance in breast cancer. Three DNA-based homologous recombination deficiency (HRD) scores (HRD-loss of heterozygosity score (LOH), HRD-telomeric allelic imbalance score (TAI), and HRD-large-scale state transition score (LST)) have been developed that are highly correlated with defects in BRCA1/2, and are associated with response to platinum therapy in triple negative breast and ovarian cancer. This study examines the frequency of BRCA1/2 defects among different breast cancer subtypes, and the ability of the HRD scores to identify breast tumors with defects in the homologous recombination DNA repair pathway. METHODS: 215 breast tumors representing all ER/HER2 subtypes were obtained from commercial vendors. Next-generation sequencing based assays were used to generate genome wide SNP profiles, BRCA1/2 mutation screening, and BRCA1 promoter methylation data. RESULTS: BRCA1/2 deleterious mutations were observed in all breast cancer subtypes. BRCA1 promoter methylation was observed almost exclusively in triple negative breast cancer. BRCA1/2 deficient tumors were identified with BRCA1/2 mutations, or BRCA1 promoter methylation, and loss of the second allele of the affected gene. All three HRD scores were highly associated with BRCA1/2 deficiency (HRD-LOH: P = 1.3 × 10(-17); HRD-TAI: P = 1.5 × 10(-19); HRD-LST: P = 3.5 × 10(-18)). A combined score (HRD-mean) was calculated using the arithmetic mean of the three scores. In multivariable analyses the HRD-mean score captured significant BRCA1/2 deficiency information not captured by the three individual scores, or by clinical variables (P values for HRD-Mean adjusted for HRD-LOH: P = 1.4 × 10(-8); HRD-TAI: P = 2.9 × 10(-7); HRD-LST: P = 2.8 × 10(-8); clinical variables: P = 1.2 × 10(-16)). CONCLUSIONS: The HRD scores showed strong correlation with BRCA1/2 deficiency regardless of breast cancer subtype. The frequency of elevated scores suggests that a significant proportion of all breast tumor subtypes may carry defects in the homologous recombination DNA repair pathway. The HRD scores can be combined to produce a more robust predictor of HRD. The combination of a robust score, and the FFPE compatible assay described in this study, may facilitate use of agents targeting homologous recombination DNA repair in the clinical setting.

Mapping of homozygous deletions in verified esophageal adenocarcinoma cell lines and xenografts.

Human esophageal adenocarcinoma (EAC) cell lines and xenografts are powerful tools in the search for genetic alterations because these models are composed of pure human cancer cell populations without admixture of normal human cells. In particular detection of homozygous deletions (HDs) is easier using these pure populations of cancer cells. Identification of HDs could potentially lead to the subsequent identification of new tumor suppressor genes (TSGs) involved in esophageal adenocarcinogenesis. Genome wide single nucleotide polymorphism (SNP) arrays were used to identify HDs in 10 verified EAC cell lines and nine EAC xenografts. In total, 61 HDs (range 1-6 per sample) were detected and confirmed by polymerase chain reaction. Besides HDs observed in common fragile genomic regions (n = 26), and gene deserts (n = 8), 27 HDs were located in gene-containing regions. HDs were noted for known TSGs, including CDKN2A, SMAD4 and CDH3/CDH1. Twenty-two new chromosomal regions were detected harboring potentially new TSGs involved in EAC carcinogenesis. Two of these regions of homozygous loss, encompassing the ITGAV and RUNX1 gene, were detected in multiple samples indicating a potential role in the carcinogenesis of EAC. To exclude culturing artifacts, these last two deletions were confirmed by fluorescent in situ hybridization in the primary tumors of which the involved cell lines and xenografts were derived. In summary, in this report we describe the identification of HDs in a series of verified EAC cell lines and xenografts. The deletions documented here are a step forward identifying the key genes involved in EAC development.

Computational method for estimating DNA copy numbers in normal samples, cancer cell lines, and solid tumors using array comparative genomic hybridization.

Genomic copy number variations are a typical feature of cancer. These variations may influence cancer outcomes as well as effectiveness of treatment. There are many computational methods developed to detect regions with deletions and amplifications without estimating actual copy numbers (CN) in these regions. We have developed a computational method capable of detecting regions with deletions and amplifications as well as estimating actual copy numbers in these regions. The method is based on determining how signal intensity from different probes is related to CN, taking into account changes in the total genome size, and incorporating into analysis contamination of the solid tumors with benign tissue. Hidden Markov Model is used to obtain the most likely CN solution. The method has been implemented for Affymetrix 500K GeneChip arrays and Agilent 244K oligonucleotide arrays. The results of CN analysis for normal cell lines, cancer cell lines, and tumor samples are presented. The method is capable of detecting copy number alterations in tumor samples with up to 80% contamination with benign tissue. Analysis of 178 cancer cell lines reveals multiple regions of common homozygous deletions and strong amplifications encompassing known tumor suppressor genes and oncogenes as well as novel cancer related genes.

Somatic mutations in BRCA1 and BRCA2 could expand the number of patients that benefit from poly (ADP ribose) polymerase inhibitors in ovarian cancer.

PURPOSE: The prevalence of BRCA(1/2) mutations in germline DNA from unselected ovarian cancer patients is 11% to 15.3%. It is important to determine the frequency of somatic BRCA(1/2) changes, given the sensitivity of BRCA-mutated cancers to poly (ADP ribose) polymerase-1 (PARP1) inhibitors and platinum analogs. PATIENTS AND METHODS: In 235 unselected ovarian cancers, BRCA(1/2) was sequenced in 235, assessed by copy number analysis in 95, and tiling arrays in 65. 113 tumors were sequenced for TP53. BRCA(1/2) transcript levels were assessed by quantitative polymerase chain reaction in 220. When available for tumors with BRCA(1/2) mutations, germline DNA was sequenced. RESULTS: Forty-four mutations (19%) in BRCA1 (n = 31)/BRCA2 (n = 13) were detected, including one homozygous BRCA1 intragenic deletion. BRCA(1/2) mutations were particularly common (23%) in high-grade serous cancers. In 28 patients with available germline DNA, nine (42.9%) of 21 and two (28.6%) of seven BRCA1 and BRCA2 mutations were demonstrated to be somatic, respectively. Five mutations not previously identified in germline DNA were more commonly somatic than germline (four of 11 v one of 17; P = .062). There was a positive association between BRCA1 and TP53 mutations (P = .012). BRCA(1/2) mutations were associated with improved progression-free survival (PFS) after platinum-based chemotherapy in univariate (P = .032; hazard ratio [HR] = 0.65; 95% CI, 0.43 to 0.98) and multivariate (P = .019) analyses. BRCA(1/2) deficiency, defined as BRCA(1/2) mutations or expression loss (in 24 [13.3%] BRCA(1/2)-wild-type cancers), was present in 67 ovarian cancers (30%) and was also significantly associated with PFS in univariate (P = .026; HR = 0.67; 95% CI, 0.47 to 0.96) and multivariate (P = .008) analyses. CONCLUSION: BRCA(1/2) somatic and germline mutations and expression loss are sufficiently common in ovarian cancer to warrant assessment for prediction of benefit in clinical trials of PARP1 inhibitors.

Evidence for epistasis between SLC6A4 and a chromosome 4 gene as risk factors in major depression.

Linkage analysis on Utah pedigrees with strong family histories of major depression including only cases with the SLC6A4 HTTLPR short allele revealed a linkage peak on chromosome 4 (maximum HLOD = 3.5). This evidence suggests epistasis between SLC6A4 and an unknown gene as risk factors for major depression.

Genome-wide association scan identifies a prostaglandin-endoperoxide synthase 2 variant involved in risk of knee osteoarthritis.

Osteoarthritis (OA), the most prevalent form of arthritis in the elderly, is characterized by the degradation of articular cartilage and has a strong genetic component. Our aim was to identify genetic variants involved in risk of knee OA in women. A pooled genome-wide association scan with the Illumina550 Duo array was performed in 255 controls and 387 cases. Twenty-eight variants with p < 1 x 10(-5) were estimated to have probabilities of being false positives

Identification of ZNF313/RNF114 as a novel psoriasis susceptibility gene.

Psoriasis is an immune-mediated skin disorder that is inherited as a multifactorial trait. Linkage studies have clearly identified a primary disease susceptibility locus lying within the major histocompatibility complex (MHC), but have generated conflicting results for other genomic regions. To overcome this difficulty, we have carried out a genome-wide association scan, where we analyzed more than 408,000 SNPs in an initial sample of 318 cases and 288 controls. Outside of the MHC, we observed a single cluster of disease-associated markers, spanning 47 kb on chromosome 20q13. The analysis of two replication data sets confirmed this association, with SNP rs495337 yielding a combined P-value of 1.4 x 10(-8) in an overall sample of 2679 cases and 2215 controls. Rs495337 maps to the SPATA2 transcript and is in absolute linkage disequilibrium with five SNPs lying in the adjacent ZNF313 gene (also known as RNF114). Real-time PCR experiments showed that, unlike SPATA2, ZNF313 is abundantly expressed in skin, T-lymphocytes and dendritic cells. Furthermore, an analysis of the expression data available from the Genevar database indicated that rs495337 is associated with increased ZNF313 transcripts levels (P = 0.003), suggesting that the disease susceptibility allele may be a ZNF313 regulatory variant tagged by rs495337. Homology searches indicated that ZNF313 is a paralogue of TRAC-1, an ubiquitin ligase regulating T-cell activation. We performed cell-free assays and confirmed that like TRAC-1, ZNF313 binds ubiquitin via an ubiquitin-interaction motif (UIM). These findings collectively identify a novel psoriasis susceptibility gene, with a putative role in the regulation of immune responses.

IL23R R381Q and ATG16L1 T300A are strongly associated with Crohn's disease in a study of New Zealand Caucasians with inflammatory bowel disease.

OBJECTIVE: Recently, separate genome-wide association analyses have identified nonsynonymous SNPs in IL23R and ATG16L1 (rs11209026; c1142G>A, R381Q, and rs2241880; c1338A>G, T300A, respectively) as strong candidate susceptibility factors for Crohn's disease (CD) in whites. The aim of our study was to test whether these SNPs are associated with CD in a population-based cohort of New Zealand Caucasian inflammatory bowel disease (IBD) patients. METHODS: Allele frequencies of rs11209026 and rs2241880 were determined in 496 CD patients, 466 ulcerative colitis (UC) patients, and 591 controls. Distribution of the relevant alleles was compared between controls and IBD patients. rs11209026 and rs2241880 genotype distributions were examined both within IBD clinical subphenotypes and CARD15 genotypes. RESULTS: rs11209026 and rs2241880 were both associated with CD (P valuers11209026=0.0026, OR 0.54, 95% CI 0.36-0.81; P valuers2241880=0.0001, OR 1.41, 95% CI 1.18-1.67). In addition, there was evidence for association of rs11209026 with UC (P value=0.037, OR 0.66, 95% CI 0.45-0.98). No significant association was observed between IL23R genotype or ATG16L1 genotype and IBD subphenotypes. IL23R was associated with CD and UC only in the absence of CARD15 mutations, whereas ATG16L1 was associated with CD in the presence and absence of CARD15 mutations. CONCLUSIONS: We replicated the previously reported associations between CD and rs11209026 and rs2241880, confirming that IL23R and ATG16L1 are susceptibility loci for CD in the New Zealand population. We also provide further evidence for association of rs11209026 with UC and a report of an additive effect between IL23R and CARD15 genotypes in CD.

Meta-analysis of genome-wide linkage studies in BMI and obesity.

OBJECTIVE: The objective was to provide an overall assessment of genetic linkage data of BMI and BMI-defined obesity using a nonparametric genome scan meta-analysis. RESEARCH METHODS AND PROCEDURES: We identified 37 published studies containing data on over 31,000 individuals from more than >10,000 families and obtained genome-wide logarithm of the odds (LOD) scores, non-parametric linkage (NPL) scores, or maximum likelihood scores (MLS). BMI was analyzed in a pooled set of all studies, as a subgroup of 10 studies that used BMI-defined obesity, and for subgroups ascertained through type 2 diabetes, hypertension, or subjects of European ancestry. RESULTS: Bins at chromosome 13q13.2- q33.1, 12q23-q24.3 achieved suggestive evidence of linkage to BMI in the pooled analysis and samples ascertained for hypertension. Nominal evidence of linkage to these regions and suggestive evidence for 11q13.3-22.3 were also observed for BMI-defined obesity. The FTO obesity gene locus at 16q12.2 also showed nominal evidence for linkage. However, overall distribution of summed rank p values <0.05 is not different from that expected by chance. The strongest evidence was obtained in the families ascertained for hypertension at 9q31.1-qter and 12p11.21-q23 (p < 0.01). CONCLUSION: Despite having substantial statistical power, we did not unequivocally implicate specific loci for BMI or obesity. This may be because genes influencing adiposity are of very small effect, with substantial genetic heterogeneity and variable dependence on environmental factors. However, the observation that the FTO gene maps to one of the highest ranking bins for obesity is interesting and, while not a validation of this approach, indicates that other potential loci identified in this study should be investigated further.

Genetic study of the melanin-concentrating hormone receptor 2 in childhood and adulthood severe obesity.

CONTEXT: The melanin-concentrating hormone receptor 2 (MCHR2) is a G protein-coupled receptor for melanin-concentrating hormone, a neuropeptide that plays an important role in feeding behaviors. MCHR2 maps on chromosome 6q16.3, in a susceptibility locus for childhood obesity. OBJECTIVE: The aim of this study was to investigate the association between MCHR2 variation and human obesity. DESIGN: Case control and family-based studies were performed. PARTICIPANTS: A total of 141 obese children and 24 nonobese adult subjects was sequenced, and case-control analyses were conducted using 628 severely obese children and 1,401 controls. RESULTS: There were 11 single nucleotide polymorphisms (SNPs) identified. We showed nominal association among -38,245 ATG A/G SNP (P = 0.03; 95% confidence interval 1.02-1.34; odds ratio 1.17), A76A T/C SNP (P = 0.03; 95% confidence interval 0.58-0.97; odds ratio 0.75), and childhood obesity. Analysis of 645 trios with childhood obesity supported further the A76A T/C association, showing an overtransmission to obese children of the at risk T allele (59.0%; P = 0.01), especially in children with most severe forms of obesity (Z score of body mass index > 4) (67.0%; P = 0.003). The A76A at risk T allele was also associated with overeating during meals (P = 0.02) in an additional group of 102 nonobese children. None of the MCHR2 variants, including the A76A SNP, showed association with adult severe obesity, although a trend for association of the T allele of this variant with food disinhibition (P = 0.06) and higher hunger (P = 0.09) was found. This variant was not associated with childhood obesity in an independent case-control study, including 1,573 subjects (P = 0.98). Moreover, the A76A SNP did not explain the linkage on the 6q locus. CONCLUSION: Our results altogether suggest that MCHR2 is not a major contributor to polygenic obesity and support a modest effect of the A76A SNP on food intake abnormalities in childhood.

TBC1D1 is a candidate for a severe obesity gene and evidence for a gene/gene interaction in obesity predisposition.

The molecular etiology of obesity predisposition is largely unknown. Here, we present evidence that genetic variation in TBC1D1 confers risk for severe obesity in females. We identified a coding variant (R125W) in TBC1D1 that segregated with the disease in 4p15-14-linked obesity pedigrees. In cases derived from pedigrees with the strongest linkage evidence, the variant was significantly associated with obesity (P=0.000007) and chromosomes carrying R125W accounted for the majority of the evidence that originally linked 4p15-14 with the disease. In addition, by selecting families that segregated R125W with obesity, we were able to generate highly significant linkage evidence for an obesity predisposition locus at 4q34-35. This result provides additional and confirming evidence that R125W affects obesity susceptibility, delimits the location of an obesity gene at 4q34-35 and identifies a gene/gene interaction that influences the risk for obesity predisposition. Finally, although the function of TBC1D1 is unknown, the protein is structurally similar to a known regulator of insulin-mediated Glut4 translocation.

Genome-wide linkage analyses of extended Utah pedigrees identifies loci that influence recurrent, early-onset major depression and anxiety disorders.

Major depressive disorder (MDD) is a common, clinically heterogeneous disorder often found comorbid with other disorders. We studied recurrent, early-onset MDD (MDD-RE) and anxiety disorders in combination to define powerful phenotypes for genetic study. We used 87 large, extended Utah pedigrees to investigate linkage to 3 phenotypes: "MDD-RE;" "MDD-RE or anxiety;" and "MDD-RE and anxiety;" where in the latter definition the disorders must appear comorbid within an individual. Pedigrees ranged in size from 2 to 6 generations and contained 3 to 42 individuals affected with MDD or anxiety (718 total). In primary analyses, we identified three regions with at least suggestive genome-wide evidence for linkage on chromosomes 3centr, 7p, and 18q. Both 7p and 18q are replication findings for related phenotypes. The best linkage evidence was for a novel locus at 3p12.3-q12.3 (LOD = 3.88, "MDD-RE or anxiety") and 18q21.33-q22.2 (LOD = 3.75, "MDD-RE and anxiety"), a well-established susceptibility locus for bipolar disorder. In our secondary sex-specific analyses, we identified two further regions of interest on chromosomes 4q and 15q. Using linked pedigrees, we localized 3centr and 18q to 9.8 and 12.2 cM, respectively, with potential for further localization with the addition of markers in specific pedigrees. Our success in replication and novel locus identification illustrates the utility of large extended pedigrees for common disorders, such as MDD. Further, it supports the hypothesis that MDD and anxiety disorders have over-lapping genetic etiologies and suggests that comorbid diagnoses may be useful in defining more genetically homogeneous forms of MDD for linkage mapping.

Predisposition locus for major depression at chromosome 12q22-12q23.2.

Major depression disorder is a common psychiatric disease with a major economic impact on society. In many cases, no effective treatment is available. The etiology of major depression is complex, but it is clear that the disease is, to a large extent, determined genetically, especially among individuals with a familial history of major depression, presumably through the involvement of multiple predisposition genes in addition to an environmental component. As a first step toward identification of chromosomal loci contributing to genetic predisposition to major depression, we have conducted a genomewide scan by using 628 microsatellite markers on 1,890 individuals from 110 Utah pedigrees with a strong family history of major depression. We identified significant linkage to major depression in males at marker D12S1300 (multipoint heterogeneity LOD score 4.6; P=.00003 after adjustment for multiple testing). With additional markers, the linkage evidence became highly significant, with the multipoint heterogeneity LOD score at marker D12S1706 increasing to 6.1 (P=.0000007 after adjustment for multiple testing). This study confirms the presence of one or more genes involved in psychiatric diseases on the q arm of chromosome 12 and provides strong evidence for the existence of a sex-specific predisposition gene to major depression at 12q22-q23.2.

A genome scan for loci influencing anti-atherogenic serum bilirubin levels.

Epidemiological studies have shown an association of decreased serum bilirubin levels with coronary artery disease. Two segregation analyses in large pedigrees have suggested a major gene responsible for high bilirubin levels occurring in about 12% of the population. Based on a recessive model from a previous segregation analysis, we performed a genome scan using 587 markers genotyped in 862 individuals from 48 Utah pedigrees to detect loci linked to high bilirubin levels. As a complementary approach, non-parametric linkage (NPL) analysis was performed. These two methods identified four regions showing evidence for linkage. The first region is on chromosome 2q34-37 with multipoint LOD and NPL scores of 3.01 and 3.22, respectively, for marker D2S1363. This region contains a previously described gene, uridine diphosphate glycosyltransferase 1, which has been associated with high bilirubin levels. A polymorphism in the promoter of this gene was recently shown to be responsible for Gilbert syndrome which is associated with mild hyperbilirubinemia. The other regions were found on chromosomes 9q21, 10q25-26, and 18q12 with maximum NPL scores of 2.39, 1.55, and 2.79, respectively. Furthermore, we investigated in these pedigrees the association between bilirubin levels and coronary artery disease. One-hundred and sixty-one male and 41 female subjects had already suffered a coronary artery disease event. Male patients showed significantly lower bilirubin concentrations than age-matched controls. This association, however, was not observed in females. These results provide evidence that loci influencing bilirubin variation exist on chromosomes 2q34-37, 9q21, 10q25-26, and 18q12 and confirms the association of low bilirubin levels with coronary artery disease in males.