Integrating genetic and network analysis to characterize genes related to mouse weight

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

Identification of susceptibility loci for complex diseases in a case-control association study using the Genetic Analysis Workshop 14 dataset

Although current methods in genetic epidemiology have been extremely successful in identifying genetic loci responsible for Mendelian traits, most common diseases do not follow simple Mendelian modes of inheritance. It is important to consider how our current methodologies function in the realm of complex diseases. The aim of this study was to determine the ability of conventional association methods to fine map a locus of interest. Six study populations were selected from 10 replicates (New York) from the Genetic Analysis Workshop 14 simulated dataset and analyzed for association between the disease trait and locus D2. Genotypes from 45 single-nucleotide polymorphisms in the telomeric region of chromosome 3 were analyzed by Pearson's chi-square tests for independence to test for association with the disease trait of interest. A significant association was detected within the region; however, it was found 3 cM from the documented location of the D2 disease locus. This result was most likely due to the method used for data simulation. In general, this study showed that conventional case-control association methods could detect disease loci responsible for the development of complex traits.     

Identification of gene expression profiles and key genes in subchondral bone of osteoarthritis using weighted gene coexpression network analysis

Osteoarthritis (OA) significantly influences the quality life of people around the world. It is urgent to find an effective way to understand the genetic etiology of OA. We used weighted gene coexpression network analysis (WGCNA) to explore the key genes involved in the subchondral bone pathological process of OA. Fifty gene expression profiles of GSE51588 were downloaded from the Gene Expression Omnibus database. The OA‐associated genes and gene ontologies were acquired from JuniorDoc. Weighted gene coexpression network analysis was used to find disease‐related networks based on 21756 gene expression correlation coefficients, hub‐genes with the highest connectivity in each module were selected, and the correlation between module eigengene and clinical traits was calculated. The genes in the traits‐related gene coexpression modules were subject to functional annotation and pathway enrichment analysis using ClusterProfiler. A total of 73 gene modules were identified, of which, 12 modules were found with high connectivity with clinical traits. Five modules were found with enriched OA‐associated genes. Moreover, 310 OA‐associated genes were found, and 34 of them were among hub‐genes in each module. Consequently, enrichment results indicated some key metabolic pathways, such as extracellular matrix (ECM)‐receptor interaction (hsa04512), focal adhesion (hsa04510), the phosphatidylinositol 3'‐kinase (PI3K)‐Akt signaling pathway (PI3K‐AKT) (hsa04151), transforming growth factor beta pathway, and Wnt pathway. We intended to identify some core genes, collagen (COL)6A3, COL6A1, ITGA11, BAMBI, and HCK, which could influence downstream signaling pathways once they were activated. In this study, we identified important genes within key coexpression modules, which associate with a pathological process of subchondral bone in OA. Functional analysis results could provide important information to understand the mechanism of OA.     

Weighted Interaction SNP Hub (WISH) network method for building genetic networks for complex diseases and traits using whole genome genotype data

Background

High-throughput genotype (HTG) data has been used primarily in genome-wide association (GWA) studies; however, GWA results explain only a limited part of the complete genetic variation of traits. In systems genetics, network approaches have been shown to be able to identify pathways and their underlying causal genes to unravel the biological and genetic background of complex diseases and traits, e.g., the Weighted Gene Co-expression Network Analysis (WGCNA) method based on microarray gene expression data. The main objective of this study was to develop a scale-free weighted genetic interaction network method using whole genome HTG data in order to detect biologically relevant pathways and potential genetic biomarkers for complex diseases and traits.

Results

We developed the Weighted Interaction SNP Hub (WISH) network method that uses HTG data to detect genome-wide interactions between single nucleotide polymorphism (SNPs) and its relationship with complex traits. Data dimensionality reduction was achieved by selecting SNPs based on its: 1) degree of genome-wide significance and 2) degree of genetic variation in a population. Network construction was based on pairwise Pearson's correlation between SNP genotypes or the epistatic interaction effect between SNP pairs. To identify modules the Topological Overlap Measure (TOM) was calculated, reflecting the degree of overlap in shared neighbours between SNP pairs. Modules, clusters of highly interconnected SNPs, were defined using a tree-cutting algorithm on the SNP dendrogram created from the dissimilarity TOM (1-TOM). Modules were selected for functional annotation based on their association with the trait of interest, defined by the Genome-wide Module Association Test (GMAT). We successfully tested the established WISH network method using simulated and real SNP interaction data and GWA study results for carcass weight in a pig resource population; this resulted in detecting modules and key functional and biological pathways related to carcass weight.

Conclusions

We developed the WISH network method which is a novel 'systems genetics' approach to study genetic networks underlying complex trait variation. The WISH network method reduces data dimensionality and statistical complexity in associating genotypes with phenotypes in GWA studies and enables researchers to identify biologically relevant pathways and potential genetic biomarkers for any complex trait of interest.

     

Phenotypic Integration Without Modularity: Testing Hypotheses About the Distribution of Pleiotropic Quantitative Trait Loci in a Continuous Space

Theories of phenotypic integration have relied heavily on the concept of modularity in order to model the ways in which traits in an organism correlate and covary. Recent investigations suggest that, while some functional and developmental processes may be morphologically and ontogenetically localized, and thus modular in a developmental sense, there is a great deal of overlap among these influences on patterns of integration in the adult form. This can result in blurry boundaries between hypothesized modules constructed to test hypotheses about phenotypic integration. This investigation tests hypotheses about the contribution of pleiotropic quantitative trait loci (QTL) to phenotypic integration in the mouse mandible without using a priori categorical hypotheses about which traits constitute a module. We ask two main questions: (1) Are the effects of pleiotropic QTL localized to highly correlated traits or more spread out among traits than one might expect by chance? (2) Does the pattern of trait influence when all pleiotropic QTL are considered together deviate from what we might expect if QTL affect traits without regard for the correlations among traits? We find that a large proportion of pleiotropic QTL affect traits that are more highly correlated than we expect by chance with the remainder having effects that are distributed as if by chance. Furthermore, the overall distribution of the effects of pleiotropic QTL differs significantly from the null distribution of no association between pleiotropic effects on traits and correlations among traits. The main modular hypothesis used by earlier studies often does not predict the distribution of sets of traits sharing a common QTL. These results suggest that there is a clear tendency for pleiotropic effects of QTL to be localized but that the localization may be best thought of as occurring in a continuous space rather being clustered in discrete modules.     

Using the expectation or the distribution of the identity by descent for mapping quantitative trait loci under the random model.

We examine the ability of four implementations of the random model to map quantitative trait loci (QTLs). The implementations use either the expectation or the distribution of the identity-by-descent value at a putative QTL and either a 2 x 1 vector of sib-pair traits or their scalar difference. When the traits of both sibs are used, there is little difference between the expectation and distribution methods, while the expectation method suffers in both precision and power when the difference between traits is used. This is consistent with the prediction that the difference between the expectation and distribution methods is inversely proportional to the amount of information available for mapping. We find, though, that the amount of information must be very low for this difference to be noticeable. This is exemplified when both marker loci are fixed. In this case, while the expectation method is powerless to detect the QTL, the distribution method can still detect the presence (but not the position) of the QTL 59% of the time (when using trait values) or 14% of the time (when using trait differences). We also note a confounding between estimates of the QTL, polygenic, and error variance. The degree of confounding is small when the vector of trait values is used but can be substantial when the expectation method and trait differences are used. We discuss this in light of the general ability of the random model to partition these components.     

Gene Tagging with Random Amplified Polymorphic DNA (RAPD) Markers for Molecular Breeding in Plants

Markers are of interest to plant breeders as a source of genetic information on crops and for use in indirect selection of traits to which the markers are linked. In the classic breeding approach, the markers were invariably the visible morphological and other phenotypic characters, and the breeders expended considerable effort and time in refining the crosses as the tight linkage or association of the desired characters with the obvious phenotypic characters was never unequivocally established. Furthermore, indirect selection for a trait using such morphological markers was not practical due to (1) a paucity of suitable markers, (2) the undesirable pleiotropic effects of many morphological markers on plant phenotype, and (3) the inability to score multiple morphological mutant traits in a single segregating population. With the advancement in molecular biology, the use of molecular markers in plant breeding has become very commonplace and has given rise to “molecular breeding”. Molecular breeding involves primarily “gene tagging”, followed by “marker-assisted selection” of desired genes or genomes. Gene tagging refers to the identification of existing DNA or the introduction of new DNA that can function as a tag or label for the gene of interest. In order for the DNA sequences to be conserved as a tag, important prerequisites exist. This review also summarizes the achievements in gene tagging that have been made over the last 7 to 8 years.     

Regression association analysis of fruit traits with molecular markers in cherries

The use of molecular markers supports the study of genetic marker–trait association of biological and agronomic interest in diverse genetic material. In this research, association between simple sequence repeat (SSR) and random amplified polymorphic DNA (RAPD) markers with fruit traits were investigated in two collections of cherries by applying multiple regression analysis (MRA). Thirty-eight SSR alleles and 135 RAPD fragments were found associated with 14 of affecting fruit traits. Some of SSR and RAPD markers were associated with more than one fruit trait in MRA. Such an association may arise due to pleiotropic effect of the linked quantitative trait locus on different traits. For example, some SSR and RAPD markers were associated with all four traits including fruit cracking, fruit firmness, total soluble solid (TSS) and fruit shape. Also, some markers had correlations with all four characters of TSS, anthocyanin, fruit skin color and fruit flesh color, indicating a significant correlation among these traits. Therefore, it is possible to use these markers along with morphological traits in cherry breeding programs for identification of suitable parents to produce mapping populations and hybrid cultivars. Also, these results could be useful in marker-assisted breeding programs when no other genetic information is available.     

The Influence of Modularity on Cranial Morphological Disparity in Carnivora and Primates (Mammalia)

Background

Although variation provides the raw material for natural selection and evolution, few empirical data exist about the factors controlling morphological variation. Because developmental constraints on variation are expected to act by influencing trait correlations, studies of modularity offer promising approaches that quantify and summarize patterns of trait relationships. Modules, highly-correlated and semi-autonomous sets of traits, are observed at many levels of biological organization, from genes to colonies. The evolutionary significance of modularity is considerable, with potential effects including constraining the variation of individual traits, circumventing pleiotropy and canalization, and facilitating the transformation of functional structures. Despite these important consequences, there has been little empirical study of how modularity influences morphological evolution on a macroevolutionary scale. Here, we conduct the first morphometric analysis of modularity and disparity in two clades of placental mammals, Primates and Carnivora, and test if trait integration within modules constrains or facilitates morphological evolution.

Principal Findings

We used both randomization methods and direct comparisons of landmark variance to compare disparity in the six cranial modules identified in previous studies. The cranial base, a highly-integrated module, showed significantly low disparity in Primates and low landmark variance in both Primates and Carnivora. The vault, zygomatic-pterygoid and orbit modules, characterized by low trait integration, displayed significantly high disparity within Carnivora. 14 of 24 results from analyses of disparity show no significant relationship between module integration and morphological disparity. Of the ten significant or marginally significant results, eight support the hypothesis that integration within modules constrains morphological evolution in the placental skull. Only the molar module, a highly-integrated and functionally important module, showed significantly high disparity in Carnivora, in support of the facilitation hypothesis.

Conclusions

This analysis of within-module disparity suggested that strong integration of traits had little influence on morphological evolution over large time scales. However, where significant results were found, the primary effect of strong integration of traits was to constrain morphological variation. Thus, within Primates and Carnivora, there was some support for the hypothesis that integration of traits within cranial modules limits morphological evolution, presumably by limiting the variation of individual traits.     

Detection of genes for ordinal traits in nuclear families and a unified approach for association studies

There is growing interest in genomewide association analysis using single-nucleotide polymorphisms (SNPs), because traditional linkage studies are not as powerful in identifying genes for common, complex diseases. Tests for linkage disequilibrium have been developed for binary and quantitative traits. However, since many human conditions and diseases are measured in an ordinal scale, methods need to be developed to investigate the association of genes and ordinal traits. Thus, in the current report we propose and derive a score test statistic that identifies genes that are associated with ordinal traits when gametic disequilibrium between a marker and trait loci exists. Through simulation, the performance of this new test is examined for both ordinal traits and quantitative traits. The proposed statistic not only accommodates and is more powerful for ordinal traits, but also has similar power to that of existing tests when the trait is quantitative. Therefore, our proposed statistic has the potential to serve as a unified approach to identifying genes that are associated with any trait, regardless of how the trait is measured. We further demonstrated the advantage of our test by revealing a significant association (P = 0.00067) between alcohol dependence and a SNP in the growth-associated protein 43.     

Gene-SCOUT: identifying genes with similar continuous trait fingerprints from phenome-wide association analyses

Large-scale phenome-wide association studies performed using densely-phenotyped cohorts such as the UK Biobank (UKB), reveal many statistically robust gene-phenotype relationships for both clinical and continuous traits. Here, we present Gene-SCOUT, a tool used to identify genes with similar continuous trait fingerprints to a gene of interest. A fingerprint reflects the continuous traits identified to be statistically associated with a gene of interest based on multiple underlying rare variant genetic architectures. Similarities between genes are evaluated by the cosine similarity measure, to capture concordant effect directionality, elucidating clusters of genes in a high dimensional space. The underlying gene-biomarker population-scale association statistics were obtained from a gene-level rare variant collapsing analysis performed on over 1500 continuous traits using 394 692 UKB participant exomes, with additional metabolomic trait associations provided through Nightingale Health''s recent study of 121 394 of these participants. We demonstrate that gene similarity estimates from Gene-SCOUT provide stronger enrichments for clinical traits compared to existing methods. Furthermore, we provide a fully interactive web-resource (http://genescout.public.cgr.astrazeneca.com) to explore the pre-calculated exome-wide similarities. This resource enables a user to examine the biological relevance of the most similar genes for Gene Ontology (GO) enrichment and UKB clinical trait enrichment statistics, as well as a detailed breakdown of the traits underpinning a given fingerprint.     

The "eyespot module" and eyespots as modules: development, evolution, and integration of a complex phenotype

Organisms are inherently modular, yet modules also evolve in response to selection for functional integration or functional specialization of traits. For serially repeated homologous traits, there is a clear expectation that selection on the function of individual traits will reduce the integration between traits and subdivide a single ancestral module. The eyespots on butterfly wings are one example of serially repeated morphological traits that share a common developmental mechanism but are subject to natural and sexual selection for divergent functions. Here, I test two hypotheses about the organization of the eyespot pattern into independent dorsal-ventral and anterior-posterior modules, using a graphical modeling technique to examine patterns of eyespot covariation among and within wing surfaces in the butterfly Bicyclus anynana. Although there is a hierarchical and complex pattern of integration among eyespots, the results show a surprising mismatch between patterns of eyespot integration and the developmental and evolutionary eyespot units identified in previous empirical studies. These results are discussed in light of the relationships between developmental, functional, and evolutionary modules, and they suggest that developmental sources of independent trait variation are often masked by developmental sources of trait integration.     

The candidate gene approach in plant genetics: a review

The candidate gene (CG) approach has been applied in plant genetics in the past decade for the characterisation and cloning of Mendelian and quantitative trait loci (QTLs). It constitutes a complementary strategy to map-based cloning and insertional mutagenesis. The goal of this paper is to present an overview of CG analyses in plant genetics. CG analysis is based on the hypothesis that known-function genes (the candidate genes) could correspond to loci controlling traits of interest. CGs refer either to cloned genes presumed to affect a given trait (`functional CGs') or to genes suggested by their close proximity on linkage maps to loci controlling the trait (`positional CGs'). In plant genetics, the most common way to identify a CG is to look for map co-segregation between CGs and loci affecting the trait. Statistical association analyses between molecular polymorphisms of the CG and variation in the trait of interest have also been carried out in a few studies. The final validation of a CG will be provided through physiological analyses, genetic transformation and/or sexual complementation. Theoretical and practical applications of validated CGs in plant genetics and breeding are discussed.     

EMMLi: A maximum likelihood approach to the analysis of modularity

Identification of phenotypic modules, semiautonomous sets of highly correlated traits, can be accomplished through exploratory (e.g., cluster analysis) or confirmatory approaches (e.g., RV coefficient analysis). Although statistically more robust, confirmatory approaches are generally unable to compare across different model structures. For example, RV coefficient analysis finds support for both two‐ and six‐module models for the therian mammalian skull. Here, we present a maximum likelihood approach that takes into account model parameterization. We compare model log‐likelihoods of trait correlation matrices using the finite‐sample corrected Akaike Information Criterion, allowing for comparison of hypotheses across different model structures. Simulations varying model complexity and within‐ and between‐module contrast demonstrate that this method correctly identifies model structure and parameters across a wide range of conditions. We further analyzed a dataset of 3‐D data, consisting of 61 landmarks from 181 macaque (Macaca fuscata) skulls, distributed among five age categories, testing 31 models, including no modularity among the landmarks and various partitions of two, three, six, and eight modules. Our results clearly support a complex six‐module model, with separate within‐ and intermodule correlations. Furthermore, this model was selected for all five age categories, demonstrating that this complex pattern of integration in the macaque skull appears early and is highly conserved throughout postnatal ontogeny. Subsampling analyses demonstrate that this method is robust to relatively low sample sizes, as is commonly encountered in rare or extinct taxa. This new approach allows for the direct comparison of models with different parameterizations, providing an important tool for the analysis of modularity across diverse systems.     

Population genetics of genomics-based crop improvement methods

Many genome-wide association studies (GWAS) in humans are concluding that, even with very large sample sizes and high marker densities, most of the genetic basis of complex traits may remain unexplained. At the same time, recent research in plant GWAS is showing much greater success with fewer resources. Both GWAS and genomic selection (GS), a method for predicting phenotypes by the use of genome-wide marker data, are receiving considerable attention among plant breeders. In this review we explore how differences in population genetic histories, as well as past selection for traits of interest, have produced trait architectures and patterns of linkage disequilibrium (LD) that frequently differ dramatically between domesticated plants and humans, making detection of quantitative trait loci (QTL) effects in crops more rewarding and less costly than in humans.     

A Bayesian Partition Method for Detecting Pleiotropic and Epistatic eQTL Modules

Studies of the relationship between DNA variation and gene expression variation, often referred to as “expression quantitative trait loci (eQTL) mapping”, have been conducted in many species and resulted in many significant findings. Because of the large number of genes and genetic markers in such analyses, it is extremely challenging to discover how a small number of eQTLs interact with each other to affect mRNA expression levels for a set of co-regulated genes. We present a Bayesian method to facilitate the task, in which co-expressed genes mapped to a common set of markers are treated as a module characterized by latent indicator variables. A Markov chain Monte Carlo algorithm is designed to search simultaneously for the module genes and their linked markers. We show by simulations that this method is more powerful for detecting true eQTLs and their target genes than traditional QTL mapping methods. We applied the procedure to a data set consisting of gene expression and genotypes for 112 segregants of S. cerevisiae. Our method identified modules containing genes mapped to previously reported eQTL hot spots, and dissected these large eQTL hot spots into several modules corresponding to possibly different biological functions or primary and secondary responses to regulatory perturbations. In addition, we identified nine modules associated with pairs of eQTLs, of which two have been previously reported. We demonstrated that one of the novel modules containing many daughter-cell expressed genes is regulated by AMN1 and BPH1. In conclusion, the Bayesian partition method which simultaneously considers all traits and all markers is more powerful for detecting both pleiotropic and epistatic effects based on both simulated and empirical data.