Modelling genotype-phenotype relationships and human disease with genetic interaction networks

Evolution of gene regulatory networks by means of selection and random genetic drift

The evolution of a population by means of genetic drift and natural selection operating on a gene regulatory network (GRN) of an individual has not been scrutinized in depth. Thus, the relative importance of various evolutionary forces and processes on shaping genetic variability in GRNs is understudied. In this study, we implemented a simulation framework, called EvoNET, that simulates forward-in-time the evolution of GRNs in a population. The fitness effect of mutations is not constant, rather fitness of each individual is evaluated on the phenotypic level, by measuring its distance from an optimal phenotype. Each individual goes through a maturation period, where its GRN may reach an equilibrium, thus deciding its phenotype. Afterwards, individuals compete to produce the next generation. We examine properties of the GRN evolution, such as robustness against the deleterious effect of mutations and the role of genetic drift. We are able to confirm previous hypotheses regarding the effect of mutations and we provide new insights on the interplay between random genetic drift and natural selection.

Identification and validation of DNA methylation-driven gene PCDHB4 as a novel tumor suppressor for glioblastoma diagnosis and prognosis

Aberrant DNA methylation is a critical regulator of gene expression in the development and progression of glioblastoma (GBM). However, the impact of methylation-driven gene PCDHB4 changes on GBM occurrence and progression remains unclear. Therefore, this study aimed to identify the PCDHB4 gene for early diagnosis and prognostic evaluation and clarify its functional role in GBM. Methylation-driven gene PCDHB4 was selected for GBM using the multi-omics integration method based on publicly available data sets. The diagnostic capabilities of PCDHB4 methylation and 5-hydroxymethylcytosines were validated in tissue and blood cell-free DNA (cfDNA) samples, respectively. Combined survival analysis of PCDHB4 methylation and immune infiltration cells evaluated the prognostic predictive performance of GBM patients. We identified that the PCDHB4 gene achieved high discriminative capabilities for GBM and normal tissues with an area under the curve value of 0.941. PCDHB4 hypermethylation was observed in cfDNA blood samples from GBM patients. Compared with GBM patients with PCDHB4 hypermethylation level, patients with PCDHB4 hypomethylation level had significantly poorer overall survival (p = 0.035). In addition, GBM patients with PCDHB4 hypermethylation and high infiltration of CD4+ T cell activation level had a favorable survival (p = 0.026). Moreover, we demonstrated that mRNA expression of PCDHB4 was downregulated in GBM tissues and upregulated in GBM cell lines with PCDHB4 demethylation, and PCDHB4 overexpression inhibited GBM cell proliferation and migration. In summary, we discovered a novel methylation-driven gene PCDHB4 for the diagnosis and prognosis of GBM and demonstrated that PCDHB4 is a tumor suppressor in vitro experiments.

Analyzing omics data by feature combinations based on kernel functions

Defining meaningful feature (molecule) combinations can enhance the study of disease diagnosis and prognosis. However, feature combinations are complex and various in biosystems, and the existing methods examine the feature cooperation in a single, fixed pattern for all feature pairs, such as linear combination. To identify the appropriate combination between two features and evaluate feature combination more comprehensively, this paper adopts kernel functions to study feature relationships and proposes a new omics data analysis method KF-[Formula: see text]-TSP. Besides linear combination, KF-[Formula: see text]-TSP also explores the nonlinear combination of features, and allows hybridizing multiple kernel functions to evaluate feature interaction from multiple views. KF-[Formula: see text]-TSP selects [Formula: see text] > 0 top-scoring pairs to build an ensemble classifier. Experimental results show that KF-[Formula: see text]-TSP with multiple kernel functions which evaluates feature combinations from multiple views is better than that with only one kernel function. Meanwhile, KF-[Formula: see text]-TSP performs better than TSP family algorithms and the previous methods based on conversion strategy in most cases. It performs similarly to the popular machine learning methods in omics data analysis, but involves fewer feature pairs. In the procedure of physiological and pathological changes, molecular interactions can be both linear and nonlinear. Hence, KF-[Formula: see text]-TSP, which can measure molecular combination from multiple perspectives, can help to mine information closely related to physiological and pathological changes and study disease mechanism.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Learning high-order interactions for polygenic risk prediction

Within the framework of precision medicine, the stratification of individual genetic susceptibility based on inherited DNA variation has paramount relevance. However, one of the most relevant pitfalls of traditional Polygenic Risk Scores (PRS) approaches is their inability to model complex high-order non-linear SNP-SNP interactions and their effect on the phenotype (e.g. epistasis). Indeed, they incur in a computational challenge as the number of possible interactions grows exponentially with the number of SNPs considered, affecting the statistical reliability of the model parameters as well. In this work, we address this issue by proposing a novel PRS approach, called High-order Interactions-aware Polygenic Risk Score (hiPRS), that incorporates high-order interactions in modeling polygenic risk. The latter combines an interaction search routine based on frequent itemsets mining and a novel interaction selection algorithm based on Mutual Information, to construct a simple and interpretable weighted model of user-specified dimensionality that can predict a given binary phenotype. Compared to traditional PRSs methods, hiPRS does not rely on GWAS summary statistics nor any external information. Moreover, hiPRS differs from Machine Learning-based approaches that can include complex interactions in that it provides a readable and interpretable model and it is able to control overfitting, even on small samples. In the present work we demonstrate through a comprehensive simulation study the superior performance of hiPRS w.r.t. state of the art methods, both in terms of scoring performance and interpretability of the resulting model. We also test hiPRS against small sample size, class imbalance and the presence of noise, showcasing its robustness to extreme experimental settings. Finally, we apply hiPRS to a case study on real data from DACHS cohort, defining an interaction-aware scoring model to predict mortality of stage II-III Colon-Rectal Cancer patients treated with oxaliplatin.

Development and validation of 3-CpG methylation prognostic signature based on different survival indicators for colorectal cancer

Abnormal DNA methylation is considered a vital hallmark to regulate gene expression and influence the development and progression of colorectal cancer (CRC). Although CRC-related methylation prognostic models have been developed, their clinical application is limited due to the lack of external validation and extension to other survival evaluation indicators. Therefore, this study aimed to develop and validate novel methylation prognostic models correlated with different survival indicators for individualized prognosis prediction for CRC patients. The prognostic-related CpG sites of methylation-driven genes screened by the MethylMix algorithm were identified and validated in The Cancer Genome Atlas (TCGA) CRC methylation data and our methylation data. The prognostic models correlated with different survival evaluation indicators (overall survival [OS] and disease-free survival [DFS]) were developed and validated in the TCGA CRC dataset (N = 376) and our independent CRC dataset (N = 227). We utilized the combination of selected 3-CpG methylation sites in three genes (DAPP1, FAM3D, and PIGR) to construct a prognostic risk-score model. In the training dataset, Kaplan-Meier survival analysis demonstrated that high-risk patients had significantly poorer survival than low-risk patients (p_OS  = .0014; p_DFS  < .001). Then, the 3-CpG methylation signature was successfully validated as an independent predictor in the testing data set (p_OS  = .016; p_DFS  = .016). A prognostic nomogram was constructed and validated. Additionally, mediation analysis revealed the direct effect of the methylation signature on CRC prognosis (p_OS  = 9.149e-06; p_DFS  = .001). In summary, our study revealed that the 3-CpG methylation signature might be a potential prognostic indicator to facilitate individualized survival prediction for CRC patients.

Single gene, two diseases, and multiple clinical presentations: Biotin-thiamine-responsive basal ganglia disease

AIM

To present seven new genetically confirmed cases of biotin-thiamin-responsive basal ganglia disease (BTBGD) with different clinical and brain magnetic resonance imaging (MRI) characteristics.

MATERIAL AND METHODS

Genetic variants, clinical presentations, brain MRI findings, treatment response, and prognosis of seven selected patients with BTBGD, diagnosed with SLC19A3 mutations were described.

RESULTS

Among seven patients diagnosed with BTBGD, two had early infantile form, four had classic childhood form, and one was asymptomatic. Four different homozygous variants were found in the SLC19A3. Two patients with early infantile form presented with encephalopathy, dystonia, and refractory seizure in the neonatal period and have different variants. Their MRI findings were similar and pathognomonic for the early infantile form. Three siblings had same variants: one presented seizure and encephalopathy at the age of 4 months, one presented seizure at 14 years, and another was asymptomatic at 20 years. Only one of them had normal MRI findings, and the others MRI findings were similar and suggestive of the classic form. Other two siblings; one of them presented with developmental delay, seizure, and dystonia at 18 months and the other presented with subacute encephalopathy and ataxia at 20 months. Their MRI findings were also similar and suggestive of the classic form.

CONCLUSION

BTBGD may present with dissimilar clinical characteristics or remain asymptomatic for a long time period even in a family or patients with same variants. Brain MRI patterns may be important for the early diagnosis of BTBGD that would save children's lives.

Introducing the first whole genomes of nationals from the United Arab Emirates

Whole Genome Sequencing (WGS) provides an in depth description of genome variation. In the era of large-scale population genome projects, the assembly of ethnic-specific genomes combined with mapping human reference genomes of underrepresented populations has improved the understanding of human diversity and disease associations. In this study, for the first time, whole genome sequences of two nationals of the United Arab Emirates (UAE) at >27X coverage are reported. The two Emirati individuals were predominantly of Central/South Asian ancestry. An in-house customized pipeline using BWA, Picard followed by the GATK tools to map the raw data from whole genome sequences of both individuals was used. A total of 3,994,521 variants (3,350,574 Single Nucleotide Polymorphisms (SNPs) and 643,947 indels) were identified for the first individual, the UAE S001 sample. A similar number of variants, 4,031,580 (3,373,501 SNPs and 658,079 indels), were identified for UAE S002. Variants that are associated with diabetes, hypertension, increased cholesterol levels, and obesity were also identified in these individuals. These Whole Genome Sequences has provided a starting point for constructing a UAE reference panel which will lead to improvements in the delivery of precision medicine, quality of life for affected individuals and a reduction in healthcare costs. The information compiled will likely lead to the identification of target genes that could potentially lead to the development of novel therapeutic modalities.

Machine Learning Methods as a Tool for Predicting Risk of Illness Applying Next-Generation Sequencing Data

Next-generation sequencing (NGS) data present an untapped potential to improve microbial risk assessment (MRA) through increased specificity and redefinition of the hazard. Most of the MRA models do not account for differences in survivability and virulence among strains. The potential of machine learning algorithms for predicting the risk/health burden at the population level while inputting large and complex NGS data was explored with Listeria monocytogenes as a case study. Listeria data consisted of a percentage similarity matrix from genome assemblies of 38 and 207 strains of clinical and food origin, respectively. Basic Local Alignment (BLAST) was used to align the assemblies against a database of 136 virulence and stress resistance genes. The outcome variable was frequency of illness, which is the percentage of reported cases associated with each strain. These frequency data were discretized into seven ordinal outcome categories and used for supervised machine learning and model selection from five ensemble algorithms. There was no significant difference in accuracy between the models, and support vector machine with linear kernel was chosen for further inference (accuracy of 89% [95% CI: 68%, 97%]). The virulence genes FAM002725, FAM002728, FAM002729, InlF, InlJ, Inlk, IisY, IisD, IisX, IisH, IisB, lmo2026, and FAM003296 were important predictors of higher frequency of illness. InlF was uniquely truncated in the sequence type 121 strains. Most important risk predictor genes occurred at highest prevalence among strains from ready-to-eat, dairy, and composite foods. We foresee that the findings and approaches described offer the potential for rethinking the current approaches in MRA.

Novel and recurrent variants in AVPR2 in 19 families with X-linked congenital nephrogenic diabetes insipidus

Congenital nephrogenic diabetes insipidus (CNDI) is characterized by the reduced ability of renal collecting duct cells to reabsorb water in response to the antidiuretic effect of vasopressin. Chronic polyuria and polydipsia are the hallmarks of the disease. Approximately 90% of all patients with CNDI have X-linked inherited disease caused by variants in the arginine vasopressin receptor 2 (AVPR2) gene. We present genetic findings in 34 individuals from 19 kindreds including one or more family members with CNDI. Coding regions of AVPR2 were sequenced bi-directionally. We identified eight novel disease-causing variants in AVPR2, p.Arg68Alafs*124, p.Ser171Arg, p.Gln174Pro, p.Trp200Arg, p.Gly201Cys, p.Gly220Arg, p.Val226Glu, and p.Gln291Pro in nine kindreds. In all three families with more than one affected individual, the novel variants segregated with the disease. We also identified eight recurrent disease-causing variants, p.Val88Met, p.Leu111Valfs*80, p.Arg113Trp, p.Tyr124*, p.Ser167Leu, p.Thr207Asn, p.Arg247Alafs*12, and p.Arg337* in ten kindreds. Our findings contribute to the growing list of AVPR2 variants causing X-linked CNDI.

CONCLUSION

Being a rapid diagnostic tool for CNDI, direct sequencing of AVPR2 should be encouraged in newborns with familial predisposition to CNDI. What is Known: • Disease-causing variants in AVPR2 cause X-linked congenital nephrogenic diabetes insipidus (CNDI). • DNA sequencing of AVPR2 is rapid, facilitates differential diagnosis, early intervention, and genetic diagnosis thus reducing morbidity in CNDI. What is New: • We identified eight novel disease-causing variants in AVPR2: p.Arg68Alafs*124, p.Ser171Arg, p.Gln174Pro, p.Trp200Arg, p.Gly201Cys, p.Gly220Arg, p.Val226Glu, and p.Gln291Pro, thereby adding to the growing list of AVPR2 disease-causing variants and emphasizing the importance of genetic testing in CNDI.

Unravelling the complexity of signalling networks in cancer: A review of the increasing role for computational modelling

Cancer induction is a highly complex process involving hundreds of different inducers but whose eventual outcome is the same. Clearly, it is essential to understand how signalling pathways and networks generated by these inducers interact to regulate cell behaviour and create the cancer phenotype. While enormous strides have been made in identifying key networking profiles, the amount of data generated far exceeds our ability to understand how it all "fits together". The number of potential interactions is astronomically large and requires novel approaches and extreme computation methods to dissect them out. However, such methodologies have high intrinsic mathematical and conceptual content which is difficult to follow. This review explains how computation modelling is progressively finding solutions and also revealing unexpected and unpredictable nano-scale molecular behaviours extremely relevant to how signalling and networking are coherently integrated. It is divided into linked sections illustrated by numerous figures from the literature describing different approaches and offering visual portrayals of networking and major conceptual advances in the field. First, the problem of signalling complexity and data collection is illustrated for only a small selection of known oncogenes. Next, new concepts from biophysics, molecular behaviours, kinetics, organisation at the nano level and predictive models are presented. These areas include: visual representations of networking, Energy Landscapes and energy transfer/dissemination (entropy); diffusion, percolation; molecular crowding; protein allostery; quinary structure and fractal distributions; energy management, metabolism and re-examination of the Warburg effect. The importance of unravelling complex network interactions is then illustrated for some widely-used drugs in cancer therapy whose interactions are very extensive. Finally, use of computational modelling to develop micro- and nano- functional models ("bottom-up" research) is highlighted. The review concludes that computational modelling is an essential part of cancer research and is vital to understanding network formation and molecular behaviours that are associated with it. Its role is increasingly essential because it is unravelling the huge complexity of cancer induction otherwise unattainable by any other approach.

Insights into Hunter syndrome from the structure of iduronate-2-sulfatase

Hunter syndrome is a rare but devastating childhood disease caused by mutations in the IDS gene encoding iduronate-2-sulfatase, a crucial enzyme in the lysosomal degradation pathway of dermatan sulfate and heparan sulfate. These complex glycosaminoglycans have important roles in cell adhesion, growth, proliferation and repair, and their degradation and recycling in the lysosome is essential for cellular maintenance. A variety of disease-causing mutations have been identified throughout the IDS gene. However, understanding the molecular basis of the disease has been impaired by the lack of structural data. Here, we present the crystal structure of human IDS with a covalently bound sulfate ion in the active site. This structure provides essential insight into multiple mechanisms by which pathogenic mutations interfere with enzyme function, and a compelling explanation for severe Hunter syndrome phenotypes. Understanding the structural consequences of disease-associated mutations will facilitate the identification of patients that may benefit from specific tailored therapies.

Hunter syndrome is a lysosomal storage disease caused by mutations in the enzyme iduronate-2-sulfatase (IDS). Here, the authors present the IDS crystal structure and give mechanistic insights into mutations that cause Hunter syndrome.

Machine learning, statistical learning and the future of biological research in psychiatry

Psychiatric research has entered the age of 'Big Data'. Datasets now routinely involve thousands of heterogeneous variables, including clinical, neuroimaging, genomic, proteomic, transcriptomic and other 'omic' measures. The analysis of these datasets is challenging, especially when the number of measurements exceeds the number of individuals, and may be further complicated by missing data for some subjects and variables that are highly correlated. Statistical learning-based models are a natural extension of classical statistical approaches but provide more effective methods to analyse very large datasets. In addition, the predictive capability of such models promises to be useful in developing decision support systems. That is, methods that can be introduced to clinical settings and guide, for example, diagnosis classification or personalized treatment. In this review, we aim to outline the potential benefits of statistical learning methods in clinical research. We first introduce the concept of Big Data in different environments. We then describe how modern statistical learning models can be used in practice on Big Datasets to extract relevant information. Finally, we discuss the strengths of using statistical learning in psychiatric studies, from both research and practical clinical points of view.

Structural and Functional Characterization of a Caenorhabditis elegans Genetic Interaction Network within Pathways

A genetic interaction (GI) is defined when the mutation of one gene modifies the phenotypic expression associated with the mutation of a second gene. Genome-wide efforts to map GIs in yeast revealed structural and functional properties of a GI network. This provided insights into the mechanisms underlying the robustness of yeast to genetic and environmental insults, and also into the link existing between genotype and phenotype. While a significant conservation of GIs and GI network structure has been reported between distant yeast species, such a conservation is not clear between unicellular and multicellular organisms. Structural and functional characterization of a GI network in these latter organisms is consequently of high interest. In this study, we present an in-depth characterization of ~1.5K GIs in the nematode Caenorhabditis elegans. We identify and characterize six distinct classes of GIs by examining a wide-range of structural and functional properties of genes and network, including co-expression, phenotypical manifestations, relationship with protein-protein interaction dense subnetworks (PDS) and pathways, molecular and biological functions, gene essentiality and pleiotropy. Our study shows that GI classes link genes within pathways and display distinctive properties, specifically towards PDS. It suggests a model in which pathways are composed of PDS-centric and PDS-independent GIs coordinating molecular machines through two specific classes of GIs involving pleiotropic and non-pleiotropic connectors. Our study provides the first in-depth characterization of a GI network within pathways of a multicellular organism. It also suggests a model to understand better how GIs control system robustness and evolution.

Network biology has focused for years on protein-protein interaction (PPI) networks, identifying nodes with central structural functions and modules associated to bioprocesses, phenotypes and diseases. Network biology field moved to a higher level of abstraction, and started characterizing a less intuitive kind of interactions, called genetic interactions (GIs) or epistasis. Mostly due to technical challenges associated to the genome-wide mapping of GIs, these studies primarily focused on unicellular organisms. They uncovered modules embedded within the structure of these networks and started characterizing their relationship with PPI-network and biological functions. We provide here the first in-depth characterization of a network composed of ~600 GIs within signaling and metabolic pathways of a multicellular organism, the nematode Caenorhabditis elegans. We characterize the structure of this network, and the function of GI classes found in this network. We also discuss how these GI classes contribute to the genomic robustness and the adaptive evolution of multicellular organisms.

A survey of disease connections for CD4+ T cell master genes and their directly linked genes

Genome-wide association studies and other genetic analyses have identified a large number of genes and variants implicating a variety of disease etiological mechanisms. It is imperative for the study of human diseases to put these genetic findings into a coherent functional context. Here we use system biology tools to examine disease connections of five master genes for CD4+ T cell subtypes (TBX21, GATA3, RORC, BCL6, and FOXP3). We compiled a list of genes functionally interacting (protein-protein interaction, or by acting in the same pathway) with the master genes, then we surveyed the disease connections, either by experimental evidence or by genetic association. Embryonic lethal genes (also known as essential genes) are over-represented in master genes and their interacting genes (55% versus 40% in other genes). Transcription factors are significantly enriched among genes interacting with the master genes (63% versus 10% in other genes). Predicted haploinsufficiency is a feature of most these genes. Disease-connected genes are enriched in this list of genes: 42% of these genes have a disease connection according to Online Mendelian Inheritance in Man (OMIM) (versus 23% in other genes), and 74% are associated with some diseases or phenotype in a Genome Wide Association Study (GWAS) (versus 43% in other genes). Seemingly, not all of the diseases connected to genes surveyed were immune related, which may indicate pleiotropic functions of the master regulator genes and associated genes.

Quantitative Fitness Analysis Identifies exo1∆ and Other Suppressors or Enhancers of Telomere Defects in Schizosaccharomyces pombe

Synthetic genetic array (SGA) has been successfully used to identify genetic interactions in S. cerevisiae and S. pombe. In S. pombe, SGA methods use either cycloheximide (C) or heat shock (HS) to select double mutants before measuring colony size as a surrogate for fitness. Quantitative Fitness Analysis (QFA) is a different method for determining fitness of microbial strains. In QFA, liquid cultures are spotted onto solid agar and growth curves determined for each spot by photography and model fitting. Here, we compared the two S. pombe SGA methods and found that the HS method was more reproducible for us. We also developed a QFA procedure for S. pombe. We used QFA to identify genetic interactions affecting two temperature sensitive, telomere associated query mutations (taz1Δ and pot1-1). We identify exo1∆ and other gene deletions as suppressors or enhancers of S. pombe telomere defects. Our study identifies known and novel gene deletions affecting the fitness of strains with telomere defects. The interactions we identify may be relevant in human cells.

Thiamine transporter-2 deficiency: outcome and treatment monitoring

Background

The clinical characteristics distinguishing treatable thiamine transporter-2 deficiency (ThTR2) due to SLC19A3 genetic defects from the other devastating causes of Leigh syndrome are sparse.

Methods

We report the clinical follow-up after thiamine and biotin supplementation in four children with ThTR2 deficiency presenting with Leigh and biotin-thiamine-responsive basal ganglia disease phenotypes. We established whole-blood thiamine reference values in 106 non-neurological affected children and monitored thiamine levels in SLC19A3 patients after the initiation of treatment. We compared our results with those of 69 patients with ThTR2 deficiency after a review of the literature.

Results

At diagnosis, the patients were aged 1 month to 17 years, and all of them showed signs of acute encephalopathy, generalized dystonia, and brain lesions affecting the dorsal striatum and medial thalami. One patient died of septicemia, while the remaining patients evidenced clinical and radiological improvements shortly after the initiation of thiamine. Upon follow-up, the patients received a combination of thiamine (10–40 mg/kg/day) and biotin (1–2 mg/kg/day) and remained stable with residual dystonia and speech difficulties. After establishing reference values for the different age groups, whole-blood thiamine quantification was a useful method for treatment monitoring.

Conclusions

ThTR2 deficiency is a reversible cause of acute dystonia and Leigh encephalopathy in the pediatric years. Brain lesions affecting the dorsal striatum and medial thalami may be useful in the differential diagnosis of other causes of Leigh syndrome. Further studies are needed to validate the therapeutic doses of thiamine and how to monitor them in these patients.

Paternal high-fat diet consumption induces common changes in the transcriptomes of retroperitoneal adipose and pancreatic islet tissues in female rat offspring

We previously showed that paternal high-fat diet (HFD) consumption programs β-cell dysfunction in female rat offspring, together with transcriptome alterations in islets. Here we investigated the retroperitoneal white adipose tissue (RpWAT) transcriptome using gene and pathway enrichment and pathway analysis to determine whether commonly affected network topologies exist between these two metabolically related tissues. In RpWAT, 5108 genes were differentially expressed due to a paternal HFD; the top 5 significantly enriched networks identified by pathway analysis in offspring of HFD fathers compared with those of fathers fed control diet were: mitochondrial and cellular response to stress, telomerase signaling, cell death and survival, cell cycle, cellular growth and proliferation, and cancer. A total of 187 adipose olfactory receptor genes were down-regulated. Interrogation against the islet transcriptome identified specific gene networks and pathways, including olfactory receptor genes that were similarly affected in both tissues (411 common genes, P<0.05). In particular, we highlight a common molecular network, cell cycle and cancer, with the same hub gene, Myc, suggesting early onset developmental changes that persist, shared responses to programmed systemic factors, or crosstalk between tissues. Thus, paternal HFD consumption triggers unique gene signatures, consistent with premature aging and chronic degenerative disorders, in both RpWAT and pancreatic islets of daughters.

Genetic interaction networks: better understand to better predict

A genetic interaction (GI) between two genes generally indicates that the phenotype of a double mutant differs from what is expected from each individual mutant. In the last decade, genome scale studies of quantitative GIs were completed using mainly synthetic genetic array technology and RNA interference in yeast and Caenorhabditis elegans. These studies raised questions regarding the functional interpretation of GIs, the relationship of genetic and molecular interaction networks, the usefulness of GI networks to infer gene function and co-functionality, the evolutionary conservation of GI, etc. While GIs have been used for decades to dissect signaling pathways in genetic models, their functional interpretations are still not trivial. The existence of a GI between two genes does not necessarily imply that these two genes code for interacting proteins or that the two genes are even expressed in the same cell. In fact, a GI only implies that the two genes share a functional relationship. These two genes may be involved in the same biological process or pathway; or they may also be involved in compensatory pathways with unrelated apparent function. Considering the powerful opportunity to better understand gene function, genetic relationship, robustness and evolution, provided by a genome-wide mapping of GIs, several in silico approaches have been employed to predict GIs in unicellular and multicellular organisms. Most of these methods used weighted data integration. In this article, we will review the later knowledge acquired on GI networks in metazoans by looking more closely into their relationship with pathways, biological processes and molecular complexes but also into their modularity and organization. We will also review the different in silico methods developed to predict GIs and will discuss how the knowledge acquired on GI networks can be used to design predictive tools with higher performances.

Predicting functional gene interactions with the hierarchical interaction score

Systems biology aims to unravel the vast network of functional interactions that govern biological systems. To date, the inference of gene interactions from large-scale 'omics data is typically achieved using correlations. We present the hierarchical interaction score (HIS) and show that the HIS outperforms commonly used methods in the inference of functional interactions between genes measured in large-scale experiments, making it a valuable statistic for systems biology.

Genetic variants and their interactions in disease risk prediction - machine learning and network perspectives

A central challenge in systems biology and medical genetics is to understand how interactions among genetic loci contribute to complex phenotypic traits and human diseases. While most studies have so far relied on statistical modeling and association testing procedures, machine learning and predictive modeling approaches are increasingly being applied to mining genotype-phenotype relationships, also among those associations that do not necessarily meet statistical significance at the level of individual variants, yet still contributing to the combined predictive power at the level of variant panels. Network-based analysis of genetic variants and their interaction partners is another emerging trend by which to explore how sub-network level features contribute to complex disease processes and related phenotypes. In this review, we describe the basic concepts and algorithms behind machine learning-based genetic feature selection approaches, their potential benefits and limitations in genome-wide setting, and how physical or genetic interaction networks could be used as a priori information for providing improved predictive power and mechanistic insights into the disease networks. These developments are geared toward explaining a part of the missing heritability, and when combined with individual genomic profiling, such systems medicine approaches may also provide a principled means for tailoring personalized treatment strategies in the future.

Quantitative genetic-interaction mapping in mammalian cells

Mapping genetic interactions (GIs) by simultaneously perturbing pairs of genes is a powerful tool for understanding complex biological phenomena. Here we describe an experimental platform for generating quantitative GI maps in mammalian cells using a combinatorial RNA interference strategy. We performed ∼11,000 pairwise knockdowns in mouse fibroblasts, focusing on 130 factors involved in chromatin regulation to create a GI map. Comparison of the GI and protein-protein interaction (PPI) data revealed that pairs of genes exhibiting positive GIs and/or similar genetic profiles were predictive of the corresponding proteins being physically associated. The mammalian GI map identified pathways and complexes but also resolved functionally distinct submodules within larger protein complexes. By integrating GI and PPI data, we created a functional map of chromatin complexes in mouse fibroblasts, revealing that the PAF complex is a central player in the mammalian chromatin landscape.

Genotype to phenotype: lessons from model organisms for human genetics

To what extent can variation in phenotypic traits such as disease risk be accurately predicted in individuals? In this Review, I highlight recent studies in model organisms that are relevant both to the challenge of accurately predicting phenotypic variation from individual genome sequences ('whole-genome reverse genetics') and for understanding why, in many cases, this may be impossible. These studies argue that only by combining genetic knowledge with in vivo measurements of biological states will it be possible to make accurate genetic predictions for individual humans.

Epigenetic epistatic interactions constrain the evolution of gene expression

Reduced activity of two genes in combination often has a more detrimental effect than expected. Such epistatic interactions not only occur when genes are mutated but also due to variation in gene expression, including among isogenic individuals in a controlled environment. We hypothesized that these 'epigenetic' epistatic interactions could place important constraints on the evolution of gene expression. Consistent with this, we show here that yeast genes with many epistatic interaction partners typically show low expression variation among isogenic individuals and low variation across different conditions. In addition, their expression tends to remain stable in response to the accumulation of mutations and only diverges slowly between strains and species. Yeast promoter architectures, the retention of gene duplicates, and the divergence of expression between humans and chimps are also consistent with selective pressure to reduce the likelihood of harmful epigenetic epistatic interactions. Based on these and previous analyses, we propose that the tight regulation of epistatic interaction network hubs makes an important contribution to the maintenance of a robust, 'canalized' phenotype. Moreover, that epigenetic epistatic interactions may contribute substantially to fitness defects when single genes are deleted.

Systems biology in animal sciences

Systems biology is a rapidly expanding field of research and is applied in a number of biological disciplines. In animal sciences, omics approaches are increasingly used, yielding vast amounts of data, but systems biology approaches to extract understanding from these data of biological processes and animal traits are not yet frequently used. This paper aims to explain what systems biology is and which areas of animal sciences could benefit from systems biology approaches. Systems biology aims to understand whole biological systems working as a unit, rather than investigating their individual components. Therefore, systems biology can be considered a holistic approach, as opposed to reductionism. The recently developed 'omics' technologies enable biological sciences to characterize the molecular components of life with ever increasing speed, yielding vast amounts of data. However, biological functions do not follow from the simple addition of the properties of system components, but rather arise from the dynamic interactions of these components. Systems biology combines statistics, bioinformatics and mathematical modeling to integrate and analyze large amounts of data in order to extract a better understanding of the biology from these huge data sets and to predict the behavior of biological systems. A 'system' approach and mathematical modeling in biological sciences are not new in itself, as they were used in biochemistry, physiology and genetics long before the name systems biology was coined. However, the present combination of mass biological data and of computational and modeling tools is unprecedented and truly represents a major paradigm shift in biology. Significant advances have been made using systems biology approaches, especially in the field of bacterial and eukaryotic cells and in human medicine. Similarly, progress is being made with 'system approaches' in animal sciences, providing exciting opportunities to predict and modulate animal traits.

Hierarchical modularity and the evolution of genetic interactomes across species

To date, cross-species comparisons of genetic interactomes have been restricted to small or functionally related gene sets, limiting our ability to infer evolutionary trends. To facilitate a more comprehensive analysis, we constructed a genome-scale epistasis map (E-MAP) for the fission yeast Schizosaccharomyces pombe, providing phenotypic signatures for ~60% of the nonessential genome. Using these signatures, we generated a catalog of 297 functional modules, and we assigned function to 144 previously uncharacterized genes, including mRNA splicing and DNA damage checkpoint factors. Comparison with an integrated genetic interactome from the budding yeast Saccharomyces cerevisiae revealed a hierarchical model for the evolution of genetic interactions, with conservation highest within protein complexes, lower within biological processes, and lowest between distinct biological processes. Despite the large evolutionary distance and extensive rewiring of individual interactions, both networks retain conserved features and display similar levels of functional crosstalk between biological processes, suggesting general design principles of genetic interactomes.

Mining for genotype-phenotype relations in Saccharomyces using partial least squares

Background

Multivariate approaches are important due to their versatility and applications in many fields as it provides decisive advantages over univariate analysis in many ways. Genome wide association studies are rapidly emerging, but approaches in hand pay less attention to multivariate relation between genotype and phenotype. We introduce a methodology based on a BLAST approach for extracting information from genomic sequences and Soft- Thresholding Partial Least Squares (ST-PLS) for mapping genotype-phenotype relations.

Results

Applying this methodology to an extensive data set for the model yeast Saccharomyces cerevisiae, we found that the relationship between genotype-phenotype involves surprisingly few genes in the sense that an overwhelmingly large fraction of the phenotypic variation can be explained by variation in less than 1% of the full gene reference set containing 5791 genes. These phenotype influencing genes were evolving 20% faster than non-influential genes and were unevenly distributed over cellular functions, with strong enrichments in functions such as cellular respiration and transposition. These genes were also enriched with known paralogs, stop codon variations and copy number variations, suggesting that such molecular adjustments have had a disproportionate influence on Saccharomyces yeasts recent adaptation to environmental changes in its ecological niche.

Conclusions

BLAST and PLS based multivariate approach derived results that adhere to the known yeast phylogeny and gene ontology and thus verify that the methodology extracts a set of fast evolving genes that capture the phylogeny of the yeast strains. The approach is worth pursuing, and future investigations should be made to improve the computations of genotype signals as well as variable selection procedure within the PLS framework.

Molecular mechanisms of epistasis within and between genes

'Disease-causing' mutations do not cause disease in all individuals. One possible important reason for this is that the outcome of a mutation can depend upon other genetic variants in a genome. These epistatic interactions between mutations occur both within and between molecules, and studies in model organisms show that they are extremely prevalent. However, epistatic interactions are still poorly understood at the molecular level, and consequently difficult to predict de novo. Here I provide an overview of our current understanding of the molecular mechanisms that can cause epistasis, and areas where more research is needed. A more complete understanding of epistasis will be vital for making accurate predictions about the phenotypes of individuals.

Synthetic biology: putting synthesis into biology

The ability to manipulate living organisms is at the heart of a range of emerging technologies that serve to address important and current problems in environment, energy, and health. However, with all its complexity and interconnectivity, biology has for many years been recalcitrant to engineering manipulations. The recent advances in synthesis, analysis, and modeling methods have finally provided the tools necessary to manipulate living systems in meaningful ways and have led to the coining of a field named synthetic biology. The scope of synthetic biology is as complicated as life itself—encompassing many branches of science and across many scales of application. New DNA synthesis and assembly techniques have made routine customization of very large DNA molecules. This in turn has allowed the incorporation of multiple genes and pathways. By coupling these with techniques that allow for the modeling and design of protein functions, scientists have now gained the tools to create completely novel biological machineries. Even the ultimate biological machinery—a self‐replicating organism—is being pursued at this moment. The aim of this article is to dissect and organize these various components of synthetic biology into a coherent picture. WIREs Syst Biol Med 2011 3 7–20 DOI: 10.1002/wsbm.104

This article is categorized under: 1

Analytical and Computational Methods > Computational Methods

2

Laboratory Methods and Technologies > Genetic/Genomic Methods

3

Laboratory Methods and Technologies > Metabolomics

Quantitative maps of genetic interactions in yeast - comparative evaluation and integrative analysis

Background

High-throughput genetic screening approaches have enabled systematic means to study how interactions among gene mutations contribute to quantitative fitness phenotypes, with the aim of providing insights into the functional wiring diagrams of genetic interaction networks on a global scale. However, it is poorly known how well these quantitative interaction measurements agree across the screening approaches, which hinders their integrated use toward improving the coverage and quality of the genetic interaction maps in yeast and other organisms.

Results

Using large-scale data matrices from epistatic miniarray profiling (E-MAP), genetic interaction mapping (GIM), and synthetic genetic array (SGA) approaches, we carried out here a systematic comparative evaluation among these quantitative maps of genetic interactions in yeast. The relatively low association between the original interaction measurements or their customized scores could be improved using a matrix-based modelling framework, which enables the use of single- and double-mutant fitness estimates and measurements, respectively, when scoring genetic interactions. Toward an integrative analysis, we show how the detections from the different screening approaches can be combined to suggest novel positive and negative interactions which are complementary to those obtained using any single screening approach alone. The matrix approximation procedure has been made available to support the design and analysis of the future screening studies.

Conclusions

We have shown here that even if the correlation between the currently available quantitative genetic interaction maps in yeast is relatively low, their comparability can be improved by means of our computational matrix approximation procedure, which will enable integrative analysis and detection of a wider spectrum of genetic interactions using data from the complementary screening approaches.

Defining the pathogenicity of creatine deficiency syndrome

This work examined nine patients with creatine deficiency syndrome (CDS): six with a creatine transport (CRTR) defect and three with a GAMT defect. Eleven nucleotide variations were detected: six in SLC6A8 and five in GAMT. These changes were analyzed at the mRNA level and specific alleles (most of which bore premature stop codons) were selected as nulls because they provoked nonsense-mediated decay activation. The impact of these CDS mutations on metabolic stress (ROS production, p38MAPK activation, aberrant proliferation and apoptosis) was analyzed in patient fibroblast cultures. Oxidative stress contributed toward the severe form of CDS, with increases seen in the intracellular ROS content and the percentage of apoptotic cells. An altered cell cycle was also seen in a number of CRTR and GAMT fibroblast cell lines (mostly those carrying null alleles). p38MAPK activation only correlated with oxidative stress in the CRTR cells. Based on intracellular creatine levels, the contribution of energy depletion toward metabolic stress was demonstrable only in selected CRTR cells. Together, these findings suggest that the apoptotic response to genotoxic damage in the present CDS cells may have been triggered by different cell signaling pathways. They also suggest that reducing oxidative stress could be helpful in treating CDS. Hum Mutat 32:1-10, 2011. © 2011 Wiley-Liss, Inc.

Finding the "dark matter" in human and yeast protein network prediction and modelling

Accurate modelling of biological systems requires a deeper and more complete knowledge about the molecular components and their functional associations than we currently have. Traditionally, new knowledge on protein associations generated by experiments has played a central role in systems modelling, in contrast to generally less trusted bio-computational predictions. However, we will not achieve realistic modelling of complex molecular systems if the current experimental designs lead to biased screenings of real protein networks and leave large, functionally important areas poorly characterised. To assess the likelihood of this, we have built comprehensive network models of the yeast and human proteomes by using a meta-statistical integration of diverse computationally predicted protein association datasets. We have compared these predicted networks against combined experimental datasets from seven biological resources at different level of statistical significance. These eukaryotic predicted networks resemble all the topological and noise features of the experimentally inferred networks in both species, and we also show that this observation is not due to random behaviour. In addition, the topology of the predicted networks contains information on true protein associations, beyond the constitutive first order binary predictions. We also observe that most of the reliable predicted protein associations are experimentally uncharacterised in our models, constituting the hidden or "dark matter" of networks by analogy to astronomical systems. Some of this dark matter shows enrichment of particular functions and contains key functional elements of protein networks, such as hubs associated with important functional areas like the regulation of Ras protein signal transduction in human cells. Thus, characterising this large and functionally important dark matter, elusive to established experimental designs, may be crucial for modelling biological systems. In any case, these predictions provide a valuable guide to these experimentally elusive regions.

"CRASH"ing with the worm: insights into L1CAM functions and mechanisms

The L1 family of cell adhesion molecules (L1CAMs) in vertebrates has long been studied for its roles in nervous system development and function. Members of this family have been associated with distinct neurological disorders that include CRASH, autism, 3p syndrome, and schizophrenia. The conservation of L1CAMs in Drosophila and Caenorhabditis elegans allows the opportunity to take advantage of these simple model organisms and their accessible genetic manipulations to dissect L1CAM functions and mechanisms of action. This review summarizes the discoveries of L1CAMs made in C. elegans, showcasing this simple model organism as a powerful system to uncover L1CAM mechanisms and roles in healthy and diseased states.

Genome-wide scoring of positive and negative epistasis through decomposition of quantitative genetic interaction fitness matrices

Recent technological developments in genetic screening approaches have offered the means to start exploring quantitative genotype-phenotype relationships on a large-scale. What remains unclear is the extent to which the quantitative genetic interaction datasets can distinguish the broad spectrum of interaction classes, as compared to existing information on mutation pairs associated with both positive and negative interactions, and whether the scoring of varying degrees of such epistatic effects could be improved by computational means. To address these questions, we introduce here a computational approach for improving the quantitative discrimination power encoded in the genetic interaction screening data. Our matrix approximation model decomposes the original double-mutant fitness matrix into separate components, representing variability across the array and query mutants, which can be utilized for estimating and correcting the single-mutant fitness effects, respectively. When applied to three large-scale quantitative interaction datasets in yeast, we could improve the accuracy of scoring various interaction classes beyond that obtained with the original fitness data, especially in synthetic genetic array (SGA) and in genetic interaction mapping (GIM) datasets. In addition to the known pairs of interactions used in the evaluation of the computational approach, a number of novel interaction pairs were also predicted, along with underlying biological mechanisms, which remained undetected by the original datasets. It was shown that the optimal choice of the scoring function depends heavily on the screening approach and on the interaction class under analysis. Moreover, a simple preprocessing of the fitness matrix could further enhance the discrimination power of the epistatic miniarray profiling (E-MAP) dataset. These systematic evaluation results provide in-depth information on the optimal analysis of the future, large-scale screening experiments. In general, the modeling framework, enabling accurate identification and classification of genetic interactions, provides a solid basis for completing and mining the genetic interaction networks in yeast and other organisms.

Predicting genetic modifier loci using functional gene networks

Most phenotypes are genetically complex, with contributions from mutations in many different genes. Mutations in more than one gene can combine synergistically to cause phenotypic change, and systematic studies in model organisms show that these genetic interactions are pervasive. However, in human association studies such nonadditive genetic interactions are very difficult to identify because of a lack of statistical power--simply put, the number of potential interactions is too vast. One approach to resolve this is to predict candidate modifier interactions between loci, and then to specifically test these for associations with the phenotype. Here, we describe a general method for predicting genetic interactions based on the use of integrated functional gene networks. We show that in both Saccharomyces cerevisiae and Caenorhabditis elegans a single high-coverage, high-quality functional network can successfully predict genetic modifiers for the majority of genes. For C. elegans we also describe the construction of a new, improved, and expanded functional network, WormNet 2. Using this network we demonstrate how it is possible to rapidly expand the number of modifier loci known for a gene, predicting and validating new genetic interactions for each of three signal transduction genes. We propose that this approach, termed network-guided modifier screening, provides a general strategy for predicting genetic interactions. This work thus suggests that a high-quality integrated human gene network will provide a powerful resource for modifier locus discovery in many different diseases.

Clinical review 2: Genetic determinants of bone density and fracture risk--state of the art and future directions

CONTEXT

Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved.

EVIDENCE ACQUISITION AND SYNTHESIS

In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward.

CONCLUSIONS

Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

FAM-MDR: a flexible family-based multifactor dimensionality reduction technique to detect epistasis using related individuals

We propose a novel multifactor dimensionality reduction method for epistasis detection in small or extended pedigrees, FAM-MDR. It combines features of the Genome-wide Rapid Association using Mixed Model And Regression approach (GRAMMAR) with Model-Based MDR (MB-MDR). We focus on continuous traits, although the method is general and can be used for outcomes of any type, including binary and censored traits. When comparing FAM-MDR with Pedigree-based Generalized MDR (PGMDR), which is a generalization of Multifactor Dimensionality Reduction (MDR) to continuous traits and related individuals, FAM-MDR was found to outperform PGMDR in terms of power, in most of the considered simulated scenarios. Additional simulations revealed that PGMDR does not appropriately deal with multiple testing and consequently gives rise to overly optimistic results. FAM-MDR adequately deals with multiple testing in epistasis screens and is in contrast rather conservative, by construction. Furthermore, simulations show that correcting for lower order (main) effects is of utmost importance when claiming epistasis. As Type 2 Diabetes Mellitus (T2DM) is a complex phenotype likely influenced by gene-gene interactions, we applied FAM-MDR to examine data on glucose area-under-the-curve (GAUC), an endophenotype of T2DM for which multiple independent genetic associations have been observed, in the Amish Family Diabetes Study (AFDS). This application reveals that FAM-MDR makes more efficient use of the available data than PGMDR and can deal with multi-generational pedigrees more easily. In conclusion, we have validated FAM-MDR and compared it to PGMDR, the current state-of-the-art MDR method for family data, using both simulations and a practical dataset. FAM-MDR is found to outperform PGMDR in that it handles the multiple testing issue more correctly, has increased power, and efficiently uses all available information.

Genes confer similar robustness to environmental, stochastic, and genetic perturbations in yeast

Gene inactivation often has little or no apparent consequence for the phenotype of an organism. This property-enetic (or mutational) robustness-is pervasive, and has important implications for disease and evolution, but is not well understood. Dating back to at least Waddington, it has been suggested that mutational robustness may be related to the requirement to withstand environmental or stochastic perturbations. Here I show that global quantitative data from yeast are largely consistent with this idea. Considering the effects of mutations in all nonessential genes shows that genes that confer robustness to environmental or stochastic change also buffer the effects of genetic change, and with similar efficacy. This means that selection during evolution for environmental or stochastic robustness (also referred to as canalization) may frequently have the side effect of increasing genetic robustness. A dynamic environment may therefore promote the evolution of phenotypic complexity. It also means that "hub" genes in genetic interaction (synthetic lethal) networks are generally genes that confer environmental resilience and phenotypic stability.

Gene methylation profiles of normal mucosa, and benign and malignant colorectal tumors identify early onset markers

It is increasingly clear that complex networks of relationships between genes and/or proteins govern neoplastic processes. Our understanding of these networks is expanded by the use of functional genomic and proteomic approaches in addition to computational modeling. Concurrently, whole-genome association scans and mutational screens of cancer genomes identify novel cancer genes. Together, these analyses have vastly increased our knowledge of cancer, in terms of both "part lists" and their functional associations. However, genetic interactions have hitherto only been studied in depth in model organisms and remain largely unknown for human systems. Here, we discuss the importance and potential benefits of identifying genetic interactions at the human genome level for creating a better understanding of cancer susceptibility and progression and developing novel effective anticancer therapies. We examine gene expression profiles in the presence and absence of co-amplification of the 8q24 and 20q13 chromosomal regions in breast tumors to illustrate the molecular consequences and complexity of genetic interactions and their role in tumorigenesis. Finally, we highlight current strategies for targeting tumor dependencies and outline potential matrix screening designs for uncovering molecular vulnerabilities in cancer cells.

Predicting quantitative genetic interactions by means of sequential matrix approximation

Despite the emerging experimental techniques for perturbing multiple genes and measuring their quantitative phenotypic effects, genetic interactions have remained extremely difficult to predict on a large scale. Using a recent high-resolution screen of genetic interactions in yeast as a case study, we investigated whether the extraction of pertinent information encoded in the quantitative phenotypic measurements could be improved by computational means. By taking advantage of the observation that most gene pairs in the genetic interaction screens have no significant interactions with each other, we developed a sequential approximation procedure which ranks the mutation pairs in order of evidence for a genetic interaction. The sequential approximations can efficiently remove background variation in the double-mutation screens and give increasingly accurate estimates of the single-mutant fitness measurements. Interestingly, these estimates not only provide predictions for genetic interactions which are consistent with those obtained using the measured fitness, but they can even significantly improve the accuracy with which one can distinguish functionally-related gene pairs from the non-interacting pairs. The computational approach, in general, enables an efficient exploration and classification of genetic interactions in other studies and systems as well.

Description of International Caenorhabditis elegans Experiment first flight (ICE-FIRST)

Traveling, living and working in space is now a reality. The number of people and length of time in space is increasing. With new horizons for exploration it becomes more important to fully understand and provide countermeasures to the effects of the space environment on the human body. In addition, space provides a unique laboratory to study how life and physiologic functions adapt from the cellular level to that of the entire organism. Caenorhabditis elegans is a genetic model organism used to study physiology on Earth. Here we provide a description of the rationale, design, methods, and space culture validation of the ICE-FIRST payload, which engaged C. elegans researchers from four nations. Here we also show C. elegans growth and development proceeds essentially normally in a chemically defined liquid medium on board the International Space Station (10.9 day round trip). By setting flight constraints first and bringing together established C. elegans researchers second, we were able to use minimal stowage space to successfully return a total of 53 independent samples, each containing more than a hundred individual animals, to investigators within one year of experiment concept. We believe that in the future, bringing together individuals with knowledge of flight experiment operations, flight hardware, space biology, and genetic model organisms should yield similarly successful payloads.

Extracting biology from high-dimensional biological data

The promise of the genome project was that a complete sequence would provide us with information that would transform biology and medicine. But the 'parts list' that has emerged from the genome project is far from the 'wiring diagram' and 'circuit logic' we need to understand the link between genotype, environment and phenotype. While genomic technologies such as DNA microarrays, proteomics and metabolomics have given us new tools and new sources of data to address these problems, a number of crucial elements remain to be addressed before we can begin to close the loop and develop a predictive quantitative biology that is the stated goal of so much of current biological research, including systems biology. Our approach to this problem has largely been one of integration, bringing together a vast wealth of information to better interpret the experimental data we are generating in genomic assays and creating publicly available databases and software tools to facilitate the work of others. Recently, we have used a similar approach to trying to understand the biological networks that underlie the phenotypic responses we observe and starting us on the road to developing a predictive biology.

Genes differentially expressed by Aspergillus flavus strains after loss of aflatoxin production by serial transfers

Aflatoxins are carcinogenic fungal secondary metabolites produced by Aspergillus flavus and other closely related species. Levels of aflatoxins in agricultural commodities are stringently regulated by many countries because of the health hazard, and thus, aflatoxins are of major concern to both producers and consumers. A cluster of genes responsible for aflatoxin biosynthesis has been identified; however, expression of these genes is a complex and poorly understood phenomenon. To better understand the molecular events that are associated with aflatoxin production, three separate nonaflatoxigenic A. flavus strains were produced through serial transfers of aflatoxigenic parental strains. The three independent aflatoxigenic/nonaflatoxigenic pairs were compared via transcription profiling by microarray analyses. Cross comparisons identified 22 features in common between the aflatoxigenic/nonaflatoxigenic pairs. Physical mapping of the 22 features using the Aspergillus oryzae genome sequence for reference identified 16 unique genes. Aflatoxin biosynthetic and regulatory gene expression levels were not significantly different between the aflatoxigenic/nonaflatoxigenic pairs, which suggests that the inability to produce aflatoxins is not due to decreased expression of known biosynthetic or regulatory genes. Of the 16 in common genes, only one gene homologous to glutathione S-transferase genes showed higher expression in the nonaflatoxigenic progeny relative to the parental strains. This gene, named hcc, was selected for over-expression in an aflatoxigenic A. flavus strain to determine if it was directly responsible for loss of aflatoxin production. Although hcc transformants showed six- to ninefold increase in expression, no discernible changes in colony morphology or aflatoxin production were detected. Possible roles of hcc and other identified genes are discussed in relation to regulation of aflatoxin biosynthesis.

A global analysis of genetic interactions in Caenorhabditis elegans

Background

Understanding gene function and genetic relationships is fundamental to our efforts to better understand biological systems. Previous studies systematically describing genetic interactions on a global scale have either focused on core biological processes in protozoans or surveyed catastrophic interactions in metazoans. Here, we describe a reliable high-throughput approach capable of revealing both weak and strong genetic interactions in the nematode Caenorhabditis elegans.

Results

We investigated interactions between 11 'query' mutants in conserved signal transduction pathways and hundreds of 'target' genes compromised by RNA interference (RNAi). Mutant-RNAi combinations that grew more slowly than controls were identified, and genetic interactions inferred through an unbiased global analysis of the interaction matrix. A network of 1,246 interactions was uncovered, establishing the largest metazoan genetic-interaction network to date. We refer to this approach as systematic genetic interaction analysis (SGI). To investigate how genetic interactions connect genes on a global scale, we superimposed the SGI network on existing networks of physical, genetic, phenotypic and coexpression interactions. We identified 56 putative functional modules within the superimposed network, one of which regulates fat accumulation and is coordinated by interactions with bar-1(ga80), which encodes a homolog of β-catenin. We also discovered that SGI interactions link distinct subnetworks on a global scale. Finally, we showed that the properties of genetic networks are conserved between C. elegans and Saccharomyces cerevisiae, but that the connectivity of interactions within the current networks is not.

Conclusions

Synthetic genetic interactions may reveal redundancy among functional modules on a global scale, which is a previously unappreciated level of organization within metazoan systems. Although the buffering between functional modules may differ between species, studying these differences may provide insight into the evolution of divergent form and function.