Nutritional Control of L1 Arrest and Recovery in Caenorhabditis elegans by Insulin-like Peptides and Signaling

Animals must coordinate development with fluctuating nutrient availability. Nutrient availability governs post-embryonic development in Caenorhabditis elegans: larvae that hatch in the absence of food do not initiate post-embryonic development but enter "L1 arrest" (or "L1 diapause") and can survive starvation for weeks, while rapidly resume normal development once get fed. Insulin-like signaling (IIS) has been shown to be a key regulator of L1 arrest and recovery. However, the C. elegans genome encodes 40 insulin-like peptides (ILPs), and it is unknown which peptides participate in nutritional control of L1 arrest and recovery. Work in other contexts has identified putative receptor agonists and antagonists, but the extent of specificity versus redundancy is unclear beyond this distinction.

We measured mRNA expression dynamics with high temporal resolution for all 40 insulin-like genes during entry into and recovery from L1 arrest. Nutrient availability influences expression of the majority of insulin-like genes, with variable dynamics suggesting complex regulation. We identified 13 candidate agonists and 8 candidate antagonists based on expression in response to nutrient availability. We selected ten candidate agonists (daf-28, ins-3, ins-4, ins-5, ins-6, ins-7, ins-9, ins-26, ins-33 and ins-35) for further characterization in L1 stage larvae. We used destabilized reporter genes to determine spatial expression patterns. Expression of candidate agonists was largely overlapping in L1 stage larvae, suggesting a role of the intestine, chemosensory neurons ASI and ASJ, and the interneuron PVT in systemic control of L1 development. Transcriptional regulation of candidate agonists was most significant in the intestine, as if nutrient uptake was a more important influence on transcription than sensory perception. Scanning in the 5' upstream promoter region of these 40 ILPs, We found that transcription factor PQM-1 and GATA putative binding sites are depleted in the promoter region of antagonists. A novel motif was also found to be over-represented in ILPs.

Phenotypic analysis of single and compound deletion mutants did not reveal effects on L1 recovery/developmental dynamics, though simultaneous disruption of ins-4 and daf-28 extended survival of L1 arrest without enhancing thermal tolerance, while overexpression of ins-4, ins-6 or daf-28 shortened L1 survival. Simultaneous disruption of several ILPs showed a temperature independent, transient dauer phenotype. These results revealed the relative redundancy and specificity among agonistic ILPs.

TGF- β and steroid hormone (SH) signaling have been reported to control the dauer formation along with IIS. Our preliminary results suggest they may also mediate the IIS control of L1 arrest and recovery, as the expression of several key components of TGF-β and SH signaling pathway genes are negatively regulated by DAF-16, and loss-of-function of these genes partially represses daf-16 null phenotype in L1 arrest, and causes a retardation in L1 development.

In summary, my dissertation study focused on the IIS, characterized the dynamics and sites of ILPs expression in response to nutrient availability, revealed the function of specific agonistic ILPs in L1 arrest, and suggested potential cross-regulation among IIS, TGF-β signaling and SH signaling in controlling L1 arrest and recovery. These findings provide insights into how post-embryonic development is governed by insulin-like signaling and nutrient availability.

Ins-4 and daf-28 function redundantly to regulate C. elegans L1 arrest

© 2014 The Authors.Caenorhabditis elegans larvae reversibly arrest development in the first larval stage in response to starvation (L1 arrest or L1 diapause). Insulin-like signaling is a critical regulator of L1 arrest. However, the C. elegans genome encodes 40 insulin-like peptides, and it is unknown which peptides participate in nutritional control of L1 development. Work in other contexts has revealed that insulin-like genes can promote development ("agonists") or developmental arrest ("antagonists"), suggesting that such agonists promote L1 development in response to feeding. We measured mRNA expression dynamics with high temporal resolution for all 40 insulin-like genes during entry into and recovery from L1 arrest. Nutrient availability influences expression of the majority of insulin-like genes, with variable dynamics suggesting complex regulation. We identified thirteen candidate agonists and eight candidate antagonists based on expression in response to nutrient availability. We selected ten candidate agonists (. daf-28, ins-3, ins-4, ins-5, ins-6, ins-7, ins-9, ins-26, ins-33 and ins-35) for further characterization in L1 stage larvae. We used destabilized reporter genes to determine spatial expression patterns. Expression of candidate agonists is largely overlapping in L1 stage larvae, suggesting a role of the intestine, chemosensory neurons ASI and ASJ, and the interneuron PVT in control of L1 development. Transcriptional regulation of candidate agonists is most significant in the intestine, as if internal nutrient status is a more important influence on transcription than sensory perception. Phenotypic analysis of single and compound deletion mutants did not reveal effects on L1 developmental dynamics, though simultaneous disruption of ins-4 and daf-28 increases survival of L1 arrest. Furthermore, overexpression of ins-4, ins-6 or daf-28 alone decreases survival and promotes cell division during starvation. These results suggest extensive functional overlap among insulin-like genes in nutritional control of L1 development while highlighting the role of ins-4, daf-28 and to a lesser extent ins-6.

Unfolded protein response genes regulated by CED-1 are required for Caenorhabditis elegans innate immunity.

The endoplasmic reticulum stress response, also known as the unfolded protein response (UPR), has been implicated in the normal physiology of immune defense and in several disorders, including diabetes, cancer, and neurodegenerative disease. Here, we show that the apoptotic receptor CED-1 and a network of PQN/ABU proteins involved in a noncanonical UPR response are required for proper defense to pathogen infection in Caenorhabditis elegans. A full-genome microarray analysis indicates that CED-1 functions to activate the expression of pqn/abu genes. We also show that ced-1 and pqn/abu genes are required for the survival of C. elegans exposed to live Salmonella enterica, and that overexpression of pqn/abu genes confers protection against pathogen-mediated killing. The results indicate that unfolded protein response genes, regulated in a CED-1-dependent manner, are involved in the C. elegans immune response to live bacteria.

Polymorphisms in the ACE and ADRB2 genes and risks of aging-associated phenotypes: the case of myocardial infarction.

Multiple functions of the beta2-adrenergic receptor (ADRB2) and angiotensin-converting enzyme (ACE) genes warrant studies of their associations with aging-related phenotypes. We focus on multimarker analyses and analyses of the effects of compound genotypes of two polymorphisms in the ADRB2 gene, rs1042713 and rs1042714, and 11 polymorphisms of the ACE gene, on the risk of such an aging-associated phenotype as myocardial infarction (MI). We used the data from a genotyped sample of the Framingham Heart Study Offspring (FHSO) cohort (n = 1500) followed for about 36 years with six examinations. The ADRB2 rs1042714 (C-->G) polymorphism and two moderately correlated (r(2) = 0.77) ACE polymorphisms, rs4363 (A-->G) and rs12449782 (A-->G), were significantly associated with risks of MI in this aging cohort in multimarker models. Predominantly linked ACE genotypes exhibited opposite effects on MI risks, e.g., the AA (rs12449782) genotype had a detrimental effect, whereas the predominantly linked AA (rs4363) genotype exhibited a protective effect. This trade-off occurs as a result of the opposite effects of rare compound genotypes of the ACE polymorphisms with a single dose of the AG heterozygote. This genetic trade-off is further augmented by the selective modulating effect of the rs1042714 ADRB2 polymorphism. The associations were not altered by adjustment for common MI risk factors. The results suggest that effects of single specific genetic variants of the ADRB2 and ACE genes on MI can be readily altered by gene-gene or/and gene-environmental interactions, especially in large heterogeneous samples. Multimarker genetic analyses should benefit studies of complex aging-associated phenotypes.

Recursive directional ligation by plasmid reconstruction allows rapid and seamless cloning of oligomeric genes.

This paper reports a new strategy, recursive directional ligation by plasmid reconstruction (PRe-RDL), to rapidly clone highly repetitive polypeptides of any sequence and specified length over a large range of molecular weights. In a single cycle of PRe-RDL, two halves of a parent plasmid, each containing a copy of an oligomer, are ligated together, thereby dimerizing the oligomer and reconstituting a functional plasmid. This process is carried out recursively to assemble an oligomeric gene with the desired number of repeats. PRe-RDL has several unique features that stem from the use of type IIs restriction endonucleases: first, PRe-RDL is a seamless cloning method that leaves no extraneous nucleotides at the ligation junction. Because it uses type IIs endonucleases to ligate the two halves of the plasmid, PRe-RDL also addresses the major limitation of RDL in that it abolishes any restriction on the gene sequence that can be oligomerized. The reconstitution of a functional plasmid only upon successful ligation in PRe-RDL also addresses two other limitations of RDL: the significant background from self-ligation of the vector observed in RDL, and the decreased efficiency of ligation due to nonproductive circularization of the insert. PRe-RDL can also be used to assemble genes that encode different sequences in a predetermined order to encode block copolymers or append leader and trailer peptide sequences to the oligomerized gene.

Cytokinesis proteins Tum and Pav have a nuclear role in Wnt regulation.

Wg/Wnt signals specify cell fates in both invertebrate and vertebrate embryos and maintain stem-cell populations in many adult tissues. Deregulation of the Wnt pathway can transform cells to a proliferative fate, leading to cancer. We have discovered that two Drosophila proteins that are crucial for cytokinesis have a second, largely independent, role in restricting activity of the Wnt pathway. The fly homolog of RacGAP1, Tumbleweed (Tum)/RacGAP50C, and its binding partner, the kinesin-like protein Pavarotti (Pav), negatively regulate Wnt activity in fly embryos and in cultured mammalian cells. Unlike many known regulators of the Wnt pathway, these molecules do not affect stabilization of Arm/beta-catenin (betacat), the principal effector molecule in Wnt signal transduction. Rather, they appear to act downstream of betacat stabilization to control target-gene transcription. Both Tum and Pav accumulate in the nuclei of interphase cells, a location that is spatially distinct from their cleavage-furrow localization during cytokinesis. We show that this nuclear localization is essential for their role in Wnt regulation. Thus, we have identified two modulators of the Wnt pathway that have shared functions in cell division, which hints at a possible link between cytokinesis and Wnt activity during tumorigenesis.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Risk of ovarian cancer and inherited variants in relapse-associated genes.

BACKGROUND: We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk. METHODS AND FINDINGS: Women with and without invasive ovarian cancer (749 cases, 1,041 controls) were genotyped at 136 single nucleotide polymorphisms (SNPs) within 13 candidate genes. Risk was estimated for each SNP and for overall variation within each gene. At the gene-level, variation within MSL1 (male-specific lethal-1 homolog) was associated with risk of serous cancer (p = 0.03); haplotypes within PRPF31 (PRP31 pre-mRNA processing factor 31 homolog) were associated with risk of invasive disease (p = 0.03). MSL1 rs7211770 was associated with decreased risk of serous disease (OR 0.81, 95% CI 0.66-0.98; p = 0.03). SNPs in MFSD7, BTN3A3, ZNF200, PTPRS, and CCND1A were inversely associated with risk (p<0.05), and there was increased risk at HEXIM1 rs1053578 (p = 0.04, OR 1.40, 95% CI 1.02-1.91). CONCLUSIONS: Tumor studies can reveal novel genes worthy of follow-up for cancer susceptibility. Here, we found that inherited markers in the gene encoding MSL1, part of a complex that modifies the histone H4, may decrease risk of invasive serous ovarian cancer.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.

Chromosomal organization of adrenergic receptor genes.

The adrenergic receptors (ARs) (subtypes alpha 1, alpha 2, beta 1, and beta 2) are a prototypic family of guanine nucleotide binding regulatory protein-coupled receptors that mediate the physiological effects of the hormone epinephrine and the neurotransmitter norepinephrine. We have previously assigned the genes for beta 2- and alpha 2-AR to human chromosomes 5 and 10, respectively. By Southern analysis of somatic cell hybrids and in situ chromosomal hybridization, we have now mapped the alpha 1-AR gene to chromosome 5q32----q34, the same position as beta 2-AR, and the beta 1-AR gene to chromosome 10q24----q26, the region where alpha 2-AR is located. In mouse, both alpha 2- and beta 1-AR genes were assigned to chromosome 19, and the alpha 1-AR locus was localized to chromosome 11. Pulsed field gel electrophoresis has shown that the alpha 1- and beta 2-AR genes in humans are within 300 kilobases (kb) and the distance between the alpha 2- and beta 1-AR genes is less than 225 kb. The proximity of these two pairs of AR genes and the sequence similarity that exists among all the ARs strongly suggest that they are evolutionarily related. Moreover, they likely arose from a common ancestral receptor gene and subsequently diverged through gene duplication and chromosomal duplication to perform their distinctive roles in mediating the physiological effects of catecholamines. The AR genes thus provide a paradigm for understanding the evolution of such structurally conserved yet functionally divergent families of receptor molecules.

cAMP stimulates transcription of the beta 2-adrenergic receptor gene in response to short-term agonist exposure.

In addition to conveying cellular responses to an effector molecule, receptors are often themselves regulated by their effectors. We have demonstrated that epinephrine modulates both the rate of transcription of the beta 2-adrenergic receptor (beta 2AR) gene and the steady-state level of beta 2AR mRNA in DDT1MF-2 cells. Short-term (30 min) exposure to epinephrine (100 nM) stimulates the rate of beta 2AR gene transcription, resulting in a 3- to 4-fold increase in steady-state beta 2AR mRNA levels. These effects are mimicked by 1 mM N6,O2'-dibutyryladenosine 3',5'-cyclic monophosphate (Bt2cAMP) or foskolin but not by phorbol esters. The half-life of the beta 2AR mRNA after addition of actinomycin D (46.7 +/- 10.2 min; mean +/- SEM; n = 5) remained unchanged after 30 min of epinephrine treatment (46.8 +/- 10.6 min; mean +/- SEM; n = 4), indicating that a change in transcription rate is the predominant factor responsible for the increase of beta 2AR mRNA. Whereas brief exposure to epinephrine or Bt2cAMP does not significantly affect the total number of cellular beta 2ARs (assessed by ligand binding), continued exposure results in a gradual decline in beta 2AR number to approximately 20% (epinephrine) or approximately 45% (Bt2cAMP) of the levels in control cells by 24 hr. Similar decreases in agonist-stimulated adenylyl cyclase activity are observed. This loss of receptors with prolonged agonist exposure is accompanied by a 50% reduction in beta 2AR mRNA. Transfection of the beta 2AR promoter region cloned onto a reporter gene (bacterial chloramphenicol acetyltransferase) allowed demonstration of a 2- to 4-fold induction of transcription by agents that elevate cAMP levels, such as forskolin or phosphodiesterase inhibitors. These results establish the presence of elements within the proximal promoter region of the beta 2AR gene responsible for the transcriptional enhancing activity of cAMP and demonstrate that beta 2AR gene expression is regulated by a type of feedback mechanism involving the second messenger cAMP.

Transcription factors MYOCD, SRF, Mesp1 and SMARCD3 enhance the cardio-inducing effect of GATA4, TBX5, and MEF2C during direct cellular reprogramming.

Transient overexpression of defined combinations of master regulator genes can effectively induce cellular reprogramming: the acquisition of an alternative predicted phenotype from a differentiated cell lineage. This can be of particular importance in cardiac regenerative medicine wherein the heart lacks the capacity to heal itself, but simultaneously contains a large pool of fibroblasts. In this study we determined the cardio-inducing capacity of ten transcription factors to actuate cellular reprogramming of mouse embryonic fibroblasts into cardiomyocyte-like cells. Overexpression of transcription factors MYOCD and SRF alone or in conjunction with Mesp1 and SMARCD3 enhanced the basal but necessary cardio-inducing effect of the previously reported GATA4, TBX5, and MEF2C. In particular, combinations of five or seven transcription factors enhanced the activation of cardiac reporter vectors, and induced an upregulation of cardiac-specific genes. Global gene expression analysis also demonstrated a significantly greater cardio-inducing effect when the transcription factors MYOCD and SRF were used. Detection of cross-striated cells was highly dependent on the cell culture conditions and was enhanced by the addition of valproic acid and JAK inhibitor. Although we detected Ca(2+) transient oscillations in the reprogrammed cells, we did not detect significant changes in resting membrane potential or spontaneously contracting cells. This study further elucidates the cardio-inducing effect of the transcriptional networks involved in cardiac cellular reprogramming, contributing to the ongoing rational design of a robust protocol required for cardiac regenerative therapies.

Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer.

BACKGROUND: We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study. METHODS/PRINCIPAL FINDINGS: The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)(per allele) = 0.66; 95% credible interval (CI) = 0.44-1.00) and rs6005835 (median OR(per allele) = 0.69; 95% CI = 0.53-0.91) in CHEK2, rs2078486 (median OR(per allele) = 1.65; 95% CI = 1.21-2.25) and rs12951053 (median OR(per allele) = 1.65; 95% CI = 1.20-2.26) in TP53, rs411697 (median OR (rare homozygote) = 0.53; 95% CI = 0.35 - 0.79) in BACH1 and rs10131 (median OR( rare homozygote) = not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study. CONCLUSIONS/SIGNIFICANCE: Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.