Search tips
Search criteria

Results 1-25 (28)

Clipboard (0)

Select a Filter Below

more »
Year of Publication
1.  The UK Adult Twin Registry (TwinsUK Resource) 
TwinsUK is a nation-wide registry of volunteer twins in the UK, with about 12,000 registered twins (83% female, equal number of monozygotic and dizygotic twins, predominantly middle-aged and older). Over the last 20 years, questionnaire and blood/urine/tissue samples have been collected on over 7,000 subjects, as well as three comprehensive phenotyping assessments in the clinical facilities of the Department of Twin Research and Genetic Epidemiology, King’s College London. The primary focus of study has been the genetic basis of healthy ageing process and complex diseases including cardiovascular, metabolic, musculoskeletal, and ophthalmologic disorders. Alongside the detailed clinical, biochemical, behavioural, and socio-economic characterisation of the study population, the major strength of TwinsUK is availability of several ‘omics’ technologies for the participants. These include genome-wide scans of single nucleotide variants, next-generation sequencing, exome sequencing, epigenetic markers (MeDIP sequencing), gene expression arrays and RNA sequencing, telomere length measures, metabolomic profiles, and gut flora microbiomics. The scientific community now can freely access parts of the phenotype data from the ‘TwinsUK Resource’ and interested researchers are encouraged to contact us via our website ( for future collaborations.
PMCID: PMC3927054  PMID: 23088889
2.  The Rate of Nonallelic Homologous Recombination in Males Is Highly Variable, Correlated between Monozygotic Twins and Independent of Age 
PLoS Genetics  2014;10(3):e1004195.
Nonallelic homologous recombination (NAHR) between highly similar duplicated sequences generates chromosomal deletions, duplications and inversions, which can cause diverse genetic disorders. Little is known about interindividual variation in NAHR rates and the factors that influence this. We estimated the rate of deletion at the CMT1A-REP NAHR hotspot in sperm DNA from 34 male donors, including 16 monozygotic (MZ) co-twins (8 twin pairs) aged 24 to 67 years old. The average NAHR rate was 3.5×10−5 with a seven-fold variation across individuals. Despite good statistical power to detect even a subtle correlation, we observed no relationship between age of unrelated individuals and the rate of NAHR in their sperm, likely reflecting the meiotic-specific origin of these events. We then estimated the heritability of deletion rate by calculating the intraclass correlation (ICC) within MZ co-twins, revealing a significant correlation between MZ co-twins (ICC = 0.784, p = 0.0039), with MZ co-twins being significantly more correlated than unrelated pairs. We showed that this heritability cannot be explained by variation in PRDM9, a known regulator of NAHR, or variation within the NAHR hotspot itself. We also did not detect any correlation between Body Mass Index (BMI), smoking status or alcohol intake and rate of NAHR. Our results suggest that other, as yet unidentified, genetic or environmental factors play a significant role in the regulation of NAHR and are responsible for the extensive variation in the population for the probability of fathering a child with a genomic disorder resulting from a pathogenic deletion.
Author Summary
Many genetic disorders are caused by deletions of specific regions of DNA in sperm or egg cells that go on to produce a child. This can occur through ectopic homologous recombination between highly similar segments of DNA at different positions within the genome. Little is known about the differences in rates of deletion between individuals or the factors that influence this. We analysed the rate of deletion at one such section of DNA in sperm DNA from 34 male donors, including 16 monozygotic co-twins. We observed a seven-fold variation in deletion rate across individuals. Deletion rate is significantly correlated between monozygote co-twins, indicating that deletion rate is heritable. This heritability cannot be explained by age, any known genetic regulator of deletion rate, Body Mass Index, smoking status or alcohol intake. Our results suggest that other, as yet unidentified, genetic or environmental factors play a significant role in the regulation of deletion. These factors are responsible for the extensive variation in the population for the probability of fathering a child with a genomic disorder resulting from a pathogenic deletion.
PMCID: PMC3945173  PMID: 24603440
3.  Cohort Profile: TwinsUK and Healthy Ageing Twin Study 
The UK's largest registry of adult twins, or TwinsUK Registry, started in 1992 and encompasses about 12 000 volunteer twins from all over the United Kingdom. More than 70% of the registered twins have filled at least one detailed health questionnaire and about half of them undergone a baseline comprehensive assessment and two follow-up clinical evaluations. The most recent follow-up visit, known as Healthy Ageing Twin Study (HATS), involved 3125 female twins aged >40 years with at least one previous clinical assessment to enable inspection of longitudinal changes in ageing traits and their genetic and environmental components. The study benefits from several state-of-the-art OMICs studies including genome-wide association, next-generation genome and transcriptome sequencing, and epigenetic and metabolomic profiles. This makes our cohort as one of the most deeply phenotyped and genotyped in the world. Several collaborative projects in the field of epidemiology of complex disorders are ongoing in our cohort and interested researchers are encouraged to get in contact for future collaborations.
PMCID: PMC3600616  PMID: 22253318
4.  The human gut and groundwater harbor non-photosynthetic bacteria belonging to a new candidate phylum sibling to Cyanobacteria 
eLife  2013;2:e01102.
Cyanobacteria were responsible for the oxygenation of the ancient atmosphere; however, the evolution of this phylum is enigmatic, as relatives have not been characterized. Here we use whole genome reconstruction of human fecal and subsurface aquifer metagenomic samples to obtain complete genomes for members of a new candidate phylum sibling to Cyanobacteria, for which we propose the designation ‘Melainabacteria’. Metabolic analysis suggests that the ancestors to both lineages were non-photosynthetic, anaerobic, motile, and obligately fermentative. Cyanobacterial light sensing may have been facilitated by regulators present in the ancestor of these lineages. The subsurface organism has the capacity for nitrogen fixation using a nitrogenase distinct from that in Cyanobacteria, suggesting nitrogen fixation evolved separately in the two lineages. We hypothesize that Cyanobacteria split from Melainabacteria prior or due to the acquisition of oxygenic photosynthesis. Melainabacteria remained in anoxic zones and differentiated by niche adaptation, including for symbiosis in the mammalian gut.
eLife digest
Microbes are ubiquitous in the world and exist in complex communities called microbiomes that have colonized many environments, including the human gut. Until modern techniques for sequencing nucleic acids became available, many of the organisms found in these microbiomes could not be studied because they could not be cultured in the laboratory. However, advances in sequencing technology have made it possible to study the evolution and properties of these microbes, including their impact on human health.
Bacteria belonging to the phylum Cyanobacteria had a significant effect on the prehistoric Earth because they were the first organisms to produce gaseous oxygen as a byproduct of photosynthesis, and thus shaped the Earth’s oxygen-rich atmosphere. Early plants took up these bacteria in a symbiotic relationship, and plastids—the organelles in plant cells that perform photosynthesis and produce oxygen–are the descendants of Cyanobacteria.
Organisms evolutionarily related to Cyanobacteria have been found in the human gut and in various aquatic sources, but these bacteria have not been studied because it has not been possible to isolate or culture them. Now, Di Rienzi, Sharon et al. have used modern sequencing techniques to obtain complete genomes for some of these bacteria, which they assign to a new phylum called Melainabacteria.
By analyzing these genomes, Di Rienzi, Sharon et al. were able to make predictions about the cell structure and metabolic abilities of Melainabacteria. Like Cyanobacteria, they have two membranes surrounding the bacterial cell; unlike Cyanobacteria, however, they have flagella that propel them through liquid or across surfaces. Most interestingly, Melainabacteria are not able to perform photosynthesis, but instead produce energy through fermentation and release hydrogen gas that can be consumed by other microorganisms.
The genome of the bacteria isolated from water reveals that it has the capacity to fix nitrogen. Cyanobacteria can also fix atmospheric nitrogen, but the protein complexes used by the two phyla are not related, which suggests that nitrogen fixation evolved after the evolutionary divergence of Cyanobacteria and Melainabacteria.
By exploring previously published datasets of bacterial communities, Di Rienzi, Sharon et al. found that Melainabacteria are common in aquatic habitats. They are also prevalent in the guts of herbivorous mammals and humans with a predominantly vegetarian diet. Melainabacteria from the human gut also synthesize several B and K vitamins, which suggests that these bacteria are beneficial to their host because in addition to aiding with the digestion of plant fibers, they are also a source of vitamins.
PMCID: PMC3787301  PMID: 24137540
Cyanobacteria; Melainabacteria; photosynthesis; nitrogen fixation; human gut; subsurface; Human; Other
5.  Mitochondrial DNA heteroplasmy in diabetes and normal adults: role of acquired and inherited mutational patterns in twins 
Human Molecular Genetics  2012;21(19):4214-4224.
Heteroplasmy, the mixture of mitochondrial genomes (mtDNA), varies among individuals and cells. Heteroplasmy levels alter the penetrance of pathological mtDNA mutations, and the susceptibility to age-related diseases such as Parkinson's disease. Although mitochondrial dysfunction occurs in age-related type 2 diabetes mellitus (T2DM), the involvement of heteroplasmy in diabetes is unclear. We hypothesized that the heteroplasmic mutational (HM) pattern may change in T2DM. To test this, we used next-generation sequencing, i.e. massive parallel sequencing (MPS), along with PCR–cloning–Sanger sequencing to analyze HM in blood and skeletal muscle DNA samples from monozygotic (MZ) twins either concordant or discordant for T2DM. Great variability was identified in the repertoires and amounts of HMs among individuals, with a tendency towards more mutations in skeletal muscle than in blood. Whereas many HMs were unique, many were either shared among twin pairs or among tissues of the same individual, regardless of their prevalence. This suggested a heritable influence on even low abundance HMs. We found no clear differences between T2DM and controls. However, we found ∼5-fold increase of HMs in non-coding sequences implying the influence of negative selection (P < 0.001). This negative selection was evident both in moderate to highly abundant heteroplasmy (>5% of the molecules per sample) and in low abundance heteroplasmy (<5% of the molecules). Although our study found no evidence supporting the involvement of HMs in the etiology of T2DM, the twin study found clear evidence of a heritable influence on the accumulation of HMs as well as the signatures of selection in heteroplasmic mutations.
PMCID: PMC3441121  PMID: 22736028
6.  Longitudinal study of variation in body mass index in middle-aged UK females 
Age  2011;34(5):1285-1294.
The importance of changing patterns of obesity in society and its implications for public health are well recognized. However, the adult life course of body mass index (BMI) changes in individuals over time is largely unknown and has mostly been extrapolated from cross-sectional studies. The present study examines individual specific variation of BMI during a 15-year follow-up period in a community-based sample of UK females. We attempted to establish whether there is a common, generalized pattern which captures variation in BMI over time. The participants of this study belong to a prospective population cohort of British women studied intensively since 1989: the Chingford Study. The sample originally consisted of 1,003 women aged 45–68 years, who were assessed annually for BMI during follow-up period. Polynomial regression models were used to assess longitudinal BMI variation. We observed a great stability in individual BMI variation during the follow-up period, reflected by high correlations between the baseline BMI and follow-up BMI 10 and 15 years later (r = 0.876, N = 810, and r = 0.824, N = 638, respectively). We also found that three different major age-related patterns in BMI could be clearly identified: no change in 30.6% in 58% it increased and in 11.4% it decreased with age. Thus, our data suggest that individual age-related changes in BMI are very different. Therefore, simply combining all individuals into groups by any other criteria (age, sex, etc.) and overlooking the distinctive patterns of BMI change may lead to biased inferences in epidemiologic and etiologic research of the future.
Electronic supplementary material
The online version of this article (doi:10.1007/s11357-011-9299-0) contains supplementary material, which is available to authorized users.
PMCID: PMC3448995  PMID: 21853263
BMI; Follow-up, curve fitting; Age-dependent patterns; Longitudinal; Weight gain
7.  Genome-wide meta-analysis points to CTC1 and ZNF676 as genes regulating telomere homeostasis in humans 
Human Molecular Genetics  2012;21(24):5385-5394.
Leukocyte telomere length (LTL) is associated with a number of common age-related diseases and is a heritable trait. Previous genome-wide association studies (GWASs) identified two loci on chromosomes 3q26.2 (TERC) and 10q24.33 (OBFC1) that are associated with the inter-individual LTL variation. We performed a meta-analysis of 9190 individuals from six independent GWAS and validated our findings in 2226 individuals from four additional studies. We confirmed previously reported associations with OBFC1 (rs9419958 P = 9.1 × 10−11) and with the telomerase RNA component TERC (rs1317082, P = 1.1 × 10−8). We also identified two novel genomic regions associated with LTL variation that map near a conserved telomere maintenance complex component 1 (CTC1; rs3027234, P = 3.6 × 10−8) on chromosome17p13.1 and zinc finger protein 676 (ZNF676; rs412658, P = 3.3 × 10−8) on 19p12. The minor allele of rs3027234 was associated with both shorter LTL and lower expression of CTC1. Our findings are consistent with the recent observations that point mutations in CTC1 cause short telomeres in both Arabidopsis and humans affected by a rare Mendelian syndrome. Overall, our results provide novel insights into the genetic architecture of inter-individual LTL variation in the general population.
PMCID: PMC3510758  PMID: 23001564
8.  Inherited Variation in Vitamin D Genes Is Associated With Predisposition to Autoimmune Disease Type 1 Diabetes 
Diabetes  2011;60(5):1624-1631.
Vitamin D deficiency (25-hydroxyvitamin D [25(OH)D] <50 nmol/L) is commonly reported in both children and adults worldwide, and growing evidence indicates that vitamin D deficiency is associated with many extraskeletal chronic disorders, including the autoimmune diseases type 1 diabetes and multiple sclerosis.
We measured 25(OH)D concentrations in 720 case and 2,610 control plasma samples and genotyped single nucleotide polymorphisms from seven vitamin D metabolism genes in 8,517 case, 10,438 control, and 1,933 family samples. We tested genetic variants influencing 25(OH)D metabolism for an association with both circulating 25(OH)D concentrations and disease status.
Type 1 diabetic patients have lower circulating levels of 25(OH)D than similarly aged subjects from the British population. Only 4.3 and 18.6% of type 1 diabetic patients reached optimal levels (≥75 nmol/L) of 25(OH)D for bone health in the winter and summer, respectively. We replicated the associations of four vitamin D metabolism genes (GC, DHCR7, CYP2R1, and CYP24A1) with 25(OH)D in control subjects. In addition to the previously reported association between type 1 diabetes and CYP27B1 (P = 1.4 × 10−4), we obtained consistent evidence of type 1 diabetes being associated with DHCR7 (P = 1.2 × 10−3) and CYP2R1 (P = 3.0 × 10−3).
Circulating levels of 25(OH)D in children and adolescents with type 1 diabetes vary seasonally and are under the same genetic control as in the general population but are much lower. Three key 25(OH)D metabolism genes show consistent evidence of association with type 1 diabetes risk, indicating a genetic etiological role for vitamin D deficiency in type 1 diabetes.
PMCID: PMC3292339  PMID: 21441443
9.  Identification of an imprinted master trans-regulator at the KLF14 locus related to multiple metabolic phenotypes 
Nature genetics  2011;43(6):561-564.
Genome-wide association studies have identified many genetic variants associated with complex traits. However, at only a minority of loci have the molecular mechanisms mediating these associations been characterized. In parallel, whilst cis-regulatory patterns of gene expression have been extensively explored, the identification of trans-regulatory effects in humans has attracted less attention. We demonstrate that the Type 2 diabetes and HDL-cholesterol associated cis-acting eQTL of the maternally-expressed transcription factor KLF14 acts as a master trans-regulator of adipose gene expression. Expression levels of genes regulated by this trans-eQTL are highly-correlated with concurrently-measured metabolic traits, and a subset of the trans-genes harbor variants directly-associated with metabolic phenotypes. This trans-eQTL network provides a mechanistic understanding of the effect of the KLF14 locus on metabolic disease risk, providing a potential model for other complex traits.
PMCID: PMC3192952  PMID: 21572415
10.  Novel genes for QTc interval. How much heritability is explained, and how much is left to find? 
Genome Medicine  2010;2(5):35.
The corrected QT (QTc) interval is a complex quantitative trait, believed to be influenced by several genetic and environmental factors. It is a strong prognostic indicator of cardiovascular mortality in patients with and without cardiac disease. More than 700 mutations have been described in 12 genes (LQT1-LQT12) involved in congenital long QT syndrome. However, the heritability (genetic contribution) of QTc interval in the general population cannot be adequately explained by these long QT syndrome genes. In order to further investigate the genetic architecture underlying QTc interval in the general population, genome-wide association studies, in which up to one million single nucleotide polymorphisms are assayed in thousands of individuals, are now being employed and have already led to the discovery of variants in seven novel loci and five loci that are known to cause congenital long or short QT syndrome. Here we show that a combined risk score using 11 of these loci explains about 10% of the heritability of QTc. Additional discovery of both common and rare variants will yield further etiological insight and accelerate clinical applications.
PMCID: PMC2887079  PMID: 20519034
11.  The effect of genome-wide association scan quality control on imputation outcome for common variants 
Imputation is an extremely valuable tool in conducting and synthesising genome-wide association studies (GWASs). Directly typed SNP quality control (QC) is thought to affect imputation quality. It is, therefore, common practise to use quality-controlled (QCed) data as an input for imputing genotypes. This study aims to determine the effect of commonly applied QC steps on imputation outcomes. We performed several iterations of imputing SNPs across chromosome 22 in a dataset consisting of 3177 samples with Illumina 610k (Illumina, San Diego, CA, USA) GWAS data, applying different QC steps each time. The imputed genotypes were compared with the directly typed genotypes. In addition, we investigated the correlation between alternatively QCed data. We also applied a series of post-imputation QC steps balancing elimination of poorly imputed SNPs and information loss. We found that the difference between the unQCed data and the fully QCed data on imputation outcome was minimal. Our study shows that imputation of common variants is generally very accurate and robust to GWAS QC, which is not a major factor affecting imputation outcome. A minority of common-frequency SNPs with particular properties cannot be accurately imputed regardless of QC stringency. These findings may not generalise to the imputation of low frequency and rare variants.
PMCID: PMC3083623  PMID: 21267008
genome-wide association study; imputation; quality control; single nucleotide polymorphism
12.  A Common Variant in the Telomerase RNA Component Is Associated with Short Telomere Length 
PLoS ONE  2010;5(9):e13048.
Telomeres shorten as cells divide. This shortening is compensated by the enzyme telomerase. We evaluated the effect of common variants in the telomerase RNA component (TERC) gene on telomere length (TL) in the population-based Health Aging and Body Composition (Health ABC) Study and in two replication samples (the TwinsUK Study and the Amish Family Osteoporosis Study, AFOS).
Five variants were identified in the TERC region by sequence analysis and only one SNP was common (rs2293607, G/A). The frequency of the G allele was 0.26 and 0.07 in white and black, respectively. Testing for association between TL and rs2293607 was performed using linear regression models or variance component analysis conditioning on relatedness among subjects.
The adjusted mean TL was significantly shorter in 665 white carriers of the G allele compared to 887 non-carriers from the Health ABC Study (4.69±0.05 kbp vs. 4.86±0.04 kbp, measured by quantitative PCR, p = 0.005). This association was replicated in another white sample from the TwinsUK Study (6.90±0.03 kbp in 301 carriers compared to 7.06±0.03 kbp in 395 non-carriers, measured by Southern blots, p = 0.009). A similar pattern of association was observed in whites from the family-based AFOS and blacks from the Health ABC cohort, although not statistically significant, possibly due to the lower allele frequency in these populations. Combined analysis using 2,953 white subjects from 3 studies showed a significant association between TL and rs2293607 (β = −0.19±0.04 kbp, p = 0.001).
Our study shows a significant association between a common variant in TERC and TL in humans, suggesting that TERC may play a role in telomere homeostasis.
PMCID: PMC2946401  PMID: 20885959
13.  Evidence that bone mineral density plays a role in degenerative disc disease: the UK Twin Spine Study 
Annals of the Rheumatic Diseases  2010;69(12):2102-2106.
Osteoarthritis (OA) and osteoporosis are often considered to lie at opposite ends of a spectrum of bone phenotypes. Lumbar degenerative disc disease (LDD) may be associated with low back pain (LBP) and is similar in many ways to OA. LDD is reported in small studies to be associated with increased spine bone mineral density (BMD). The present work aimed to confirm this association in a large population sample using MRI and explore the relationship further, in particular to determine whether it is mediated genetically.
A population based sample (N=908, age range 32–74 years) of UK female twins having MRI of the lumbar spine was used in this study. LDD traits and summary measures and their relationship with BMD at the lumbar spine and hip were examined using multivariate multiple regression and maximum likelihood based variance decomposition.
There was a significant positive correlation between LDD and BMD at the lumbar spine and hip, which remained significant after adjustment for confounders. Both traits were highly heritable and the associations between them were mediated genetically.
A clear, significant and independent association of BMD at hip and lumbar spine with LDD was found which is, in part, genetically mediated. The association with the non-axial site, the hip, is of particular interest and suggests a systemic bone effect. This should encourage the search for pleiotropic genes to help in the understanding of the bone–cartilage relationship. Moreover, genetic variants identified could provide novel therapeutic targets in the management of LBP.
PMCID: PMC3002767  PMID: 20570838
14.  Digital Quantification of Human Eye Color Highlights Genetic Association of Three New Loci 
PLoS Genetics  2010;6(5):e1000934.
Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.
Author Summary
We measured human eye color to hue and saturation values from high-resolution, digital, full-eye photographs of several thousand Dutch Europeans. This quantitative approach, which is extremely cost-effective, portable, and time efficient, revealed that human eye color varies along more dimensions than the one represented by the blue-green-brown categories studied previously. Our work represents the first genome-wide study of quantitative human eye color. We clearly identified 3 new loci, LYST, 17q25.3, TTC3/DSCR9, in contributing to the natural and subtle eye color variation along multiple dimensions, providing new leads towards a more detailed understanding of the genetic basis of human eye color. Our quantitative prediction model explained over 50% of eye color variance, representing the highest accuracy achieved so far in genomic prediction of human complex and quantitative traits, with relevance for future forensic applications.
PMCID: PMC2865509  PMID: 20463881
15.  Association of FTO gene variants with body composition in UK twins 
Annals of Human Genetics  2012;76(5):333-341.
The association of FTO gene variants with body mass index (BMI) and other obesity characteristics is well established. However, uncertainties remain whether the association is present only in young populations and whether it is attributable to body fat mass specifically. We aimed to clarify these two questions in a large sample (N= 4,523 individuals) of middle-aged and older (range 40–80 years) British female twins. The women were assessed for BMI, waist and hip circumference, total lean (LBM) and fat (FBM) body mass. Since the majority of FTO association signals have been reported in a haploblock bordering 52,355–52,408 kb (on chromosome 16q12.2), we examined five genotyped and 43 imputed SNPs mapped to this block. Canonical correlation and other association analyses showed significant and consistent association between the selected SNP and studied body composition phenotypes, with p-values reaching p= 0.000004. Of particular interest, in addition to the expected significant associations between FTO variants and FBM, we also identified significant associations with LBM. These results suggest that the association between FTO variants and body composition phenotypes is present across a wide range of ages, and that FTO appears primarily to affect the amount of body soft tissue, influencing both fat and lean mass.
PMCID: PMC3532623  PMID: 22817777
BMI; lean and fat body mass; waist and hip circumference; SNP; association analysis
16.  Efficacy and safety of strontium ranelate in the treatment of knee osteoarthritis: results of a double-blind, randomised placebo-controlled trial 
Annals of the Rheumatic Diseases  2012;72(2):179-186.
Strontium ranelate is currently used for osteoporosis. The international, double-blind, randomised, placebo-controlled Strontium ranelate Efficacy in Knee OsteoarthrItis triAl evaluated its effect on radiological progression of knee osteoarthritis.
Patients with knee osteoarthritis (Kellgren and Lawrence grade 2 or 3, and joint space width (JSW) 2.5–5 mm) were randomly allocated to strontium ranelate 1 g/day (n=558), 2 g/day (n=566) or placebo (n=559). The primary endpoint was radiographical change in JSW (medial tibiofemoral compartment) over 3 years versus placebo. Secondary endpoints included radiological progression, Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) score, and knee pain. The trial is registered (ISRCTN41323372).
The intention-to-treat population included 1371 patients. Treatment with strontium ranelate was associated with smaller degradations in JSW than placebo (1 g/day: −0.23 (SD 0.56) mm; 2 g/day: −0.27 (SD 0.63) mm; placebo: −0.37 (SD 0.59) mm); treatment-placebo differences were 0.14 (SE 0.04), 95% CI 0.05 to 0.23, p<0.001 for 1 g/day and 0.10 (SE 0.04), 95% CI 0.02 to 0.19, p=0.018 for 2 g/day. Fewer radiological progressors were observed with strontium ranelate (p<0.001 and p=0.012 for 1 and 2 g/day). There were greater reductions in total WOMAC score (p=0.045), pain subscore (p=0.028), physical function subscore (p=0.099) and knee pain (p=0.065) with strontium ranelate 2 g/day. Strontium ranelate was well tolerated.
Treatment with strontium ranelate 1 and 2 g/day is associated with a significant effect on structure in patients with knee osteoarthritis, and a beneficial effect on symptoms for strontium ranelate 2 g/day.
PMCID: PMC3599139  PMID: 23117245
Knee Osteoarthritis; Osteoarthritis; Outcomes research
17.  Genome Wide Association Analysis of a Founder Population Identified TAF3 as a Gene for MCHC in Humans 
PLoS ONE  2013;8(7):e69206.
The red blood cell related traits are highly heritable but their genetics are poorly defined. Only 5–10% of the total observed variance is explained by the genetic loci found to date, suggesting that additional loci should be searched using approaches alternative to large meta analysis. GWAS (Genome Wide Association Study) for red blood cell traits in a founder population cohort from Northern Italy identified a new locus for mean corpuscular hemoglobin concentration (MCHC) in the TAF3 gene. The association was replicated in two cohorts (rs1887582, P = 4.25E–09). TAF3 encodes a transcription cofactor that participates in core promoter recognition complex, and is involved in zebrafish and mouse erythropoiesis. We show here that TAF3 is required for transcription of the SPTA1 gene, encoding alpha spectrin, one of the proteins that link the plasma membrane to the actin cytoskeleton. Mutations in SPTA1 are responsible for hereditary spherocytosis, a monogenic disorder of MCHC, as well as for the normal MCHC level. Based on our results, we propose that TAF3 is required for normal erythropoiesis in human and that it might have a role in controlling the ratio between hemoglobin (Hb) and cell volume and in the dynamics of RBC maturation in healthy individuals. Finally, TAF3 represents a potential candidate or a modifier gene for disorders of red cell membrane.
PMCID: PMC3729833  PMID: 23935956
18.  Congruence of Additive and Non-Additive Effects on Gene Expression Estimated from Pedigree and SNP Data 
PLoS Genetics  2013;9(5):e1003502.
There is increasing evidence that heritable variation in gene expression underlies genetic variation in susceptibility to disease. Therefore, a comprehensive understanding of the similarity between relatives for transcript variation is warranted—in particular, dissection of phenotypic variation into additive and non-additive genetic factors and shared environmental effects. We conducted a gene expression study in blood samples of 862 individuals from 312 nuclear families containing MZ or DZ twin pairs using both pedigree and genotype information. From a pedigree analysis we show that the vast majority of genetic variation across 17,994 probes is additive, although non-additive genetic variation is identified for 960 transcripts. For 180 of the 960 transcripts with non-additive genetic variation, we identify expression quantitative trait loci (eQTL) with dominance effects in a sample of 339 unrelated individuals and replicate 31% of these associations in an independent sample of 139 unrelated individuals. Over-dominance was detected and replicated for a trans association between rs12313805 and ETV6, located 4MB apart on chromosome 12. Surprisingly, only 17 probes exhibit significant levels of common environmental effects, suggesting that environmental and lifestyle factors common to a family do not affect expression variation for most transcripts, at least those measured in blood. Consistent with the genetic architecture of common diseases, gene expression is predominantly additive, but a minority of transcripts display non-additive effects.
Author Summary
Gene expression levels are known to influence common disease susceptibility in humans, with GWAS significant SNPs frequently found in regulatory regions. The expression levels of most genes are influenced by genetic variants, often located close to the gene itself. Expression Quantitative Trait Loci (eQTL) mapping studies have been very successful in identifying SNPs associated with expression levels; however, little is currently known about the extent of additive and non-additive genetic variance and the role of common environment on gene expression. Here we report a comprehensive study of the sources of genetic and non-genetic variation for gene expression levels using both pedigree and genotype information. We show that the majority of transcripts exhibit only additive genetic variance with congruence from independent methods using pedigree and genotype approaches. However, there are a small number of probes whose expression levels are influenced by non-additive genetic variance. For some of these probes we identify SNPs acting in a dominant and over-dominant manner that replicate in an independent sample. Surprisingly, only 17 probes exhibit significant levels of common environmental effects, suggesting that environmental and lifestyle factors common to a family do not affect expression variation for most transcripts, at least those measured in blood.
PMCID: PMC3656157  PMID: 23696747
19.  First all-in-one diagnostic tool for DNA intelligence: genome-wide inference of biogeographic ancestry, appearance, relatedness, and sex with the Identitas v1 Forensic Chip 
When a forensic DNA sample cannot be associated directly with a previously genotyped reference sample by standard short tandem repeat profiling, the investigation required for identifying perpetrators, victims, or missing persons can be both costly and time consuming. Here, we describe the outcome of a collaborative study using the Identitas Version 1 (v1) Forensic Chip, the first commercially available all-in-one tool dedicated to the concept of developing intelligence leads based on DNA. The chip allows parallel interrogation of 201,173 genome-wide autosomal, X-chromosomal, Y-chromosomal, and mitochondrial single nucleotide polymorphisms for inference of biogeographic ancestry, appearance, relatedness, and sex. The first assessment of the chip’s performance was carried out on 3,196 blinded DNA samples of varying quantities and qualities, covering a wide range of biogeographic origin and eye/hair coloration as well as variation in relatedness and sex. Overall, 95 % of the samples (N = 3,034) passed quality checks with an overall genotype call rate >90 % on variable numbers of available recorded trait information. Predictions of sex, direct match, and first to third degree relatedness were highly accurate. Chip-based predictions of biparental continental ancestry were on average ~94 % correct (further support provided by separately inferred patrilineal and matrilineal ancestry). Predictions of eye color were 85 % correct for brown and 70 % correct for blue eyes, and predictions of hair color were 72 % for brown, 63 % for blond, 58 % for black, and 48 % for red hair. From the 5 % of samples (N = 162) with <90 % call rate, 56 % yielded correct continental ancestry predictions while 7 % yielded sufficient genotypes to allow hair and eye color prediction. Our results demonstrate that the Identitas v1 Forensic Chip holds great promise for a wide range of applications including criminal investigations, missing person investigations, and for national security purposes.
Electronic supplementary material
The online version of this article (doi:10.1007/s00414-012-0788-1) contains supplementary material, which is available to authorized users.
PMCID: PMC3631519  PMID: 23149900
DNA intelligence; Forensic DNA phenotyping; SNP; Prediction; Relatedness; Kinship; Ancestry; Eye color; Hair color; Sex
20.  A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans 
PLoS Genetics  2012;8(9):e1002932.
Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.
Author Summary
Monozygotic twins look more alike than dizygotic twins or other siblings, and siblings in turn look more alike than unrelated individuals, indicating that human facial morphology has a strong genetic component. We quantitatively assessed human facial shape phenotypes based on statistical shape analyses of facial landmarks obtained from three-dimensional magnetic resonance images of the head. These phenotypes turned out to be highly promising for studying the genetic basis of human facial variation in that they showed high heritability in our twin data. A subsequent genome-wide association study (GWAS) identified five candidate genes affecting facial shape in Europeans: PRDM16, PAX3, TP63, C5orf50, and COL17A1. In addition, our data suggest that genetic variants associated with NSCL/P also influence normal facial shape variation. Overall, this study provides novel and confirmatory links between common DNA variants and normal variation in human facial morphology. Our results also suggest that the high heritability of facial phenotypes seems to be explained by a large number of DNA variants with relatively small individual effect size, a phenomenon well known for other complex human traits, such as adult body height.
PMCID: PMC3441666  PMID: 23028347
21.  Extent, Causes, and Consequences of Small RNA Expression Variation in Human Adipose Tissue 
PLoS Genetics  2012;8(5):e1002704.
Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.
Author Summary
Genetic information is transmitted to the cell only through RNA molecules. A special class of RNAs is comprised of the small (up to 30 nucleotide) ones, known to be potent regulators of various cellular processes. At the same time, they have not been as widely studied as messenger RNAs—we do not know how much variation in their sequence and expression level occurs naturally in human populations or how this variability influences other traits. We measured small RNA levels and genetic variability in fat tissue from 131 individuals by high-throughput sequencing. We could associate the expression levels with genetic background of the individuals, as well as changes in metabolic traits. Surprisingly, we found no large scale influence of small RNA variation on mRNA levels, their main regulatory target. Overall, our study is the first to give a quantitative picture of the naturally occurring variation in these important regulatory molecules in human fat tissue.
PMCID: PMC3349731  PMID: 22589741
22.  Genetic polymorphisms of innate immunity-related inflammatory pathways and their association with factors related to type 2 diabetes 
BMC Medical Genetics  2011;12:95.
Type 2 diabetes mellitus (T2DM) has been linked to a state of pre-clinical chronic inflammation resulting from abnormalities in the innate immune pathway. Serum levels of pro-inflammatory cytokines and acute-phase proteins, collectively known as 'inflammatory network', are elevated in the pre-, or early, stages of T2DM and increase with disease progression. Genetic variation can affect the innate immune response to certain environmental factors, and may, therefore, determine an individual's lifetime risk of disease.
We conducted a cross-sectional study in 6,720 subjects from the TwinsUK Registry to evaluate the association between 18 single nucleotide polymorphisms (SNPs) in five genes (TLR4, IL1A, IL6, TNFA, and CRP) along the innate immunity-related inflammatory pathway and biomarkers of predisposition to T2DM [fasting insulin and glucose, HDL- and LDL- cholesterols, triglycerides (TGs), amyloid-A, sensitive C-reactive protein (sCRP) and vitamin D binding protein (VDBP) and body mass index (BMI)].
Of 18 the SNPs examined for their association with nine metabolic phenotypes of interest, six were significantly associated with five metabolic phenotypes (Bonferroni correction, P ≤ 0.0027). Fasting insulin was associated with SNPs in IL6 and TNFA, serum HDL-C with variants of TNFA and CRP and serum sCRP level with SNPs in CRP. Cross-correlation analysis among the different metabolic factors related to risk of T2DM showed several significant associations. For example, BMI was directly correlated with glucose (r = 0.11), insulin (r = 0.15), sCRP (r = 0.23), LDL-C (r = 0.067) and TGs (r = 0.18) but inversely with HDL-C (r = -0.14). sCRP was also positively correlated (P < 0.0001) with insulin (r = 0.17), amyloid-A (r = 0.39), TGs (r = 0.26), and VDBP (r = 0.36) but inversely with HDL-C (r = -0.12).
Genetic variants in the innate immunity pathway and its related inflammatory cascade is associated with some metabolic risk factors for T2DM; an observation that may provide a rationale for further studying their role as biomarkers for disease early risk prediction.
PMCID: PMC3161932  PMID: 21756351
23.  The Architecture of Gene Regulatory Variation across Multiple Human Tissues: The MuTHER Study 
PLoS Genetics  2011;7(2):e1002003.
While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.
Author Summary
Regulation of gene expression is a fundamental cellular process determining a large proportion of the phenotypic variance. Previous studies have identified genetic loci influencing gene expression levels (eQTLs), but the complexity of their tissue-specific properties has not yet been well-characterized. In this study, we perform cis-eQTL analysis in a unique matched co-twin design for three human tissues derived simultaneously from the same set of individuals. The study design allows validation of the substantial discoveries we make in each tissue. We explore in depth the tissue-dependent features of regulatory variants and estimate the proportions of shared and specific effects. We use continuous measures of eQTL sharing to circumvent the statistical power limitations of comparing direct overlap of eQTLs in multiple tissues. In this framework, we demonstrate that 30% of eQTLs are shared among tissues, while 29% are exclusively tissue-specific. Furthermore, we show that the fold change in expression between eQTL genotypic classes differs between tissues. Even among shared eQTLs, we report a substantial proportion (10%–20%) of significant tissue differences in magnitude of these effects. The complexities we highlight here are essential for understanding the impact of regulatory variants on complex traits.
PMCID: PMC3033383  PMID: 21304890
24.  An Integration of Genome-Wide Association Study and Gene Expression Profiling to Prioritize the Discovery of Novel Susceptibility Loci for Osteoporosis-Related Traits 
PLoS Genetics  2010;6(6):e1000977.
Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.6×10−8), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6×10−13; SOX6, p = 6.4×10−10) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.
Author Summary
BMD and hip geometry are two major predictors of osteoporotic fractures, the most severe consequence of osteoporosis in elderly persons. We performed sex-specific genome-wide association studies (GWAS) for BMD at the lumbar spine and femor neck skeletal sites as well as hip geometric indices (NSA, NL, and NW) in the Framingham Osteoporosis Study and then replicated the top findings in two independent studies. Three novel loci were significant: in women, including chromosome 1p13.2 (RAP1A) for NW; in men, 2q11.2 (TBC1D8) for NSA and 18q11.2 (OSBPL1A) for NW. We confirmed a previously reported region on 8q24.12 (TNFRSF11B/OPG) for lumbar spine BMD in women. In addition, we integrated GWAS signals with eQTL in several tissues and publicly available expression signature profiling in cellular and whole-animal models, and prioritized 16 candidate genes/loci based on their potential involvement in skeletal metabolism. Among three prioritized loci (GPR177, SOX6, and CASR genes) associated with BMD in women, GPR177 and SOX6 have been successfully replicated later in a large-scale meta-analysis, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of using expression profiling to support the candidacy of suggestive GWAS signals that may contain important genes of interest.
PMCID: PMC2883588  PMID: 20548944
25.  Common Genetic Variants near the Brittle Cornea Syndrome Locus ZNF469 Influence the Blinding Disease Risk Factor Central Corneal Thickness 
PLoS Genetics  2010;6(5):e1000947.
Central corneal thickness (CCT), one of the most highly heritable human traits (h2 typically>0.9), is important for the diagnosis of glaucoma and a potential risk factor for glaucoma susceptibility. We conducted genome-wide association studies in five cohorts from Australia and the United Kingdom (total N = 5058). Three cohorts were based on individually genotyped twin collections, with the remaining two cohorts genotyped on pooled samples from singletons with extreme trait values. The pooled sample findings were validated by individual genotyping the pooled samples together with additional samples also within extreme quantiles. We describe methods for efficient combined analysis of the results from these different study designs. We have identified and replicated quantitative trait loci on chromosomes 13 and 16 for association with CCT. The locus on chromosome 13 (nearest gene FOXO1) had an overall meta-analysis p-value for all the individually genotyped samples of 4.6×10−10. The locus on chromosome 16 was associated with CCT with p = 8.95×10−11. The nearest gene to the associated chromosome 16 SNPs was ZNF469, a locus recently implicated in Brittle Cornea Syndrome (BCS), a very rare disorder characterized by abnormal thin corneas. Our findings suggest that in addition to rare variants in ZNF469 underlying CCT variation in BCS patients, more common variants near this gene may contribute to CCT variation in the general population.
Author Summary
Central corneal thickness (CCT) is an important eye measurement. It has been considered as a prognosticator for the development of glaucoma, with a thin cornea potentially increasing the risk of developing a subtype known as open-angle glaucoma. CCT is highly heritable, yet its genetic determinants are poorly characterized. We have revealed two loci near gene FOXO1 and ZNF469 associated with CCT in this multi-stage genome-wide association study examining over 5,000 samples. It is of particular interest that, while rare mutations in ZNF469 cause Brittle Cornea Syndrome, more common variants near this gene also contribute to CCT variation in the general population. Furthermore, given the relation between CCT and glaucoma, results from our CCT studies will implement the search for the disease-susceptibility genes of glaucoma.
PMCID: PMC2869325  PMID: 20485516

Results 1-25 (28)