Methodology for Good Machine Learning with Multi-Omics Data.

In 2020, Novartis Pharmaceuticals Corporation and the U.S. Food and Drug Administration (FDA) started a 4-year scientific collaboration to approach complex new data modalities and advanced analytics. The scientific question was to find novel radio-genomics-based prognostic and predictive factors for HR+/HER- metastatic breast cancer under a Research Collaboration Agreement. This collaboration has been providing valuable insights to help successfully implement future scientific projects, particularly using artificial intelligence and machine learning. This tutorial aims to provide tangible guidelines for a multi-omics project that includes multidisciplinary expert teams, spanning across different institutions. We cover key ideas, such as "maintaining effective communication" and "following good data science practices," followed by the four steps of exploratory projects, namely (1) plan, (2) design, (3) develop, and (4) disseminate. We break each step into smaller concepts with strategies for implementation and provide illustrations from our collaboration to further give the readers actionable guidance.

Intrinsic Subtype and Overall Survival of Patients with Advanced HR+/HER2- Breast Cancer Treated with Ribociclib and ET: Correlative Analysis of MONALEESA-2, -3, -7.

PURPOSE: The MONALEESA-2, -3, -7 trials demonstrated statistically significant and clinically meaningful progression-free survival and overall survival (OS) benefits with ribociclib plus endocrine therapy (ET) versus ET alone in hormone receptor-positive, HER2-negative (HR+/HER2-) advanced breast cancer (ABC). Understanding the association of intrinsic subtypes with survival outcomes could potentially guide treatment decisions. Here, we evaluated the association of intrinsic subtypes with OS in MONALEESA-2, -3, -7. EXPERIMENTAL DESIGN: Tumor samples from MONALEESA-2, -3, -7 underwent PAM50-based subtyping. The relationship between subtypes and OS was assessed using univariable and multivariable Cox proportional hazards models. Multivariable models were adjusted for clinical prognostic factors. RESULTS: Overall, 990 tumors (among 2,066 patients) from ribociclib (n = 580) and placebo (n = 410) arms were profiled. Subtype distribution was luminal A, 54.5%; luminal B, 28.0%; HER2-enriched (HER2E) 14.6%; and basal-like, 2.8%; and was consistent across treatment arms. The luminal A subtype had the best OS outcomes in both arms, while basal-like had the worst. Patients with HER2E (HR, 0.60; P = 0.018), luminal B (HR, 0.69; P = 0.023), and luminal A (HR, 0.75; P = 0.021) subtypes derived OS benefit with ribociclib. Patients with basal-like subtype did not derive benefit from ribociclib (HR, 1.92; P = 0.137); however, patient numbers were small (n = 28). CONCLUSIONS: The prognostic value of intrinsic subtypes for OS was confirmed in this pooled analysis of the MONALEESA trials (largest dataset in HR+/HER2- ABC). While basal-like subtype did not benefit, a consistent OS benefit was observed with ribociclib added to ET across luminal and HER2E subtypes.

Clonal hematopoiesis detection in patients with cancer using cell-free DNA sequencing.

In the context of cancer, clonal hematopoiesis of indeterminate potential (CHIP) is associated with the development of therapy-related myeloid neoplasms and shorter overall survival. Cell-free DNA (cfDNA) sequencing is becoming widely adopted for genomic screening of patients with cancer but has not been used extensively to determine CHIP status because of a requirement for matched blood and tumor sequencing. We present an accurate classification approach to determine the CH status from cfDNA sequencing alone, applying our model to 4324 oncology clinical cfDNA samples. Using this method, we determined that 30.3% of patients in this cohort have evidence of CH, and the incidence of CH varies by tumor type. Matched RNA sequencing data show evidence of increased inflammation, especially neutrophil activation, within the tumors and tumor microenvironments of patients with CH. In addition, patients with CH had evidence of neutrophil activation systemically, pointing to a potential mechanism of action for the worse outcomes associated with CH status. Neutrophil activation may be one of many mechanisms, however, because patients with estrogen receptor-positive breast cancer harboring TET2 frameshift mutations had worse outcomes but similar neutrophil frequencies to patients without CH. Together, these data show the feasibility of detecting CH through cfDNA sequencing alone and an application of this method, demonstrating increased inflammation in patients with CH both systemically and in the tumor microenvironment.

Updated Overall Survival of Ribociclib plus Endocrine Therapy versus Endocrine Therapy Alone in Pre- and Perimenopausal Patients with HR+/HER2- Advanced Breast Cancer in MONALEESA-7: A Phase III Randomized Clinical Trial.

PURPOSE: Ribociclib plus endocrine therapy (ET) demonstrated a statistically significant progression-free survival and overall survival (OS) benefit in the phase III MONALEESA-7 trial of pre-/perimenopausal patients with hormone receptor (HR)-positive (HR+), HER2-negative (HER2-) advanced breast cancer (ABC). The median OS was not reached in the ribociclib arm in the protocol-specified final analysis; we hence performed an exploratory OS and additional outcomes analysis with an extended follow-up (median, 53.5 months). PATIENTS AND METHODS: Patients were randomized to receive ET [goserelin plus nonsteroidal aromatase inhibitor (NSAI) or tamoxifen] with ribociclib or placebo. OS was evaluated with a stratified Cox proportional hazard model and summarized with Kaplan-Meier methods. RESULTS: The intent-to-treat population included 672 patients. Median OS was 58.7 months with ribociclib versus 48.0 months with placebo [hazard ratio = 0.76; 95% confidence interval (CI), 0.61-0.96]. Kaplan-Meier estimated OS at 48 months was 60% and 50% with ribociclib and placebo, respectively. Subgroup analyses were generally consistent with the OS benefit, including patients who received NSAI and patients aged less than 40 years. Subsequent antineoplastic therapies following discontinuation were balanced between the ribociclib (77%) and placebo (78%) groups. Use of cyclin-dependent kinase 4/6 inhibitors after discontinuation was higher with placebo (26%) versus ribociclib (13%). Time to first chemotherapy was significantly delayed with ribociclib versus placebo. No drug-drug interactions were observed between ribociclib and either NSAI. CONCLUSIONS: Ribociclib plus ET continued to show significantly longer OS than ET alone in pre-/perimenopausal patients, including patients aged less than 40 years, with HR+/HER2- ABC with 53.5 months of median follow-up (ClinicalTrials.gov, NCT02278120).

Genomic Profiling of Premenopausal HR+ and HER2- Metastatic Breast Cancer by Circulating Tumor DNA and Association of Genetic Alterations With Therapeutic Response to Endocrine Therapy and Ribociclib.

PURPOSE: This analysis evaluated the genomic landscape of premenopausal patients with hormone receptor-positive and human epidermal growth factor receptor 2-negative advanced breast cancer and the association of genetic alterations with response to ribociclib in the phase III MONALEESA-7 trial. METHODS: Premenopausal patients were randomly assigned 1:1 to receive endocrine therapy plus ribociclib or placebo. Plasma collected at baseline was sequenced using targeted next-generation sequencing for approximately 600 relevant cancer genes. The association of circulating tumor DNA alterations with progression-free survival (PFS) was evaluated to identify biomarkers of response and resistance to ribociclib. RESULTS: Baseline circulating tumor DNA was sequenced in 565 patients; 489 had evidence of ≥ 1 alteration. The most frequent alterations included PIK3CA (28%), TP53 (19%), CCND1 (10%), MYC (8%), GATA3 (8%), receptor tyrosine kinases (17%), and the Chr8p11.23 locus (12%). A treatment benefit of ribociclib was seen with wild-type (hazard ratio [HR] 0.45 [95% CI, 0.33 to 0.62]) and altered (HR 0.57 [95% CI, 0.36 to 0.9]) PIK3CA. Overall, patients with altered CCND1 had shorter PFS regardless of treatment, suggesting CCND1 as a potential prognostic biomarker. Benefit with ribociclib was seen in patients with altered (HR 0.21 [95% CI, 0.08 to 0.54]) or wild-type (HR 0.52 [95% CI, 0.39 to 0.68]) CCND1, but greater benefit was observed with altered, suggesting predictive potential of CCND1. Alterations in TP53, MYC, Chr8p11.23 locus, and receptor tyrosine kinases were associated with worse PFS, but ribociclib benefit was independent of alteration status. CONCLUSION: In this study-to our knowledge, the first large study of premenopausal patients with hormone receptor-positive and human epidermal growth factor receptor 2-negative advanced breast cancer-multiple genomic alterations were associated with poor outcome. A PFS benefit of ribociclib was observed regardless of gene alteration status, although in this exploratory analysis, a magnitude of benefits varied by alteration.

Correlative Biomarker Analysis of Intrinsic Subtypes and Efficacy Across the MONALEESA Phase III Studies.

PURPOSE: The prognostic and predictive value of intrinsic subtypes in hormone receptor-positive, human epidermal growth factor receptor 2-negative advanced breast cancer treated with endocrine therapy and ribociclib (RIB) is currently unknown. We evaluated the association of intrinsic subtypes with progression-free survival (PFS) in the MONALEESA trials. METHODS: A retrospective and exploratory PAM50-based analysis of tumor samples from the phase III MONALEESA-2, MONALEESA-3, and MONALEESA-7 trials was undertaken. The prognostic relationship of PAM50-based subtypes with PFS and risk of disease progression by subtype and treatment were evaluated using a multivariable Cox proportional hazards model, adjusting for age, prior chemotherapy, performance status, visceral disease, bone-only metastases, histological grade, number of metastatic sites, prior endocrine therapy, and de novo metastatic disease. RESULTS: Overall, 1,160 tumors from the RIB (n = 672) and placebo (n = 488) cohorts were robustly profiled. Subtype distribution was luminal A (LumA), 46.7%; luminal B (LumB), 24.0%; normal-like, 14.0%; HER2-enriched (HER2E), 12.7%; and basal-like, 2.6% and was generally consistent across treatment arms and trials. The associations between subtypes and PFS were statistically significant in both arms (P < .001). The risks of disease progression for LumB, HER2E, and basal-like subtypes were 1.44, 2.31, and 3.96 times higher compared with those for LumA, respectively. All subtypes except basal-like demonstrated significant PFS benefit with RIB. HER2E (hazard ratio [HR], 0.39; P < .0001), LumB (HR, 0.52; P < .0001), LumA (HR, 0.63; P = .0007), and normal-like (HR, 0.47; P = .0005) subtypes derived benefit from RIB. Patients with basal-like subtype (n = 30) did not derive benefit from RIB (HR, 1.15; P = .77). CONCLUSION: In this retrospective exploratory analysis of hormone receptor-positive and human epidermal growth factor receptor 2-negative advanced breast cancer, each intrinsic subtype exhibited a consistent PFS benefit with RIB, except for basal-like.

Aberrant FGFR signaling mediates resistance to CDK4/6 inhibitors in ER+ breast cancer.

Using an ORF kinome screen in MCF-7 cells treated with the CDK4/6 inhibitor ribociclib plus fulvestrant, we identified FGFR1 as a mechanism of drug resistance. FGFR1-amplified/ER+ breast cancer cells and MCF-7 cells transduced with FGFR1 were resistant to fulvestrant ± ribociclib or palbociclib. This resistance was abrogated by treatment with the FGFR tyrosine kinase inhibitor (TKI) lucitanib. Addition of the FGFR TKI erdafitinib to palbociclib/fulvestrant induced complete responses of FGFR1-amplified/ER+ patient-derived-xenografts. Next generation sequencing of circulating tumor DNA (ctDNA) in 34 patients after progression on CDK4/6 inhibitors identified FGFR1/2 amplification or activating mutations in 14/34 (41%) post-progression specimens. Finally, ctDNA from patients enrolled in MONALEESA-2, the registration trial of ribociclib, showed that patients with FGFR1 amplification exhibited a shorter progression-free survival compared to patients with wild type FGFR1. Thus, we propose breast cancers with FGFR pathway alterations should be considered for trials using combinations of ER, CDK4/6 and FGFR antagonists.

Molecular Alterations and Buparlisib Efficacy in Patients with Squamous Cell Carcinoma of the Head and Neck: Biomarker Analysis from BERIL-1.

Purpose: The preplanned exploratory analysis of the BERIL-1 trial presented here aimed to identify biomarkers of response to the combination of buparlisib and paclitaxel.Patients and Methods: BERIL-1 was a multicenter, randomized, double-blind, placebo-controlled phase II study. Patients with recurrent or metastatic squamous cell carcinoma of the head and neck (SCCHN) progressing on/after one previous platinum-based chemotherapy regimen in the recurrent or metastatic setting were treated with either buparlisib plus paclitaxel or placebo plus paclitaxel. Archival tumor tissue and ctDNA samples were analyzed for molecular alterations and immune infiltration using next-generation sequencing or immunohistochemistry.Results: Biomarker analyses were performed in randomized patients (n = 158) with available biomarker data. The most frequently (>5%) mutated genes were TP53, FAT1, TET2, KMT2D, PIK3CA, NOTCH1, NFE2L2, NOTCH2, CCND1, and CDKN2A Patients with SCCHN tumors (from various primary sites) having HPV-negative status (HR = 0.51), TP53 alterations (HR = 0.55) or low mutational load (HR = 0.57) derived overall survival (OS) benefit with the combination of buparlisib and paclitaxel. OS benefit with this combination was also increased in patients with presence of intratumoral TILs ≥10% (HR = 0.51), stromal TILs ≥15% (HR = 0.53), intratumoral CD8-positive cells ≥5% (HR = 0.45), stromal CD8-positive cells ≥10% (HR = 0.47), or CD8-positive cells in invasive margins >25% (HR = 0.37). A trend for improved progression-free survival with the combination of buparlisib and paclitaxel was also observed in these patients.Conclusions: The BERIL-1 biomarker analyses showed that patients with TP53 alterations, HPV-negative status, low mutational load, or high infiltration of TILs or CD8-positive cells derived survival benefit with the combination of buparlisib and paclitaxel. Clin Cancer Res; 24(11); 2505-16. ©2018 AACR.

Association of MTOR Mutations With Developmental Brain Disorders, Including Megalencephaly, Focal Cortical Dysplasia, and Pigmentary Mosaicism.

IMPORTANCE: Focal cortical dysplasia (FCD), hemimegalencephaly, and megalencephaly constitute a spectrum of malformations of cortical development with shared neuropathologic features. These disorders are associated with significant childhood morbidity and mortality. OBJECTIVE: To identify the underlying molecular cause of FCD, hemimegalencephaly, and diffuse megalencephaly. DESIGN, SETTING, AND PARTICIPANTS: Patients with FCD, hemimegalencephaly, or megalencephaly (mean age, 11.7 years; range, 2-32 years) were recruited from Pediatric Hospital A. Meyer, the University of Hong Kong, and Seattle Children's Research Institute from June 2012 to June 2014. Whole-exome sequencing (WES) was performed on 8 children with FCD or hemimegalencephaly using standard-depth (50-60X) sequencing in peripheral samples (blood, saliva, or skin) from the affected child and their parents and deep (150-180X) sequencing in affected brain tissue. Targeted sequencing and WES were used to screen 93 children with molecularly unexplained diffuse or focal brain overgrowth. Histopathologic and functional assays of phosphatidylinositol 3-kinase-AKT (serine/threonine kinase)-mammalian target of rapamycin (mTOR) pathway activity in resected brain tissue and cultured neurons were performed to validate mutations. MAIN OUTCOMES AND MEASURES: Whole-exome sequencing and targeted sequencing identified variants associated with this spectrum of developmental brain disorders. RESULTS: Low-level mosaic mutations of MTOR were identified in brain tissue in 4 children with FCD type 2a with alternative allele fractions ranging from 0.012 to 0.086. Intermediate-level mosaic mutation of MTOR (p.Thr1977Ile) was also identified in 3 unrelated children with diffuse megalencephaly and pigmentary mosaicism in skin. Finally, a constitutional de novo mutation of MTOR (p.Glu1799Lys) was identified in 3 unrelated children with diffuse megalencephaly and intellectual disability. Molecular and functional analysis in 2 children with FCD2a from whom multiple affected brain tissue samples were available revealed a mutation gradient with an epicenter in the most epileptogenic area. When expressed in cultured neurons, all MTOR mutations identified here drive constitutive activation of mTOR complex 1 and enlarged neuronal size. CONCLUSIONS AND RELEVANCE: In this study, mutations of MTOR were associated with a spectrum of brain overgrowth phenotypes extending from FCD type 2a to diffuse megalencephaly, distinguished by different mutations and levels of mosaicism. These mutations may be sufficient to cause cellular hypertrophy in cultured neurons and may provide a demonstration of the pattern of mosaicism in brain and substantiate the link between mosaic mutations of MTOR and pigmentary mosaicism in skin.

Genetic association analysis of 300 genes identifies a risk haplotype in SLC18A2 for post-traumatic stress disorder in two independent samples.

The genetic architecture of post-traumatic stress disorder (PTSD) remains poorly understood with the vast majority of genetic association studies reporting on single candidate genes. We conducted a large genetic study in trauma-exposed European-American women (N=2538; 845 PTSD cases, 1693 controls) by testing 3742 SNPs across more than 300 genes and conducting polygenic analyses using results from the Psychiatric Genome-Wide Association Studies Consortium (PGC). We tested the association between each SNP and two measures of PTSD, a severity score and diagnosis. We found a significant association between PTSD (diagnosis) and SNPs (top SNP: rs363276, odds ratio (OR)=1.4, p=2.1E-05) in SLC18A2 (vesicular monoamine transporter 2). A haplotype analysis of 9 SNPs in SLC18A2, including rs363276, identified a risk haplotype (CGGCGGAAG, p=0.0046), and the same risk haplotype was associated with PTSD in an independent cohort of trauma-exposed African-Americans (p=0.049; N=748, men and women). SLC18A2 is involved in transporting monoamines to synaptic vesicles and has been implicated in a number of neuropsychiatric disorders including major depression. Eight genes previously associated with PTSD had SNPs with nominally significant associations (p<0.05). The polygenic analyses suggested that there are SNPs in common between PTSD severity and bipolar disorder. Our data are consistent with a genetic architecture for PTSD that is highly polygenic, influenced by numerous SNPs with weak effects, and may overlap with mood disorders. Genome-wide studies with very large samples sizes are needed to detect these types of effects.

A polygenic burden of rare disruptive mutations in schizophrenia.

Schizophrenia is a common disease with a complex aetiology, probably involving multiple and heterogeneous genetic factors. Here, by analysing the exome sequences of 2,536 schizophrenia cases and 2,543 controls, we demonstrate a polygenic burden primarily arising from rare (less than 1 in 10,000), disruptive mutations distributed across many genes. Particularly enriched gene sets include the voltage-gated calcium ion channel and the signalling complex formed by the activity-regulated cytoskeleton-associated scaffold protein (ARC) of the postsynaptic density, sets previously implicated by genome-wide association and copy-number variation studies. Similar to reports in autism, targets of the fragile X mental retardation protein (FMRP, product of FMR1) are enriched for case mutations. No individual gene-based test achieves significance after correction for multiple testing and we do not detect any alleles of moderately low frequency (approximately 0.5 to 1 per cent) and moderately large effect. Taken together, these data suggest that population-based exome sequencing can discover risk alleles and complements established gene-mapping paradigms in neuropsychiatric disease.

Interaction between genetic variants and exposure to Hurricane Katrina on post-traumatic stress and post-traumatic growth: a prospective analysis of low income adults.

BACKGROUND: There is considerable variation in psychological reactions to natural disasters, with responses ranging from relatively mild and transitory symptoms to severe and persistent posttraumatic stress (PTS). Some survivors also report post-traumatic growth (PTG), or positive psychological changes due to the experience and processing of the disaster and its aftermath. Gene-environment interaction (GxE) studies could offer new insight into the factors underlying variability in post-disaster psychological responses. However, few studies have explored GxE in a disaster context. METHODS: We examined whether ten common variants in seven genes (BDNF, CACNA1C, CRHR1, FKBP5, OXTR, RGS2, SLC6A4) modified associations between Hurricane Katrina exposure and PTS and PTG. Data were from a prospective study of 205 low-income non-Hispanic Black parents residing in New Orleans prior to and following Hurricane Katrina. RESULTS: We found a significant association (after correction) between RGS2 (rs4606; p=0.0044) and PTG, which was mainly driven by a cross-over GxE (p=0.006), rather than a main genetic effect (p=0.071). The G (minor allele) was associated with lower PTG scores for low levels of Hurricane exposure and higher PTG scores for moderate and high levels of exposure. We also found a nominally significant association between variation in FKBP5 (rs1306780, p=0.0113) and PTG, though this result did not survive correction for multiple testing. LIMITATIONS: Although the inclusion of low-income non-Hispanic Black parents allowed us to examine GxE among a highly vulnerable group, our findings may not generalize to other populations or groups experiencing other natural disasters. Moreover, not all participants invited to participate in the genetic study provided saliva. CONCLUSIONS: To our knowledge, this is the first study to identify GxE in the context of post-traumatic growth. Future studies are needed to clarify the role of GxE in PTS and PTG and post-disaster psychological responses, especially among vulnerable populations.

Genome-wide association study implicates a novel RNA gene, the lincRNA AC068718.1, as a risk factor for post-traumatic stress disorder in women.

Posttraumatic stress disorder (PTSD) is a common and debilitating mental disorder with a particularly high burden for women. Emerging evidence suggests PTSD may be more heritable among women and evidence from animal models and human correlational studies suggest connections between sex-linked biology and PTSD vulnerability, which may extend to the disorder's genetic architecture. We conducted a genome-wide association study (GWAS) of PTSD in a primarily African American sample of women from the Detroit Neighborhood Health Study (DNHS) and tested for replication in an independent cohort of primarily European American women from the Nurses Health Study II (NHSII). We genotyped 413 DNHS women - 94 PTSD cases and 319 controls exposed to at least one traumatic event - on the Illumina HumanOmniExpress BeadChip for >700,000 markers and tested 578 PTSD cases and 1963 controls from NHSII for replication. We performed a network-based analysis integrating data from GWAS-derived independent regions of association and the Reactome database of functional interactions. We found genome-wide significant association for one marker mapping to a novel RNA gene, lincRNA AC068718.1, for which we found suggestive evidence of replication in NHSII. Our network-based analysis indicates that our top GWAS results were enriched for pathways related to telomere maintenance and immune function. Our findings implicate a novel RNA gene, lincRNA AC068718.1, as risk factor for PTSD in women and add to emerging evidence that non-coding RNA genes may play a crucial role in shaping the landscape of gene regulation with putative pathological effects that lead to phenotypic differences.

The ankyrin-3 gene is associated with posttraumatic stress disorder and externalizing comorbidity.

BACKGROUND: The ankyrin 3 gene (ANK3) produces the ankyrin G protein that plays an integral role in regulating neuronal activity. Previous studies have linked ANK3 to bipolar disorder and schizophrenia. A recent mouse study suggests that ANK3 may regulate behavioral disinhibition and stress reactivity. This led us to hypothesize that ANK3 might also be associated with stress-related psychopathology such as posttraumatic stress disorder (PTSD), as well as disorders of the externalizing spectrum such as antisocial personality disorder and substance-related disorders that are etiologically linked to impulsivity and temperamental disinhibition. METHODS: We examined the possibility of association between ANK3 SNPs and both PTSD and externalizing (defined by a factor score representing a composite of adult antisociality and substance abuse) in a cohort of white non-Hispanic combat veterans and their intimate partners (n=554). Initially, we focused on rs9804190-a SNP previously reported to be associated with bipolar disorder, schizophrenia, and ankyrin G expression in brain. Then we examined 358 additional ANK3 SNPs utilizing a multiple-testing correction. RESULTS: rs9804190 was associated with both externalizing and PTSD (p=0.028 and p=0.042 respectively). Analysis of other ANK3 SNPs identified several that were more strongly associated with either trait. The most significant association with externalizing was observed at rs1049862 (p=0.00040, pcorrected=0.60). The most significant association with PTSD (p=0.00060, pcorrected=0.045) was found with three SNPs in complete linkage disequilibrium (LD)-rs28932171, rs11599164, and rs17208576. CONCLUSIONS: These findings support a role of ANK3 in risk of stress-related and externalizing disorders, beyond its previous associations with bipolar disorder and schizophrenia.

Pleiotropy in complex traits: challenges and strategies.

Genome-wide association studies have identified many variants that each affects multiple traits, particularly across autoimmune diseases, cancers and neuropsychiatric disorders, suggesting that pleiotropic effects on human complex traits may be widespread. However, systematic detection of such effects is challenging and requires new methodologies and frameworks for interpreting cross-phenotype results. In this Review, we discuss the evidence for pleiotropy in contemporary genetic mapping studies, new and established analytical approaches to identifying pleiotropic effects, sources of spurious cross-phenotype effects and study design considerations. We also outline the molecular and clinical implications of such findings and discuss future directions of research.

Fetal hemoglobin in sickle cell anemia: genetic studies of the Arab-Indian haplotype.

Sickle cell anemia is common in the Middle East and India where the HbS gene is sometimes associated with the Arab-Indian (AI) ß-globin gene (HBB) cluster haplotype. In this haplotype of sickle cell anemia, fetal hemoglobin (HbF) levels are 3-4 fold higher than those found in patients with HbS haplotypes of African origin. Little is known about the genetic elements that modulate HbF in AI haplotype patients. We therefore studied Saudi HbS homozygotes with the AI haplotype (mean HbF 19.2±7.0%, range 3.6 to 39.6%) and employed targeted genotyping of polymorphic sites to explore cis- and trans- acting elements associated with high HbF expression. We also described sequences which appear to be unique to the AI haplotype for which future functional studies are needed to further define their role in HbF modulation. All cases, regardless of HbF concentration, were homozygous for AI haplotype-specific elements cis to HBB. SNPs in BCL11A and HBS1L-MYB that were associated with HbF in other populations explained only 8.8% of the variation in HbF. KLF1 polymorphisms associated previously with high HbF were not present in the 44 patients tested. More than 90% of the HbF variance in sickle cell patients with the AI haplotype remains unexplained by the genetic loci that we studied. The dispersion of HbF levels among AI haplotype patients suggests that other genetic elements modulate the effects of the known cis- and trans-acting regulators. These regulatory elements, which remain to be discovered, might be specific in the Saudi and some other populations where HbF levels are especially high.

A genome-wide association study of total bilirubin and cholelithiasis risk in sickle cell anemia.

Serum bilirubin levels have been associated with polymorphisms in the UGT1A1 promoter in normal populations and in patients with hemolytic anemias, including sickle cell anemia. When hemolysis occurs circulating heme increases, leading to elevated bilirubin levels and an increased incidence of cholelithiasis. We performed the first genome-wide association study (GWAS) of bilirubin levels and cholelithiasis risk in a discovery cohort of 1,117 sickle cell anemia patients. We found 15 single nucleotide polymorphisms (SNPs) associated with total bilirubin levels at the genome-wide significance level (p value <5 × 10(-8)). SNPs in UGT1A1, UGT1A3, UGT1A6, UGT1A8 and UGT1A10, different isoforms within the UGT1A locus, were identified (most significant rs887829, p = 9.08 × 10(-25)). All of these associations were validated in 4 independent sets of sickle cell anemia patients. We tested the association of the 15 SNPs with cholelithiasis in the discovery cohort and found a significant association (most significant p value 1.15 × 10(-4)). These results confirm that the UGT1A region is the major regulator of bilirubin metabolism in African Americans with sickle cell anemia, similar to what is observed in other ethnicities.

Naïve Bayesian Classifier and Genetic Risk Score for Genetic Risk Prediction of a Categorical Trait: Not so Different after all!

One of the most popular modeling approaches to genetic risk prediction is to use a summary of risk alleles in the form of an unweighted or a weighted genetic risk score, with weights that relate to the odds for the phenotype in carriers of the individual alleles. Recent contributions have proposed the use of Bayesian classification rules using Naïve Bayes classifiers. We examine the relation between the two approaches for genetic risk prediction and show that the methods are mathematically related. In addition, we study the properties of the two approaches and describe how they can be generalized to include various models of inheritance.

Human longevity and common variations in the LMNA gene: a meta-analysis.

A mutation in the LMNA gene is responsible for the most dramatic form of premature aging, Hutchinson-Gilford progeria syndrome (HGPS). Several recent studies have suggested that protein products of this gene might have a role in normal physiological cellular senescence. To explore further LMNA's possible role in normal aging, we genotyped 16 SNPs over a span of 75.4 kb of the LMNA gene on a sample of long-lived individuals (LLI) (US Caucasians with age ≥ 95 years, N=873) and genetically matched younger controls (N=443). We tested all common nonredundant haplotypes (frequency ≥ 0.05) based on subgroups of these 16 SNPs for association with longevity. The most significant haplotype, based on four SNPs, remained significant after adjustment for multiple testing (OR=1.56, P=2.5 × 10(-5) , multiple-testing-adjusted P=0.0045). To attempt to replicate these results, we genotyped 3619 subjects from four independent samples of LLI and control subjects from (i) the New England Centenarian Study (NECS) (N=738), (ii) the Southern Italian Centenarian Study (SICS) (N=905), (iii) France (N=1103), and (iv) the Einstein Ashkenazi Longevity Study (N= 702). We replicated the association with the most significant haplotype from our initial analysis in the NECS sample (OR=1.60, P=0.0023), but not in the other three samples (P > 0.15). In a meta-analysis combining all five samples, the best haplotype remained significantly associated with longevity after adjustment for multiple testing in the initial and follow-up samples (OR=1.18, P=7.5 × 10(-4) , multiple-testing-adjusted P=0.037). These results suggest that LMNA variants may play a role in human lifespan.