Search

Interactions between proteins are highly conserved across species. As a result, the molecular basis of multiple diseases affecting humans can be studied in model organisms that offer many alternative experimental opportunities. One such organism-Caenorhabditis elegans-has been used to produce much molecular quantitative genetics and systems biology data over the past decade. We present WormQTL(HD) (Human Disease), a database that quantitatively and systematically links expression Quantitative Trait Loci (eQTL) findings in C. elegans to gene-disease associations in man. WormQTL(HD), available online at http://www.wormqtl-hd.org, is a user-friendly set of tools to reveal functionally coherent, evolutionary conserved gene networks. These can be used to predict novel gene-to-gene associations and the functions of genes underlying the disease of interest. We created a new database that links C. elegans eQTL data sets to human diseases (34 337 gene-disease associations from OMIM, DGA, GWAS Central and NHGRI GWAS Catalogue) based on overlapping sets of orthologous genes associated to phenotypes in these two species. We utilized QTL results, high-throughput molecular phenotypes, classical phenotypes and genotype data covering different developmental stages and environments from WormQTL database. All software is available as open source, built on MOLGENIS and xQTL workbench.

Here, we present WormQTL (http://www.wormqtl.org), an easily accessible database enabling search, comparative analysis and meta-analysis of all data on variation in Caenorhabditis spp. Over the past decade, Caenorhabditis elegans has become instrumental for molecular quantitative genetics and the systems biology of natural variation. These efforts have resulted in a valuable amount of phenotypic, high-throughput molecular and genotypic data across different developmental worm stages and environments in hundreds of C. elegans strains. WormQTL provides a workbench of analysis tools for genotype-phenotype linkage and association mapping based on but not limited to R/qtl (http://www.rqtl.org). All data can be uploaded and downloaded using simple delimited text or Excel formats and are accessible via a public web user interface for biologists and R statistic and web service interfaces for bioinformaticians, based on open source MOLGENIS and xQTL workbench software. WormQTL welcomes data submissions from other worm researchers.

Quantitative trait loci (QTL) hotspots (genomic locations affecting many traits) are a common feature in genetical genomics studies and are biologically interesting since they may harbor critical regulators. Therefore, statistical procedures to assess the significance of hotspots are of key importance. One approach, randomly allocating observed QTL across the genomic locations separately by trait, implicitly assumes all traits are uncorrelated. Recently, an empirical test for QTL hotspots was proposed on the basis of the number of traits that exceed a predetermined LOD value, such as the standard permutation LOD threshold. The permutation null distribution of the maximum number of traits across all genomic locations preserves the correlation structure among the phenotypes, avoiding the detection of spurious hotspots due to nongenetic correlation induced by uncontrolled environmental factors and unmeasured variables. However, by considering only the number of traits above a threshold, without accounting for the magnitude of the LOD scores, relevant information is lost. In particular, biologically interesting hotspots composed of a moderate to small number of traits with strong LOD scores may be neglected as nonsignificant. In this article we propose a quantile-based permutation approach that simultaneously accounts for the number and the LOD scores of traits within the hotspots. By considering a sliding scale of mapping thresholds, our method can assess the statistical significance of both small and large hotspots. Although the proposed approach can be applied to any type of heritable high-volume "omic" data set, we restrict our attention to expression (e)QTL analysis. We assess and compare the performances of these three methods in simulations and we illustrate how our approach can effectively assess the significance of moderate and small hotspots with strong LOD scores in a yeast expression data set.

Summary: xQTL workbench is a scalable web platform for the mapping of quantitative trait loci (QTLs) at multiple levels: for example gene expression (eQTL), protein abundance (pQTL), metabolite abundance (mQTL) and phenotype (phQTL) data. Popular QTL mapping methods for model organism and human populations are accessible via the web user interface. Large calculations scale easily on to multi-core computers, clusters and Cloud. All data involved can be uploaded and queried online: markers, genotypes, microarrays, NGS, LC-MS, GC-MS, NMR, etc. When new data types come available, xQTL workbench is quickly customized using the Molgenis software generator., Availability: xQTL workbench runs on all common platforms, including Linux, Mac OS X and Windows. An online demo system, installation guide, tutorials, software and source code are available under the LGPL3 license from http://www.xqtl.org., Contact: m.a.swertz@rug.nl.

During a meeting of the SYSGENET working group 'Bioinformatics', currently available software tools and databases for systems genetics in mice were reviewed and the needs for future developments discussed. The group evaluated interoperability and performed initial feasibility studies. To aid future compatibility of software and exchange of already developed software modules, a strong recommendation was made by the group to integrate HAPPY and R/qtl analysis toolboxes, GeneNetwork and XGAP database platforms, and TIQS and xQTL processing platforms. R should be used as the principal computer language for QTL data analysis in all platforms and a 'cloud' should be used for software dissemination to the community. Furthermore, the working group recommended that all data models and software source code should be made visible in public repositories to allow a coordinated effort on the use of common data structures and file formats.

Motivation: Sample mix-ups can arise during sample collection, handling, genotyping or data management. It is unclear how often sample mix-ups occur in genome-wide studies, as there currently are no post hoc methods that can identify these mix-ups in unrelated samples. We have therefore developed an algorithm (MixupMapper) that can both detect and correct sample mix-ups in genome-wide studies that study gene expression levels., Results: We applied MixupMapper to five publicly available human genetical genomics datasets. On average, 3% of all analyzed samples had been assigned incorrect expression phenotypes: in one of the datasets 23% of the samples had incorrect expression phenotypes. The consequences of sample mix-ups are substantial: when we corrected these sample mix-ups, we identified on average 15% more significant cis-expression quantitative trait loci (cis-eQTLs). In one dataset, we identified three times as many significant cis-eQTLs after correction. Furthermore, we show through simulations that sample mix-ups can lead to an underestimation of the explained heritability of complex traits in genome-wide association datasets., Availability and Implementation: MixupMapper is freely available at http://www.genenetwork.nl/mixupmapper/

We extensively phenotyped 6000 Arabidopsis plants with experimentally perturbed DNA methylomes as well as a diverse panel of natural accessions in a common garden. We found that alterations in DNA methylation not only caused heritable phenotypic diversity but also produced heritability patterns closely resembling those of the natural accessions. Our findings indicate that epigenetically induced and naturally occurring variation in complex traits share part of their polygenic architecture and may offer complementary adaptation routes in ecological settings.

Genes can be silenced with short-interfering RNA molecules (siRNA). siRNAs are widely used to identify gene functions and have high potential for therapeutic treatments. It is critical that the siRNA specifically targets the expression of the gene of interest but has no off-target effects on other genes. Although siRNAs were initially considered to be exclusively active on mature mRNAs in the cytoplasm, additional studies have shown that siRNAs are present in the nucleus as well, suggesting that pre-mRNA sequences containing introns and other untranslated regions can also be targeted. In this study, we investigated the extent to which off-targets may occur in Drosophila melanogaster by looking at mature mRNA sequences and pre-mature RNA sequences separately. First, an in silico approach revealed that, based on sequence similarity, numerous off-targets are predicted to occur in RNAi experiments. Second, existing microarray data were used to investigate a possible effect of the predicted off-targets based on analysis of in vitro data. We found that the occurrence of off-targets in both mature and pre-mature RNA sequences in RNAi experiments can be extensive and significant. Possibilities are discussed how to minimize off-target effects.

Motivation: R/qtl is free and powerful software for mapping and exploring quantitative trait loci (QTL). R/qtl provides a fully comprehensive range of methods for a wide range of experimental cross types. We recently added multiple QTL mapping (MQM) to R/qtl. MQM adds higher statistical power to detect and disentangle the effects of multiple linked and unlinked QTL compared with many other methods. MQM for R/qtl adds many new features including improved handling of missing data, analysis of 10,000 s of molecular traits, permutation for determining significance thresholds for QTL and QTL hot spots, and visualizations for cis-trans and QTL interaction effects. MQM for R/qtl is the first free and open source implementation of MQM that is multi-platform, scalable and suitable for automated procedures and large genetical genomics datasets., Availability: R/qtl is free and open source multi-platform software for the statistical language R, and is made available under the GPLv3 license. R/qtl can be installed from http://www.rqtl.org/. R/qtl queries should be directed at the mailing list, see http://www.rqtl.org/list/., Contact: kbroman@biostat.wisc.edu.

Alternative splicing is considered a major mechanism for creating multicellular diversity from a limited repertoire of genes. Here, we performed the first study of genetic variation controlling alternative splicing patterns by comprehensively identifying quantitative trait loci affecting the differential expression of transcript isoforms in a large recombinant inbred population of Caenorhabditis elegans, using a new generation of whole-genome very-high-density oligonucleotide microarrays. Using 60 experimental lines, we were able to detect 435 genes with substantial heritable variation, of which 36% were regulated at a distance (in trans). Nonetheless, we find only a very small number of examples of heritable variation in alternative splicing (22 transcripts), and most of these genes colocalize with the associated genomic loci. Our findings suggest that the regulatory mechanism of alternative splicing in C. elegans is robust toward genetic variation at the genome-wide scale, which is in striking contrast to earlier observations in humans.

Motivation: ChIP-chip and ChIP-seq technologies provide genome-wide measurements of various types of chromatin marks at an unprecedented resolution. With ChIP samples collected from different tissue types and/or individuals, we can now begin to characterize stochastic or systematic changes in epigenetic patterns during development (intra-individual) or at the population level (inter-individual). This requires statistical methods that permit a simultaneous comparison of multiple ChIP samples on a global as well as locus-specific scale. Current analytical approaches are mainly geared toward single sample investigations, and therefore have limited applicability in this comparative setting. This shortcoming presents a bottleneck in biological interpretations of multiple sample data., Results: To address this limitation, we introduce a parametric classification approach for the simultaneous analysis of two (or more) ChIP samples. We consider several competing models that reflect alternative biological assumptions about the global distribution of the data. Inferences about locus-specific and genome-wide chromatin differences are reached through the estimation of multivariate mixtures. Parameter estimates are obtained using an incremental version of the Expectation-Maximization algorithm (IEM). We demonstrate efficient scalability and application to three very diverse ChIP-chip and ChIP-seq experiments. The proposed approach is evaluated against several published ChIP-chip and ChIP-seq software packages. We recommend its use as a first-pass algorithm to identify candidate regions in the epigenome, possibly followed by some type of second-pass algorithm to fine-tune detected peaks in accordance with biological or technological criteria., Availability: R source code is available at http://gbic.biol.rug.nl/supplementary/2009/ChromatinProfiles/. Access to Chip-seq data: GEO repository GSE17937.

Combining global gene-expression profiling and genetic analysis of natural allelic variation (genetical genomics) has great potential in dissecting the genetic pathways underlying complex phenotypes. Efficient use of microarrays is paramount in experimental design as the cost of conducting this type of study is high. For those organisms where recombinant inbred lines are available for mapping, the "distant pair design" maximizes the number of informative contrasts over all marker loci. Here, we describe an extension of this design, named the "optimal pair design," for use with F2 crosses between outbred lines. The performance of this design is investigated by simulation and compared to several other two-color microarray designs. We show that, for a given number of microarrays, the optimal pair design outperforms all other designs considered for detection of expression quantitative trait loci (eQTL) with additive effects by linkage analysis. We also discuss the suitability of this design for outbred crosses in organisms with large genomes and for detection of dominance.

Motivation: Affymetrix arrays use multiple probes per gene to measure mRNA abundances. Standard software takes averages over probes. Important information may be lost if polymorphisms in the mRNA affect the hybridization of individual probes., Results: We present custom software to analyze genetical genomics experiments in human, mouse and other organisms: (i) an R package providing functions for QTL analysis at the individual probe level and (ii) Perl scripts providing custom tracks in the UCSC Genome Browser to check for sequence polymorphisms in probe regions., Availability: http://gbic.biol.rug.nl/supplementary.

Whole-genome profiling of gene expression in a segregating population has the potential to identify the regulatory consequences of natural allelic variation. Costs of such studies are high and require that resources--microarrays and population--are used as efficiently as possible. We show that current studies can be improved significantly by a new design for two-color microarrays. Our "distant pair design" profiles twice as many individuals as there are arrays, cohybridizes individuals with dissimilar genomes, gives more weight to known regulatory loci if wished, and therewith maximizes the power for decomposing expression variation into regulatory factors. It can also exploit a large population (larger than twice the number of available microarrays) as a useful resource to select the most dissimilar pairs of individuals from. Our approach identifies more regulatory factors than alternative strategies do in computer simulations for realistic genome sizes, and similar promising results are obtained in an application on Arabidopsis thaliana. Our results will aid the design and analysis of future studies on gene expression and will help to shed more light on gene regulatory networks.

Short-oligonucleotide arrays typically contain multiple probes per gene. In genetical genomics applications a statistical model for the individual probe signals can help in separating "true" differential mRNA expression from "ghost" effects caused by polymorphisms, misdesigned probes, and batch effects. It can also help in detecting alternative splicing, start, or termination.

Gene expression can be studied at a genome-wide scale with the aid of modern microarray technologies. Expression profiling of tens to hundreds of individuals in a genetic population can reveal the consequences of genetic variation. In this paper it is argued that the design and analysis of such a study is not a matter of simply applying the existing and more-or-less standard computational tools for microarrays to a new type of experimental data. It is shown how to fully exploit the power of genetics through optimal experimental design and analysis for two major microarray technologies, cDNA two-colour arrays and Affymetrix short oligonucleotide arrays.

Motivation: Genomic research laboratories need adequate infrastructure to support management of their data production and research workflow. But what makes infrastructure adequate? A lack of appropriate criteria makes any decision on buying or developing a system difficult. Here, we report on the decision process for the case of a molecular genetics group establishing a microarray laboratory., Results: Five typical requirements for experimental genomics database systems were identified: (i) evolution ability to keep up with the fast developing genomics field; (ii) a suitable data model to deal with local diversity; (iii) suitable storage of data files in the system; (iv) easy exchange with other software; and (v) low maintenance costs. The computer scientists and the researchers of the local microarray laboratory considered alternative solutions for these five requirements and chose the following options: (i) use of automatic code generation; (ii) a customized data model based on standards; (iii) storage of datasets as black boxes instead of decomposing them in database tables; (iv) loosely linking to other programs for improved flexibility; and (v) a low-maintenance web-based user interface. Our team evaluated existing microarray databases and then decided to build a new system, Molecular Genetics Information System (MOLGENIS), implemented using code generation in a period of three months. This case can provide valuable insights and lessons to both software developers and a user community embarking on large-scale genomic projects., Availability: http://www.molgenis.nl

Epistasis is a common and important phenomenon, as indicated by results from a number of recent experiments. Unfortunately, the discovery of epistatic quantitative trait loci (QTL) is difficult since one must search for multiple QTL simultaneously in two or more dimensions. Such a multidimensional search necessitates many statistical tests, and a high statistical threshold must be adopted to avoid false positives. Furthermore, the large number of (interaction) parameters in comparison with the number of observations results in a serious danger of overfitting and overinterpretation of the data. In this article we present a new statistical framework for mapping epistasis in inbred line crosses. It is based on reducing the high dimensionality of the problem in two ways. First, epistatic QTL are mapped in a one-dimensional genome scan for high interactions between QTL and the genetic background. Second, the dimension of the search is bounded by penalized likelihood methods. We use simulated backcross data to illustrate the new approach.

The chromatin loop model predicts that genes within the same chromatin domain exhibit coordinated regulation. We here present the first direct experimental support for this model in plants. Two reporter genes, the E. coli beta-glucuronidase gene and the firefly luciferase gene, driven by different promoters, were placed between copies of the chicken lysozyme A element, a member of the matrix-associated region (MAR) group of chromatin boundary elements, and introduced in tobacco (Nicotiana tabacum). In plants carrying A elements, quantitative enzyme activities and mRNA levels of both genes show high correlations compared to control plants. The A element thus creates an artificial chromatin domain that yields coordinated expression. Surprisingly, enzyme activities correlated poorly with their respective mRNA levels. We hypothesize that this indicates the occurrence of "error pipelines" in data generation: systematic errors of a given analytical method will point in the same direction and cancel out in correlation analysis, resulting in better correlations. In combining different methods of analysis, however, such errors do not cancel out and as a result relevant correlations can be masked. Such error pipelines will have to be taken into account when different types of (e.g., whole-genome) data sets are combined in quantitative analyses.

Fruit flavor is a result of a complex mixture of numerous compounds. The formation of these compounds is closely correlated with the metabolic changes occurring during fruit maturation. Here, we describe the use of DNA microarrays and appropriate statistical analyses to dissect a complex developmental process. In doing so, we have identified a novel strawberry alcohol acyltransferase (SAAT) gene that plays a crucial role in flavor biogenesis in ripening fruit. Volatile esters are quantitatively and qualitatively the most important compounds providing fruity odors. Biochemical evidence for involvement of the SAAT gene in formation of fruity esters is provided by characterizing the recombinant protein expressed in Escherichia coli. The SAAT enzyme showed maximum activity with aliphatic medium-chain alcohols, whose corresponding esters are major components of strawberry volatiles. The enzyme was capable of utilizing short- and medium-chain, branched, and aromatic acyl-CoA molecules as cosubstrates. The results suggest that the formation of volatile esters in fruit is subject to the availability of acyl-CoA molecules and alcohol substrates and is dictated by the temporal expression pattern of the SAAT gene(s) and substrate specificity of the SAAT enzyme(s).

A mixture model approach is presented for the mapping of one or more quantitative trait loci (QTLs) in complex populations. In order to exploit the full power of complete linkage maps the simultaneous likelihood of phenotype and a multilocus (all markers and putative QTLs) genotype is computed. Maximum likelihood estimation in our mixture models is implemented via an Expectation-Maximization algorithm: exact, stochastic or Monte Carlo EM by using a simple and flexible Gibbs sampler. Parameters include allele frequencies of markers and QTLs, discrete or normal effects of biallelic or multiallelic QTLs, and homogeneous or heterogeneous residual variances. As an illustration a dairy cattle data set consisting of twenty half-sib families has been reanalyzed. We discuss the potential which our and other approaches have for realistic multiple-QTL analyses in complex populations.

Six transgenic tobacco lines, each homozygous for the beta-glucuronidase (GUS) gene at a different locus, and wild type were selfed and intercrossed to evaluate GUS activity in all possible hemizygous, homozygous and dihybrid combinations of GUS alleles. The transgenic lines are characterized by their GUS activity (two low, three intermediate, one high), T-DNA complexity (four single-copy, two more complex single-locus) and the presence of the chicken lysozyme matrix-associated region (MAR) around the full T-DNA (two lines). Gene action and interaction was analyzed by weighted linear regression with parameters for additivity, dominance and epistasis. The analysis showed that each of the four single-copy lines acted fully additively. In contrast, the two complex single-locus lines showed classical single-locus overdominance and were epistatic dominant over all other GUS alleles. The latter is manifested in severe suppression of GUS activity in dihybrid lines, irrespective of the presence of MAR elements around the GUS gene. Such elements apparently do not protect against epistatic dominance. The quantitative data suggested that the epistatic dominance and overdominance are based on the same molecular mechanism. Our approach of a genetic analysis of quantitative variation in well-characterized transgenic lines provides a powerful tool to gain insight into complex plant traits.

In this paper we address the mapping of multiple quantitative trait loci (QTLs) in line crosses for which the genetic data are highly incomplete. Such complicated situations occur, for instance, when dominant markers are used or when unequally informative markers are used in experiments with outbred populations. We describe a general and flexible Monte Carlo expectation-maximization (Monte Carlo EM) algorithm for fitting multiple-QTL models to such data. Implementation of this algorithm is straightforward in standard statistical software, but computation may take much time. The method may be generalized to cope with more complex models for animal and human pedigrees. A practical example is presented, where a three-QTL model is adopted in an outbreeding situation with dominant markers. The example is concerned with the linkage between randomly amplified polymorphic DNA (RAPD) markers and QTLs for partial resistance to Fusarium oxysporum in lily.

To study the role of matrix-associated regions (MARs) in establishing independent chromatin domains in plants, two transgenes were cloned between chicken lysozyme A elements. These transgenes were the neomycin phosphotransferase (NPTII) gene under control of the nopaline synthase (nos) promoter and the [beta]-glucuronidase (GUS) gene controlled by the double cauliflower mosaic virus (dCaMV) 35S promoter. The A elements are supposed to establish an artificial chromatin domain upon integration into the plant DNA, resulting in an independent unit of transcriptional regulation. Such a domain is thought to be characterized by a correlated and position-independent, hence copy number-dependent, expression of the genes within the domain. The presence of MARs resulted in a higher relative transformation efficiency, demonstrating MAR influence on NPTII gene expression. However, variation in NPTII gene expression was not significantly reduced. The selection bias for NPTII gene expression during transformation could not fully account for the lack of reduction in variation of NPTII gene expression. Topological interactions between the promoter and A element may interfere with the A element as a domain boundary. In contrast, the GUS gene on the same putative chromatin domain showed a highly significant reduction in variation of gene expression, as expected from previous results. Surprisingly, no copy number-dependent GUS gene expression was observed: all plants showed approximately the same GUS activity. We concluded, therefore, that dCaMV 35S-GUS gene expression in mature tobacco plants is regulated by some form of dosage compensation.

Although the interval mapping method is widely used for mapping quantitative trait loci (QTLs), it is not very well suited for mapping multiple QTLs. Here, we present the results of a computer simulation to study the application of exact and approximate models for multiple QTLs. In particular, we focus on an automatic two-stage procedure in which in the first stage "important" markers are selected in multiple regression on markers. In the second stage a QTL is moved along the chromosomes by using the preselected markers as cofactors, except for the markers flanking the interval under study. A refined procedure for cases with large numbers of marker cofactors is described. Our approach will be called MQM mapping, where MQM is an acronym for "multiple-QTL models" as well as for "marker-QTL-marker." Our simulation work demonstrates the great advantage of MQM mapping compared to interval mapping in reducing the chance of a type I error (i.e., a QTL is indicated at a location where actually no QTL is present) and in reducing the chance of a type II error (i.e., a QTL is not detected).

A very general method is described for multiple linear regression of a quantitative phenotype on genotype [putative quantitative trait loci (QTLs) and markers] in segregating generations obtained from line crosses. The method exploits two features, (a) the use of additional parental and F1 data, which fixes the joint QTL effects and the environmental error, and (b) the use of markers as cofactors, which reduces the genetic background noise. As a result, a significant increase of QTL detection power is achieved in comparison with conventional QTL mapping. The core of the method is the completion of any missing genotypic (QTL and marker) observations, which is embedded in a general and simple expectation maximization (EM) algorithm to obtain maximum likelihood estimates of the model parameters. The method is described in detail for the analysis of an F2 generation. Because of the generality of the approach, it is easily applicable to other generations, such as backcross progenies and recombinant inbred lines. An example is presented in which multiple QTLs for plant height in tomato are mapped in an F2 progeny, using additional data from the parents and their F1 progeny.

Matrix-associated regions may be useful for studying the role of chromatin architecture in transgene activity of transformed plants. The chicken lysozyme A element was shown to have specific affinity for tobacco nuclear matrices, and its influence on the variability of transgene expression in tobacco plants was studied. T-DNA constructs in which this element flanked either the [beta]-glucuronidase (GUS) reporter gene or both reporter and selection gene were introduced in tobacco. The variation in GUS gene activity was reduced significantly among mature first-generation transgenic plants carrying the A element. Average GUS activity became somewhat higher, but the maximum attainable level of gene expression was similar for all three constructs. Transient gene expression assays showed that the A element did not contain general enhancer functions; therefore, its presence seemed to prevent the lower levels of transgene expression. The strongest reduction in variability was found in plants transformed with the construct carrying the A elements at the borders of the T-DNA. In this population, expression levels became copy number dependent. The presence of two A elements in the T-DNA did not interfere with meiosis.

The interval mapping method is widely used for the mapping of quantitative trait loci (QTLs) in segregating generations derived from crosses between inbred lines. The efficiency of detecting and the accuracy of mapping multiple QTLs by using genetic markers are much increased by employing multiple QTL models instead of the single QTL models (and no QTL models) used in interval mapping. However, the computational work involved with multiple QTL models is considerable when the number of QTLs is large. In this paper it is proposed to combine multiple linear regression methods with conventional interval mapping. This is achieved by fitting one QTL at a time in a given interval and simultaneously using (part of) the markers as cofactors to eliminate the effects of additional QTLs. It is shown that the proposed method combines the easy computation of the single QTL interval mapping method with much of the efficiency and accuracy of multiple QTL models.