Strategies for performing genotype–phenotype association studies in nonhuman primates

Anthropoid primate models offer opportunities to study genetic influence on alcohol consumption and alcohol-related intermediate phenotypes in socially and behaviorally complex animal models that are closely related to humans, and in which functionally equivalent or orthologous genetic variants exist. This review will discuss the methods commonly used for performing candidate gene-based studies in rhesus macaques in order to model how functional genetic variation moderates risk for human psychiatric disorders. Various in silico and in vitro approaches to identifying functional genetic variants for performance of these studies will be discussed. Next, I will provide examples of how this approach can be used for performing candidate gene-based studies and for examining gene by environment interactions. Finally, these approaches will then be placed in the context of how function-guided studies can inform us of genetic variants that may be under selection across species, demonstrating how functional genetic variants that may have conferred selective advantage at some point in the evolutionary history of humans could increase risk for addictive disorders in modern society.     

The role of phenotype in gene discovery in the whole genome sequencing era

As whole genome sequence becomes a routine component of gene discovery studies in humans, we will have an exhaustive catalog of genetic variation and the challenge becomes understanding the phenotypic consequences of these variants. Statistical genetic methods and analytical approaches that are concerned with optimizing phenotypes for gene discovery for complex traits offer two general categories of advantages. They may increase power to localize genes of interest and also aid in interpreting associations between genetic variants and disease outcomes by suggesting potential mechanisms and pathways through which genes may affect outcomes. Such phenotype optimization approaches include use of allied phenotypes such as symptoms or ages of onset to reduce genetic heterogeneity within a set of cases, study of quantitative risk factors or endophenotypes, joint analyses of related phenotypes, and derivation of new phenotypes designed to extract independent measures underlying the correlations among a set of related phenotypes through approaches such as principal components. New opportunities are also presented by technological advances that permit efficient collection of hundreds or thousands of phenotypes on an individual, including phenotypes more proximal to the level of gene action such as levels of gene expression, microRNAs, or metabolic and proteomic profiles.     

Behavioural syndromes in fishes: a review with implications for ecology and fisheries management

This review examines the contribution of research on fishes to the growing field of behavioural syndromes. Current knowledge of behavioural syndromes in fishes is reviewed with respect to five main axes of animal personality: (1) shyness-boldness, (2) exploration-avoidance, (3) activity, (4) aggressiveness and (5) sociability. Compared with other taxa, research on fishes has played a leading role in describing the shy-bold personality axis and has made innovative contributions to the study of the sociability dimension by incorporating social network theory. Fishes are virtually the only major taxon in which behavioural correlations have been compared between populations. This research has guided the field in examining how variation in selection regime may shape personality. Recent research on fishes has also made important strides in understanding genetic and neuroendocrine bases for behavioural syndromes using approaches involving artificial selection, genetic mapping, candidate gene and functional genomics. This work has illustrated consistent individual variation in highly complex neuroendocrine and gene expression pathways. In contrast, relatively little work on fishes has examined the ontogenetic stability of behavioural syndromes or their fitness consequences. Finally, adopting a behavioural syndrome framework in fisheries management issues including artificial propagation, habitat restoration and invasive species, may promote restoration success. Few studies, however, have examined the ecological relevance of behavioural syndromes in the field. Knowledge of how behavioural syndromes play out in the wild will be crucial to incorporating such a framework into management practices.     

Genomewide analysis of inherited variation associated with phosphorylation of PI3K/AKT/mTOR signaling proteins

While there exists a wealth of information about genetic influences on gene expression, less is known about how inherited variation influences the expression and post-translational modifications of proteins, especially those involved in intracellular signaling. The PI3K/AKT/mTOR signaling pathway contains several such proteins that have been implicated in a number of diseases, including a variety of cancers and some psychiatric disorders. To assess whether the activation of this pathway is influenced by genetic factors, we measured phosphorylated and total levels of three key proteins in the pathway (AKT1, p70S6K, 4E-BP1) by ELISA in 122 lymphoblastoid cell lines from 14 families. Interestingly, the phenotypes with the highest proportion of genetic influence were the ratios of phosphorylated to total protein for two of the pathway members: AKT1 and p70S6K. Genomewide linkage analysis suggested several loci of interest for these phenotypes, including a linkage peak for the AKT1 phenotype that contained the AKT1 gene on chromosome 14. Linkage peaks for the phosphorylated:total protein ratios of AKT1 and p70S6K also overlapped on chromosome 3. We selected and genotyped candidate genes from under the linkage peaks, and several statistically significant associations were found. One polymorphism in HSP90AA1 was associated with the ratio of phosphorylated to total AKT1, and polymorphisms in RAF1 and GRM7 were associated with the ratio of phosphorylated to total p70S6K. These findings, representing the first genomewide search for variants influencing human protein phosphorylation, provide useful information about the PI3K/AKT/mTOR pathway and serve as a valuable proof of concept for studies integrating human genomics and proteomics.     

Genome-wide Introgression Lines and their Use in Genetic and Molecular Dissection of Complex Phenotypes in Rice (Oryza sativa L.)

Tremendous efforts have been taken worldwide to develop genome-wide genetic stocks for rice functional genomic (FG) research since the rice genome was completely sequenced. To facilitate FG research of complex polygenic phenotypes in rice, we report the development of over 20 000 introgression lines (ILs) in three elite rice genetic backgrounds for a wide range of complex traits, including resistances/tolerances to many biotic and abiotic stresses, morpho-agronomic traits, physiological traits, etc., by selective introgression. ILs within each genetic background are phenotypically similar to their recurrent parent but each carries one or a few traits introgressed from a known donor. Together, these ILs contain a significant portion of loci affecting the selected complex phenotypes at which allelic diversity exists in the primary gene pool of rice. A forward genetics strategy was proposed and demonstrated with examples on how to use these ILs for large-scale FG research. Complementary to the genome-wide insertional mutants, these ILs opens a new way for highly efficient discovery, candidate gene identification and cloning of important QTLs for specific phenotypes based on convergent evidence from QTL position, expression profiling, functional and molecular diversity analyses of candidate genes, highlights the importance of genetic networks underlying complex phenotypes in rice that may ultimately lead to more complete understanding of the genetic and molecular bases of quantitative trait variation in rice. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s11103-005-8519-3     

Multiple regression methods show great potential for rare variant association tests

The investigation of associations between rare genetic variants and diseases or phenotypes has two goals. Firstly, the identification of which genes or genomic regions are associated, and secondly, discrimination of associated variants from background noise within each region. Over the last few years, many new methods have been developed which associate genomic regions with phenotypes. However, classical methods for high-dimensional data have received little attention. Here we investigate whether several classical statistical methods for high-dimensional data: ridge regression (RR), principal components regression (PCR), partial least squares regression (PLS), a sparse version of PLS (SPLS), and the LASSO are able to detect associations with rare genetic variants. These approaches have been extensively used in statistics to identify the true associations in data sets containing many predictor variables. Using genetic variants identified in three genes that were Sanger sequenced in 1998 individuals, we simulated continuous phenotypes under several different models, and we show that these feature selection and feature extraction methods can substantially outperform several popular methods for rare variant analysis. Furthermore, these approaches can identify which variants are contributing most to the model fit, and therefore both goals of rare variant analysis can be achieved simultaneously with the use of regression regularization methods. These methods are briefly illustrated with an analysis of adiponectin levels and variants in the ADIPOQ gene.     

Increased 90-kDa ribosomal S6 kinase (Rsk) activity is protective against mutant huntingtin toxicity

Uncovering the underlying genetic component of any disease is key to the understanding of its pathophysiology and may open new avenues for development of therapeutic strategies and biomarkers. In the past several years, there has been an explosion of genome-wide association studies (GWAS) resulting in the discovery of novel candidate genes conferring risk for complex diseases, including neurodegenerative diseases. Despite this success, there still remains a substantial genetic component for many complex traits and conditions that is unexplained by the GWAS findings. Additionally, in many cases, the mechanism of action of the newly discovered disease risk variants is not inherently obvious. Furthermore, a genetic region with multiple genes may be identified via GWAS, making it difficult to discern the true disease risk gene. Several alternative approaches are proposed to overcome these potential shortcomings of GWAS, including the use of quantitative, biologically relevant phenotypes. Gene expression levels represent an important class of endophenotypes. Genetic linkage and association studies that utilize gene expression levels as endophenotypes determined that the expression levels of many genes are under genetic influence. This led to the postulate that there may exist many genetic variants that confer disease risk via modifying gene expression levels. Results from the handful of genetic studies which assess gene expression level endophenotypes in conjunction with disease risk suggest that this combined phenotype approach may both increase the power for gene discovery and lead to an enhanced understanding of their mode of action. This review summarizes the evidence in support of gene expression levels as promising endophenotypes in the discovery and characterization of novel candidate genes for complex diseases, which may also represent a novel approach in the genetic studies of Alzheimer's and other neurodegenerative diseases.     

Animal models of anxiety and their molecular dissection

Genetic approaches to psychiatric illness need appropriate animal models both for investigating how genetic variants give rise to behavioural disorder and for identifying genes that may be important in human conditions. Yet the relevance of many animal models to psychiatric illness is often not clear. Here I discuss how genetic approaches can be used to validate animal models of anxiety, an approach which is applicable to other animal models. One drawback of genetic validation is the difficulty inherent in identifying the molecular variants that influence the phenotype. I review genetic approaches that have the potential to overcome this problem.     

Adaptation to climate through flowering phenology: a case study in Medicago truncatula

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness‐related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome‐wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.     

Gene–environment interactions at the FKBP5 locus: sensitive periods,mechanisms and pleiotropism

Psychiatric phenotypes are multifactorial and polygenic, resulting from the complex interplay of genes and environmental factors that act cumulatively throughout an organism's lifetime. Adverse life events are strong predictors of risk for a number of psychiatric disorders and a number of studies have focused on gene–environment interactions (GxEs) occurring at genetic loci involved in the stress response. Such a locus that has received increasing attention is the gene encoding FK506 binding protein 51 (FKBP5), a heat shock protein 90 cochaperone of the steroid receptor complex that among other functions regulates sensitivity of the glucocorticoid receptor. Interactions between FKBP5 gene variants and life stressors alter the risk not only for mood and anxiety disorders, but also for a number of other disease phenotypes. In this review, we will focus on molecular and system‐wide mechanisms of this GxE with the aim of establishing a framework that explains GxE interactions. We will also discuss how an understanding of the biological effects of this GxE may lead to novel therapeutic approaches .     

Polymorphisms in vitamin D receptor, osteopontin, insulin-like growth factor 1 and insulin, and their associations with bone, egg and growth traits in a layer--broiler cross in chickens

Bone strength traits in chickens are gaining importance due to economic losses and welfare concerns associated with bone fractures and other abnormalities. A chicken F2 resource population was generated from layer and broiler genetic lines, and traits relating to bone strength, egg production, egg quality and growth rate were measured in approximately 500 F2 hens. Four biological candidate genes (vitamin D receptor, VDR; insulin, INS; insulin-like growth factor 1, IGF1; and osteopontin, SPP1) were selected for investigation. Single nucleotide polymorphisms (SNPs) were identified for each candidate gene by comparing sequences between grandparent lines. Polymerase chain reaction restriction-fragment length polymorphism or SNaPshot assays were developed to genotype the F2 population and to evaluate associations between each SNP genotype and multiple phenotypes. Significant associations (P < 0.0125) were found between VDR and bone mineral content of the humerus at 35 weeks of age; between IGF1 and SPP1 and 5-week body weight; and between INS and 55-week body weight.     

Network‐based global inference of human disease genes

Deciphering the genetic basis of human diseases is an important goal of biomedical research. On the basis of the assumption that phenotypically similar diseases are caused by functionally related genes, we propose a computational framework that integrates human protein–protein interactions, disease phenotype similarities, and known gene–phenotype associations to capture the complex relationships between phenotypes and genotypes. We develop a tool named CIPHER to predict and prioritize disease genes, and we show that the global concordance between the human protein network and the phenotype network reliably predicts disease genes. Our method is applicable to genetically uncharacterized phenotypes, effective in the genome‐wide scan of disease genes, and also extendable to explore gene cooperativity in complex diseases. The predicted genetic landscape of over 1000 human phenotypes, which reveals the global modular organization of phenotype–genotype relationships. The genome‐wide prioritization of candidate genes for over 5000 human phenotypes, including those with under‐characterized disease loci or even those lacking known association, is publicly released to facilitate future discovery of disease genes.     

Post‐GWAS in Psychiatric Genetics: A Developmental Perspective on the “Other” Next Steps

As psychiatric genetics enters an era where gene identification is finally yielding robust, replicable genetic associations and polygenic risk scores, it is important to consider next steps and delineate how that knowledge will be applied to ultimately ameliorate suffering associated with substance use and psychiatric disorders. Much of the post‐genome‐wide association study discussion has focused on the potential of genetic information to elucidate the underlying biology and use this information for the development of more effective pharmaceutical treatments. In this review we focus on additional areas of research that should follow gene identification. By taking genetic findings into longitudinal, developmental studies, we can map the pathways by which genetic risk manifests across development, elucidating the early behavioral manifestations of risk, and studying how various environments and interventions moderate that risk across developmental stages. The delineation of risk across development will advance our understanding of mechanism, sex differences and risk and resilience processes in different racial/ethnic groups. Here, we review how the extant twin study literature can be used to guide these efforts. Together, these new lines of research will enable us to develop more informed, tailored prevention and intervention efforts.     

New connections across pathways and cellular processes: industrialized mutant screening reveals novel associations between diverse phenotypes in Arabidopsis

In traditional mutant screening approaches, genetic variants are tested for one or a small number of phenotypes. Once bona fide variants are identified, they are typically subjected to a limited number of secondary phenotypic screens. Although this approach is excellent at finding genes involved in specific biological processes, the lack of wide and systematic interrogation of phenotype limits the ability to detect broader syndromes and connections between genes and phenotypes. It could also prevent detection of the primary phenotype of a mutant. As part of a systems biology approach to understand plastid function, large numbers of Arabidopsis thaliana homozygous T-DNA lines are being screened with parallel morphological, physiological, and chemical phenotypic assays (www.plastid.msu.edu). To refine our approaches and validate the use of this high-throughput screening approach for understanding gene function and functional networks, approximately 100 wild-type plants and 13 known mutants representing a variety of phenotypes were analyzed by a broad range of assays including metabolite profiling, morphological analysis, and chlorophyll fluorescence kinetics. Data analysis using a variety of statistical approaches showed that such industrial approaches can reliably identify plant mutant phenotypes. More significantly, the study uncovered previously unreported phenotypes for these well-characterized mutants and unexpected associations between different physiological processes, demonstrating that this approach has strong advantages over traditional mutant screening approaches. Analysis of wild-type plants revealed hundreds of statistically robust phenotypic correlations, including metabolites that are not known to share direct biosynthetic origins, raising the possibility that these metabolic pathways have closer relationships than is commonly suspected.     

Haplotype structure,adaptive history and associations with exploratory behaviour of the DRD4 gene region in four great tit (Parus major) populations

The assessment of genetic architecture and selection history in genes for behavioural traits is fundamental to our understanding of how these traits evolve. The dopamine receptor D4 (DRD4) gene is a prime candidate for explaining genetic variation in novelty seeking behaviour, a commonly assayed personality trait in animals. Previously, we showed that a single nucleotide polymorphism in exon 3 of this gene is associated with exploratory behaviour in at least one of four Western European great tit (Parus major) populations. These heterogeneous association results were explained by potential variable linkage disequilibrium (LD) patterns between this marker and the causal variant or by other genetic or environmental differences among the populations. Different adaptive histories are further hypothesized to have contributed to these population differences. Here, we genotyped 98 polymorphisms of the complete DRD4 gene including the flanking regions for 595 individuals of the four populations. We show that the LD structure, specifically around the original exon 3 SNP is conserved across the four populations and does not explain the heterogeneous association results. Study‐wide significant associations with exploratory behaviour were detected in more than one haplotype block around exon 2, 3 and 4 in two of the four tested populations with different allele effect models. This indicates genetic heterogeneity in the association between multiple DRD4 polymorphisms and exploratory behaviour across populations. The association signals were in or close to regions with signatures of positive selection. We therefore hypothesize that variation in exploratory and other dopamine‐related behaviour evolves locally by occasional adaptive shifts in the frequency of underlying genetic variants.     

Female influence on pre- and post-copulatory sexual selection and its genetic basis in Drosophila melanogaster

Genetic variation among females is likely to influence the outcome of both pre- and post-copulatory sexual selection in Drosophila melanogaster. Here we use association testing to survey natural variation in 10 candidate female genes for their effects on female reproduction. Females from 91 chromosome two substitution lines were scored for phenotypes affecting pre- and post-copulatory sexual selection such as mating and remating rate, propensity to use sperm from the second male to mate, and measures of fertility. There were significant genetic contributions to phenotypic variation for all the traits measured. Resequencing of the 10 candidate genes in the 91 lines yielded 68 non-synonymous polymorphisms which were tested for associations with the measured phenotypes. Twelve significant associations (markerwise P<0.01) were identified. Polymorphisms in the putative serine protease homolog CG9897 and the putative odorant binding protein CG11797 associated with female propensity to remate and met an experimentwise significance of P<0.05. Several other associations, including those impacting both fertility and female remating rate suggest that sperm storage might be an important factor mitigating female influence on sexual selection.     

On the genetic architecture of rapidly adapting and convergent life history traits in guppies

The genetic basis of traits shapes and constrains how adaptation proceeds in nature; rapid adaptation can proceed using stores of polygenic standing genetic variation or hard selective sweeps, and increasing polygenicity fuels genetic redundancy, reducing gene re-use (genetic convergence). Guppy life history traits evolve rapidly and convergently among natural high- and low-predation environments in northern Trinidad. This system has been studied extensively at the phenotypic level, but little is known about the underlying genetic architecture. Here, we use four independent F2 QTL crosses to examine the genetic basis of seven (five female, two male) guppy life history phenotypes and discuss how these genetic architectures may facilitate or constrain rapid adaptation and convergence. We use RAD-sequencing data (16,539 SNPs) from 370 male and 267 female F2 individuals. We perform linkage mapping, estimates of genome-wide and per-chromosome heritability (multi-locus associations), and QTL mapping (single-locus associations). Our results are consistent with architectures of many loci of small-effect for male age and size at maturity and female interbrood period. Male trait associations are clustered on specific chromosomes, but female interbrood period exhibits a weak genome-wide signal suggesting a potentially highly polygenic component. Offspring weight and female size at maturity are also associated with a single significant QTL each. These results suggest rapid, repeatable phenotypic evolution of guppies may be facilitated by polygenic trait architectures, but subsequent genetic redundancy may limit gene re-use across populations, in agreement with an absence of strong signatures of genetic convergence from recent analyses of wild guppies.Subject terms: Evolutionary genetics, Quantitative trait     

Integrating candidate gene and quantitative genetic approaches to understand variation in timing of breeding in wild tit populations

Two commonly used techniques for estimating the effect of genes on traits in wild populations are the candidate gene approach and quantitative genetic analyses. However, whether these two approaches measure the same underlying processes remains unresolved. Here, we use these two methods to test whether they are alternative or complementary approaches to understanding genetic variation in the timing of reproduction - a key trait involved in adaptation to climate change - in wild tit populations. Our analyses of the candidate gene Clock show weak correlates with timing variables in blue tits, but no association in great tits, confirming earlier results. Quantitative genetic analyses revealed very low levels of both direct (female) and indirect (male) additive genetic variation in timing traits for both species, in contrast to previous studies on these traits, and much lower than generally assumed. Hence, neither method suggests strong genetic effects on the timing of breeding in birds, and further work should seek to assess the generality of these conclusions. We discuss how differences in the genetic control of traits, species life-history and confounding environmental variables may determine how useful integrating these two techniques is to understand the phenotypic variation in wild populations.     

Mouse models for psychiatric disorders

Genes involved in psychiatric disorders are difficult to identify, and those that have been proposed so far remain ambiguous. As it is unrealistic to expect the development of, say, a ‘schizophrenic’ or ‘autistic’ mouse, mice are unlikely to have the same role in gene identification in psychiatry as circling mice did in the discovery of human deafness genes. However, many psychiatric disorders are associated with intermediate phenotypes that can be modeled and studied in mice, including physiological or anatomical brain changes and behavioral traits. Mouse models help to evaluate the effect of a human candidate gene mutation on an intermediate trait, and to identify new candidate genes. Once a gene or pathway has been identified, mice are also used to study the interplay of different genes in that system.