Phenotype mining is a novel approach for elucidating the genetic basis of complex phenotypic variation. It involves a search of rich phenotype databases for measures correlated with genetic variation, as identified in genome-wide genotyping or sequencing studies. An initial implementation of phenotype mining in a prospective unselected population cohort, the Northern Finland 1966 Birth Cohort (NFBC1966), identifies neurodevelopment-related traits—intellectual deficits, poor school performance and hearing abnormalities—which are more frequent among individuals with large (>500 kb) deletions than among other cohort members. Observation of extensive shared single nucleotide polymorphism haplotypes around deletions suggests an opportunity to expand phenotype mining from cohort samples to the populations from which they derive.

Vervet monkeys are a frequently studied animal model in neuroscience research. Although equally distantly related to humans, the ancestors of vervets diverged from those of macaques and baboons more than eleven million years ago, antedating the divergence of the ancestors of humans, chimpanzees and gorillas. To facilitate anatomic localization in the vervet brain, two linked on-line electronic atlases are described, one based on registered MRI scans from hundreds of vervets (http://www.loni.ucla.edu/Research/Atlases/Data/vervet/vervetmratlas/vervetmratlas.html) and the other based on a high-resolution cryomacrotome study of a single vervet (http://www.loni.ucla.edu/Research/Atlases/Data/vervet/vervetatlas/vervetatlas.html). The averaged MRI atlas is also available as a volume in Neuroimaging Informatics Technology Initiative format. In the cryomacrotome atlas, various sulcal and subcortical structures have been anatomically labeled and surface rendered views are provided along the primary planes of section. Both atlases simultaneously provide views in all three primary planes of section, rapid navigation by clicking on the displayed images, and stereotaxic coordinates in the averaged MRI atlas space. Despite the extended time period since their divergence, the major sulcal and subcortical landmarks in vervets are highly conserved relative to those described in macaques.

We carried out a genome-wide association study of hemoglobin levels in 16,001 individuals of European and Indian Asian ancestry. The most closely associated SNP (rs855791) results in nonsynonymous (V736A) change in the serine protease domain of TMPRSS6 and a blood hemoglobin concentration 0.13 (95% CI 0.09–0.17) g/dl lower per copy of allele A (P = 1.6 × 10−13). Our findings suggest that TMPRSS6, a regulator of hepcidin synthesis and iron handling, is crucial in hemoglobin level maintenance.

Two independent lines of evidence support the localization of a schizophrenia susceptibility locus to the proximal long arm of chromosome 5. A partial trisomy of chromosome 5 (5q11.2–q13.3) cosegregates with the disorder in a Canadian family of Chinese descent, and DNA markers from proximal 5q cosegregate with schizophrenia (plus related disorders) in families of British and Icelandic descent. We constructed a human:hamster hybrid cell line (HHW 1064) whose only human complement is a chromosome 5 that is missing the trisomic region associated with schizophrenia. In combination with a “matched” cell hybrid (HHW 105) containing an intact chromosome 5, we physically mapped DNA markers relative to the trisomy. “Schizophrenia-linked” DNA markers p105–153Ra (D5S39) and p105-599Ha (D5S76) map within the trisomy and proximal to the 5q11.2 breakpoint, respectively. The hybrid cell lines HHW 105 and HHW 1064 together provide a means to identify and generate syntenic DNA markers to further investigate the location of a schizophrenia locus.

Although genome-wide association studies (GWASs) have identified numerous loci associated with complex traits, imprecise modeling of the genetic relatedness within study samples may cause substantial inflation of test statistics and possibly spurious associations. Variance component approaches, such as efficient mixed-model association (EMMA), can correct for a wide range of sample structures by explicitly accounting for pairwise relatedness between individuals, using high-density markers to model the phenotype distribution; but such approaches are computationally impractical. We report here a variance component approach implemented in publicly available software, EMMA eXpedited (EMMAX), that reduces the computational time for analyzing large GWAS data sets from years to hours. We apply this method to two human GWAS data sets, performing association analysis for ten quantitative traits from the Northern Finland Birth Cohort and seven common diseases from the Wellcome Trust Case Control Consortium. We find that EMMAX outperforms both principal component analysis and genomic control in correcting for sample structure.

Phenomics is an emerging transdiscipline dedicated to the systematic study of phenotypes on a genome-wide scale. New methods for high-throughput genotyping have changed the priority for biomedical research to phenotyping, but the human phenome is vast and its dimensionality remains unknown. Phenomics research strategies capable of linking genetic variation to public health concerns need to prioritize development of mechanistic frameworks that relate neural systems functioning to human behavior. New approaches to phenotype definition will benefit from crossing neuropsychiatric syndromal boundaries, and defining phenotypic features across multiple levels of expression from proteome to syndrome. The demand for high throughput phenotyping may stimulate a migration from conventional laboratory to web-based assessment of behavior, and this offers the promise of dynamic phenotyping –the iterative refinement of phenotype assays based on prior genotype-phenotype associations. Phenotypes that can be studied across species may provide greatest traction, particularly given rapid development in transgenic modeling. Phenomics research demands vertically integrated research teams, novel analytic strategies and informatics infrastructure to help manage complexity. The Consortium for Neuropsychiatric Phenomics at UCLA has been supported by the NIH Roadmap Initiative to illustrate these principles, and is developing applications that may help investigators assemble, visualize, and ultimately test multi-level phenomics hypotheses. As the transdiscipline of phenomics matures, and work is extended to large-scale international collaborations, there is promise that systematic new knowledgebases will help fulfill the promise of personalized medicine and the rational diagnosis and treatment of neuropsychiatric syndromes.

We investigated the strength of evidence for association of the 5-HTTLPR polymorphism and the personality trait of Harm Avoidance. We used new primary data from a large sample of adults drawn from the Finnish population. We also applied meta-analytic techniques to synthesize existing published data. The large number studies of the 5-HTTLPR polymorphism allowed us to apply a formal test of publication bias, as well as formally investigate the impact of potential moderating factors such as measurement instrument.

Univariate ANOVA of primary data (n = 3,872), with 5-HTTLPR genotype as a between-groups factor, indicated no evidence of association with Harm Avoidance (p = 0.99). Meta-analysis indicated no evidence of significant association of 5-HTTLPR with Harm Avoidance (d = 0.02, p = 0.37), or EPQ Neuroticism (d = 0.01, p = 0.71), although there was evidence of association with NEO Neuroticism (d = 0.18, p < 0.001).

Our analyses indicate that the 5-HTTLPR variant is not associated with Harm Avoidance. Together with our previous analyses of a large sample of participants with extreme Neuroticism scores (defined by the EPQ), we have data that excludes a meaningful genetic effect of the 5-HTTLPR on two measures of anxiety-related personality traits. There remains the possibility that the variant influences the NEO personality questionnaire measure of Neuroticism. However, a large, well-powered primary study is required to test this hypothesis directly and adequately.

Background

Although asthma is highly prevalent among certain Hispanic subgroups, genetic determinants of asthma and asthma‐related traits have not been conclusively identified in Hispanic populations. A study was undertaken to identify genomic regions containing susceptibility loci for pulmonary function and bronchodilator responsiveness (BDR) in Costa Ricans.

Methods

Eight extended pedigrees were ascertained through schoolchildren with asthma in the Central Valley of Costa Rica. Short tandem repeat (STR) markers were genotyped throughout the genome at an average spacing of 8.2 cM. Multipoint variance component linkage analyses of forced expiratory volume in 1 second (FEV1) and FEV1/ forced vital capacity (FVC; both pre‐bronchodilator and post‐bronchodilator) and BDR were performed in these eight families (pre‐bronchodilator spirometry, n = 640; post‐bronchodilator spirometry and BDR, n = 624). Nine additional STR markers were genotyped on chromosome 7. Secondary analyses were repeated after stratification by cigarette smoking.

Results

Among all subjects, the highest logarithm of the odds of linkage (LOD) score for FEV1 (post‐bronchodilator) was found on chromosome 7q34–35 (LOD = 2.45, including the additional markers). The highest LOD scores for FEV1/FVC (pre‐bronchodilator) and BDR were found on chromosomes 2q (LOD = 1.53) and 9p (LOD = 1.53), respectively. Among former and current smokers there was near‐significant evidence of linkage to FEV1/FVC (post‐bronchodilator) on chromosome 5p (LOD = 3.27) and suggestive evidence of linkage to FEV1 on chromosomes 3q (pre‐bronchodilator, LOD = 2.74) and 4q (post‐bronchodilator, LOD = 2.66).

Conclusions

In eight families of children with asthma in Costa Rica, there is suggestive evidence of linkage to FEV1 on chromosome 7q34–35. In these families, FEV1/FVC may be influenced by an interaction between cigarette smoking and a locus (loci) on chromosome 5p.

Genetic approaches to dissecting complex traits in animal models increasingly use transcript levels as a molecular phenotype and as validation for predictions of gene function. A recent study in BMC Biology using these approaches shows the complexity of the genetic contribution to aggressive behavior in Drosophila.

Now that genome-wide association studies (GWAS) are dominating the landscape of genetic research on neuropsychiatric syndromes, investigators are being faced with complexity on an unprecedented scale. It is now clear that phenomics, the systematic study of phenotypes on a genome-wide scale, comprises a rate-limiting step on the road to genomic discovery. To gain traction on the myriad paths leading from genomic variation to syndromal manifestations, informatics strategies must be deployed to navigate increasingly broad domains of knowledge and help researchers find the most important signals. The success of the Gene Ontology project suggests the potential benefits of developing schemata to represent higher levels of phenotypic expression. Challenges in cognitive ontology development include the lack of formal definitions of key concepts and relations among entities, the inconsistent use of terminology across investigators and time, and the fact that relations among cognitive concepts are not likely to be well represented by simple hierarchical “tree” structures. Because cognitive concept labels are labile, there is a need to represent empirical findings at the cognitive test indicator level. This level of description has greater consistency, and benefits from operational definitions of its concepts and relations to quantitative data. Considering cognitive test indicators as the foundation of cognitive ontologies carries several implications, including the likely utility of cognitive task taxonomies. The concept of cognitive “test speciation” is introduced to mark the evolution of paradigms sufficiently unique that their results cannot be “mated” productively with others in meta-analysis. Several projects have been initiated to develop cognitive ontologies at the Consortium for Neuropsychiatric Phenomics (www.phenomics.ucla.edu), in hope that these ultimately will enable more effective collaboration, and facilitate connections of information about cognitive phenotypes to other levels of biological knowledge. Several free web applications are available already to support examination and visualization of cognitive concepts in the literature (PubGraph, PubAtlas, PubBrain) and to aid collaborative development of cognitive ontologies (Phenowiki and the Cognitive Atlas). It is hoped that these tools will help formalize inference about cognitive concepts in behavioral and neuroimaging studies, and facilitate discovery of the genetic bases of both healthy cognition and cognitive disorders.

Genome-wide gene expression studies may provide substantial insight into gene activities and biological pathways differing between tissues and individuals. We investigated such gene expression variation by analyzing expression profiles in brain tissues derived from eight different brain regions and from blood in 12 monkeys from a biomedically important non-human primate model, the vervet (Chlorocebus aethiops sabaeus). We characterized brain regional differences in gene expression, focusing on transcripts for which inter-individual variation of expression in brain correlates well with variation in blood from the same individuals. Using stringent criteria, we identified 29 transcripts whose expression is measurable, stable, replicable, variable between individuals, relevant to brain function and heritable. Polymorphisms identified in probe regions could, in a minority of transcripts, confound the interpretation of the observed inter-individual variation. The high heritability of levels of these transcripts in a large vervet pedigree validated our approach of focusing on transcripts that showed higher inter-individual compared with intra-individual variation. These selected transcripts are candidate expression Quantitative Trait Loci, differentially regulating transcript levels in the brain among individuals. Given the high degree of conservation of tissue expression profiles between vervets and humans, our findings may facilitate the understanding of regional and individual transcriptional variation and its genetic mechanisms in humans. The approach employed here—utilizing higher quality tissue and more precise dissection of brain regions than is usually possible in humans—may therefore provide a powerful means to investigate variation in gene expression relevant to complex brain related traits, including human neuropsychiatric diseases.

Rationale: The basis for gender influences on allergen-specific IgEs is unclear.

Objectives: To perform regular and sex-stratified genomewide linkage analyses of IgE to each of three allergens (Ascaris lumbricoides, Blatella germanica [German cockroach]), and Dermatophagoides pteronyssinus [dust mite]) and to conduct an association study of a candidate gene in a linked genomic region.

Methods: Genomewide linkage analyses of allergen-specific IgEs were conducted in 653 members of eight large families of Costa Rican children with asthma. An analysis of the association between single-nucleotide polymorphisms in thymic stromal lymphopoietin (TSLP) and IgE measurements was conducted in 417 parent–child trios in Costa Rica. Significant results were replicated in 470 families of white children in the Childhood Asthma Management Program (CAMP).

Measurements and Main Results: Among all subjects, there was suggestive evidence of linkage (LOD ⩾ 2.72) to IgE to Ascaris (on chromosome 7q) and IgE to dust mite (on chromosomes 7p and 12q). In a sex-stratified analysis, there was significant evidence of linkage to IgE to cockroach on chromosome 5q23 (peak LOD, 4.14 at 127 cM) in female subjects. TSLP is located within the 1.5 LOD-unit support interval for this linkage peak and has female-specific effects on lung disease in mice. In a sex-stratified analysis, the T allele of single-nucleotide polymorphism rs2289276 in TSLP was associated with reductions in IgE to cockroach (in Costa Rican girls) and total IgE (in girls in Costa Rica and in CAMP; P value for sex-by-genotype interaction, <0.01 in both studies).

Conclusions: Consistent with findings in murine models, a variant in TSLP may have female-specific effects on allergic phenotypes.

Recent studies have identified a small number of genomic rearrangements that occur frequently in the general population. Bioinformatics tools are now available for systematic genome-wide surveys of higher-order structures predisposing to such common variations in genomic architecture. Segmental duplications (SDs) constitute up to 5 per cent of the genome and play an important role in generating additional rearrangements and in disease aetiology. We conducted a genome-wide database search for a form of SD, palindromic segmental duplications (PSDs), which consist of paired, inverted duplications, and which predispose to inversions, duplications and deletions. The survey was complemented by a search for SDs in tandem orientation (TSDs) that can mediate duplications and deletions but not inversions. We found more than 230 distinct loci with higher-order genomic structure that can mediate genomic variation, of these about 180 contained a PSD. A number of these sites were previously identified as harbouring common inversions or as being associated with specific genomic diseases characterised by duplication, deletions or inversions. Most of the regions, however, were previously unidentified; their characterisation should identify further common rearrangements and may indicate localisations for additional genomic disorders. The widespread distribution of complex chromosomal architecture suggests a potentially high degree of plasticity of the human genome and could uncover another level of genetic variation within human populations.

Background

Investigation of the environmental influences on human behavioral phenotypes is important for our understanding of the causation of psychiatric disorders. However, there are complexities associated with the assessment of environmental influences on behavior.

Methods/Principal Findings

We conducted a series of analyses using a prospective, longitudinal study of a nationally representative birth cohort from Finland (the Northern Finland 1966 Birth Cohort). Participants included a total of 3,761 male and female cohort members who were living in Finland at the age of 16 years and who had complete temperament scores. Our initial analyses (Wessman et al., in press) provide evidence in support of four stable and robust temperament clusters. Using these temperament clusters, as well as independent temperament dimensions for comparison, we conducted a data-driven analysis to assess the influence of a broad set of life course measures, assessed pre-natally, in infancy, and during adolescence, on adult temperament.

Results

Measures of early environment, neurobehavioral development, and adolescent behavior significantly predict adult temperament, classified by both cluster membership and temperament dimensions. Specifically, our results suggest that a relatively consistent set of life course measures are associated with adult temperament profiles, including maternal education, characteristics of the family’s location and residence, adolescent academic performance, and adolescent smoking.

Conclusions

Our finding that a consistent set of life course measures predict temperament clusters indicate that these clusters represent distinct developmental temperament trajectories and that information about a subset of life course measures has implications for adult health outcomes.

Asymmetry is a prominent feature of human brains with important functional consequences. Many asymmetric traits show population bias, but little is known about the genetic and environmental sources contributing to inter-individual variance. Anatomic asymmetry has been observed in Old World monkeys, but the evidence for the direction and extent of asymmetry is equivocal and only one study has estimated the genetic contributions to inter-individual variance. In this study we characterize a range of qualitative and quantitative asymmetry measures in structural brain MRIs acquired from an extended pedigree of Old World vervet monkeys (n = 357), and implement variance component methods to estimate the proportion of trait variance attributable to genetic and environmental sources. Four of six asymmetry measures show pedigree-level bias and one of the traits has a significant heritability estimate of about 30%. We also found that environmental variables more significantly influence the width of the right compared to the left prefrontal lobe.

Recent genome-wide association (GWA) studies described 95 loci controlling serum lipid levels. These common variants explain ∼25% of the heritability of the phenotypes. To date, no unbiased screen for gene–environment interactions for circulating lipids has been reported. We screened for variants that modify the relationship between known epidemiological risk factors and circulating lipid levels in a meta-analysis of genome-wide association (GWA) data from 18 population-based cohorts with European ancestry (maximum N = 32,225). We collected 8 further cohorts (N = 17,102) for replication, and rs6448771 on 4p15 demonstrated genome-wide significant interaction with waist-to-hip-ratio (WHR) on total cholesterol (TC) with a combined P-value of 4.79×10−9. There were two potential candidate genes in the region, PCDH7 and CCKAR, with differential expression levels for rs6448771 genotypes in adipose tissue. The effect of WHR on TC was strongest for individuals carrying two copies of G allele, for whom a one standard deviation (sd) difference in WHR corresponds to 0.19 sd difference in TC concentration, while for A allele homozygous the difference was 0.12 sd. Our findings may open up possibilities for targeted intervention strategies for people characterized by specific genomic profiles. However, more refined measures of both body-fat distribution and metabolic measures are needed to understand how their joint dynamics are modified by the newly found locus.

Author Summary

Circulating serum lipids contribute greatly to the global health by affecting the risk for cardiovascular diseases. Serum lipid levels are partly inherited, and already 95 loci affecting high- and low-density lipoprotein cholesterol, total cholesterol, and triglycerides have been found. Serum lipids are also known to be affected by multiple epidemiological risk factors like body composition, lifestyle, and sex. It has been hypothesized that there are loci modifying the effects between risk factors and serum lipids, but to date only candidate gene studies for interactions have been reported. We conducted a genome-wide screen with meta-analysis approach to identify loci having interactions with epidemiological risk factors on serum lipids with over 30,000 population-based samples. When combining results from our initial datasets and 8 additional replication cohorts (maximum N = 17,102), we found a genome-wide significant locus in chromosome 4p15 with a joint P-value of 4.79×10−9 modifying the effect of waist-to-hip ratio on total cholesterol. In the area surrounding this genetic variant, there were two genes having association between the genotypes and the gene expression in adipose tissue, and we also found enrichment of association in genes belonging to lipid metabolism related functions.

Background

C-reactive protein (CRP) is a heritable marker of chronic inflammation that is strongly associated with cardiovascular disease. We aimed to identify genetic variants that are associated with CRP levels.

Methods and Results

We performed a genome wide association (GWA) analysis of CRP in 66,185 participants from 15 population-based studies. We sought replication for the genome wide significant and suggestive loci in a replication panel comprising 16,540 individuals from ten independent studies. We found 18 genome-wide significant loci and we provided evidence of replication for eight of them. Our results confirm seven previously known loci and introduce 11 novel loci that are implicated in pathways related to the metabolic syndrome (APOC1, HNF1A, LEPR, GCKR, HNF4A, and PTPN2), immune system (CRP, IL6R, NLRP3, IL1F10, and IRF1), or that reside in regions previously not known to play a role in chronic inflammation (PPP1R3B, SALL1, PABPC4, ASCL1, RORA, and BCL7B). We found significant interaction of body mass index (BMI) with LEPR (p<2.9×10−6). A weighted genetic risk score that was developed to summarize the effect of risk alleles was strongly associated with CRP levels and explained approximately 5% of the trait variance; however, there was no evidence for these genetic variants explaining the association of CRP with coronary heart disease.

Conclusion

We identified 18 loci that were associated with CRP levels. Our study highlights immune response and metabolic regulatory pathways involved in the regulation of chronic inflammation.

Background & Aims

Progressive familial intrahepatic cholestasis (PFIC) with normal serum levels of gamma-glutamyltranspeptidase can result from mutations in ATP8B1 (encoding familial intrahepatic cholestasis 1 [FIC1]) or ABCB11 (encoding bile salt export pump [BSEP]). We evaluated clinical and laboratory features of disease in patients diagnosed with PFIC, who carried mutations in ATP8B1 (FIC1 deficiency) or ABCB11 (BSEP deficiency). Our goal was to identify features that distinguish presentation and course of these 2 disorders, thus facilitating diagnosis and elucidating the differing consequences of ATP8B1 and ABCB11 mutations.

Methods

A retrospective multi-center study was conducted, using questionnaires and chart review. Available clinical and biochemical data from 145 PFIC patients with mutations in either ATP8B1 (61 “FIC1 patients”) or ABCB11 (84 “BSEP patients”) were evaluated.

Results

At presentation, serum aminotransferase and bile salt levels were higher in BSEP patients; serum alkaline phosphatase values were higher, and serum albumin values were lower, in FIC1 patients. Elevated white blood cell counts, and giant or multinucleate cells at liver biopsy, were more common in BSEP patients. BSEP patients more often had gallstones and portal hypertension. Diarrhea, pancreatic disease, rickets, pneumonia, abnormal sweat tests, hearing impairment, and poor growth were more common in FIC1 patients. Among BSEP patients, the course of disease was less rapidly progressive in patients bearing the D482G mutation.

Conclusions

Severe forms of FIC1 and BSEP deficiency differed. BSEP patients manifested more severe hepatobiliary disease, while FIC1 patients showed greater evidence of extrahepatic disease.

Context:

Plasma levels of C-reactive protein (CRP) are independently associated with risk of coronary heart disease, but whether CRP is causally associated with coronary heart disease or merely a marker of underlying atherosclerosis is uncertain.

Objective:

To investigate association of genetic loci with CRP levels and risk of coronary heart disease.

Design, setting and participants:

We first carried out a genome-wide association (n=17,967) and replication study (n=14,747) to identify genetic loci associated with plasma CRP concentrations. Data collection took place between 1989 and 2008 and genotyping between 2003 and 2008. We carried out a Mendelian randomisation study of the most closely associated SNP in the CRP locus and published data on other CRP variants involving a total of 28,112 cases and 100,823 controls, to investigate the association of CRP variants with coronary heart disease. We compared our finding with that predicted from meta-analysis of observational studies of CRP levels and risk of coronary heart disease. For the other loci associated with CRP levels, we selected the most closely associated SNP for testing against coronary heart disease among 14,365 cases and 32,069 controls.

Main outcome measure:

Risk of coronary heart disease.

Results:

Polymorphisms in five genetic loci were strongly associated with CRP levels (% difference per minor allele): SNP rs6700896 in LEPR (−14.7% [95% Confidence Interval {CI}], −17.5 – −11.9, P=1.6×10−21), rs4537545 in IL6R (−10.8% [95% CI, −13.8 – −7.7], P=5.1×10−11), rs7553007 in CRP locus (−20.7% [95% CI, −23.5 – −17.9], P=3.3×10−38), rs1183910 in HNF1A (−13.6% [95% CI, −16.4 – −10.6], P=1.2×10−17) and rs4420638 in APOE-CI-CII (−21.8% [95% CI, −25.4 – −18.1], P=2.1×10−25). Association of SNP rs7553007 in the CRP locus with coronary heart disease gave odds ratio (OR) 0.98 (95% CI, 0.94 – 1.01) per 20% lower CRP. Our Mendelian randomisation study of variants in the CRP locus showed no association with coronary heart disease: OR 1.00 (95% CI, 0.97 – 1.02) per 20% lower CRP, compared with OR 0.94 (95% CI, 0.94 – 0.95) predicted from meta-analysis of the observational studies of CRP levels and coronary heart disease (Z-score −3.45, P<.001). SNPs rs6700896 in LEPR (OR 1.06 [95% CI, 1.02 – 1.09] per minor allele), rs4537545 in IL6R (OR 0.94 [95% CI, 0.91 – 0.97]) and rs4420638 in the APOE-CI-CII cluster (OR 1.16 [95% CI, 1.12 – 1.21]) were all associated with risk of coronary heart disease.

Conclusions:

The lack of concordance between the effect on coronary heart disease risk of CRP genotypes and CRP levels argues against a causal association of CRP with coronary heart disease.

The area and volume of brain structural features, as assessed by high-resolution 3D magnetic resonance imaging (MRI), are among the most heritable measures relating to the human central nervous system. We have conducted MRI scanning of all available monkeys over 2 years of age (n=357) from the extended multigenerational pedigree of the Vervet Research Colony (VRC). Using a combination of automated and manual segmentation we have quantified several correlated but distinct brain structural phenotypes. The estimated heritabilities (h2) for these measures in the VRC are higher than those reported previously for such features in humans or in other non human primates (NHP): total brain volume (h2=0.99, standard error (se)=0.06), cerebral volume (h2=0.98, se=0.06), cerebellar volume (h2=0.86, se=0.09), hippocampal volume (h2=0.95, se=0.07) and corpus callosum cross-sectional areas (h2=0.87, se=0.07). These findings indicate that, in the controlled environment and with the inbreeding structure of the VRC, additive genetic factors account for almost all of the observed variance in brain structure, and suggest the potential of the VRC for genetic mapping of quantitative trait loci (QTL) underlying such variance.

Genome-wide association studies (GWAS) of longitudinal birth cohorts enable joint investigation of environmental and genetic influences on complex traits. We report GWAS results for nine quantitative metabolic traits (triglycerides, high-density lipoprotein, low-density lipoprotein, glucose, insulin, C-reactive protein, body mass index, and systolic and diastolic blood pressure) in the Northern Finland Birth Cohort 1966 (NFBC1966), drawn from the most genetically isolated Finnish regions. We replicate most previously reported associations for these traits and identify nine new associations, several of which highlight genes with metabolic functions: high-density lipoprotein with NR1H3 (LXRA), low-density lipoprotein with AR and FADS1-FADS2, glucose with MTNR1B, and insulin with PANK1. Two of these new associations emerged after adjustment of results for body mass index. Gene-environment interaction analyses suggested additional associations, which will require validation in larger samples. The currently identified loci, together with quantified environmental exposures, explain little of the trait variation in NFBC1966. The association observed between low-density lipoprotein and an infrequent variant in AR suggests the potential of such a cohort for identifying associations with both common, low-impact and rarer, high-impact quantitative trait loci.

Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797-22,562 persons, aged 18-104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 × 10-8), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 × 10-11; LDL, P = 2.6 × 10-10), TMEM57 (TC, P = 5.4 × 10-10), CTCF-PRMT8 region (HDL, P = 8.3 × 10-16), DNAH11 (LDL, P = 6.1 × 10-9), FADS3-FADS2 (TC, P = 1.5 × 10-10; LDL, P = 4.4 × 10-13) and MADD-FOLH1 region (HDL, P = 6 × 10-11). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors.

Recent genome-wide association (GWA) studies have identified dozens of common variants associated with adult height. However, it is unknown how these variants influence height growth during childhood. We derived peak height velocity in infancy (PHV1) and puberty (PHV2) and timing of pubertal height growth spurt from parametric growth curves fitted to longitudinal height growth data to test their association with known height variants. The study consisted of N = 3,538 singletons from the prospective Northern Finland Birth Cohort 1966 with genotype data and frequent height measurements (on average 20 measurements per person) from 0–20 years. Twenty-six of the 48 variants tested associated with adult height (p<0.05, adjusted for sex and principal components) in this sample, all in the same direction as in previous GWA scans. Seven SNPs in or near the genes HHIP, DLEU7, UQCC, SF3B4/SV2A, LCORL, and HIST1H1D associated with PHV1 and five SNPs in or near SOCS2, SF3B4/SV2A, C17orf67, CABLES1, and DOT1L with PHV2 (p<0.05). We formally tested variants for interaction with age (infancy versus puberty) and found biologically meaningful evidence for an age-dependent effect for the SNP in SOCS2 (p = 0.0030) and for the SNP in HHIP (p = 0.045). We did not have similar prior evidence for the association between height variants and timing of pubertal height growth spurt as we had for PHVs, and none of the associations were statistically significant after correction for multiple testing. The fact that in this sample, less than half of the variants associated with adult height had a measurable effect on PHV1 or PHV2 is likely to reflect limited power to detect these associations in this dataset. Our study is the first genetic association analysis on longitudinal height growth in a prospective cohort from birth to adulthood and gives grounding for future research on the genetic regulation of human height during different periods of growth.

Author Summary

Family studies have shown that adult height is largely genetically determined. Identification of common genetic factors has been expedited with recent advances in genotyping techniques. However, factors regulating childhood height growth remain unclear. We investigated genetic variants of adult height for associations with peak height velocity in infancy (PHV1) and puberty (PHV2) and timing of pubertal growth spurt in a population based sample of 3,538 Finns born in 1966. Most variants studied associated with adult height in this sample. Of the 48 genetic variants tested, seven of them associated with PHV1 and five with PHV2. However, only one of these associated with both, and we found suggestive evidence for differential effects at different stages of growth for some of the variants. In this sample, less than half of the variants associated with adult height had a measurable effect on PHV1 or PHV2. However, these differences may reflect lower statistical power to detect associations with height velocities compared to adult height. This study provides a foundation for further biological investigation into the genes acting at each stage of height growth.

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.

Author Summary

The history of Latin America has entailed a complex process of population mixture between Native and recent immigrants across a vast geographic region. Few details are known about this process or about how it has shaped the genetic makeup of contemporary Latin American populations. To perform a broad exploration of the genetic diversity of Latin America we carried out genome-wide analyses in 13 mestizo populations sampled from 7 countries across the region. We observe a marked variation in ancestry both within and between mestizo populations. This variation in ancestry correlates with pre-Columbian Native population density in the areas examined and with recent patterns of demographic growth of the sites sampled. We also find evidence that the mixture at the origin of these populations involved mainly immigrant European men and Native and African women. Finally, mestizo populations show a differentiated Amerindian genetic background, consistent with a predominantly local Native ancestry. Mestizos thus still reveal the genetic imprint of the pre-Columbian Native American population diversification. Our study helps delineate the genetic landscape of Latin America and has a number of implications for gene identification analyses in populations from the region.