From the operating table to global science: How a near-death experience sparked my passion for life.

Kevin Christopher Brown Jr. is a winner of the fourth annual Rising Black Scientists Awards for a scholar in the life and health sciences. We asked emerging Black scientists to tell us about their scientific vision and goals, experiences that sparked their interest in science, how they want to contribute to a more inclusive scientific community, and how these all fit together on their journey. This is his story.

Colorectal Cancer Cells Enter a Diapause-like DTP State to Survive Chemotherapy.

Cancer cells enter a reversible drug-tolerant persister (DTP) state to evade death from chemotherapy and targeted agents. It is increasingly appreciated that DTPs are important drivers of therapy failure and tumor relapse. We combined cellular barcoding and mathematical modeling in patient-derived colorectal cancer models to identify and characterize DTPs in response to chemotherapy. Barcode analysis revealed no loss of clonal complexity of tumors that entered the DTP state and recurred following treatment cessation. Our data fit a mathematical model where all cancer cells, and not a small subpopulation, possess an equipotent capacity to become DTPs. Mechanistically, we determined that DTPs display remarkable transcriptional and functional similarities to diapause, a reversible state of suspended embryonic development triggered by unfavorable environmental conditions. Our study provides insight into how cancer cells use a developmentally conserved mechanism to drive the DTP state, pointing to novel therapeutic opportunities to target DTPs.

Children develop robust and sustained cross-reactive spike-specific immune responses to SARS-CoV-2 infection.

SARS-CoV-2 infection is generally mild or asymptomatic in children but a biological basis for this outcome is unclear. Here we compare antibody and cellular immunity in children (aged 3-11 years) and adults. Antibody responses against spike protein were high in children and seroconversion boosted responses against seasonal Beta-coronaviruses through cross-recognition of the S2 domain. Neutralization of viral variants was comparable between children and adults. Spike-specific T cell responses were more than twice as high in children and were also detected in many seronegative children, indicating pre-existing cross-reactive responses to seasonal coronaviruses. Importantly, children retained antibody and cellular responses 6 months after infection, whereas relative waning occurred in adults. Spike-specific responses were also broadly stable beyond 12 months. Therefore, children generate robust, cross-reactive and sustained immune responses to SARS-CoV-2 with focused specificity for the spike protein. These findings provide insight into the relative clinical protection that occurs in most children and might help to guide the design of pediatric vaccination regimens.

Robust SARS-CoV-2-specific T cell immunity is maintained at 6 months following primary infection.

The immune response to SARS-CoV-2 is critical in controlling disease, but there is concern that waning immunity may predispose to reinfection. We analyzed the magnitude and phenotype of the SARS-CoV-2-specific T cell response in 100 donors at 6 months following infection. T cell responses were present by ELISPOT and/or intracellular cytokine staining analysis in all donors and characterized by predominant CD4+ T cell responses with strong interleukin (IL)-2 cytokine expression. Median T cell responses were 50% higher in donors who had experienced a symptomatic infection, indicating that the severity of primary infection establishes a 'set point' for cellular immunity. T cell responses to spike and nucleoprotein/membrane proteins were correlated with peak antibody levels. Furthermore, higher levels of nucleoprotein-specific T cells were associated with preservation of nucleoprotein-specific antibody level although no such correlation was observed in relation to spike-specific responses. In conclusion, our data are reassuring that functional SARS-CoV-2-specific T cell responses are retained at 6 months following infection.

Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance.

Large-scale genomic studies have identified multiple somatic aberrations in breast cancer, including copy number alterations and point mutations. Still, identifying causal variants and emergent vulnerabilities that arise as a consequence of genetic alterations remain major challenges. We performed whole-genome small hairpin RNA (shRNA) "dropout screens" on 77 breast cancer cell lines. Using a hierarchical linear regression algorithm to score our screen results and integrate them with accompanying detailed genetic and proteomic information, we identify vulnerabilities in breast cancer, including candidate "drivers," and reveal general functional genomic properties of cancer cells. Comparisons of gene essentiality with drug sensitivity data suggest potential resistance mechanisms, effects of existing anti-cancer drugs, and opportunities for combination therapy. Finally, we demonstrate the utility of this large dataset by identifying BRD4 as a potential target in luminal breast cancer and PIK3CA mutations as a resistance determinant for BET-inhibitors.

High-Resolution CRISPR Screens Reveal Fitness Genes and Genotype-Specific Cancer Liabilities.

The ability to perturb genes in human cells is crucial for elucidating gene function and holds great potential for finding therapeutic targets for diseases such as cancer. To extend the catalog of human core and context-dependent fitness genes, we have developed a high-complexity second-generation genome-scale CRISPR-Cas9 gRNA library and applied it to fitness screens in five human cell lines. Using an improved Bayesian analytical approach, we consistently discover 5-fold more fitness genes than were previously observed. We present a list of 1,580 human core fitness genes and describe their general properties. Moreover, we demonstrate that context-dependent fitness genes accurately recapitulate pathway-specific genetic vulnerabilities induced by known oncogenes and reveal cell-type-specific dependencies for specific receptor tyrosine kinases, even in oncogenic KRAS backgrounds. Thus, rigorous identification of human cell line fitness genes using a high-complexity CRISPR-Cas9 library affords a high-resolution view of the genetic vulnerabilities of a cell.

Systematic exploration of dynamic splicing networks reveals conserved multistage regulators of neurogenesis.

Alternative splicing (AS) is a critical regulatory layer; yet, factors controlling functionally coordinated splicing programs during developmental transitions are poorly understood. Here, we employ a screening strategy to identify factors controlling dynamic splicing events important for mammalian neurogenesis. Among previously unknown regulators, Rbm38 acts widely to negatively control neural AS, in part through interactions mediated by the established repressor of splicing, Ptbp1. Puf60, a ubiquitous factor, is surprisingly found to promote neural splicing patterns. This activity requires a conserved, neural-differential exon that remodels Puf60 co-factor interactions. Ablation of this exon rewires distinct AS networks in embryonic stem cells and at different stages of mouse neurogenesis. Single-cell transcriptome analyses further reveal distinct roles for Rbm38 and Puf60 isoforms in establishing neuronal identity. Our results describe important roles for previously unknown regulators of neurogenesis and establish how an alternative exon in a widely expressed splicing factor orchestrates temporal control over cell differentiation.

Functional genomic landscape of cancer-intrinsic evasion of killing by T cells.

The genetic circuits that allow cancer cells to evade destruction by the host immune system remain poorly understood1-3. Here, to identify a phenotypically robust core set of genes and pathways that enable cancer cells to evade killing mediated by cytotoxic T lymphocytes (CTLs), we performed genome-wide CRISPR screens across a panel of genetically diverse mouse cancer cell lines that were cultured in the presence of CTLs. We identify a core set of 182 genes across these mouse cancer models, the individual perturbation of which increases either the sensitivity or the resistance of cancer cells to CTL-mediated toxicity. Systematic exploration of our dataset using genetic co-similarity reveals the hierarchical and coordinated manner in which genes and pathways act in cancer cells to orchestrate their evasion of CTLs, and shows that discrete functional modules that control the interferon response and tumour necrosis factor (TNF)-induced cytotoxicity are dominant sub-phenotypes. Our data establish a central role for genes that were previously identified as negative regulators of the type-II interferon response (for example, Ptpn2, Socs1 and Adar1) in mediating CTL evasion, and show that the lipid-droplet-related gene Fitm2 is required for maintaining cell fitness after exposure to interferon-γ (IFNγ). In addition, we identify the autophagy pathway as a conserved mediator of the evasion of CTLs by cancer cells, and show that this pathway is required to resist cytotoxicity induced by the cytokines IFNγ and TNF. Through the mapping of cytokine- and CTL-based genetic interactions, together with in vivo CRISPR screens, we show how the pleiotropic effects of autophagy control cancer-cell-intrinsic evasion of killing by CTLs and we highlight the importance of these effects within the tumour microenvironment. Collectively, these data expand our knowledge of the genetic circuits that are involved in the evasion of the immune system by cancer cells, and highlight genetic interactions that contribute to phenotypes associated with escape from killing by CTLs.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity.

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

Duration of Protection against Mild and Severe Disease by Covid-19 Vaccines.

BACKGROUND: Vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (Covid-19), have been used since December 2020 in the United Kingdom. Real-world data have shown the vaccines to be highly effective against Covid-19 and related severe disease and death. Vaccine effectiveness may wane over time since the receipt of the second dose of the ChAdOx1-S (ChAdOx1 nCoV-19) and BNT162b2 vaccines. METHODS: We used a test-negative case-control design to estimate vaccine effectiveness against symptomatic Covid-19 and related hospitalization and death in England. Effectiveness of the ChAdOx1-S and BNT162b2 vaccines was assessed according to participant age and status with regard to coexisting conditions and over time since receipt of the second vaccine dose to investigate waning of effectiveness separately for the B.1.1.7 (alpha) and B.1.617.2 (delta) variants. RESULTS: Vaccine effectiveness against symptomatic Covid-19 with the delta variant peaked in the early weeks after receipt of the second dose and then decreased by 20 weeks to 44.3% (95% confidence interval [CI], 43.2 to 45.4) with the ChAdOx1-S vaccine and to 66.3% (95% CI, 65.7 to 66.9) with the BNT162b2 vaccine. Waning of vaccine effectiveness was greater in persons 65 years of age or older than in those 40 to 64 years of age. At 20 weeks or more after vaccination, vaccine effectiveness decreased less against both hospitalization, to 80.0% (95% CI, 76.8 to 82.7) with the ChAdOx1-S vaccine and 91.7% (95% CI, 90.2 to 93.0) with the BNT162b2 vaccine, and death, to 84.8% (95% CI, 76.2 to 90.3) and 91.9% (95% CI, 88.5 to 94.3), respectively. Greater waning in vaccine effectiveness against hospitalization was observed in persons 65 years of age or older in a clinically extremely vulnerable group and in persons 40 to 64 years of age with underlying medical conditions than in healthy adults. CONCLUSIONS: We observed limited waning in vaccine effectiveness against Covid-19-related hospitalization and death at 20 weeks or more after vaccination with two doses of the ChAdOx1-S or BNT162b2 vaccine. Waning was greater in older adults and in those in a clinical risk group.

Covid-19 Vaccine Effectiveness against the Omicron (B.1.1.529) Variant.

BACKGROUND: A rapid increase in coronavirus disease 2019 (Covid-19) cases due to the omicron (B.1.1.529) variant of severe acute respiratory syndrome coronavirus 2 in highly vaccinated populations has aroused concerns about the effectiveness of current vaccines. METHODS: We used a test-negative case-control design to estimate vaccine effectiveness against symptomatic disease caused by the omicron and delta (B.1.617.2) variants in England. Vaccine effectiveness was calculated after primary immunization with two doses of BNT162b2 (Pfizer-BioNTech), ChAdOx1 nCoV-19 (AstraZeneca), or mRNA-1273 (Moderna) vaccine and after a booster dose of BNT162b2, ChAdOx1 nCoV-19, or mRNA-1273. RESULTS: Between November 27, 2021, and January 12, 2022, a total of 886,774 eligible persons infected with the omicron variant, 204,154 eligible persons infected with the delta variant, and 1,572,621 eligible test-negative controls were identified. At all time points investigated and for all combinations of primary course and booster vaccines, vaccine effectiveness against symptomatic disease was higher for the delta variant than for the omicron variant. No effect against the omicron variant was noted from 20 weeks after two ChAdOx1 nCoV-19 doses, whereas vaccine effectiveness after two BNT162b2 doses was 65.5% (95% confidence interval [CI], 63.9 to 67.0) at 2 to 4 weeks, dropping to 8.8% (95% CI, 7.0 to 10.5) at 25 or more weeks. Among ChAdOx1 nCoV-19 primary course recipients, vaccine effectiveness increased to 62.4% (95% CI, 61.8 to 63.0) at 2 to 4 weeks after a BNT162b2 booster before decreasing to 39.6% (95% CI, 38.0 to 41.1) at 10 or more weeks. Among BNT162b2 primary course recipients, vaccine effectiveness increased to 67.2% (95% CI, 66.5 to 67.8) at 2 to 4 weeks after a BNT162b2 booster before declining to 45.7% (95% CI, 44.7 to 46.7) at 10 or more weeks. Vaccine effectiveness after a ChAdOx1 nCoV-19 primary course increased to 70.1% (95% CI, 69.5 to 70.7) at 2 to 4 weeks after an mRNA-1273 booster and decreased to 60.9% (95% CI, 59.7 to 62.1) at 5 to 9 weeks. After a BNT162b2 primary course, the mRNA-1273 booster increased vaccine effectiveness to 73.9% (95% CI, 73.1 to 74.6) at 2 to 4 weeks; vaccine effectiveness fell to 64.4% (95% CI, 62.6 to 66.1) at 5 to 9 weeks. CONCLUSIONS: Primary immunization with two doses of ChAdOx1 nCoV-19 or BNT162b2 vaccine provided limited protection against symptomatic disease caused by the omicron variant. A BNT162b2 or mRNA-1273 booster after either the ChAdOx1 nCoV-19 or BNT162b2 primary course substantially increased protection, but that protection waned over time. (Funded by the U.K. Health Security Agency.).

Meaningful Endpoints for Idiopathic Pulmonary Fibrosis (IPF) Clinical Trials: Emphasis on 'Feels, Functions, Survives'. Report of a Collaborative Discussion in a Symposium with Direct Engagement from Representatives of Patients, Investigators, the National Institutes of Health, a Patient Advocacy Organization, and a Regulatory Agency.

Background: Idiopathic pulmonary fibrosis (IPF) carries significant mortality and unpredictable progression, with limited therapeutic options. Designing trials with patient-meaningful endpoints, enhancing the reliability and interpretability of results, and streamlining the regulatory approval process are of critical importance to advancing clinical care in IPF. Methods: A landmark in-person symposium in June 2023 assembled 43 participants from the US and internationally, including patients with IPF, investigators, and regulatory representatives, to discuss the immediate future of IPF clinical trial endpoints. Patient advocates were central to discussions, which evaluated endpoints according to regulatory standards and the FDA's 'feels, functions, survives' criteria. Results: Three themes emerged: 1) consensus on endpoints mirroring the lived experiences of patients with IPF; 2) consideration of replacing forced vital capacity (FVC) as the primary endpoint, potentially by composite endpoints that include 'feels, functions, survives' measures or FVC as components; 3) support for simplified, user-friendly patient-reported outcomes (PROs) as either components of primary composite endpoints or key secondary endpoints, supplemented by functional tests as secondary endpoints and novel biomarkers as supportive measures (FDA Guidance for Industry (Multiple Endpoints in Clinical Trials) available at: https://www.fda.gov/media/162416/download). Conclusions: This report, detailing the proceedings of this pivotal symposium, suggests a potential turning point in designing future IPF clinical trials more attuned to outcomes meaningful to patients, and documents the collective agreement across multidisciplinary stakeholders on the importance of anchoring IPF trial endpoints on real patient experiences-namely, how they feel, function, and survive. There is considerable optimism that clinical care in IPF will progress through trials focused on patient-centric insights, ultimately guiding transformative treatment strategies to enhance patients' quality of life and survival.

Effectiveness of mRNA COVID-19 Monovalent and Bivalent Vaccine Booster Doses Against Omicron Severe Outcomes Among Adults Aged ≥50 Years in Ontario, Canada: A Canadian Immunization Research Network Study.

We estimated the effectiveness of booster doses of monovalent and bivalent mRNA COVID-19 vaccines against Omicron-associated severe outcomes among adults aged ≥50 years in Ontario, Canada. Monovalent and bivalent mRNA COVID-19 booster doses provided similar strong initial protection against severe outcomes. Uncertainty remains around waning of protection from these vaccines.

The Effect of COVID-19 Vaccination on Outpatient Antibiotic Prescribing in Older Adults: A Self-Controlled Risk-Interval Study.

BACKGROUND: Coronavirus disease 2019 (COVID-19) vaccination has been associated with reduced outpatient antibiotic prescribing among older adults with laboratory-confirmed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We assessed the impact of COVID-19 vaccination on outpatient antibiotic prescribing in the broader population of older adults, regardless of SARS-CoV-2 infection status. METHODS: We included adults aged ≥65 years who received their first, second, and/or third COVID-19 vaccine dose from December 2020 to December 2022. We used a self-controlled risk-interval design and included cases who received an antibiotic prescription 2-6 weeks before vaccination (pre-vaccination or control interval) or after vaccination (post-vaccination or risk interval). We used conditional logistic regression to estimate the odds of being prescribed (1) any antibiotic, (2) a typical "respiratory" infection antibiotic, or (3) a typical "urinary tract" infection antibiotic (negative control) in the post-vaccination interval versus the pre-vaccination interval. We accounted for temporal changes in antibiotic prescribing using background monthly antibiotic prescribing counts. RESULTS: 469 923 vaccine doses met inclusion criteria. The odds of receiving any antibiotic or a respiratory antibiotic prescription were lower in the post-vaccination versus pre-vaccination interval (aOR, .973; 95% CI, .968-.978; aOR, .961; 95% CI, .953-.968, respectively). There was no association between vaccination and urinary antibiotic prescriptions (aOR, .996; 95% CI, .987-1.006). Periods with high (>10%) versus low (<5%) SARS-CoV-2 test positivity demonstrated greater reductions in antibiotic prescribing (aOR, .875; 95% CI, .845-.905; aOR, .996; 95% CI, .989-1.003, respectively). CONCLUSIONS: COVID-19 vaccination was associated with reduced outpatient antibiotic prescribing in older adults, especially during periods of high SARS-CoV-2 circulation.

Association between infectious diseases consultation and mortality in hospitalized patients with Gram-negative bloodstream infection: a retrospective population-wide cohort study.

OBJECTIVES: Data supporting routine infectious diseases (ID) consultation in Gram-negative bloodstream infection (GN-BSI) are limited. We evaluated the association between ID consultation and mortality in patients with GN-BSI in a retrospective population-wide cohort study in Ontario using linked health administrative databases. METHODS: Hospitalized adult patients with GN-BSI between April 2017 and December 2021 were included. The primary outcome was time to all-cause mortality censored at 30 days, analyzed using a mixed effects Cox proportional hazards model with hospital as a random effect. ID consultation 1-10 days after the first positive blood culture was treated as a time-varying exposure. RESULTS: Of 30,159 patients with GN-BSI across 53 hospitals, 11,013 (36.5%) received ID consultation. Median prevalence of ID consultation for patients with GN-BSI across hospitals was 35.0% with wide variability (range 2.7-76.1%, interquartile range 19.6-41.1%). 1041 (9.5%) patients who received ID consultation died within 30 days, compared to 1797 (9.4%) patients without ID consultation. In the fully-adjusted multivariable model, ID consultation was associated with mortality benefit (adjusted HR 0.82, 95% CI 0.77-0.88, p < 0.0001; translating to absolute risk reduction of -3.8% or NNT of 27). Exploratory subgroup analyses of the primary outcome showed that ID consultation could have greater benefit in patients with high-risk features (nosocomial infection, polymicrobial or non-Enterobacterales infection, antimicrobial resistance, or non-urinary tract source). CONCLUSIONS: Early ID consultation was associated with reduced mortality in patients with GN-BSI. If resources permit, routine ID consultation for this patient population should be considered to improve patient outcomes.

Rare germline deleterious variants increase susceptibility for lung cancer.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci.

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1.

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

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

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