The Circadian Rhythm Gene <i>Arntl2</i> Is a Metastasis Susceptibility Gene for Estrogen Receptor-Negative Breast Cancer

<div><p>Breast cancer mortality is primarily due to metastasis rather than primary tumors, yet relatively little is understood regarding the etiology of metastatic breast cancer. Previously, using a mouse genetics approach, we demonstrated that inherited germline polymorphisms contribute to metastatic disease, and that these single nucleotide polymorphisms (SNPs) could be used to predict outcome in breast cancer patients. In this study, a backcross between a highly metastatic (FVB/NJ) and low metastatic (MOLF/EiJ) mouse strain identified <i>Arntl2</i>, a gene encoding a circadian rhythm transcription factor, as a metastasis susceptibility gene associated with progression, specifically in estrogen receptor-negative breast cancer patients. Integrated whole genome sequence analysis with DNase hypersensitivity sites reveals SNPs in the predicted promoter of <i>Arntl2</i>. Using CRISPR/Cas9-mediated substitution of the MOLF promoter, we demonstrate that the SNPs regulate <i>Arntl2</i> transcription and affect metastatic burden. Finally, analysis of SNPs associated with <i>ARNTL2</i> expression in human breast cancer patients revealed reproducible associations of <i>ARNTL2</i> expression quantitative trait loci (eQTL) SNPs with disease-free survival, consistent with the mouse studies.</p></div

Human eQTL SNPs associated with ARNTL2 expression that are prognostic for ER- breast cancer.

<p>Human eQTL SNPs associated with ARNTL2 expression that are prognostic for ER- breast cancer.</p

Distal chromosome 6 is a metastasis modifier.

<p>(A) Survival plot of FVB/NJ x (MOLF/EiJ x MMTV-PyMT) mice. p-value calculated by Log-Rank test. FVB = FVB/NJ x MMTV-PyMT; MOLF = MOLF/EiJ x MMTV-PyMT. (B) Surface metastasis count of FVB/NJ x (MOLF/EiJ x MMTV-PyMT) mice. p-value calculated by Mann-Whitney test. FVB = FVB/NJ x MMTV-PyMT; MOLF = MOLF/EiJ x MMTV-PyMT. (C) Depiction of breeding scheme. (D) LOD score plot for metastasis susceptibility, latency and tumor burden loci in the backcross. The y-axis represents the LOD score observed across the genome. The mouse chromosomes are presented head-to-tail on the x-axis. The lower horizontal dashed line represents the suggestive LOD score for a quantitative trait locus after correction for multiple testing by permutation testing. The upper dashed line is the significance threshold after permutation testing. (E) Approximate location of metastasis-associated eQTL candidate genes on distal chromosome 6. (F) The Gene expression-based Outcome for Breast Cancer Online (GOBO) database was queried for the eQTL genes on distal chromosome 6. Distant metastasis-free survival (DMFS) plotted as Kaplan-Meier curves for patients with ER+ and ER- tumors expressing high (blue), intermediate (red), or low (gray) levels of the eQTL genes.</p

Tumor-autonomous effect of Arntl2.

<p>(A) <i>Arntl2</i>^+/+ and <i>Arntl2</i>^-/- mice were bred to MMTV-PyMT+ mice to generate spontaneous primary tumors and lung metastases. At the time of euthanasia, autochthonous PyMT-induced tumors were weighed. N.S. = not significant. (B) Surface lung metastasis count of spontaneous metastases from (A). P-value calculated by two-tailed Mann-Whitney test (C) Primary tumor weight of orthotopically injected 4T1 cells into F1 hybrids from <i>Arntl2</i>^+/+ and Arntl2^+/- mice crossed with BALB/cJ mice. N.S. = not significant. (D) Surface lung metastasis count of the injected mice in (C). N.S. = not significant.</p

SNPs within the <i>Arntl2</i> promoter modulate chromatin binding.

<p>(A) Differential binding of MOLF and FVB promoter DNA probes to chromatin proteins. Relative densitometry compared to beads only is reported above each sample. Densitometry was calculated using ImageJ. (B) Relative mRNA expression of <i>Stkl3</i> in wildtype and heterozygous CRISPR 6DT1 cells. (C) Working model depicting differential chromatin binding of MOLF and FVB to regions of the nucleus. LADs = laminar-associated domains (repressive state).</p

Genome editing of the <i>Arntl2</i> promoter.

<p>(A) Haplotypes for the <i>Arntl2</i> putative promoter from the FVB, MOLF, PWK and WSB inbred strains. Brain RNAseq data shows that wild-derived PWK/PhJ carries the same haplotype as MOLF/EiJ, and WSB/EiJ shares 9 of the 10 FVB/NJ SNPs at this location. (B) Relative expression of <i>Arntl2</i> in brain RNA from the PWK and WSB strains. (C) Confirmation of genomically-edited 6DT1 cells. (D) Relative mRNA expression of <i>Arntl2</i> in CRISPR-engineered 6DT1 cells as measured by qPCR. (E) Primary tumor weight of injected CRISPR-engineered WT and heterozygous substituted cells in FVB mice. N.S. = not significant. (F) Surface lung metastasis count on injected cells in (E). P-value calculated by two-tailed Mann-Whitney test.</p

Effect of Arntl2 expression <i>in vivo</i>.

<p>(A) Relative <i>Arntl2</i> mRNA expression in 4T1 cells with control (scramble) and <i>Arntl2</i> shRNA cell lines. (B) Protein expression of ARNTL2 in transiently transfected HEK293 cells with shRNA and overexpression vector expressing myc-tag. Actin is the loading control. (C) Weight of primary tumor from orthotopically injected BALB/cJ mice with 4T1 control and C7 and C11 <i>Arntl2</i> shRNA constructs. N.S. = not significant. (D) Count of surface pulmonary nodules from (C). P-value calculated by two-tailed Mann-Whitney test. (E) Primary tumor weight of <i>Arntl2</i> overexpressing 4T1 cells. N.S. = not significant. (F) Lung metastases from (E), p-value calculated by one-tailed Mann-Whitney test.</p

Association of Urinary N7-(1-hydroxyl-3-buten-1-yl) Guanine (EB-GII) Adducts and Butadiene-Mercapturic Acids with Lung Cancer Development in Cigarette Smokers

Approximately 10% of smokers will develop lung cancer. Sensitive predictive biomarkers are needed to identify susceptible individuals. 1,3-Butadiene (BD) is among the most abundant tobacco smoke carcinogens. BD is metabolically activated to 3,4-epoxy-1-butene (EB), which is detoxified via the glutathione conjugation/mercapturic acid pathway to form monohydroxybutenyl mercapturic acid (MHBMA) and dihydroxybutyl mercapturic acid (DHBMA). Alternatively, EB can react with guanine nucleobases of DNA to form N7-(1-hydroxyl-3-buten-1-yl) guanine (EB-GII) adducts. We employed isotope dilution LC/ESI-HRMS/MS methodologies to quantify MHBMA, DHBMA, and EB-GII in urine of smokers who developed lung cancer (N = 260) and matched smoking controls (N = 259) from the Southern Community Cohort (white and African American). The concentrations of all three biomarkers were significantly higher in smokers that subsequently developed lung cancer as compared to matched smoker controls after adjusting for age, sex, and race/ethnicity (p p p = 0.0007 for DHBMA). The odds ratio (OR) for lung cancer development was 1.63 for MHBMA, 1.37 for DHBMA, and 1.97 for EB-GII, with a higher OR in African American subjects than in whites. The association of urinary EB-GII, MHBMA, and DHBMA with lung cancer status did not remain upon adjustment for total nicotine equivalents. These findings reveal that urinary MHBMA, DHBMA, and EB-GII are directly correlated with the BD dose delivered via smoking and are associated with lung cancer risk

Effects of <i>BCL2</i> gene variants on associations between risk scores for established risk factors and endometrial cancer risk, the Shanghai Endometrial Cancer Genetics Study (SECGS).

*<p>Adjusted for age (continuous) and education.</p

Association of genetic variants in BCL-2-family genes with endometrial cancer risk among 1,028 cases and 1,922 controls, the Shanghai Endometrial Cancer Genetics Study (SECGS).

a<p>Major/minor alleles.</p>b<p>Minor allele frequency in controls.</p>c<p>Adjusted for age and education. SNPs associated with endometrial cancer risk (<i>P</i> value<0.05) are represented in bold. <i>P</i> values were not corrected for multiple tests.</p><p><i>Notes</i>:</p><p>-<i>AA</i> indicates homozygotes for the major allele, who were used as the reference group (not shown in the table), <i>Aa</i> indicates heterozygotes and <i>aa</i> indicates homozygotes for the minor allele.</p><p>-No deviation from Hardy-Weinberg Equilibrium (HWE) was observed for any given SNP among controls.</p><p>-5UTR, 5′ untranslated region; 3′ UTR, 3′ untranslated region.</p><p>-Significant SNPs are presented in bold.</p