A Type II-B Cas9 nuclease with minimized off-targets and reduced chromosomal translocations in vivo.

Streptococcus pyogenes Cas9 (SpCas9) and derived enzymes are widely used as genome editors, but their promiscuous nuclease activity often induces undesired mutations and chromosomal rearrangements. Several strategies for mapping off-target effects have emerged, but they suffer from limited sensitivity. To increase the detection sensitivity, we develop an off-target assessment workflow that uses Duplex Sequencing. The strategy increases sensitivity by one order of magnitude, identifying previously unknown SpCas9's off-target mutations in the humanized PCSK9 mouse model. To reduce off-target risks, we perform a bioinformatic search and identify a high-fidelity Cas9 variant of the II-B subfamily from Parasutterella secunda (PsCas9). PsCas9 shows improved specificity as compared to SpCas9 across multiple tested sites, both in vitro and in vivo, including the PCSK9 site. In the future, while PsCas9 will offer an alternative to SpCas9 for research and clinical use, the Duplex Sequencing workflow will enable a more sensitive assessment of Cas9 editing outcomes.

APOL1 promotes endothelial cell activation beyond the glomerulus.

Apolipoprotein L1 (APOL1) high-risk genotypes are associated with increased risk of chronic kidney disease (CKD) in people of West African ancestry. Given the importance of endothelial cells (ECs) in CKD, we hypothesized that APOL1 high-risk genotypes may contribute to disease via EC-intrinsic activation and dysfunction. Single cell RNA sequencing (scRNA-seq) analysis of the Kidney Precision Medicine Project dataset revealed APOL1 expression in ECs from various renal vascular compartments. Utilizing two public transcriptomic datasets of kidney tissue from African Americans with CKD and a dataset of APOL1-expressing transgenic mice, we identified an EC activation signature; specifically, increased intercellular adhesion molecule 1 (ICAM-1) expression and enrichment in leukocyte migration pathways. In vitro, APOL1 expression in ECs derived from genetically modified human induced pluripotent stem cells and glomerular ECs triggered changes in ICAM-1 and platelet endothelial cell adhesion molecule 1 (PECAM-1) leading to an increase in monocyte attachment. Overall, our data suggest the involvement of APOL1 as an inducer of EC activation in multiple renal vascular beds with potential effects beyond the glomerular vasculature.

Hepatic patatin-like phospholipase domain-containing 3 levels are increased in I148M risk allele carriers and correlate with NAFLD in humans.

In nonalcoholic fatty liver disease (NAFLD) the patatin-like phospholipase domain-containing 3 (PNPLA3) rs738409 variant is a contributor. In mice, the Pnpla3 148M variant accumulates on lipid droplets and probably leads to sequestration of a lipase cofactor leading to impaired mobilization of triglycerides. To advance our understanding of the localization and abundance of PNPLA3 protein in humans, we used liver biopsies from patients with NAFLD to investigate the link to NAFLD and the PNPLA3 148M genotype. We experimentally qualified an antibody against human PNPLA3. Hepatic PNPLA3 protein fractional area and localization were determined by immunohistochemistry in biopsies from a well-characterized NAFLD cohort of 67 patients. Potential differences in hepatic PNPLA3 protein levels among patients related to degree of steatosis, lobular inflammation, ballooning, and fibrosis, and PNPLA3 I148M gene variants were assessed. Immunohistochemistry staining in biopsies from patients with NAFLD showed that hepatic PNPLA3 protein was predominantly localized to the membranes of small and large lipid droplets in hepatocytes. PNPLA3 protein levels correlated strongly with steatosis grade (p = 0.000027) and were also significantly higher in patients with lobular inflammation (p = 0.009), ballooning (p = 0.022), and significant fibrosis (stage 2-4, p = 0.014). In addition, PNPLA3 levels were higher in PNPLA3 rs738409 148M (CG, GG) risk allele carriers compared to 148I (CC) nonrisk allele carriers (p = 0.0029). Conclusion: PNPLA3 protein levels were associated with increased hepatic lipid content and disease severity in patients with NAFLD and were higher in PNPLA3 rs738409 (148M) risk allele carriers. Our hypothesis that increased hepatic levels of PNPLA3 may be part of the pathophysiological mechanism of NAFLD is supported.

The long noncoding RNA TUNAR modulates Wnt signaling and regulates human ß-cell proliferation.

Many long noncoding RNAs (lncRNAs) are enriched in pancreatic islets and several lncRNAs are linked to type 2 diabetes (T2D). Although they have emerged as potential players in ß-cell biology and T2D, little is known about their functions and mechanisms in human ß-cells. We identified an islet-enriched lncRNA, TUNAR (TCL1 upstream neural differentiation-associated RNA), which was upregulated in ß-cells of patients with T2D and promoted human ß-cell proliferation via fine-tuning of the Wnt pathway. TUNAR was upregulated following Wnt agonism by a glycogen synthase kinase-3 (GSK3) inhibitor in human ß-cells. Reciprocally, TUNAR repressed a Wnt antagonist Dickkopf-related protein 3 (DKK3) and stimulated Wnt pathway signaling. DKK3 was aberrantly expressed in ß-cells of patients with T2D and displayed a synchronized regulatory pattern with TUNAR at the single cell level. Mechanistically, DKK3 expression was suppressed by the repressive histone modifier enhancer of zeste homolog 2 (EZH2). TUNAR interacted with EZH2 in ß-cells and facilitated EZH2-mediated suppression of DKK3. These findings reveal a novel cell-specific epigenetic mechanism via islet-enriched lncRNA that fine-tunes the Wnt pathway and subsequently human ß-cell proliferation.NEW & NOTEWORTHY The discovery that long noncoding RNA TUNAR regulates ß-cell proliferation may be important in designing new treatments for diabetes.

A Phenotypic Screening Approach Using Human Treg Cells Identified Regulators of Forkhead Box p3 Expression.

Regulatory T (Treg) cells, expressing the transcription factor forkhead box p3 (FOXP3), are the key cells regulating peripheral autoreactive T lymphocytes by suppressing effector T cells. FOXP3+ Treg cells play essential roles controlling immune responses in autoimmune diseases and cancer. Several clinical approaches (e.g., polyclonal expansion of Treg cells with anti-CD3 and anti-CD28 coated beads in the presence of drugs) are under evaluation. However, expression of FOXP3, recognized as the master regulator of Treg cells, in induced Treg cells have been shown to be instable, and molecular targets involved in regulating FOXP3 expression and Treg cell function have not been well-defined. Thus, new targets directly regulating FOXP3 expression and the expression of its downstream genes (e.g., cytotoxic T-lymphocyte-associated protein 4 (CTLA4)) have the potential to stabilize the Treg cell phenotype and function. This report describes the development of an automated medium-throughput 384-well plate flow cytometry phenotypic assay meauring the protein expression of FOXP3 and CTLA4 in human Treg cells. Screening a library of 4213 structurally diverse compounds allowed us to identify a variety of compounds regulating FOXP3 and CTLA4 expression. Further evaluation of these and related small molecules, followed by confirmation using siRNA-mediated gene knockdown, revealed three targets: euchromatic histone-lysine N-methyltransferase (EHMT2) and glycogen synthase kinase 3 alpha/beta (GSK3α/ß) as potent positive regulators of FOXP3 expression, and bromodomain and extra-terminal domain (BET) inhibitors as negative regulators of FOXP3 and CTLA4 expression. These targets have potential implications for establishing novel therapies for autoimmune diseases and cancer.

In vivo genome and base editing of a human PCSK9 knock-in hypercholesterolemic mouse model.

BACKGROUND: Plasma concentration of low-density lipoprotein (LDL) cholesterol is a well-established risk factor for cardiovascular disease. Inhibition of proprotein convertase subtilisin/kexin type 9 (PCSK9), which regulates cholesterol homeostasis, has recently emerged as an approach to reduce cholesterol levels. The development of humanized animal models is an important step to validate and study human drug targets, and use of genome and base editing has been proposed as a mean to target disease alleles. RESULTS: To address the lack of validated models to test the safety and efficacy of techniques to target human PCSK9, we generated a liver-specific human PCSK9 knock-in mouse model (hPCSK9-KI). We showed that plasma concentrations of total cholesterol were higher in hPCSK9-KI than in wildtype mice and increased with age. Treatment with evolocumab, a monoclonal antibody that targets human PCSK9, reduced cholesterol levels in hPCSK9-KI but not in wildtype mice, showing that the hypercholesterolemic phenotype was driven by overexpression of human PCSK9. CRISPR-Cas9-mediated genome editing of human PCSK9 reduced plasma levels of human and not mouse PCSK9, and in parallel reduced plasma concentrations of total cholesterol; genome editing of mouse Pcsk9 did not reduce cholesterol levels. Base editing using a guide RNA that targeted human and mouse PCSK9 reduced plasma levels of human and mouse PCSK9 and total cholesterol. In our mouse model, base editing was more precise than genome editing, and no off-target editing nor chromosomal translocations were identified. CONCLUSIONS: Here, we describe a humanized mouse model with liver-specific expression of human PCSK9 and a human-like hypercholesterolemia phenotype, and demonstrate that this mouse can be used to evaluate antibody and gene editing-based (genome and base editing) therapies to modulate the expression of human PCSK9 and reduce cholesterol levels. We predict that this mouse model will be used in the future to understand the efficacy and safety of novel therapeutic approaches for hypercholesterolemia.

Therapeutic Genome Editing With CRISPR/Cas9 in a Humanized Mouse Model Ameliorates α1-antitrypsin Deficiency Phenotype.

α1-antitrypsin (AAT) is a circulating serine protease inhibitor secreted from the liver and important in preventing proteolytic neutrophil elastase associated tissue damage, primarily in lungs. In humans, AAT is encoded by the SERPINA1 (hSERPINA1) gene in which a point mutation (commonly referred to as PiZ) causes aggregation of the miss-folded protein in hepatocytes resulting in subsequent liver damage. In an attempt to rescue the pathologic liver phenotype of a mouse model of human AAT deficiency (AATD), we used adenovirus to deliver Cas9 and a guide-RNA (gRNA) molecule targeting hSERPINA1. Our single dose therapeutic gene editing approach completely reverted the phenotype associated with the PiZ mutation, including circulating transaminase and human AAT (hAAT) protein levels, liver fibrosis and protein aggregation. Furthermore, liver histology was significantly improved regarding inflammation and overall morphology in hSERPINA1 gene edited PiZ mice. Genomic analysis confirmed significant disruption to the hSERPINA1 transgene resulting in a reduction of hAAT protein levels and quantitative mRNA analysis showed a reduction in fibrosis and hepatocyte proliferation as a result of editing. Our findings indicate that therapeutic gene editing in hepatocytes is possible in an AATD mouse model.

Selective Inhibition of FOXO1 Activator/Repressor Balance Modulates Hepatic Glucose Handling.

Insulin resistance is a hallmark of diabetes and an unmet clinical need. Insulin inhibits hepatic glucose production and promotes lipogenesis by suppressing FOXO1-dependent activation of G6pase and inhibition of glucokinase, respectively. The tight coupling of these events poses a dual conundrum: mechanistically, as the FOXO1 corepressor of glucokinase is unknown, and clinically, as inhibition of glucose production is predicted to increase lipogenesis. Here, we report that SIN3A is the insulin-sensitive FOXO1 corepressor of glucokinase. Genetic ablation of SIN3A abolishes nutrient regulation of glucokinase without affecting other FOXO1 target genes and lowers glycemia without concurrent steatosis. To extend this work, we executed a small-molecule screen and discovered selective inhibitors of FOXO-dependent glucose production devoid of lipogenic activity in hepatocytes. In addition to identifying a novel mode of insulin action, these data raise the possibility of developing selective modulators of unliganded transcription factors to dial out adverse effects of insulin sensitizers.

Phenotypic Screen for Cardiac Regeneration Identifies Molecules with Differential Activity in Human Epicardium-Derived Cells versus Cardiac Fibroblasts.

Activation and proliferation of resident cardiac progenitor cells has therapeutic potential to repair the heart after injury. However, research has been impeded by a lack of well-defined and characterized cell sources and difficulties in translation to screening platforms. Here, we describe the development, validation, and use of a 384-well phenotypic assay in primary human epicardium-derived cells (EPDCs) to identify compounds that induce proliferation while maintaining the progenitor phenotype. Using this assay, we screened 7400 structurally diverse compounds where greater than 90% are biologically annotated and known to modulate a broad range of biological targets. From the primary screen, we identified and validated hits and expanded upon the lead molecules of interest. A counterscreen was developed in human cardiac fibroblasts to filter out compounds with a general proliferative effect, after which the activity of selected molecules was confirmed across multiple EPDC donors. To further examine the mechanism of action of compounds with annotated targets, we performed knockdown experiments to understand whether a single known target was responsible for the proliferative effect, confirming results with protein expression and activity assays. Here, we were able to show that the annotated targets of compounds of interest were not responsible for the proliferative effect, which highlights potential differences in cell types and signaling pathways and possible polypharmacology. These studies demonstrate the feasibility of using relevant human primary cells in a phenotypic screen to identify compounds as novel biological tools and starting points for drug discovery projects, and we disclose the first small molecules to proliferate human primary EPDCs.

Genome-wide RNAi analysis reveals that simultaneous inhibition of specific mevalonate pathway genes potentiates tumor cell death.

The mevalonate (MVA) pathway is often dysregulated or overexpressed in many cancers suggesting tumor dependency on this classic metabolic pathway. Statins, which target the rate-limiting enzyme of this pathway, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), are promising agents currently being evaluated in clinical trials for anti-cancer efficacy. To uncover novel targets that potentiate statin-induced apoptosis when knocked down, we carried out a pooled genome-wide short hairpin RNA (shRNA) screen. Genes of the MVA pathway were amongst the top-scoring targets, including sterol regulatory element binding transcription factor 2 (SREBP2), 3-hydroxy-3-methylglutaryl-coenzyme A synthase 1 (HMGCS1) and geranylgeranyl diphosphate synthase 1 (GGPS1). Each gene was independently validated and shown to significantly sensitize A549 cells to statin-induced apoptosis when knocked down. SREBP2 knockdown in lung and breast cancer cells completely abrogated the fluvastatin-induced upregulation of sterol-responsive genes HMGCR and HMGCS1. Knockdown of SREBP2 alone did not affect three-dimensional growth of lung and breast cancer cells, yet in combination with fluvastatin cell growth was disrupted. Taken together, these results show that directly targeting multiple levels of the MVA pathway, including blocking the sterol-feedback loop initiated by statin treatment, is an effective and targetable anti-tumor strategy.

A novel role of the checkpoint kinase ATR in leptin signaling.

In a world with increasing incidences of obesity, it becomes critical to understand the detailed regulation of appetite. To identify novel regulators of the signaling mediated by one of the key hormones of energy homeostasis, leptin, we screened a set of compounds for their effect on the downstream Signal Transducer and Activator of Transcription 3 (STAT3) signaling. Interestingly, cells exposed to inhibitors of the Ataxia Telangiectasia and RAD3-related protein ATR increased their leptin dependent STAT3 activity. This was due to failure of the cells to induce the negative feedback mediator Suppressor of Cytokine Signaling 3 (SOCS3), suggesting that ATR has a previously unknown role in the negative feedback regulation of leptin signaling. This is an important finding not only because it sheds light on additional genes involved in leptin signaling, but also because it brings forward a new potential therapeutic intervention point for increasing leptin signaling in obese individuals.

A comprehensive platform for highly multiplexed mammalian functional genetic screens.

BACKGROUND: Genome-wide screening in human and mouse cells using RNA interference and open reading frame over-expression libraries is rapidly becoming a viable experimental approach for many research labs. There are a variety of gene expression modulation libraries commercially available, however, detailed and validated protocols as well as the reagents necessary for deconvolving genome-scale gene screens using these libraries are lacking. As a solution, we designed a comprehensive platform for highly multiplexed functional genetic screens in human, mouse and yeast cells using popular, commercially available gene modulation libraries. The Gene Modulation Array Platform (GMAP) is a single microarray-based detection solution for deconvolution of loss and gain-of-function pooled screens. RESULTS: Experiments with specially constructed lentiviral-based plasmid pools containing ~78,000 shRNAs demonstrated that the GMAP is capable of deconvolving genome-wide shRNA "dropout" screens. Further experiments with a larger, ~90,000 shRNA pool demonstrate that equivalent results are obtained from plasmid pools and from genomic DNA derived from lentivirus infected cells. Parallel testing of large shRNA pools using GMAP and next-generation sequencing methods revealed that the two methods provide valid and complementary approaches to deconvolution of genome-wide shRNA screens. Additional experiments demonstrated that GMAP is equivalent to similar microarray-based products when used for deconvolution of open reading frame over-expression screens. CONCLUSION: Herein, we demonstrate four major applications for the GMAP resource, including deconvolution of pooled RNAi screens in cells with at least 90,000 distinct shRNAs. We also provide detailed methodologies for pooled shRNA screen readout using GMAP and compare next-generation sequencing to GMAP (i.e. microarray) based deconvolution methods.

Exploring gene function and drug action using chemogenomic dosage assays.

In this chapter, we describe a series of genome-wide, cell-based assays that provide a solid basis for understanding drug-gene interactions, gene function, and for defining the mechanism of action of small molecules. Each of these assays takes advantage of the ability to grow complex pools competitively and to use high-density microarrays that report the results of such screens. The assays described here take advantage of alterations in gene dosage of Saccharomyces cerevisiae, and include HIP (haploinsufficiency profiling), HOP (homozygous profiling), and MSP (multicopy suppression profiling) as genetic tools to understand gene function and drug mechanism. The common experimental theme is that, in each assay, strains are pooled and screened in parallel to investigate the relative contribution of each gene product to sensitivity or resistance to a drug or environmental perturbation across the genome in a single assay. Further, the compendium of results from these screens can inform large-scale network analysis of genetic function, gene-gene interactions, and mechanism of drug action.

Systematic analysis of genome-wide fitness data in yeast reveals novel gene function and drug action.

We systematically analyzed the relationships between gene fitness profiles (co-fitness) and drug inhibition profiles (co-inhibition) from several hundred chemogenomic screens in yeast. Co-fitness predicted gene functions distinct from those derived from other assays and identified conditionally dependent protein complexes. Co-inhibitory compounds were weakly correlated by structure and therapeutic class. We developed an algorithm predicting protein targets of chemical compounds and verified its accuracy with experimental testing. Fitness data provide a novel, systems-level perspective on the cell.

Chemical-genetic profiling of imidazo[1,2-a]pyridines and -pyrimidines reveals target pathways conserved between yeast and human cells.

Small molecules have been shown to be potent and selective probes to understand cell physiology. Here, we show that imidazo[1,2-a]pyridines and imidazo[1,2-a]pyrimidines compose a class of compounds that target essential, conserved cellular processes. Using validated chemogenomic assays in Saccharomyces cerevisiae, we discovered that two closely related compounds, an imidazo[1,2-a]pyridine and -pyrimidine that differ by a single atom, have distinctly different mechanisms of action in vivo. 2-phenyl-3-nitroso-imidazo[1,2-a]pyridine was toxic to yeast strains with defects in electron transport and mitochondrial functions and caused mitochondrial fragmentation, suggesting that compound 13 acts by disrupting mitochondria. By contrast, 2-phenyl-3-nitroso-imidazo[1,2-a]pyrimidine acted as a DNA poison, causing damage to the nuclear DNA and inducing mutagenesis. We compared compound 15 to known chemotherapeutics and found resistance required intact DNA repair pathways. Thus, subtle changes in the structure of imidazo-pyridines and -pyrimidines dramatically alter both the intracellular targeting of these compounds and their effects in vivo. Of particular interest, these different modes of action were evident in experiments on human cells, suggesting that chemical-genetic profiles obtained in yeast are recapitulated in cultured cells, indicating that our observations in yeast can: (1) be leveraged to determine mechanism of action in mammalian cells and (2) suggest novel structure-activity relationships.

Off-target effects of psychoactive drugs revealed by genome-wide assays in yeast.

To better understand off-target effects of widely prescribed psychoactive drugs, we performed a comprehensive series of chemogenomic screens using the budding yeast Saccharomyces cerevisiae as a model system. Because the known human targets of these drugs do not exist in yeast, we could employ the yeast gene deletion collections and parallel fitness profiling to explore potential off-target effects in a genome-wide manner. Among 214 tested, documented psychoactive drugs, we identified 81 compounds that inhibited wild-type yeast growth and were thus selected for genome-wide fitness profiling. Many of these drugs had a propensity to affect multiple cellular functions. The sensitivity profiles of half of the analyzed drugs were enriched for core cellular processes such as secretion, protein folding, RNA processing, and chromatin structure. Interestingly, fluoxetine (Prozac) interfered with establishment of cell polarity, cyproheptadine (Periactin) targeted essential genes with chromatin-remodeling roles, while paroxetine (Paxil) interfered with essential RNA metabolism genes, suggesting potential secondary drug targets. We also found that the more recently developed atypical antipsychotic clozapine (Clozaril) had no fewer off-target effects in yeast than the typical antipsychotics haloperidol (Haldol) and pimozide (Orap). Our results suggest that model organism pharmacogenetic studies provide a rational foundation for understanding the off-target effects of clinically important psychoactive agents and suggest a rational means both for devising compound derivatives with fewer side effects and for tailoring drug treatment to individual patient genotypes.

PROPHECY--a yeast phenome database, update 2006.

Connecting genotype to phenotype is fundamental in biomedical research and in our understanding of disease. Phenomics--the large-scale quantitative phenotypic analysis of genotypes on a genome-wide scale--connects automated data generation with the development of novel tools for phenotype data integration, mining and visualization. Our yeast phenomics database PROPHECY is available at http://prophecy.lundberg.gu.se. Via phenotyping of 984 heterozygous diploids for all essential genes the genotypes analysed and presented in PROPHECY have been extended and now include all genes in the yeast genome. Further, phenotypic data from gene overexpression of 574 membrane spanning proteins has recently been included. To facilitate the interpretation of quantitative phenotypic data we have developed a new phenotype display option, the Comparative Growth Curve Display, where growth curve differences for a large number of mutants compared with the wild type are easily revealed. In addition, PROPHECY now offers a more informative and intuitive first-sight display of its phenotypic data via its new summary page. We have also extended the arsenal of data analysis tools to include dynamic visualization of phenotypes along individual chromosomes. PROPHECY is an initiative to enhance the growing field of phenome bioinformatics.

Genetic pleiotropy in Saccharomyces cerevisiae quantified by high-resolution phenotypic profiling.

Genetic pleiotropy, the ability of a mutation in a single gene to give rise to multiple phenotypic outcomes, constitutes an important but incompletely understood biological phenomenon. We used a high-resolution and high-precision phenotypic profiling approach to quantify the fitness contribution of genes on the five smallest yeast chromosomes during different forms of environmental stress, selected to probe a wide diversity of physiological features. We found that the extent of pleiotropy is much higher than previously claimed; 17% of the yeast genes were pleiotropic whereof one-fifth were hyper-pleiotropic. Pleiotropic genes preferentially participate in functions related to determination of protein fate, cell growth and morphogenesis, signal transduction and transcription. Contrary to what has earlier been proposed we did not find experimental evidence for slower evolutionary rate of pleiotropic genes/proteins. We also refute the existence of phenotypic islands along chromosomes but report on a remarkable loss both of pleiotropy and of phenotypic penetrance towards chromosomal ends. Thus, the here reported features of pleiotropy both have implications on our understanding of evolutionary processes as well as the mechanisms underlying disease.

PROPHECY--a database for high-resolution phenomics.

The rapid recent evolution of the field phenomics--the genome-wide study of gene dispensability by quantitative analysis of phenotypes--has resulted in an increasing demand for new data analysis and visualization tools. Following the introduction of a novel approach for precise, genome-wide quantification of gene dispensability in Saccharomyces cerevisiae we here announce a public resource for mining, filtering and visualizing phenotypic data--the PROPHECY database. PROPHECY is designed to allow easy and flexible access to physiologically relevant quantitative data for the growth behaviour of mutant strains in the yeast deletion collection during conditions of environmental challenges. PROPHECY is publicly accessible at http://prophecy.lundberg.gu.se.

High-resolution yeast phenomics resolves different physiological features in the saline response.

We present a methodology for gene functional prediction based on extraction of physiologically relevant growth variables from all viable haploid yeast knockout mutants. This quantitative phenomics approach, here applied to saline cultivation, identified marginal but functionally important phenotypes and allowed the precise determination of time to adapt to an environmental challenge, rate of growth, and efficiency of growth. We identified approximately 500 salt-sensitive gene deletions, the majority of which were previously uncharacterized and displayed salt sensitivity for only one of the three physiological features. We also report a high correlation to protein-protein interaction data; in particular, several salt-sensitive subcellular networks indicating functional modules were revealed. In contrast, no correlation was found between gene dispensability and gene expression. It is proposed that high-resolution phenomics will be instrumental in systemwide descriptions of intragenomic functional networks.