Her2 amplification, Rel-A, and Bach1 can influence APOBEC3A expression in breast cancer cells.

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells.

Genetic inhibitors of APOBEC3B-induced mutagenesis.

The cytidine deaminases APOBEC3A (A3A) and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we used a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. We also determined whether each deletion was epistatic with Ung1 loss, which indicated whether the encoded factors participate in the homologous recombination (HR)-dependent bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single-stranded DNA (ssDNA). We found that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics, we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, (2) defective CTF18-RFC complex function, and (3) defective HR-mediated bypass of APOBEC-induced lesions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancers display three- to fourfold more APOBEC-induced mutations. Mirroring our results in yeast, Rev1-mediated C-to-G substitutions are mainly responsible for increased APOBEC-signature mutations in BRCA1/2-deficient tumors, and these mutations associate with lagging strand synthesis during replication. These results identify important factors that influence DNA replication dynamics and likely the abundance of APOBEC-induced mutation during tumor progression. They also highlight a novel role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

Genetic modifiers of APOBEC-induced mutagenesis.

The cytidine deaminases APOBEC3A and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we utilized a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. Also, we determined whether each deletion was epistatic with UNG1 loss, which indicated whether the encoded factors participate in the error-free bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single stranded DNA (ssDNA). Additionally, we determined that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, which results in extremely high A3B-induced mutagenesis, (2) defective CTF18-RFC complex function, which results in high levels of A3B induced mutations specifically on the leading strand template that synergistically increase with loss of UNG1, and (3) defective HR-mediated bypass of APOBEC-induced lesions, which were epistatic with Ung1 loss and result from increased Rev1-mediated C-to-G substitutions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancer tumors display 3- to 4-fold more APOBEC-induced mutations. Mirroring our results in yeast, for BRCA1/2 deficient tumors Rev1-mediated C-to-G substitutions are solely responsible for increased APOBEC-signature mutations and these mutations occur on the lagging strand during DNA replication. Together these results identify important factors that influence the dynamics of DNA replication and likely the abundance of APOBEC-induced mutation during tumor progression as well as a novel mechanistic role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

Transcriptional Changes following Cellular Knockdown of the Schizophrenia Risk Gene SETD1A Are Enriched for Common Variant Association with the Disorder.

Loss of function mutations in SETD1A are the first experiment-wide significant findings to emerge from exome sequencing studies of schizophrenia. Although SETD1A is known to encode a histone methyltransferase, the consequences of reduced S ETD1A activity on gene expression in neural cells have, to date, been unknown. To explore transcriptional changes through which genetic perturbation of SETD1A could confer risk for schizophrenia, we have performed genome-wide gene expression profiling of a commonly used human neuroblastoma cell line in which SETD1A expression has been experimentally reduced using RNA interference (RNAi). We identified 1,031 gene expression changes that were significant in two separate RNAi conditions compared with control, including effects on genes of known neurodevelopmental importance such as DCX and DLX5. Genes that were differentially expressed following SETD1A knockdown were enriched for annotation to metabolic pathways, peptidase regulator activity and integrin-mediated regulation of cell adhesion. Moreover, differentially expressed genes were enriched for common variant association with schizophrenia, suggesting a degree of molecular convergence between this rare schizophrenia risk factor and susceptibility variants for the disorder operating more generally.

The Psychiatric Risk Gene NT5C2 Regulates Adenosine Monophosphate-Activated Protein Kinase Signaling and Protein Translation in Human Neural Progenitor Cells.

BACKGROUND: The 5'-nucleotidase, cytosolic II gene (NT5C2, cN-II) is associated with disorders characterized by psychiatric and psychomotor disturbances. Common psychiatric risk alleles at the NT5C2 locus reduce expression of this gene in the fetal and adult brain, but downstream biological risk mechanisms remain elusive. METHODS: Distribution of the NT5C2 protein in the human dorsolateral prefrontal cortex and cortical human neural progenitor cells (hNPCs) was determined using immunostaining, publicly available expression data, and reverse transcriptase quantitative polymerase chain reaction. Phosphorylation quantification of adenosine monophosphate-activated protein kinase (AMPK) alpha (Thr172) and ribosomal protein S6 (Ser235/Ser236) was performed using Western blotting to infer the degree of activation of AMPK signaling and the rate of protein translation. Knockdowns were induced in hNPCs and Drosophila melanogaster using RNA interference. Transcriptomic profiling of hNPCs was performed using microarrays, and motility behavior was assessed in flies using the climbing assay. RESULTS: Expression of NT5C2 was higher during neurodevelopment and was neuronally enriched in the adult human cortex. Knockdown in hNPCs affected AMPK signaling, a major nutrient-sensing mechanism involved in energy homeostasis, and protein translation. Transcriptional changes implicated in protein translation were observed in knockdown hNPCs, and expression changes to genes related to AMPK signaling and protein translation were confirmed using reverse transcriptase quantitative polymerase chain reaction. The knockdown in Drosophila was associated with drastic climbing impairment. CONCLUSIONS: We provide an extensive neurobiological characterization of the psychiatric risk gene NT5C2, describing its previously unknown role in the regulation of AMPK signaling and protein translation in neural stem cells and its association with Drosophila melanogaster motility behavior.

CRISPR-Cas9 genome editing in human cells occurs via the Fanconi anemia pathway.

CRISPR-Cas genome editing creates targeted DNA double-strand breaks (DSBs) that are processed by cellular repair pathways, including the incorporation of exogenous DNA via single-strand template repair (SSTR). To determine the genetic basis of SSTR in human cells, we developed a coupled inhibition-cutting system capable of interrogating multiple editing outcomes in the context of thousands of individual gene knockdowns. We found that human Cas9-induced SSTR requires the Fanconi anemia (FA) pathway, which is normally implicated in interstrand cross-link repair. The FA pathway does not directly impact error-prone, non-homologous end joining, but instead diverts repair toward SSTR. Furthermore, FANCD2 protein localizes to Cas9-induced DSBs, indicating a direct role in regulating genome editing. Since FA is itself a genetic disease, these data imply that patient genotype and/or transcriptome may impact the effectiveness of gene editing treatments and that treatments biased toward FA repair pathways could have therapeutic value.

Pleiotropic Effects of Trait-Associated Genetic Variation on DNA Methylation: Utility for Refining GWAS Loci.

Most genetic variants identified in genome-wide association studies (GWASs) of complex traits are thought to act by affecting gene regulation rather than directly altering the protein product. As a consequence, the actual genes involved in disease are not necessarily the most proximal to the associated variants. By integrating data from GWAS analyses with those from genetic studies of regulatory variation, it is possible to identify variants pleiotropically associated with both a complex trait and measures of gene regulation. In this study, we used summary-data-based Mendelian randomization (SMR), a method developed to identify variants pleiotropically associated with both complex traits and gene expression, to identify variants associated with complex traits and DNA methylation. We used large DNA methylation quantitative trait locus (mQTL) datasets generated from two different tissues (blood and fetal brain) to prioritize genes for >40 complex traits with robust GWAS data and found considerable overlap with the results of SMR analyses performed with expression QTL (eQTL) data. We identified multiple examples of variable DNA methylation associated with GWAS variants for a range of complex traits, demonstrating the utility of this approach for refining genetic association signals.

Transcriptional consequences of schizophrenia candidate miR-137 manipulation in human neural progenitor cells.

MIR137, transcribed as the microRNA miR-137, is one of the leading candidate schizophrenia susceptibility genes to arise from large genome-wide association studies (GWAS) of the disorder. Recent data suggest that miR-137 modulates the expression of other schizophrenia susceptibility genes. Although bioinformatic resources are available with which to predict genes regulated by individual microRNA, there has been a lack of empirical data on genome-wide gene expression changes following miR-137 manipulation. We have therefore performed a genome-wide assessment of transcriptional changes in a human neural progenitor cell line after miR-137 over-expression and inhibition in order to elucidate molecular pathways by which genetic perturbation of miR-137 could promote susceptibility to schizophrenia. Bioinformatically-predicted miR-137 targets showed a small but highly significant down-regulation following miR-137 over-expression. Genes that were significantly down-regulated in association with miR-137 over-expression were enriched for involvement in neuronal differentiation. Differentially expressed genes that were confirmed by qPCR included others at genome-wide significant risk loci for schizophrenia (MAD1L1 and DPYD) and BDNF. These data point to molecular pathways through which genetic variation at the MIR137 locus could confer risk for schizophrenia.

Convergent evidence that oligodendrocyte lineage transcription factor 2 (OLIG2) and interacting genes influence susceptibility to schizophrenia.

Abnormal oligodendrocyte function has been postulated as a primary etiological event in schizophrenia. Oligodendrocyte lineage transcription factor 2 (OLIG2) encodes a transcription factor central to oligodendrocyte development. Analysis of OLIG2 in a case-control sample (n = approximately 1,400) in the U.K. revealed several SNPs to be associated with schizophrenia (minimum P = 0.0001, gene-wide P = 0.0009). To obtain independent support for this association, we sought evidence for genetic interaction between OLIG2 and three genes of relevance to oligodendrocyte function for which we have reported evidence for association with schizophrenia: CNP, NRG1, and ERBB4. We found interaction effects on disease risk between OLIG2 and CNP (minimum P = 0.0001, corrected P = 0.008) for interaction with ERBB4 (minimum P = 0.002, corrected P = 0.04) but no evidence for interaction with NRG1. To investigate the biological plausibility of the interactions, we sought correlations between the expression of the genes. The results were similar to those of the genetic interaction analysis. OLIG2 expression significantly correlated in cerebral cortex with CNP (P < 10(-7)) and ERBB4 (P = 0.002, corrected P = 0.038) but not NRG1. In mouse striatum, Olig2 and Cnp expression also was correlated, and linkage analysis for trans-effects on gene expression suggests that each locus regulates the other's expression. Our data provide strong convergent evidence that variation in OLIG2 confers susceptibility to schizophrenia alone and as part of a network of genes implicated in oligodendrocyte function.

Evidence that interaction between neuregulin 1 and its receptor erbB4 increases susceptibility to schizophrenia.

There is now strong evidence that Neuregulin 1 (NRG1) is a susceptibility gene for schizophrenia. NRG1 mediates some of its effects through the tyrosine kinase receptor erbB4, and analysis of gene knock-out animals suggests that the functional interaction of NRG1 and erbB4 mediates behaviors that may model some aspects of the schizophrenia phenotype in mice. Given these findings, we have sought evidence for association between schizophrenia and erbB4. Mutation screening of erbB4 in 14 DSMIV schizophrenics revealed 15 SNPs, none of which were nonsynonymous. Analysis of the allele frequencies of each SNP in pools of 368 DSMIV schizophrenics and 368 controls provided modest evidence for association with two of the SNPs, although individual genotyping in an extended sample of 680 cases did not confirm this. However, we did find evidence for a significant interaction between the NRG1 "Icelandic" schizophrenia risk haplotype and erbB4 (P = 0.019). The NRG1 and erbB4 interacting marker was further genotyped in an independent sample of 290 cases and 634 controls from Dublin. Interaction between NRG1 and erbB4 remained significant in the combined sample of 970 cases and 1,341 controls, OR = 2.98 (CI: 1.16-7.64), P = 0.01, although it only showed a trend in the Dublin sample alone (P = 0.11, two tailed). Our data require independent replication, but tentatively suggest that NRG1 may mediate its effects on schizophrenia susceptibility through functional interaction with erbB4, and that genetic interaction between variants at the two loci increases susceptibility to schizophrenia.