A genome-wide association study of monozygotic twin-pairs suggests a locus related to variability of serum high-density lipoprotein cholesterol

Genetics of environmental sensitivity to psychiatric and neurodevelopmental phenotypes: evidence from GWAS of monozygotic twins

Individual sensitivity to environmental exposures may be genetically influenced. This genotype-by-environment interplay implies differences in phenotypic variance across genotypes. However, environmental sensitivity genetic variants have proven challenging to detect. GWAS of monozygotic twin differences is a family-based variance analysis method, which is more robust to systemic biases that impact population-based methods. We combined data from up to 21,792 monozygotic twins (10,896 pairs) from 11 studies to conduct the largest GWAS meta-analysis of monozygotic phenotypic differences in children and adolescents/adults for seven psychiatric and neurodevelopmental phenotypes: attention deficit hyperactivity disorder (ADHD) symptoms, autistic traits, anxiety and depression symptoms, psychotic-like experiences, neuroticism, and wellbeing. The SNP-heritability of variance in these phenotypes were estimated (h2: 0% to 18%), but were imprecise. We identified a total of 13 genome-wide significant associations (SNP, gene, and gene-set), including genes related to stress-reactivity for depression, growth factor-related genes for autistic traits and catecholamine uptake-related genes for psychotic-like experiences. Monozygotic twins are an important new source of evidence about the genetics of environmental sensitivity.

Twenty-Five Novel Loci for Carotid Intima-Media Thickness: A Genome-Wide Association Study in >45 000 Individuals and Meta-Analysis of >100 000 Individuals

OBJECTIVE

Carotid artery intima-media thickness (cIMT) is a widely accepted marker of subclinical atherosclerosis. Twenty susceptibility loci for cIMT were previously identified and the identification of additional susceptibility loci furthers our knowledge on the genetic architecture underlying atherosclerosis.

APPROACH AND RESULTS

We performed 3 genome-wide association studies in 45 185 participants from the UK Biobank study who underwent cIMT measurements and had data on minimum, mean, and maximum thickness. We replicated 15 known loci and identified 20 novel loci associated with cIMT at P<5×10-8. Seven novel loci (ZNF385D, ADAMTS9, EDNRA, HAND2, MYOCD, ITCH/EDEM2/MMP24, and MRTFA) were identified in all 3 phenotypes. An additional new locus (LOXL1) was identified in the meta-analysis of the 3 phenotypes. Sex interaction analysis revealed sex differences in 7 loci including a novel locus (SYNE3) in males. Meta-analysis of UK Biobank data with a previous meta-analysis led to identification of three novel loci (APOB, FIP1L1, and LOXL4). Transcriptome-wide association analyses implicated additional genes ARHGAP42, NDRG4, and KANK2. Gene set analysis showed an enrichment in extracellular organization and the PDGF (platelet-derived growth factor) signaling pathway. We found positive genetic correlations of cIMT with coronary artery disease rg=0.21 (P=1.4×10-7), peripheral artery disease rg=0.45 (P=5.3×10-5), and systolic blood pressure rg=0.30 (P=4.0×10-18). A negative genetic correlation between average of maximum cIMT and high-density lipoprotein was found rg=-0.12 (P=7.0×10-4).

CONCLUSIONS

Genome-wide association meta-analyses in >100 000 individuals identified 25 novel loci associated with cIMT providing insights into genes and tissue-specific regulatory mechanisms of proatherosclerotic processes. We found evidence for shared biological mechanisms with cardiovascular diseases.

A Nutrigenetic Update on CETP Gene–Diet Interactions on Lipid-Related Outcomes

Purpose of Review

An abnormal lipid profile is considered a main risk factor for cardiovascular diseases and evidence suggests that single nucleotide polymorphisms (SNPs) in the cholesteryl ester transfer protein (CETP) gene contribute to variations in lipid levels in response to dietary intake. The objective of this review was to identify and discuss nutrigenetic studies assessing the interactions between CETP SNPs and dietary factors on blood lipids.

Recent Findings

Relevant articles were obtained through a literature search of PubMed and Google Scholar through to July 2021. An article was included if it examined an interaction between CETP SNPs and dietary factors on blood lipids. From 49 eligible nutrigenetic studies, 27 studies reported significant interactions between 8 CETP SNPs and 17 dietary factors on blood lipids in 18 ethnicities. The discrepancies in the study findings could be attributed to genetic heterogeneity, and differences in sample size, study design, lifestyle and measurement of dietary intake. The most extensively studied ethnicities were those of Caucasian populations and majority of the studies reported an interaction with dietary fat intake. The rs708272 (TaqIB) was the most widely studied CETP SNP, where ‘B1’ allele was associated with higher CETP activity, resulting in lower high-density lipoprotein cholesterol and higher serum triglycerides under the influence of high dietary fat intake.

Summary

Overall, the findings suggest that CETP SNPs might alter blood lipid profiles by modifying responses to diet, but further large studies in multiple ethnic groups are warranted to identify individuals at risk of adverse lipid response to diet.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11883-022-00987-y.

Genetics of Biochemical Phenotypes

Biomarkers diagnose, predict or assess the risk of disease, and studies of the effects of genetic variation on biomarker phenotypes in the general population complement studies on patients diagnosed with disease. This paper traces the evolution of studies on biomarker genetics over the past 40 years through examples drawn from the work of Professor Martin and his colleagues.

Cohorts

Cohort studies are essential for conducting large studies of multiple exposures and outcomes in humans. Recently, the ability to combine data from multiple cohorts in, for example, meta-analyses, and the willingness in the genetics community to collaborate to enable replication studies has led to many new insights into the genetic and environmental determinants of human health and behaviors. The contribution of Professor Nicholas Martin to the development of cohort studies, particularly of twin and twin-family studies, over a period of several decades is reviewed. He has contributed to the development and use of both Australian and international resources. The contributions of Australian twin studies to genomewide association projects are multiple, and across multiple domains, from biomarkers, lifestyle and behavior to disorders and disease.

The Swedish Twin Registry: Content and Management as a Research Infrastructure

The Swedish Twin Registry functions as research infrastructure containing information on 216,258 twins born between 1886 and 2015, of whom 86,199 pairs have zygosity determined by DNA, an intrapair similarity algorithm, or being of opposite sex. In essence, practically all twins alive and currently 9 years or older have been invited for participation and donation of DNA on which genomewide single nucleotide polymorphisms array genotyping has been performed. Content, management and alternatives for future improvements are discussed.

The Danish Twin Registry: An Updated Overview

The Danish Twin Registry (DTR) was established in the 1950s, when twins born from 1870 to 1910 were ascertained, and has since been extended to include twins from birth cohorts until 2009. The DTR currently comprises of more than 175,000 twins from the 140 birth cohorts. This makes the DTR the oldest nationwide twin register and among the largest in the world. The combination of data from several surveys, including biological samples and repeated measurements on the same individuals, and data from Danish national registers provides a unique resource for a wide range of twin studies. This article provides an updated overview of the data in the DTR: First, we provide a summary of the establishment of the register, the different ascertainment methods and the twins included; then follows an overview of major surveys conducted in the DTR since 1994 and a description of the DTR biobank, including a description of the molecular data created so far; finally, a short description is given of the linkage to Danish national registers at Statistics Denmark and some recent examples of studies using the various data resources in the DTR are highlighted.

Generalizability and applicability of results obtained from populations of European descent regarding the effect direction and size of HDL-C level-associated genetic variants to the Hungarian general and Roma populations

OBJECTIVES

Large-scale association studies that mainly involve European populations identified many genetic loci related to high-density lipoprotein cholesterol (HDL-C) levels, one of the most important indicators of the risk for cardiovascular diseases. The question with intense speculation of whether the effect estimates obtained from European populations for different HDL-C level-related SNPs are applicable to the Roma ethnicity, the largest minority group in Europe with a South Asian origin, was addressed in the present study.

DESIGN

The associations between 21 SNPs (in the genes LIPC(G), CETP, GALNT2, HMGCP, ABCA1, KCTD10 and WWOX) and HDL-C levels were examined separately in adults of the Hungarian general (N = 1542) and Roma (N = 646) populations by linear regression. Individual effects (direction and size) of single SNPs on HDL-C levels were computed and compared between the study groups and with data published in the literature.

RESULTS

Significant associations between SNPs and HDL-C levels were more frequently found in general subjects than in Roma subjects (11 SNPs in general vs. 4 SNPs in Roma). The CETP gene variants rs1532624, rs708272 and rs7499892 consistently showed significant associations with HDL-C levels across the study groups (p ˂ 0.05), indicating a possible causal variant(s) in this region. Although nominally significant differences in effect size were found for three SNPs (rs693 in gene APOB, rs9989419 in gene CETP, and rs2548861 in gene WWOX) by comparing the general and Roma populations, most of these SNPs did not have a significant effect on HDL-C levels. The β coefficients for SNPs in the Roma population were found to be identical both in direction and magnitude to the effect obtained previously in large-scale studies on European populations.

CONCLUSIONS

The effect of the vast majority of the SNPs on HDL-C levels could be replicated in the Hungarian general and Roma populations, which indicates that the effect size measurements obtained from the literature can be used for risk estimation for both populations.

Association of BDNF rs11030104 SNP and serum lipid levels in two Chinese ethnic groups

The correlation between the BDNF rs11030104 single nucleotide polymorphism (SNP) and serum lipid levels has been understudied. The present study was conducted to detect the association of the BDNF rs11030104 SNP and several environmental factors with serum lipid levels in the Jing and Han nationalities. Genotypes of the BDNF rs11030104 SNP in 709 unrelated subjects of Han and 706 unrelated participants of Jing populations were determined by polymerase chain reaction and restriction fragment length polymorphism combined with gel electrophoresis, and further verified by direct sequencing. There was no significant difference in either genotypic or allelic frequencies between the Han and Jing populations. The genotypic and allelic frequencies of the SNP in Jing but not in Han populations were different between male and female subgroups (P<0.05 for each). The levels of serum total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) in the Jing population were different among the genotypes, the G allele carriers had lower TC and LDL-C levels than the G allele non-carriers. Subgroup analyses showed that the differences in serum TC and LDL-C levels among the genotypes were observed in the Jing males but not in females. Serum lipid profiles were also significantly associated with some environmental factors in the Han and Jing populations, or in male and female subgroups of the two ethnic groups (P<0.05 for all). Our study exhibited a correlation between the BDNF rs11030104 SNP and serum TC and LDL-C levels in the Jing males. These results indicate that there may be a racial/ethnic- and/or sex-specific association of the BDNF rs11030104 SNP and serum lipid parameters.

Heritability and Genome-Wide Association Study of Plasma Cholesterol in Chinese Adult Twins

Dyslipidemia represents a strong and independent risk factor for cardiovascular disease. Plasma cholesterol, such as total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), and high density lipoprotein cholesterol (HDL-C), is the common indicator of diagnosing dyslipidemia. Here based on 382 Chinese twin pairs, we explored the magnitude of genetic impact on TC, HDL-C, LDL-C variation and further searched for genetic susceptibility loci for them using genome-wide association study (GWAS). The ACE model was the best fit model with additive genetic parameter (A) accounting for 26.6%, common or shared environmental parameter (C) accounting for 47.8%, unique/non-shared environmental parameter (E) accounting for 25.6% for the variance in HDL-C. The AE model was the best fit model for TC (A: 61.4%; E: 38.6%) and LDL-C (A: 65.5%; E: 34.5%). While no SNPs reached the genome-wide significance level (P < 5 × 10^-8), 8, 14, 9 SNPs exceeded the suggestive significance level (P < 1 × 10^-5) for TC, HDL-C, LDL-C, respectively. The promising genetic regions for TC, HDL-C, LDL-C were on chromosome 11 around rs7107698, chromosome 5 around rs12518218, chromosome 2 around rs10490120, respectively. Gene-based analysis found 1038, 1033 and 1090 genes nominally associated with TC, HDL-C, LDL-C (P < 0.05), especially FAF1, KLKB1 for TC, KLKB1 for HDL-C, and NTRK1, FAF1, SNTB2 for LDL-C, respectively. The number of common related genes among TC, HDL-C and LDL-C was 71, including FAF1, KLKB1, etc. Pathway enrichment analysis discovered known related pathways-zinc transporters, metal ion SLC transporters for TC, cell adhesion molecules CAMs, IL-6 signaling for HDL, FC epsilon RI signaling pathway, NFAT pathway for LDL, respectively. In conclusion, the TC and LDL-C level is moderately heritable and the HDL-C level is lowly heritable in Chinese population. The genomic loci, functional genes and pathways are identified to account for the heritability of plasma cholesterol level. Our findings provide important insights into plasma cholesterol molecular physiology and expect future research to replicate and validate our results.

The SRGAP2 SNPs, their haplotypes and G × E interactions on serum lipid traits

Maonan nationality is a relatively conservative and isolated minority in China. Little is known about the association of the Slit-Robo Rho GTPase activating protein 2 gene (SRGAP2) single nucleotide polymorphisms (SNPs) and serum lipid levels in the Chinese populations. This study was performed to clarify the association of the SRGAP2 rs2483058 and rs2580520 SNPs and their haplotypes with serum lipid traits in the Maonan and Han populations. Genotyping of the 2 SNPs was performed in 2444 unrelated subjects (Han, 1210 and Maonan, 1234) by polymerase chain reaction and restriction fragment length polymorphism combined with gel electrophoresis, and then confirmed by direct sequencing. The allelic (rs2483058) and genotypic (rs2483058 and rs2580520) frequencies were different between the two ethnic groups. Four haplotypes were identified in our populations, and the rs2483058G-rs2580520C haplotype was the commonest one. The rs2483058C-rs2580520G haplotype was associated with an increased risk of dyslipidemia, and showed consistent association with serum total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), apolipoprotein (Apo) A1 levels, and the ApoA1/ApoB ratio. These results indicated that the SRGAP2 SNPs and their haplotypes were associated with serum lipid levels. Their haplotypes can explain much more serum lipid variation than any single SNP alone, especially for serum TC, HDL-C and ApoA1 levels.

MixFit: Methodology for Computing Ancestry-Related Genetic Scores at the Individual Level and Its Application to the Estonian and Finnish Population Studies

Ancestry information at the individual level can be a valuable resource for personalized medicine, medical, demographical and history research, as well as for tracing back personal history. We report a new method for quantitatively determining personal genetic ancestry based on genome-wide data. Numerical ancestry component scores are assigned to individuals based on comparisons with reference populations. These comparisons are conducted with an existing analytical pipeline making use of genotype phasing, similarity matrix computation and our addition-multidimensional best fitting by MixFit. The method is demonstrated by studying Estonian and Finnish populations in geographical context. We show the main differences in the genetic composition of these otherwise close European populations and how they have influenced each other. The components of our analytical pipeline are freely available computer programs and scripts one of which was developed in house (available at: www.geenivaramu.ee/en/tools/mixfit).

The Concordance and Heritability of Type 2 Diabetes in 34,166 Twin Pairs From International Twin Registers: The Discordant Twin (DISCOTWIN) Consortium

Twin pairs discordant for disease may help elucidate the epigenetic mechanisms and causal environmental factors in disease development and progression. To obtain the numbers of pairs, especially monozygotic (MZ) twin pairs, necessary for in-depth studies while also allowing for replication, twin studies worldwide need to pool their resources. The Discordant Twin (DISCOTWIN) consortium was established for this goal. Here, we describe the DISCOTWIN Consortium and present an analysis of type 2 diabetes (T2D) data in nearly 35,000 twin pairs. Seven twin cohorts from Europe (Denmark, Finland, Norway, the Netherlands, Spain, Sweden, and the United Kingdom) and one from Australia investigated the rate of discordance for T2D in same-sex twin pairs aged 45 years and older. Data were available for 34,166 same-sex twin pairs, of which 13,970 were MZ, with T2D diagnosis based on self-reported diagnosis and medication use, fasting glucose and insulin measures, or medical records. The prevalence of T2D ranged from 2.6% to 12.3% across the cohorts depending on age, body mass index (BMI), and national diabetes prevalence. T2D discordance rate was lower for MZ (5.1%, range 2.9-11.2%) than for same-sex dizygotic (DZ) (8.0%, range 4.9-13.5%) pairs. Across DISCOTWIN, 720 discordant MZ pairs were identified. Except for the oldest of the Danish cohorts (mean age 79), heritability estimates based on contingency tables were moderate to high (0.47-0.77). From a meta-analysis of all data, the heritability was estimated at 72% (95% confidence interval 61-78%). This study demonstrated high T2D prevalence and high heritability for T2D liability across twin cohorts. Therefore, the number of discordant MZ pairs for T2D is limited. By combining national resources, the DISCOTWIN Consortium maximizes the number of discordant MZ pairs needed for in-depth genotyping, multi-omics, and phenotyping studies, which may provide unique insights into the pathways linking genes to the development of many diseases.

A Genome-Wide Test of the Differential Susceptibility Hypothesis Reveals a Genetic Predictor of Differential Response to Psychological Treatments for Child Anxiety Disorders

BACKGROUND

The differential susceptibly hypothesis suggests that certain genetic variants moderate the effects of both negative and positive environments on mental health and may therefore be important predictors of response to psychological treatments. Nevertheless, the identification of such variants has so far been limited to preselected candidate genes. In this study we extended the differential susceptibility hypothesis from a candidate gene to a genome-wide approach to test whether a polygenic score of environmental sensitivity predicted response to cognitive behavioural therapy (CBT) in children with anxiety disorders.

METHODS

We identified variants associated with environmental sensitivity using a novel method in which within-pair variability in emotional problems in 1,026 monozygotic twin pairs was examined as a function of the pairs' genotype. We created a polygenic score of environmental sensitivity based on the whole-genome findings and tested the score as a moderator of parenting on emotional problems in 1,406 children and response to individual, group and brief parent-led CBT in 973 children with anxiety disorders.

RESULTS

The polygenic score significantly moderated the effects of parenting on emotional problems and the effects of treatment. Individuals with a high score responded significantly better to individual CBT than group CBT or brief parent-led CBT (remission rates: 70.9, 55.5 and 41.6%, respectively).

CONCLUSIONS

Pending successful replication, our results should be considered exploratory. Nevertheless, if replicated, they suggest that individuals with the greatest environmental sensitivity may be more likely to develop emotional problems in adverse environments but also benefit more from the most intensive types of treatment.

Genetic studies of plasma analytes identify novel potential biomarkers for several complex traits

Genome-wide association studies of 146 plasma protein levels in 818 individuals revealed 56 genome-wide significant associations (28 novel) with 47 analytes. Loci associated with plasma levels of 39 proteins tested have been previously associated with various complex traits such as heart disease, inflammatory bowel disease, Type 2 diabetes and multiple sclerosis. These data suggest that these plasma protein levels may constitute informative endophenotypes for these complex traits. We found three potential pleiotropic genes: ABO for plasma SELE and ACE levels, FUT2 for CA19-9 and CEA plasma levels and APOE for ApoE and CRP levels. We also found multiple independent signals in loci associated with plasma levels of ApoH, CA19-9, FetuinA, IL6r and LPa. Our study highlights the power of biological traits for genetic studies to identify genetic variants influencing clinically relevant traits, potential pleiotropic effects and complex disease associations in the same locus.

Gene-Environment Interplay in Physical, Psychological, and Cognitive Domains in Mid to Late Adulthood: Is APOE a Variability Gene?

Despite emerging interest in gene-environment interaction (GxE) effects, there is a dearth of studies evaluating its potential relevance apart from specific hypothesized environments and biometrical variance trends. Using a monozygotic within-pair approach, we evaluated evidence of G×E for body mass index (BMI), depressive symptoms, and cognition (verbal, spatial, attention, working memory, perceptual speed) in twin studies from four countries. We also evaluated whether APOE is a 'variability gene' across these measures and whether it partly represents the 'G' in G×E effects. In all three domains, G×E effects were pervasive across country and gender, with small-to-moderate effects. Age-cohort trends were generally stable for BMI and depressive symptoms; however, they were variable-with both increasing and decreasing age-cohort trends-for different cognitive measures. Results also suggested that APOE may represent a 'variability gene' for depressive symptoms and spatial reasoning, but not for BMI or other cognitive measures. Hence, additional genes are salient beyond APOE.

Association of common variants in TOMM40/APOE/APOC1 region with human longevity in a Chinese population

Apolipoprotein E (APOE), translocase of outer mitochondrial membrane 40 homolog (TOMM40) and apolipoprotein C-I (APOC1) may extend lifespan by marked delay or escape from age-related diseases. This study aimed to elucidate the association of human longevity with genetic variations in TOMM40/APOE/APOC1 region in a Chinese population. Ten tag single-nucleotide polymorphisms (SNPs) in the TOMM40/APOE/APOC1 region were successfully genotyped in 616 unrelated long-lived individuals and 846 younger controls. Of the 10 SNPs, rs7254892 in 5' upstream of TOMM40 showed significant association with human longevity (G/A-A/A vs G/G: odds ratio (OR)=1.59, 95% confidence interval (CI)=1.20-2.09, P=0.0011, Bonferroni corrected P (Pc)=0.033). The haplotype analysis suggested that individuals carrying the haplotype A-A-A-A-T-A-T-G-C-A (rs7254892-rs157580-rs2075649-rs2075650-rs157582-rs8106922-rs1160985-rs405697-rs439401-rs445925) tended to have longer lifespan than those carrying the most common haplotype G-G-A-A-C-A-C-A-T-G (OR=1.59, 95% CI=1.19-2.12, P=0.0018, Pc=0.0216). These findings indicated that variants in TOMM40/APOE/APOC1 region might be associated with human longevity. Further studies are needed to identify the causal genetic variants influencing human longevity.

Determination of Zygosity in Adult Chinese Twins Using the 450K Methylation Array versus Questionnaire Data

Previous studies have shown that both single nucleotide polymorphisms (SNPs) and questionnaires-based method can be used for twin zygosity determination, but few validation studies have been conducted using Chinese populations. In the current study, we recruited 192 same sex Chinese adult twin pairs to evaluate the validity of using genetic markers-based method and questionnaire-based method in zygosity determination. We considered the relatedness analysis based on more than 0.6 million SNPs genotyping as the golden standards for zygosity determination. After quality control, qualified twins were left for relatedness analysis based on identical by descent calculation. Then those same sex twin pairs were included in the zygosity questionnaire validation analysis. Logistic regression model was applied to assess the discriminant ability of age, sex and the three questions in zygosity determination. Leave one out cross-validation was used as a measurement of internal validation. The results of zygosity determination based on 65 SNPs in 450k methylation array were all consistent with genotyping. Age, gender, questions of appearance confused by strangers and previously perceived zygosity consisted of the most predictable model with a consistency rate of 0.8698, cross validation predictive error of 0.1347. For twin studies with genotyping and\or 450k methylation array, there would be no need to conduct other zygosity testing for the sake of costs consideration.

Gene-gene and gene-environment interactions detected by transcriptome sequence analysis in twins

Understanding the genetic architecture of gene expression is an intermediate step in understanding the genetic architecture of complex diseases. RNA sequencing technologies have improved the quantification of gene expression and allow measurement of allele-specific expression (ASE). ASE is hypothesized to result from the direct effect of cis regulatory variants, but a proper estimation of the causes of ASE has not been performed thus far. In this study, we take advantage of a sample of twins to measure the relative contributions of genetic and environmental effects to ASE, and we find substantial effects from gene × gene (G×G) and gene × environment (G×E) interactions. We propose a model where ASE requires genetic variability in cis, a difference in the sequence of both alleles, but where the magnitude of the ASE effect depends on trans genetic and environmental factors that interact with the cis genetic variants.

Genome-wide family-based linkage analysis of exome chip variants and cardiometabolic risk

Linkage analysis of complex traits has had limited success in identifying trait-influencing loci. Recently, coding variants have been implicated as the basis for some biomedical associations. We tested whether coding variants are the basis for linkage peaks of complex traits in 42 African-American (n = 596) and 90 Hispanic (n = 1,414) families in the Insulin Resistance Atherosclerosis Family Study (IRASFS) using Illumina HumanExome Beadchips. A total of 92,157 variants in African Americans (34%) and 81,559 (31%) in Hispanics were polymorphic and tested using two-point linkage and association analyses with 37 cardiometabolic phenotypes. In African Americans 77 LOD scores greater than 3 were observed. The highest LOD score was 4.91 with the APOE SNP rs7412 (MAF = 0.13) with plasma apolipoprotein B (ApoB). This SNP was associated with ApoB (P-value = 4 × 10(-19)) and accounted for 16.2% of the variance in African Americans. In Hispanic families, 104 LOD scores were greater than 3. The strongest evidence of linkage (LOD = 4.29) was with rs5882 (MAF = 0.46) in CETP with HDL. CETP variants were strongly associated with HDL (0.00049 < P-value <4.6 × 10(-12)), accounting for up to 4.5% of the variance. These loci have previously been shown to have effects on the biomedical traits evaluated here. Thus, evidence of strong linkage in this genome wide survey of primarily coding variants was uncommon. Loci with strong evidence of linkage was characterized by large contributions to the variance, and, in these cases, are common variants. Less compelling evidence of linkage and association was observed with additional loci that may require larger family sets to confirm.

IgA measurements in over 12 000 Swedish twins reveal sex differential heritability and regulatory locus near CD30L

In a broad attempt to improve the understanding of the genetic regulation of serum IgA levels, the heritability was estimated in over 12 000 Swedish twins, and a genome-wide association study was conducted in a subsample of 9617. Using the classical twin model the heritability was found to be significantly larger among females (61%) compared with males (21%), while contribution from shared environment (20%) was only seen for males. By modeling the genetic relationship matrix with IgA levels, we estimate that a substantial proportion (31%) of variance in IgA levels can ultimately be explained by the investigated SNPs. The genome-wide association study revealed significant association to two loci: (i) rs6928791 located on chromosome 6, 22 kb upstream of the gene SAM and SH3 domain containing 1 (SASH1) and (ii) rs13300483 on chromosome 9, situated 12 kb downstream the CD30 ligand (CD30L) encoding gene. The association to rs13300483 was replicated in two additional independent Swedish materials. The heritability of IgA levels is moderate and can partly be attributable to common variation in the CD30L locus.

Genetic determinants influencing human serum metabolome among African Americans

Phenotypes proximal to gene action generally reflect larger genetic effect sizes than those that are distant. The human metabolome, a result of multiple cellular and biological processes, are functional intermediate phenotypes proximal to gene action. Here, we present a genome-wide association study of 308 untargeted metabolite levels among African Americans from the Atherosclerosis Risk in Communities (ARIC) Study. Nineteen significant common variant-metabolite associations were identified, including 13 novel loci (p<1.6×10⁻¹⁰). These loci were associated with 7–50% of the difference in metabolite levels per allele, and the variance explained ranged from 4% to 20%. Fourteen genes were identified within the nineteen loci, and four of them contained non-synonymous substitutions in four enzyme-encoding genes (KLKB1, SIAE, CPS1, and NAT8); the other significant loci consist of eight other enzyme-encoding genes (ACE, GATM, ACY3, ACSM2B, THEM4, ADH4, UGT1A, TREH), a transporter gene (SLC6A13) and a polycystin protein gene (PKD2L1). In addition, four potential disease-associated paths were identified, including two direct longitudinal predictive relationships: NAT8 with N-acetylornithine, N-acetyl-1-methylhistidine and incident chronic kidney disease, and TREH with trehalose and incident diabetes. These results highlight the value of using endophenotypes proximal to gene function to discover new insights into biology and disease pathology.

Most contemporary GWAS studies have achieved increased power by increasing the size of the discovery sample to tens of thousands of individuals. An alternative approach for detecting the effects of novel loci is to measure phenotypes that more immediately reflect the effects of gene function. The metabolome consists of a collection of small molecules resulting from a variety of cellular and biologic processes, which can be considered intermediate phenotypes proximal to gene function. Here, we report a genome-wide association study identifying nineteen genetic loci influencing untargeted metabolomes traits among African Americans in the Atherosclerosis Risk in Communities (ARIC) Study. Fourteen genes mapped within nineteen loci, including twelve enzyme-encoding genes (KLKB1, SIAE, CPS1, NAT8, ACE, GATM, ACY3, ACSM2B, THEM4, ADH4, UGT1A and TREH), a transporter gene (SLC6A13) and a polycystin protein gene (PKD2L1). In addition, four potential disease-associated paths were identified, including two direct longitudinal predictive relationships: NAT8 with N-acetylornithine, N-acetyl-1-methylhistidine and incident chronic kidney disease, and TREH with trehalose and incident diabetes. These results highlight the value of using phenotypes proximal to gene function to promote novel gene discovery.

THOC5: a novel gene involved in HDL-cholesterol metabolism

Although numerous genes are known to regulate serum lipid traits, identified variants explain only a small proportion of the expected heritability. We intended to identify further genetic variants associated with lipid phenotypes in a self-contained population of Sorbs in Germany. We performed a genome-wide association study (GWAS) on LDL-cholesterol, HDL-cholesterol (HDL-C), and triglyceride (TG) levels in 839 Sorbs. All single-nucleotide polymorphisms with a P value <0.01 were subjected to a meta-analysis, including an independent Swedish cohort (Diabetes Genetics Initiative; n = ∼3,100). Novel association signals with the strongest effects were subjected to replication studies in an additional German cohort (Berlin, n = 2,031). In the initial GWAS in the Sorbs, we identified 14 loci associated with lipid phenotypes reaching P values <10⁻⁵ and confirmed significant effects for 18 previously reported loci. The combined meta-analysis of the three study cohorts (n(HDL) = 6041; n(LDL) = 5,995; n(TG) = 6,087) revealed a novel association for a variant in THOC5 (rs8135828) with serum HDL-C levels (P = 1.78 × 10⁻⁷; Z-score = -5.221). Consistently, the variant was also associated with circulating APOA1 levels in Sorbs. The small interfering RNA-mediated mRNA silencing of THOC5 in HepG2 cells resulted in lower mRNA levels of APOA1, SCARB1, and ABCG8 (all P < 0.05). We propose THOC5 to be a novel gene involved in the regulation of serum HDL-C levels.