Mapping Systemic Inflammation and Antibody Responses in Multisystem Inflammatory Syndrome in Children (MIS-C).

Initially, children were thought to be spared from disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, a month into the epidemic, a novel multisystem inflammatory syndrome in children (MIS-C) emerged. Herein, we report on the immune profiles of nine MIS-C cases. All MIS-C patients had evidence of prior SARS-CoV-2 exposure, mounting an antibody response with intact neutralization capability. Cytokine profiling identified elevated signatures of inflammation (IL-18 and IL-6), lymphocytic and myeloid chemotaxis and activation (CCL3, CCL4, and CDCP1), and mucosal immune dysregulation (IL-17A, CCL20, and CCL28). Immunophenotyping of peripheral blood revealed reductions of non-classical monocytes, and subsets of NK and T lymphocytes, suggesting extravasation to affected tissues. Finally, profiling the autoantigen reactivity of MIS-C plasma revealed both known disease-associated autoantibodies (anti-La) and novel candidates that recognize endothelial, gastrointestinal, and immune-cell antigens. All patients were treated with anti-IL-6R antibody and/or IVIG, which led to rapid disease resolution.

Inherited cancer predisposing mutations in patients with therapy-related myeloid neoplasms.

Some patients with therapy-related myeloid neoplasms (t-MN) may have unsuspected inherited cancer predisposition syndrome (CPS). We propose a set of clinical criteria to identify t-MN patients with high risk of CPS (HR-CPS). Among 225 t-MN patients with an antecedent non-myeloid malignancy, our clinical criteria identified 52 (23%) HR-CPS patients. Germline whole-exome sequencing identified pathogenic or likely pathogenic variants in 10 of 27 HR-CPS patients compared to 0 of 9 low-risk CPS patients (37% vs. 0%, p = 0.04). These simple clinical criteria identify t-MN patients most likely to benefit from genetic testing for inherited CPS.

Characterization of COVID-19 disease in pediatric oncology patients: The New York-New Jersey regional experience.

PURPOSE: Pediatric oncology patients undergoing active chemotherapy are suspected to be at a high risk for severe disease secondary to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection; however, data to support this are lacking. We aim to describe the characteristics of coronavirus disease 2019 (COVID-19) in this population and also its impact on pediatric cancer care in the New York region during the peak of the pandemic. PATIENTS AND METHODS: This multicenter, retrospective study included 13 institutions. Clinical and laboratory information on 98 patients ≤21 years of age receiving active anticancer therapy, who tested positive for SARS-CoV-2 by nasopharyngeal swab polymerase chain reaction (PCR), was collected. RESULTS: Of the 578 pediatric oncology patients tested for COVID-19, 98 were positive, of whom 73 were symptomatic. Most experienced mild disease, 28 required inpatient management, 25 needed oxygen support, and seven required mechanical ventilation. There is a slightly higher risk of severe disease in males and obese patients, though not statistically significant. Persistent lymphopenia was noted in severe cases. Delays in cancer therapy occurred in 67% of SARS-CoV-2-positive patients. Of four deaths, none were solely attributable to COVID-19. The impact of the pandemic on pediatric oncology care was significant, with 54% of institutions reporting delays in chemotherapy, 46% delays in surgery, and 30% delays in transplant. CONCLUSION: In this large multi-institutional cohort, we observed that mortality and morbidity from COVID-19 amongst pediatric oncology patients were low overall, but higher than reported in general pediatrics. Certain subgroups might be at higher risk of severe disease. Delays in cancer care due to SARS-CoV-2 remain a concern.

Genetic overlap between autoimmune diseases and non-Hodgkin lymphoma subtypes.

Epidemiologic studies show an increased risk of non-Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen-stimulation in one disease leads to another. Here we assess shared genetic risk in genome-wide-association-studies (GWAS). Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B-cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta-analysis. Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta-analysis (combing AD with NHL) implicated genes of apoptosis and telomere length. This GWAS-based analysis four NHL subtypes and three ADs revealed few weakly-associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Exome chip analyses identify genes affecting mortality after HLA-matched unrelated-donor blood and marrow transplantation.

Although survival outcomes have significantly improved, up to 40% of patients die within 1 year of HLA-matched unrelated-donor blood and marrow transplantation (BMT). To identify non-HLA genetic contributors to mortality after BMT, we performed the first exome-wide association study in the DISCOVeRY-BMT cohorts using the Illumina HumanExome BeadChip. This study includes 2473 patients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome and 2221 10/10 HLA-matched donors treated from 2000 to 2011. Single-variant and gene-level analyses were performed on overall survival (OS), transplantation-related mortality (TRM), and disease-related mortality (DRM). Genotype mismatches between recipients and donors in a rare nonsynonymous variant of testis-expressed gene TEX38 significantly increased risk of TRM, which was more dramatic when either the recipient or donor was female. Using the SKAT-O test to evaluate gene-level effects, variant genotypes of OR51D1 in recipients were significantly associated with OS and TRM. In donors, 4 (ALPP, EMID1, SLC44A5, LRP1), 1 (HHAT), and 2 genes (LYZL4, NT5E) were significantly associated with OS, TRM, and DRM, respectively. Inspection of NT5E crystal structures showed 4 of the associated variants affected the enzyme structure and likely decreased the catalytic efficiency of the enzyme. Further confirmation of these findings and additional functional studies may provide individualized risk prediction and prognosis, as well as alternative donor selection strategies.

Comparing ethnicity-based and expanded carrier screening methods at a single fertility center reveals significant differences in carrier rates and carrier couple rates.

PURPOSE: To evaluate the efficiency of expanded carrier screening (ECS) compared with ethnicity-based screening in identifying carriers. METHODS: A total of 4232 infertility patients underwent ECS from a single genetic testing laboratory at our center between June 2013 and July 2015. Self-reported ethnicity was recorded. Carrier rates based on ECS were calculated. In addition, carrier status was determined for two other screening panels: ethnicity-based guidelines or the ECS panel recommended by the American College of Obstetricians and Gynecologists (ACOG) using ECS results. Carrier rate and carrier couple rates were compared in the overall study population and in each self-reported ethnicity. RESULTS: The ECS panel used to screen the patient population identified 1243 carriers (29.4%). For the same population, ethnicity-based screening and the ACOG panel would have identified 359 (8.5%) and 659 carriers (15.6%), respectively, representing statistically significant differences. Differences in identifying carriers across self-reported ethnicities varied. In 15 couples (1.2%), both partners carried pathogenic variants for the same genes, 47% of whom would have been missed had screening been ethnicity-based. CONCLUSION: We propose that all reproductive-aged women should be offered ECS. Carrier couple rates would likely increase further with expansion of the panel, playing a pivotal role in preventing genetic disease in fertility clinics.

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant.

Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.

Choosing an expanded carrier screening panel: comparing two panels at a single fertility centre.

RESEARCH QUESTION: What are the factors contributing to similarities and differences in carrier rates between two expanded carrier screening (ECS) panels? DESIGN: Retrospective cross-sectional study. A total of 7700 infertility patients who underwent ECS from one of two genetic testing laboratories (Panel A or Panel B) using a genotyping microarray were included in the study. Individuals presenting to the Centre between June 2013 and July 2015 underwent screening via Panel A. Those presenting between August 2015 and April 2017 underwent screening via Panel B. Self-reported ethnicity was recorded. Panel content, carrier rates for the overall study population and for comparable self-reported ethnicities, carrier couple rates, and the top 10 identified disorders were compared. RESULTS: Of 4232 individuals screened by Panel A, 1243 were identified as carriers (29.4%). Panel B identified 1503 carriers among the 3468 (43.3%) participants (P < 0.0001). Carrier couple rate also varied between panels (1.2% versus 3.1%; P = 0.0017). A total of 311 disorders covering 2746 mutations were observed across the two ECS panels, with 372 (13.5%) shared mutations. Carrier rates did not differ for the shared mutations overall and across ethnicities. Significant differences were observed when comparing unique content in the overall population (P < 2 .2 × 10-16) and across ethnicities (P < 2.2 × 10-16 to 0.0010). CONCLUSIONS: Carrier rates in the overall population and across ethnicities vary widely based on panel content, and highlight the need to expand panel content as well as incorporate preconception carrier screening coupled with genetic counselling into routine assisted reproduction practice.

The APC tumor suppressor is required for epithelial cell polarization and three-dimensional morphogenesis.

The Adenomatous Polyposis Coli (APC) tumor suppressor has been previously implicated in the control of apical-basal polarity; yet, the consequence of APC loss-of-function in epithelial polarization and morphogenesis has not been characterized. To test the hypothesis that APC is required for the establishment of normal epithelial polarity and morphogenesis programs, we generated APC-knockdown epithelial cell lines. APC depletion resulted in loss of polarity and multi-layering on permeable supports, and enlarged, filled spheroids with disrupted polarity in 3D culture. Importantly, these effects of APC knockdown were independent of Wnt/ß-catenin signaling, but were rescued with either full-length or a carboxy (c)-terminal segment of APC. Moreover, we identified a gene expression signature associated with APC knockdown that points to several candidates known to regulate cell-cell and cell-matrix communication. Analysis of epithelial tissues from mice and humans carrying heterozygous APC mutations further supports the importance of APC as a regulator of epithelial behavior and tissue architecture. These data also suggest that the initiation of epithelial-derived tumors as a result of APC mutation or gene silencing may be driven by loss of polarity and dysmorphogenesis.

The genetics of breast cancer risk in the post-genome era: thoughts on study design to move past BRCA and towards clinical relevance.

More than 12 % of women will be diagnosed with breast cancer in their lifetime. Although there have been tremendous advances in elucidating genetic risk factors underlying both familial and sporadic breast cancer, much of the genetic contribution to breast cancer etiology remains unknown. The discovery of BRCA1 and BRCA2 over 20 years ago remains the seminal event in the field and has paved the way for the discovery of other high-penetrance susceptibility genes by linkage analysis. The advent of genome-wide association studies made possible the next wave of discoveries, in which over 80 low-penetrance and moderate-penetrance variants were identified. Although these studies were highly successful at discovering variants associated with both familial and sporadic breast cancer, the variants identified to date explain only 50 % of the heritability of breast cancer. In this review, we look back at the investigative strategies that have led to our current understanding of breast cancer genetics, consider the challenges of performing association studies in heterogeneous complex diseases such as breast cancer, and look ahead toward the types of study designs that may lead to the identification of the genetic variation accounting for the remaining missing heritability.

Systematic computation with functional gene-sets among leukemic and hematopoietic stem cells reveals a favorable prognostic signature for acute myeloid leukemia.

BACKGROUND: Genes that regulate stem cell function are suspected to exert adverse effects on prognosis in malignancy. However, diverse cancer stem cell signatures are difficult for physicians to interpret and apply clinically. To connect the transcriptome and stem cell biology, with potential clinical applications, we propose a novel computational "gene-to-function, snapshot-to-dynamics, and biology-to-clinic" framework to uncover core functional gene-sets signatures. This framework incorporates three function-centric gene-set analysis strategies: a meta-analysis of both microarray and RNA-seq data, novel dynamic network mechanism (DNM) identification, and a personalized prognostic indicator analysis. This work uses complex disease acute myeloid leukemia (AML) as a research platform. RESULTS: We introduced an adjustable "soft threshold" to a functional gene-set algorithm and found that two different analysis methods identified distinct gene-set signatures from the same samples. We identified a 30-gene cluster that characterizes leukemic stem cell (LSC)-depleted cells and a 25-gene cluster that characterizes LSC-enriched cells in parallel; both mark favorable-prognosis in AML. Genes within each signature significantly share common biological processes and/or molecular functions (empirical p = 6e-5 and 0.03 respectively). The 25-gene signature reflects the abnormal development of stem cells in AML, such as AURKA over-expression. We subsequently determined that the clinical relevance of both signatures is independent of known clinical risk classifications in 214 patients with cytogenetically normal AML. We successfully validated the prognosis of both signatures in two independent cohorts of 91 and 242 patients respectively (log-rank p < 0.0015 and 0.05; empirical p < 0.015 and 0.08). CONCLUSION: The proposed algorithms and computational framework will harness systems biology research because they efficiently translate gene-sets (rather than single genes) into biological discoveries about AML and other complex diseases.

Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci.

BACKGROUND: Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with diseases of the colon including inflammatory bowel diseases (IBD) and colorectal cancer (CRC). However, the functional role of many of these SNPs is largely unknown and tissue-specific resources are lacking. Expression quantitative trait loci (eQTL) mapping identifies target genes of disease-associated SNPs. This study provides a comprehensive eQTL map of distal colonic samples obtained from 40 healthy African Americans and demonstrates their relevance for GWAS of colonic diseases. RESULTS: 8.4 million imputed SNPs were tested for their associations with 16,252 expression probes representing 12,363 unique genes. 1,941 significant cis-eQTL, corresponding to 122 independent signals, were identified at a false discovery rate (FDR) of 0.01. Overall, among colon cis-eQTL, there was significant enrichment for GWAS variants for IBD (Crohn's disease [CD] and ulcerative colitis [UC]) and CRC as well as type 2 diabetes and body mass index. ERAP2, ADCY3, INPP5E, UBA7, SFMBT1, NXPE1 and REXO2 were identified as target genes for IBD-associated variants. The CRC-associated eQTL rs3802842 was associated with the expression of C11orf93 (COLCA2). Enrichment of colon eQTL near transcription start sites and for active histone marks was demonstrated, and eQTL with high population differentiation were identified. CONCLUSIONS: Through the comprehensive study of eQTL in the human colon, this study identified novel target genes for IBD- and CRC-associated genetic variants. Moreover, bioinformatic characterization of colon eQTL provides a tissue-specific tool to improve understanding of biological differences in diseases between different ethnic groups.

Establishment of Definitions and Review Process for Consistent Adjudication of Cause-specific Mortality after Allogeneic Unrelated-donor Hematopoietic Cell Transplantation.

Clinical trials commonly use adjudication committees to refine endpoints, but observational research or genome-wide association studies rarely do. Our goals were to establish definitions of cause-specific death after unrelated-donor allogeneic hematopoietic cell transplantation (URD-HCT), to estimate discordance between reported and adjudicated cause-specific death, and to identify factors contributing to inconsistency in cause-specific death determination. A consensus panel adjudicated cause-specific death in 1484 patients who died within 1 year after HCT, derived from 3532 acute leukemia or myelodysplasia patients after URD-HCT from 2000 to 2011 reported by 151 US transplant centers to the Center for International Blood and Marrow Transplant Research. Deaths were classified as disease-related or transplant-related. The panel agreed with >99% of deaths reported by centers as disease-related and 80% reported as transplant-related. Year of transplant (cohort effect) and disease status significantly influenced agreement between the panel and centers. Sensitivity analysis of deaths < 100 days post-transplant yielded the lowest agreement between the panel and centers for myelodysplastic syndrome patients. Standard predefined criteria for adjudicating cause-specific death led to consistent application to similar clinical scenarios and clearer delineation of cause-specific death categories. Other studies of competing events such as cancer-specific versus treatment-related mortality would benefit from our results. Our detailed algorithm should result in more consistent reporting of cause-specific death by centers.

Evolving risk of therapy-related acute myeloid leukemia following cancer chemotherapy among adults in the United States, 1975-2008.

Therapy-related acute myeloid leukemia (tAML) is a rare but highly fatal complication of cytotoxic chemotherapy. Despite major changes in cancer treatment, data describing tAML risks over time are sparse. Among 426068 adults initially treated with chemotherapy for first primary malignancy (9 US population-based cancer registries, 1975-2008), we identified 801 tAML cases, 4.70 times more than expected in the general population (P < .001). Over time, tAML risks increased after chemotherapy for non-Hodgkin lymphoma (n = 158; Poisson regression Ptrend < .001), declined for ovarian cancer (n = 72; Ptrend < .001), myeloma (n = 62; Ptrend = .02), and possibly lung cancer (n = 65; Ptrend = .18), and were significantly heterogeneous for breast cancer (n = 223; Phomogeneity = .005) and Hodgkin lymphoma (n = 58; Phomogeneity = .007). tAML risks varied significantly by age at first cancer and latency and were nonsignificantly heightened with radiotherapy for lung, breast, and ovarian cancers. We identified newly emerging elevated tAML risks in patients treated with chemotherapy since 2000 for esophageal, cervical, prostate, and possibly anal cancers; and since the 1990s for bone/joint and endometrial cancers. Using long-term, population-based data, we observed significant variation in tAML risk with time, consistent with changing treatment practices and differential leukemogenicity of specific therapies. tAML risks should be weighed against the benefits of chemotherapy, particularly for new agents and new indications for standard agents.

A genome-wide meta-analysis of nodular sclerosing Hodgkin lymphoma identifies risk loci at 6p21.32.

Nodular sclerosing Hodgkin lymphoma (NSHL) is a distinct, highly heritable Hodgkin lymphoma subtype. We undertook a genome-wide meta-analysis of 393 European-origin adolescent/young adult NSHL patients and 3315 controls using the Illumina Human610-Quad Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We identified 3 single nucleotide polymorphisms (SNPs) on chromosome 6p21.32 that were significantly associated with NSHL risk: rs9268542 (P = 5.35 × 10(-10)), rs204999 (P = 1.44 × 10(-9)), and rs2858870 (P = 1.69 × 10(-8)). We also confirmed a previously reported association in the same region, rs6903608 (P = 3.52 × 10(-10)). rs204999 and rs2858870 were weakly correlated (r(2) = 0.257), and the remaining pairs of SNPs were not correlated (r(2) < 0.1). In an independent set of 113 NSHL cases and 214 controls, 2 SNPs were significantly associated with NSHL and a third showed a comparable odds ratio (OR). These SNPs are found on 2 haplotypes associated with NSHL risk (rs204999-rs9268528-rs9268542-rs6903608-rs2858870; AGGCT, OR = 1.7, P = 1.71 × 10(-6); GAATC, OR = 0.4, P = 1.16 × 10(-4)). All individuals with the GAATC haplotype also carried the HLA class II DRB1*0701 allele. In a separate analysis, the DRB1*0701 allele was associated with a decreased risk of NSHL (OR = 0.5, 95% confidence interval = 0.4, 0.7). These data support the importance of the HLA class II region in NSHL etiology.

An imbalance between proliferation and differentiation underlies the development of microRNA-defective pineoblastoma.

Mutations in the microRNA processing genes DICER1 and DROSHA drive several cancers that resemble embryonic progenitors. To understand how microRNAs regulate tumorigenesis, we ablated Drosha or Dicer1 in the developing pineal gland to emulate the pathogenesis of pineoblastoma, a brain tumor that resembles undifferentiated precursors of the pineal gland. Accordingly, these mice develop pineal tumors marked by loss of microRNAs, including the let-7/miR-98-5p family, and de-repression of microRNA target genes. Pineal tumors driven by loss of Drosha or Dicer1 mimic tumors driven by Rb1 loss, as they exhibit upregulation of S-phase genes and homeobox transcription factors that regulate pineal development. Blocking proliferation of these tumors facilitates expression of pinealocyte maturation markers, with a concomitant reduction in embryonic markers. Select embryonic markers remain elevated, however, as the microRNAs that normally repress these target genes remain absent. One such microRNA target gene is the oncofetal transcription factor Plagl2 , which regulates expression of pro-growth genes, and inhibiting their signaling impairs tumor growth. Thus, we demonstrate that tumors driven by loss of microRNA processing may be therapeutically targeted by inhibiting downstream drivers of proliferation.

Multiple Myeloma Risk and Outcomes Are Associated with Pathogenic Germline Variants in DNA Repair Genes.

First-degree relatives of patients with multiple myeloma are at increased risk for the disease, but the contribution of pathogenic germline variants (PGV) in hereditary cancer genes to multiple myeloma risk and outcomes is not well characterized. To address this, we analyzed germline exomes in two independent cohorts of 895 and 786 patients with multiple myeloma. PGVs were identified in 8.6% of the Discovery cohort and 11.5% of the Replication cohort, with a notable presence of high- or moderate-penetrance PGVs (associated with autosomal dominant cancer predisposition) in DNA repair genes (3.6% and 4.1%, respectively). PGVs in BRCA1 (OR = 3.9, FDR < 0.01) and BRCA2 (OR = 7.0, FDR < 0.001) were significantly enriched in patients with multiple myeloma when compared with 134,187 healthy controls. Five of the eight BRCA2 PGV carriers exhibited tumor-specific copy number loss in BRCA2, suggesting somatic loss of heterozygosity. PGVs associated with autosomal dominant cancer predisposition were associated with younger age at diagnosis, personal or familial cancer history, and longer progression-free survival after upfront high-dose melphalan and autologous stem-cell transplantation (P < 0.01). Significance: Our findings suggest up to 10% of patients with multiple myeloma may have an unsuspected cancer predisposition syndrome. Given familial implications and favorable outcomes with high-dose melphalan and autologous stem-cell transplantation in high-penetrance PGV carriers, genetic testing should be considered for young or newly diagnosed patients with a personal or family cancer history. See related commentary by Walker, p. 375.

Natural selection and functional genetic variation in the p53 pathway.

The allele frequencies of two functional single-nucleotide polymorphisms (SNPs) in the p53 pathway, the MDM2 SNP309 and TP53 Arg72Pro, vary dramatically among populations. That the frequencies of the TP53 SNP follow a clinal distribution may suggest that selective pressure from environmental variables correlated with latitude contributed to these observed population differences. Recently, winter temperature and UV radiation were found to be significantly correlated with the TP53 and the MDM2 SNPs, respectively, in East Asians; whether these correlations are more extreme than expected based upon nonselective factors such as patterns of human migration remains unclear. Here, we genotyped these two SNPs in 971 unrelated individuals from 52 unique populations worldwide and tested for correlations with both latitude and a number of climate-related environmental variables on a global scale, controlling for these neutral processes. The TP53 SNP was associated with a significant selection signal for a few climate variables, such as short-wave radiation flux in the winter, but these signals were no longer significant after correction for multiple tests. The MDM2 SNP did not exhibit a significant signal with any climate variable. Therefore, these SNPs are unlikely to be under selective pressure driven by these variables. Thus, these data underscore the need to incorporate population history when assessing signatures of selection.