Systems biology approaches that are based on the genetics of gene expression have been fruitful in identifying genetic regulatory loci related to complex traits. We use microarray and genetic marker data from an F2 mouse intercross to examine the large-scale organization of the gene co-expression network in liver, and annotate several gene modules in terms of 22 physiological traits. We identify chromosomal loci (referred to as module quantitative trait loci, mQTL) that perturb the modules and describe a novel approach that integrates network properties with genetic marker information to model gene/trait relationships. Specifically, using the mQTL and the intramodular connectivity of a body weight–related module, we describe which factors determine the relationship between gene expression profiles and weight. Our approach results in the identification of genetic targets that influence gene modules (pathways) that are related to the clinical phenotypes of interest.

Targeted mutagenesis is an essential tool of reverse genetics that could be used experimentally to investigate basic plant biology or modify crop plants for improvement of important agricultural traits. Although targeted mutagenesis is routine in several model organisms including yeast and mouse, efficient and widely usable methods to generate targeted modifications in plant genes are not currently available. In this study we investigated the efficacy of a targeted-mutagenesis approach based on zinc-finger nucleases (ZFNs). In this procedure, ZFNs are used to generate double-strand breaks at specific genomic sites, and subsequent repair produces mutations at the break site. To determine whether ZFNs can cleave and induce mutations at specific sites within higher plant genomes, we introduced a construct carrying both a ZFN gene, driven by a heat-shock promoter, and its target into the Arabidopsis genome. Induction of ZFN expression by heat shock during seedling development resulted in mutations at the ZFN recognition sequence at frequencies as high as 0.2 mutations per target. Of 106 ZFN-induced mutations characterized, 83 (78%) were simple deletions of 1-52 bp (median of 4 bp), 14 (13%) were simple insertions of 1-4 bp, and 9 (8%) were deletions accompanied by insertions. In 10% of induced individuals, mutants were present in the subsequent generation, thus demonstrating efficient transmission of the ZFNinduced mutations. These data indicate that ZFNs can form the basis of a highly efficient method for targeted mutagenesis of plant genes.gene targeting ͉ nonhomologous end joining A major focus of plant biotechnology is genetic modification and improvement of crop plants. With this aim, large-scale genome and͞or EST sequencing projects are underway for many important plant species including rice, maize, wheat, soybean, and tomato (www.ncbi.nlm.nih.gov͞genomes͞PLANTS͞ PlantList.html). The enormous amount of genome-sequence information becoming available has intensified the need for methods that can use this sequence information to generate targeted modifications in plant genes. Targeted-mutagenesis methods could be used experimentally to investigate plant gene function or for genetic modification of important crop plants. Such methods are especially important for species, including most crops, that lack readily available mutant collections (1, 2). Furthermore, targeted mutagenesis could facilitate development of genetically modified crops lacking transgenic DNA, including genes conferring resistance to antibiotics. To date, efficient and widely usable methods for targeted modification of higher plant genomes are not available.The most widely used targeted-mutagenesis strategy is gene targeting (GT) by homologous recombination (3-6). Efficient GT procedures have been available for Ͼ20 years in yeast (7) and mouse (8). In these systems, DNA fragments are introduced into cultured cells, and repair by homologous recombination causes the introduced DNA to be incorporated into the homologous locus. Typically, GT events ...

SUMMARY To identify therapeutic targets for Glioblastoma (GBM), we performed genome-wide CRISPR-Cas9 "knockout" (KO) screens in patient-derived GBM stem-like cells (GSCs) and human neural stem/progenitors (NSCs), non-neoplastic stem cell controls, for genes required for their in vitro growth. Surprisingly, the vast majority GSC-lethal hits were found outside of molecular networks commonly altered in GBM and GSCs (e.g., oncogenic drivers). In vitro and in vivo validation of GSC-specific targets revealed several strong hits, including the wee1-like kinase, PKMYT1/Myt1. Mechanistic studies demonstrated that PKMYT1 acts redundantly with WEE1 to inhibit Cyclin B-CDK1 activity via CDK1-Y15 phosphorylation and to promote timely completion of mitosis in NSCs. However, in GSCs, this redundancy is lost, likely as a result of oncogenic signaling, causing GBM-specific lethality.

We hypothesized that a common SNP in the 3' untranslated region of the upstream transcription factor 1 (USF1), rs3737787, may affect lipid traits by influencing gene expression levels, and we investigated this possibility utilizing the Mexican population, which has a high predisposition to dyslipidemia. We first associated rs3737787 genotypes in Mexican Familial Combined Hyperlipidemia (FCHL) case/control fat biopsies, with global expression patterns. To identify sets of co-expressed genes co-regulated by similar factors such as transcription factors, genetic variants, or environmental effects, we utilized weighted gene co-expression network analysis (WGCNA). Through WGCNA in the Mexican FCHL fat biopsies we identified two significant Triglyceride (TG)-associated co-expression modules. One of these modules was also associated with FCHL, the other FCHL component traits, and rs3737787 genotypes. This USF1-regulated FCHL-associated (URFA) module was enriched for genes involved in lipid metabolic processes. Using systems genetics procedures we identified 18 causal candidate genes in the URFA module. The FCHL causal candidate gene fatty acid desaturase 3 (FADS3) was associated with TGs in a recent Caucasian genome-wide significant association study and we replicated this association in Mexican FCHL families. Based on a USF1-regulated FCHL-associated co-expression module and SNP rs3737787, we identify a set of causal candidate genes for FCHL-related traits. We then provide evidence from two independent datasets supporting FADS3 as a causal gene for FCHL and elevated TGs in Mexicans.

Genes regulated by the same miRNA can be discovered by virtue of their coexpression at the transcriptional level and the presence of a conserved miRNA-binding site in their 39 UTRs. Using this principle we have integrated the three best performing and complementary algorithms into a framework for inference of regulation by miRNAs (FIRM) from sets of coexpressed genes. We demonstrate the utility of FIRM by inferring a cancer-miRNA regulatory network through the analysis of 2240 gene coexpression signatures from 46 cancers. By analyzing this network for functional enrichment of known hallmarks of cancer we have discovered a subset of 13 miRNAs that regulate oncogenic processes across diverse cancers. We have performed experiments to test predictions from this miRNA-regulatory network to demonstrate that miRNAs of the miR-29 family (miR-29a, miR-29b, and miR-29c) regulate specific genes associated with tissue invasion and metastasis in lung adenocarcinoma. Further, we highlight the specificity of using FIRM inferences to identify miRNAregulated genes by experimentally validating that miR-767-5p, which partially shares the miR-29 seed sequence, regulates only a subset of miR-29 targets. By providing mechanistic linkage between miRNA dysregulation in cancer, their binding sites in the 39UTRs of specific sets of coexpressed genes, and their associations with known hallmarks of cancer, FIRM, and the inferred cancer miRNA-regulatory network will serve as a powerful public resource for discovery of potential cancer biomarkers.

Gene by environment (GxE) interactions are clearly important in many human diseases, but they have proven to be difficult to study on a molecular level. We report genetic analysis of thousands of transcript abundance traits in human primary endothelial cell (EC) lines in response to proinflammatory oxidized phospholipids implicated in cardiovascular disease. Of the 59 most regulated transcripts, approximately one-third showed evidence of GxE interactions. The interactions resulted primarily from effects of distal-, trans-acting loci, but a striking example of a local-GxE interaction was also observed for FGD6. Some of the distal interactions were validated by siRNA knockdown experiments, including a locus involved in the regulation of multiple transcripts involved in the ER stress pathway. Our findings add to the understanding of the overall architecture of complex human traits and are consistent with the possibility that GxE interactions are responsible, in part, for the failure of association studies to more fully explain common disease variation.

Systems biology approaches that are based on the genetics of gene expression have been fruitful in identifying genetic regulatory loci related to complex traits. We use microarray and genetic marker data from an F2 mouse intercross to examine the large-scale organization of the gene co-expression network in liver, and annotate several gene modules in terms of 22 physiological traits. We identify chromosomal loci (referred to as module quantitative trait loci, mQTL) that perturb the modules and describe a novel approach that integrates network properties with genetic marker information to model gene/trait relationships. Specifically, using the mQTL and the intramodular connectivity of a body weight-related module, we describe which factors determine the relationship between gene expression profiles and weight. Our approach results in the identification of genetic targets that influence gene modules (pathways) that are related to the clinical phenotypes of interest.

The regulation of host-pathogen interactions during Mycobacterium tuberculosis (Mtb) infection remains unresolved. MicroRNAs (miRNAs) are important regulators of the immune system, and so we used a systems biology approach to construct an miRNA regulatory network activated in macrophages during Mtb infection. Our network comprises 77 putative miRNAs that are associated with temporal gene expression signatures in macrophages early after Mtb infection. In this study, we demonstrate a dual role for one of these regulators, miR-155. On the one hand, miR-155 maintains the survival of Mtbinfected macrophages, thereby providing a niche favoring bacterial replication; on the other hand, miR-155 promotes the survival and function of Mtb-specific T cells, enabling an effective adaptive immune response. MiR-155-induced cell survival is mediated through the SH2 domain-containing inositol 5-phosphatase 1 (SHIP1)/protein kinase B (Akt) pathway. Thus, dual regulation of the same cell survival pathway in innate and adaptive immune cells leads to vastly different outcomes with respect to bacterial containment.