Protein coding variation in the J:ARC and J:DO outbred laboratory mouse stocks provides a molecular basis for distinct research applications

Outbred laboratory mice (Mus musculus) are readily available and have high fecundity, making them a popular choice in biomedical research, especially toxicological and pharmacological applications. Direct high throughput genome sequencing (HTS) of these widely used research animals is an important genetic quality control measure that enhances research reproducibility. HTS data have been used to confirm the common origin of outbred stocks and to molecularly define distinct outbred populations. But these data have also revealed unexpected population structure and homozygosity in some populations; genetic features that emerge when outbred stocks are not properly maintained. We used exome sequencing to discover and interrogate protein-coding variation in a newly established population of Swiss-derived outbred stock (J:ARC) that is closely related to other, commonly used CD-1 outbred populations. We used these data to describe the genetic architecture of the J:ARC population including heterozygosity, minor allele frequency, LD decay, and we defined novel, protein-coding sequence variation. These data reveal the expected genetic architecture for a properly maintained outbred stock and provide a basis for the on-going genetic quality control. We also compared these data to protein-coding variation found in a multiparent outbred stock, the Diversity Outbred (J:DO). We found that the more recently derived, multiparent outbred stock has significantly higher interindividual variability, greater overall genetic variation, higher heterozygosity, and fewer novel variants than the Swiss-derived J:ARC stock. However, among the novel variants found in the J:DO stock, significantly more are predicted to be protein-damaging. The fact that individuals from this population can tolerate a higher load of potentially damaging variants highlights the buffering effects of allelic diversity and the differing selective pressures in these stocks. While both outbred stocks offer significant individual heterozygosity, our data provide a molecular basis for their intended applications, where the J:DO are best suited for studies requiring maximum, population-level genetic diversity and power for mapping, while the J:ARC are best suited as a general-purpose outbred stock with robust fecundity, relatively low allelic diversity, and less potential for extreme phenotypic variability.

Association of the IGF1 gene with fasting insulin levels

Insulin-like growth factor 1 (IGF-I) has been associated with insulin resistance. Genome-wide association studies (GWASs) of fasting insulin (FI) identified single-nucleotide variants (SNVs) near the IGF1 gene, raising two hypotheses: (1) these associations are mediated by IGF-I levels and (2) these noncoding variants either tag other functional variants in the region or are directly functional. In our study, analyses including 5141 individuals from population-based cohorts suggest that FI associations near IGF1 are not mediated by IGF-I. Analyses of targeted sequencing data in 3539 individuals reveal a large number of novel rare variants at the IGF1 locus and show a FI association with a subset of rare nonsynonymous variants (PSKAT=5.7 × 10(-4)). Conditional analyses suggest that this association is partly explained by the GWAS signal and the presence of a residual independent rare variant effect (Pconditional=0.019). Annotation using ENCODE data suggests that the GWAS variants may have a direct functional role in insulin biology. In conclusion, our study provides insight into variation present at the IGF1 locus and into the genetic architecture underlying FI levels, suggesting that FI associations of SNVs near IGF1 are not mediated by IGF-I and suggesting a role for both rare nonsynonymous and common functional variants in insulin biology.

A Dominant Mutation in Human RAD51 Reveals Its Function in DNA Interstrand Crosslink Repair Independent of Homologous Recombination

Repair of DNA interstrand crosslinks requires action of multiple DNA repair pathways, including homologous recombination. Here, we report a de novo heterozygous T131P mutation in RAD51/FANCR, the key recombinase essential for homologous recombination, in a patient with Fanconi anemia-like phenotype. In vitro, RAD51-T131P displays DNA-independent ATPase activity, no DNA pairing capacity, and a co-dominant-negative effect on RAD51 recombinase function. However, the patient cells are homologous recombination proficient due to the low ratio of mutant to wild-type RAD51 in cells. Instead, patient cells are sensitive to crosslinking agents and display hyperphosphorylation of Replication Protein A due to increased activity of DNA2 and WRN at the DNA interstrand crosslinks. Thus, proper RAD51 function is important during DNA interstrand crosslink repair outside of homologous recombination. Our study provides a molecular basis for how RAD51 and its associated factors may operate in a homologous recombination-independent manner to maintain genomic integrity.

Association of levels of fasting glucose and insulin with rare variants at the chromosome 11p11.2-MADD locus: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study

BACKGROUND

Common variation at the 11p11.2 locus, encompassing MADD, ACP2, NR1H3, MYBPC3, and SPI1, has been associated in genome-wide association studies with fasting glucose and insulin (FI). In the Cohorts for Heart and Aging Research in Genomic Epidemiology Targeted Sequencing Study, we sequenced 5 gene regions at 11p11.2 to identify rare, potentially functional variants influencing fasting glucose or FI levels.

METHODS AND RESULTS

Sequencing (mean depth, 38×) across 16.1 kb in 3566 individuals without diabetes mellitus identified 653 variants, 79.9% of which were rare (minor allele frequency <1%) and novel. We analyzed rare variants in 5 gene regions with FI or fasting glucose using the sequence kernel association test. At NR1H3, 53 rare variants were jointly associated with FI (P=2.73×10(-3)); of these, 7 were predicted to have regulatory function and showed association with FI (P=1.28×10(-3)). Conditioning on 2 previously associated variants at MADD (rs7944584, rs10838687) did not attenuate this association, suggesting that there are >2 independent signals at 11p11.2. One predicted regulatory variant, chr11:47227430 (hg18; minor allele frequency=0.00068), contributed 20.6% to the overall sequence kernel association test score at NR1H3, lies in intron 2 of NR1H3, and is a predicted binding site for forkhead box A1 (FOXA1), a transcription factor associated with insulin regulation. In human HepG2 hepatoma cells, the rare chr11:47227430 A allele disrupted FOXA1 binding and reduced FOXA1-dependent transcriptional activity.

CONCLUSIONS

Sequencing at 11p11.2-NR1H3 identified rare variation associated with FI. One variant, chr11:47227430, seems to be functional, with the rare A allele reducing transcription factor FOXA1 binding and FOXA1-dependent transcriptional activity.

Genetic determinants of age-related macular degeneration in diverse populations from the PAGE study

PURPOSE

Substantial progress has been made in identifying susceptibility variants for AMD in European populations; however, few studies have been conducted to understand the role these variants play in AMD risk in diverse populations. The present study aims to examine AMD risk across diverse populations in known and suspected AMD complement factor and lipid-related loci.

METHODS

Targeted genotyping was performed across study sites for AMD and lipid trait-associated single nucleotide polymorphism (SNPs). Genetic association tests were performed at individual sites and then meta-analyzed using logistic regression assuming an additive genetic model stratified by self-described race/ethnicity. Participants included cases with early or late AMD and controls with no signs of AMD as determined by fundus photography. Populations included in this study were European Americans, African Americans, Mexican Americans, and Singaporeans from the Population Architecture using Genomics and Epidemiology (PAGE) study.

RESULTS

Index variants of AMD, rs1061170 (CFH) and rs10490924 (ARMS2), were associated with AMD at P=3.05×10(-8) and P=6.36×10(-6), respectively, in European Americans. In general, none of the major AMD index variants generalized to our non-European populations with the exception of rs10490924 in Mexican Americans at an uncorrected P value<0.05. Four lipid-associated SNPS (LPL rs328, TRIB1 rs6987702, CETP rs1800775, and KCTD10/MVK rs2338104) were associated with AMD in African Americans and Mexican Americans (P<0.05), but these associations did not survive strict corrections for multiple testing.

CONCLUSIONS

While most associations did not generalize in the non-European populations, variants within lipid-related genes were found to be associated with AMD. This study highlights the need for larger well-powered studies in non-European populations.

Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study

INTRODUCTION

We have previously shown that a tag single nucleotide polymorphism (rs10235235), which maps to the CYP3A locus (7q22.1), was associated with a reduction in premenopausal urinary estrone glucuronide levels and a modest reduction in risk of breast cancer in women age ≤50 years.

METHODS

We further investigated the association of rs10235235 with breast cancer risk in a large case control study of 47,346 cases and 47,570 controls from 52 studies participating in the Breast Cancer Association Consortium. Genotyping of rs10235235 was conducted using a custom Illumina Infinium array. Stratified analyses were conducted to determine whether this association was modified by age at diagnosis, ethnicity, age at menarche or tumor characteristics.

RESULTS

We confirmed the association of rs10235235 with breast cancer risk for women of European ancestry but found no evidence that this association differed with age at diagnosis. Heterozygote and homozygote odds ratios (ORs) were OR = 0.98 (95% CI 0.94, 1.01; P = 0.2) and OR = 0.80 (95% CI 0.69, 0.93; P = 0.004), respectively (P(trend) = 0.02). There was no evidence of effect modification by tumor characteristics. rs10235235 was, however, associated with age at menarche in controls (P(trend) = 0.005) but not cases (P(trend) = 0.97). Consequently the association between rs10235235 and breast cancer risk differed according to age at menarche (P(het) = 0.02); the rare allele of rs10235235 was associated with a reduction in breast cancer risk for women who had their menarche age ≥15 years (OR(het) = 0.84, 95% CI 0.75, 0.94; OR(hom) = 0.81, 95% CI 0.51, 1.30; P(trend) = 0.002) but not for those who had their menarche age ≤11 years (OR(het) = 1.06, 95% CI 0.95, 1.19, OR(hom) = 1.07, 95% CI 0.67, 1.72; P(trend) = 0.29).

CONCLUSIONS

To our knowledge rs10235235 is the first single nucleotide polymorphism to be associated with both breast cancer risk and age at menarche consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk. These associations are likely mediated via an effect on circulating hormone levels.

Large-scale genotyping identifies 41 new loci associated with breast cancer risk

Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ∼9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10(-8)). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.

Seven new loci associated with age-related macular degeneration

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Insights into the genetic architecture of early stage age-related macular degeneration: a genome-wide association study meta-analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10(-8)). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4-1.88, P = 2.7 × 10(-10), and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29-1.68, P = 4.9 × 10(-9)). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10(-4); tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Association of genetic variants on 8p21 and 4q12 with age-related macular degeneration in Asian populations

PURPOSE

To evaluate the association of genetic variants at chromosomes 8p21 and 4q12 with the risk of developing AMD and its two main subtypes, choroidal neovascular membrane (CNV) and polypoidal choroidal vasculopathy (PCV), in Asian populations.

METHODS

The study population comprised 2360 patients with neovascular AMD (1013 typical AMD-CNV and 1282 PCV), and 3598 controls from four independent cohorts, two of Japanese (n = 4859) and two of Chinese (n = 1099) ethnicity. We performed a meta-analysis in case-control studies of two reported single nucleotide polymorphisms (SNPs) (rs13278062 at TNFRSF10A-LOC389641 on 8p21 and rs1713985 at REST-C4orf14-POLR2B-IGFBP7 on 4q12) by using logistic regression analysis adjusted for age and sex. Subgroup analysis by CNV and PCV subtypes were performed to evaluate the significance of these two variants.

RESULTS

The reported association between rs13278062 at 8p21 and neovascular AMD was replicated in this population (P = 1.12 × 10(-4), odds ratio [OR] = 0.79, 95% confidence interval [CI] = 0.70-0.89). However, there was no association of rs1713985 at 4q12 with neovascular AMD, or its two subtypes, typical AMD-CNV and PCV (all P > 0.05). The study sample size had a statistical power of greater than 99% to detect an association of a risk allele with AMD with an OR of 1.30, as reported in the original study of rs1713985 and AMD.

CONCLUSIONS

The present results did not replicate the reported association between rs1713985 at 4q12 and neovascular AMD. However, we confirmed the association between rs13278062 at 8p21 and neovascular AMD in Asian populations.

Sex-Limited Genome-Wide Linkage Scan for Body Mass Index in an Unselected Sample of 933 Australian Twin Families

Genes involved in pathways regulating body weight may operate differently in men and women. To determine whether sex-limited genes influence the obesity-related phenotype body mass index (BMI), we have conducted a general non- scalar sex-limited genome-wide linkage scan using variance components analysis in Mx (Neale, 2002). BMI measurements and genotypic data were available for 2053 Australian female and male adult twins and their siblings from 933 families. Clinical measures of BMI were available for 64.4% of these individuals, while only self-reported measures were available for the remaining participants. The mean age of participants was 39.0 years of age (SD 12.1 years). The use of a sex-limited linkage model identified areas on the genome where quantitative trait loci (QTL) effects differ between the sexes, particularly on chromosome 8 and 20, providing us with evidence that some of the genes responsible for BMI may have different effects in men and women. Our highest linkage peak was observed at 12q24 (–log10p = 3.02), which was near the recommended threshold for suggestive linkage (–log10p = 3.13). Previous studies have found evidence for a quantitative trait locus on 12q24 affecting BMI in a wide range of populations, and candidate genes for non- insulin-dependent diabetes mellitus, a consequence of obesity, have also been mapped to this region. We also identified many peaks near a –log10p of 2 (threshold for replicating an existing finding) in many areas across the genome that are within regions previously identified by other studies, as well as in locations that harbor genes known to influence weight regulation.

Variance Components Models for Physical Activity With Age as Modifier: A Comparative Twin Study in Seven Countries

Physical activity is influenced by genetic factors whose expression may change with age. We employed an extension to the classical twin model that allows a modifier variable, age, to interact with the effects of the latent genetic and environmental factors. The model was applied to self-reported data from twins aged 19 to 50 from seven countries that collaborated in the GenomEUtwin project: Australia, Denmark, Finland, Norway, Netherlands, Sweden and United Kingdom. Results confirmed the importance of genetic influences on physical activity in all countries and showed an age-related decrease in heritability for 4 countries. In the other three countries age did not interact with heritability but those samples were smaller or had a more restricted age range. Effects of shared environment were absent, except in older Swedish participants. The study confirms the importance of taking age effects into account when exploring the genetic and environmental contribution to physical activity. It also suggests that the power of genome-wide association studies to identify the genetic variants contributing to physical activity may be larger in young adult cohorts.

RET mutational spectrum in Hirschsprung disease: evaluation of 601 Chinese patients

Rare (RVs) and common variants of the RET gene contribute to Hirschsprung disease (HSCR; congenital aganglionosis). While RET common variants are strongly associated with the commonest manifestation of the disease (males; short-segment aganglionosis; sporadic), rare coding sequence (CDS) variants are more frequently found in the lesser common and more severe forms of the disease (females; long/total colonic aganglionosis; familial).

Here we present the screening for RVs in the RET CDS and intron/exon boundaries of 601 Chinese HSCR patients, the largest number of patients ever reported. We identified 61 different heterozygous RVs (50 novel) distributed among 100 patients (16.64%). Those include 14 silent, 29 missense, 5 nonsense, 4 frame-shifts, and one in-frame amino-acid deletion in the CDS, two splice-site deletions, 4 nucleotide substitutions and a 22-bp deletion in the intron/exon boundaries and 1 single-nucleotide substitution in the 5′ untranslated region. Exonic variants were mainly clustered in RET the extracellular domain. RET RVs were more frequent among patients with the most severe phenotype (24% vs. 15% in short-HSCR). Phasing RVs with the RET HSCR-associated haplotype suggests that RVs do not underlie the undisputable association of RET common variants with HSCR. None of the variants were found in 250 Chinese controls.

Identification of four novel variants that influence central corneal thickness in multi-ethnic Asian populations

Central corneal thickness (CCT) is a highly heritable trait. Genes that significantly influence CCT can be candidate genes for common disorders in which CCT has been implicated, such as primary open-angle glaucoma (POAG) and keratoconus. Because the genetic factors controlling CCT in different Asian populations are unclear, we have built on previous work conducted on Singaporean Indians and Malays and extended our hypothesis to individuals of Chinese descent. We have followed up on all suggestive signals of association with CCT (P < 10(-4)) from the previously reported meta-analysis comprising Indians and Malays in a sample of Chinese individuals (n= 2681). In the combined sample (n= 7711), strong evidence of association was observed at four novel loci: IBTK on chromosome 6q14.1; CHSY1 on chromosome 15q26.3; and intergenic regions on chromosomes 7q11.2 and 9p23 (8.01 × 10(-11) < λ(GC) corrected P(meta) < 8.72 × 10(-8)). These four new loci explain an additional 4.3% of the total CCT variance across the sample cohorts over and above that of previously identified loci. We also extend on a previous finding at a fifth locus (AKAP13) where a new single-nucleotide polymorphism (rs1821481, P(meta) = 9.99 × 10(-9)) was found to be significantly more informative compared with the previously reported rs6496932 (P(meta) = 3.64 × 10(-5)). Performing association analysis in Asians may lead to the discovery of ethnic-specific genes that control CCT, offering further mechanistic insights into the regulation of CCT. In addition, it may also provide several candidate genes for interrogation for POAG, keratoconus and possible racial/ethnic variations.

Genome-wide association studies in Asians confirm the involvement of ATOH7 and TGFBR3, and further identify CARD10 as a novel locus influencing optic disc area

Damage to the optic nerve (e.g. from glaucoma) has an adverse and often irreversible impact on vision. Earlier studies have suggested that the size of the optic nerve head could be governed by hereditary factors. We conducted a genome-wide association study (GWAS) on 4445 Singaporean individuals (n = 2132 of Indian and n = 2313 of Malay ancestry, respectively), with replication in Rotterdam, the Netherlands (n = 9326 individuals of Caucasian ancestry) using the most widely reported parameter for optic disc traits, the optic disc area. We identified a novel locus on chromosome 22q13.1, CARD10, which strongly associates with optic disc area in both Singaporean cohorts as well as in the Rotterdam Study (RS; rs9607469, per-allele change in optic disc area = 0.051 mm(2); P(meta) = 2.73×10(-12)) and confirmed the association between CDC7/TGFBR3 (lead single nucleotide polymorphism (SNP) rs1192415, P(meta) = 7.57×10(-17)) and ATOH7 (lead SNP rs7916697, P(meta) = 2.00 × 10(-15)) and optic disc area in Asians. This is the first Asian-based GWAS on optic disc area, identifying a novel locus for the optic disc area, but also confirming the results found in Caucasian persons suggesting that there are general genetic determinants applicable to the size of the optic disc across different ethnicities.

Collagen-related genes influence the glaucoma risk factor, central corneal thickness

Central corneal thickness (CCT) is a risk factor of glaucoma, the most common cause of irreversible blindness worldwide. The identification of genetic determinants affecting CCT in the normal population will provide insights into the mechanisms underlying the association between CCT and glaucoma, as well as the pathogenesis of glaucoma itself. We conducted two genome-wide association studies for CCT in 5080 individuals drawn from two ethnic populations in Singapore (2538 Indian and 2542 Malays) and identified novel genetic loci significantly associated with CCT (COL8A2 rs96067, p(meta) = 5.40 × 10⁻¹³, interval of RXRA-COL5A1 rs1536478, p(meta) = 3.05 × 10⁻⁹). We confirmed the involvement of a previously reported gene for CCT and brittle cornea syndrome (ZNF469) [rs9938149 (p(meta) = 1.63 × 10⁻¹⁶) and rs12447690 (p(meta) = 1.92 × 10⁻¹⁴)]. Evidence of association exceeding the formal threshold for genome-wide significance was observed at rs7044529, an SNP located within COL5A1 when data from this study (n = 5080, P = 0.0012) were considered together with all published data (reflecting an additional 7349 individuals, p(Fisher) = 1.5 × 10⁻⁹). These findings implicate the involvement of collagen genes influencing CCT and thus, possibly the pathogenesis of glaucoma.

Haplotype analysis reveals a possible founder effect of RET mutation R114H for Hirschsprung's disease in the Chinese population

Background

Hirschsprung's disease (HSCR) is a congenital disorder associated with the lack of intramural ganglion cells in the myenteric and sub-mucosal plexuses along varying segments of the gastrointestinal tract. The RET gene is the major gene implicated in this gastrointestinal disease. A highly recurrent mutation in RET (RET^R114H) has recently been identified in ∼6–7% of the Chinese HSCR patients which, to date, has not been found in Caucasian patients or controls nor in Chinese controls. Due to the high frequency of RET^R114H in this population, we sought to investigate whether this mutation may be a founder HSCR mutation in the Chinese population.

Methodology and Principal Findings

To test whether all RET^R114 were originated from a single mutational event, we predicted the approximate age of RET^R114H by applying a Bayesian method to RET SNPs genotyped in 430 Chinese HSCR patients (of whom 25 individuals had the mutation) to be between 4–23 generations old depending on growth rate. We reasoned that if RET^R114H was a founder mutation then those with the mutation would share a haplotype on which the mutation resides. Including SNPs spanning 509.31 kb across RET from a recently obtained 500 K genome-wide dataset for a subset of 181 patients (14 RET^R114H patients), we applied haplotype estimation methods to determine whether there were any segments shared between patients with RET^R114H that are not present in those without the mutation or controls. Analysis yielded a 250.2 kb (51 SNP) shared segment over the RET gene (and downstream) in only those patients with the mutation with no similar segments found among other patients.

Conclusions

This suggests that RET^R114H is a founder mutation for HSCR in the Chinese population.

Multicenter dizygotic twin cohort study confirms two linkage susceptibility loci for body mass index at 3q29 and 7q36 and identifies three further potential novel loci

OBJECTIVE

To identify common loci and potential genetic variants affecting body mass index (BMI, kg m(-2)) in study populations originating from Europe.

DESIGN

We combined genome-wide linkage scans of six cohorts from Australia, Denmark, Finland, the Netherlands, Sweden and the United Kingdom with an approximately 10-cM microsatellite marker map. Variance components linkage analysis was carried out with age, sex and country of origin as covariates.

SUBJECTS

The GenomEUtwin consortium consists of twin cohorts from eight countries (Australia, Denmark, the Netherlands, Finland, Italy, Norway, Sweden and the United Kingdom) with a total data collection of more than 500,000 monozygotic and dizygotic (DZ) twin pairs. Variance due to early-life events and the environment is reduced within twin pairs, which makes DZ pairs highly valuable for linkage studies of complex traits. This study totaled 4401 European-originated twin families (10,535 individuals) from six countries (Australia, Denmark, the Netherlands, Finland, Sweden and the United Kingdom).

RESULTS

We found suggestive evidence for a quantitative trait locus on 3q29 and 7q36 in the combined sample of DZ twins (multipoint logarithm of odds score (MLOD) 2.6 and 2.4, respectively). Two individual cohorts showed strong evidence independently for three additional loci: 16q23 (MLOD=3.7) and 2p24 (MLOD=3.4) in the Dutch cohort and 20q13 (MLOD=3.2) in the Finnish cohort.

CONCLUSION

Linkage analysis of the combined data in this large twin cohort study provided evidence for suggestive linkage to BMI. In addition, two cohorts independently provided significant evidence of linkage to three new loci. The results of our study suggest a smaller environmental variance between DZ twins than full siblings, with a corresponding increase in heritability for BMI as well as an increase in linkage signal in well-replicated regions. The results are consistent with the possibility of locus heterogeneity for some genomic regions, and indicate a lack of major common quantitative trait locus variants affecting BMI in European populations.

Genetic variation in female BMI increases with number of children born but failure to replicate association between GNbeta3 variants and increased BMI in parous females

Objective:As post-pregnancy weight retention in women is a risk factor for obesity later in life, we assessed the changes in magnitude of genetic and environmental variation in BMI due to parity in an Australian female sample, comprising 11,915 female twins and their sisters.Results:Total variance of BMI increased with parity, primarily driven by an increase in nonadditive genetic variance. This finding was of particular interest given Gutersohn et al's (2000) report of recessive association between post-partum weight retention and the 825T allele of theGNβ3gene (rs5443) at 12p13.31. Hence, we attempted to replicate this association and test an additional three SNPs also located near or onGNβ3(rs10744716 (upstream), rs5442 (exon 10), rs5446 (exon 11)) in a sample of 3131 females and 1693 males from 2585 twin families. No association was found between these SNPs among females, even when allowing for a genotype by parity effect. However, a significant association was observed among males for rs10744716 (χ22= 10.22,p= 0.006; effect size = 0.47kg/m2), representing a novel association between a region upstream ofGNβ3and male BMI.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE

To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity.

DESIGN

We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins.

RESULTS

The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either.

CONCLUSIONS

In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Genetic influences on the difference in variability of height, weight and body mass index between Caucasian and East Asian adolescent twins

OBJECTIVE

Twin studies are useful for investigating the causes of trait variation between as well as within a population. The goals of the present study were two-fold: First, we aimed to compare the total phenotypic, genetic and environmental variances of height, weight and BMI between Caucasians and East Asians using twins. Secondly, we intended to estimate the extent to which genetic and environmental factors contribute to differences in variability of height, weight and BMI between Caucasians and East Asians.

DESIGN

Height and weight data from 3735 Caucasian and 1584 East Asian twin pairs (age: 13-15 years) from Australia, China, Finland, Japan, the Netherlands, South Korea, Taiwan and the United States were used for analyses. Maximum likelihood twin correlations and variance components model-fitting analyses were conducted to fulfill the goals of the present study.

RESULTS

The absolute genetic variances for height, weight and BMI were consistently greater in Caucasians than in East Asians with corresponding differences in total variances for all three body measures. In all 80 to 100% of the differences in total variances of height, weight and BMI between the two population groups were associated with genetic differences.

CONCLUSION

Height, weight and BMI were more variable in Caucasian than in East Asian adolescents. Genetic variances for these three body measures were also larger in Caucasians than in East Asians. Variance components model-fitting analyses indicated that genetic factors contributed to the difference in variability of height, weight and BMI between the two population groups. Association studies for these body measures should take account of our findings of differences in genetic variances between the two population groups.

Within-family outliers: segregating alleles or environmental effects? A linkage analysis of height from 5815 sibling pairs

Most information in linkage analysis for quantitative traits comes from pairs of relatives that are phenotypically most discordant or concordant. Confounding this, within-family outliers from non-genetic causes may create false positives and negatives. We investigated the influence of within-family outliers empirically, using one of the largest genome-wide linkage scans for height. The subjects were drawn from Australian twin cohorts consisting of 8447 individuals in 2861 families, providing a total of 5815 possible pairs of siblings in sibships. A variance component linkage analysis was performed, either including or excluding the within-family outliers. Using the entire dataset, the largest LOD scores were on chromosome 15q (LOD 2.3) and 11q (1.5). Excluding within-family outliers increased the LOD score for most regions, but the LOD score on chromosome 15 decreased from 2.3 to 1.2, suggesting that the outliers may create false negatives and false positives, although rare alleles of large effect may also be an explanation. Several regions suggestive of linkage to height were found after removing the outliers, including 1q23.1 (2.0), 3q22.1 (1.9) and 5q32 (2.3). We conclude that the investigation of the effect of within-family outliers, which is usually neglected, should be a standard quality control measure in linkage analysis for complex traits and may reduce the noise for the search of common variants of modest effect size as well as help identify rare variants of large effect and clinical significance. We suggest that the effect of within-family outliers deserves further investigation via theoretical and simulation studies.

Sex differences in genetic variation in weight: a longitudinal study of body mass index in adolescent twins

Genes that influence a phenotype earlier in life may differ from those influencing the same phenotype later, particularly during significant development periods such as puberty, when it is known that new genetic and environmental influences may become important. In the present study, body mass index (BMI) data were collected from 470 monozygotic twin pairs and 673 dizygotic twin pairs longitudinally at ages 12, 14 and 16, roughly straddling puberty. In order to examine whether there are qualitative and quantitative differences in genetic and environmental influences affecting BMI in males and females, during development, a general sex-limitation simplex model (which represents the longitudinal time series of the data) was fitted to the repeated measurements of BMI. The ADE simplex model provided the best fit to the adolescent data, with disparity in the magnitude of additive genetic influences between sexes, but no differences in the non-additive genetic (epistasis or dominance) or environmental influences. Results found may reflect many genetic and environmental influences during puberty, including the possible complex interaction between genes involved in the biological mechanism of weight regulation and the development of likely peer pressured activities such as severe exercise and diet regimes. Although, over 1,000 pairs of twins were used, this study still lacked the power to properly discriminate between additive and non-additive genetic variance.

Genetic influences on exercise participation in 37,051 twin pairs from seven countries

BACKGROUND

A sedentary lifestyle remains a major threat to health in contemporary societies. To get more insight in the relative contribution of genetic and environmental influences on individual differences in exercise participation, twin samples from seven countries participating in the GenomEUtwin project were used.

METHODOLOGY

Self-reported data on leisure time exercise behavior from Australia, Denmark, Finland, Norway, The Netherlands, Sweden and United Kingdom were used to create a comparable index of exercise participation in each country (60 minutes weekly at a minimum intensity of four metabolic equivalents).

PRINCIPAL FINDINGS

Modest geographical variation in exercise participation was revealed in 85,198 subjects, aged 19-40 years. Modeling of monozygotic and dizygotic twin resemblance showed that genetic effects play an important role in explaining individual differences in exercise participation in each country. Shared environmental effects played no role except for Norwegian males. Heritability of exercise participation in males and females was similar and ranged from 48% to 71% (excluding Norwegian males).

CONCLUSIONS

Genetic variation is important in individual exercise behavior and may involve genes influencing the acute mood effects of exercise, high exercise ability, high weight loss ability, and personality. This collaborative study suggests that attempts to find genes influencing exercise participation can pool exercise data across multiple countries and different instruments.

Bias, precision and heritability of self-reported and clinically measured height in Australian twins

Many studies of quantitative and disease traits in human genetics rely upon self-reported measures. Such measures are based on questionnaires or interviews and are often cheaper and more readily available than alternatives. However, the precision and potential bias cannot usually be assessed. Here we report a detailed quantitative genetic analysis of stature. We characterise the degree of measurement error by utilising a large sample of Australian twin pairs (857 MZ, 815 DZ) with both clinical and self-reported measures of height. Self-report height measurements are shown to be more variable than clinical measures. This has led to lowered estimates of heritability in many previous studies of stature. In our twin sample the heritability estimate for clinical height exceeded 90%. Repeated measures analysis shows that 2-3 times as many self-report measures are required to recover heritability estimates similar to those obtained from clinical measures. Bivariate genetic repeated measures analysis of self-report and clinical height measures showed an additive genetic correlation >0.98. We show that the accuracy of self-report height is upwardly biased in older individuals and in individuals of short stature. By comparing clinical and self-report measures we also showed that there was a genetic component to females systematically reporting their height incorrectly; this phenomenon appeared to not be present in males. The results from the measurement error analysis were subsequently used to assess the effects of error on the power to detect linkage in a genome scan. Moderate reduction in error (through the use of accurate clinical or multiple self-report measures) increased the effective sample size by 22%; elimination of measurement error led to increases in effective sample size of 41%.

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within-family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

Quantitative geneticists attempt to understand variation between individuals within a population for traits such as height in humans and the number of bristles in fruit flies. This has been traditionally done by partitioning the variation in underlying sources due to genetic and environmental factors, using the observed amount of variation between and within families. A problem with this approach is that one can never be sure that the estimates are correct, because nature and nurture can be confounded without one knowing it. The authors got around this problem by comparing the similarity between relatives as a function of the exact proportion of genes that they have in common, looking only within families. Using this approach, the authors estimated the amount of total variation for height in humans that is due to genetic factors from 3,375 sibling pairs. For each pair, the authors estimated the proportion of genes that they share from DNA markers. It was found that about 80% of the total variation can be explained by genetic factors, close to results that are obtained from classical studies. This study provides the first validation of an estimate of genetic variation by using a source of information that is free from nature–nurture assumptions.

Sex-limited genome-wide linkage scan for body mass index in an unselected sample of 933 Australian twin families

Genes involved in pathways regulating body weight may operate differently in men and women. To determine whether sex-limited genes influence the obesity-related phenotype body mass index (BMI), we have conducted a general nonscalar sex-limited genome-wide linkage scan using variance components analysis in Mx (Neale, 2002). BMI measurements and genotypic data were available for 2053 Australian female and male adult twins and their siblings from 933 families. Clinical measures of BMI were available for 64.4% of these individuals, while only self-reported measures were available for the remaining participants. The mean age of participants was 39.0 years of age (SD 12.1 years). The use of a sex-limited linkage model identified areas on the genome where quantitative trait loci (QTL) effects differ between the sexes, particularly on chromosome 8 and 20, providing us with evidence that some of the genes responsible for BMI may have different effects in men and women. Our highest linkage peak was observed at 12q24 (-log10p = 3.02), which was near the recommended threshold for suggestive linkage (-log10p = 3.13). Previous studies have found evidence for a quantitative trait locus on 12q24 affecting BMI in a wide range of populations, and candidate genes for noninsulin-dependent diabetes mellitus, a consequence of obesity, have also been mapped to this region. We also identified many peaks near a -log10p of 2 (threshold for replicating an existing finding) in many areas across the genome that are within regions previously identified by other studies, as well as in locations that harbor genes known to influence weight regulation.

Sex differences in heritability of BMI: a comparative study of results from twin studies in eight countries

Body mass index (BMI), a simple anthropometric measure, is the most frequently used measure of adiposity and has been instrumental in documenting the worldwide increase in the prevalence of obesity witnessed during the last decades. Although this increase in overweight and obesity is thought to be mainly due to environmental changes, i.e., sedentary lifestyles and high caloric diets, consistent evidence from twin studies demonstrates high heritability and the importance of genetic differences for normal variation in BMI. We analysed self-reported data on BMI from approximately 37,000 complete twin pairs (including opposite sex pairs) aged 20-29 and 30-39 from eight different twin registries participating in the GenomEUtwin project. Quantitative genetic analyses were conducted and sex differences were explored. Variation in BMI was greater for women than for men, and in both sexes was primarily explained by additive genetic variance in all countries. Sex differences in the variance components were consistently significant. Results from analyses of opposite sex pairs also showed evidence of sex-specific genetic effects suggesting there may be some differences between men and women in the genetic factors that influence variation in BMI. These results encourage the continued search for genes of importance to the body composition and the development of obesity. Furthermore, they suggest that strategies to identify predisposing genes may benefit from taking into account potential sex specific effects.

Heritability of adult body height: a comparative study of twin cohorts in eight countries

A major component of variation in body height is due to genetic differences, but environmental factors have a substantial contributory effect. In this study we aimed to analyse whether the genetic architecture of body height varies between affluent western societies. We analysed twin data from eight countries comprising 30,111 complete twin pairs by using the univariate genetic model of the Mx statistical package. Body height and zygosity were self-reported in seven populations and measured directly in one population. We found that there was substantial variation in mean body height between countries; body height was least in Italy (177 cm in men and 163 cm in women) and greatest in the Netherlands (184 cm and 171 cm, respectively). In men there was no corresponding variation in heritability of body height, heritability estimates ranging from 0.87 to 0.93 in populations under an additive genes/unique environment (AE) model. Among women the heritability estimates were generally lower than among men with greater variation between countries, ranging from 0.68 to 0.84 when an additive genes/shared environment/unique environment (ACE) model was used. In four populations where an AE model fit equally well or better, heritability ranged from 0.89 to 0.93. This difference between the sexes was mainly due to the effect of the shared environmental component of variance, which appears to be more important among women than among men in our study populations. Our results indicate that, in general, there are only minor differences in the genetic architecture of height between affluent Caucasian populations, especially among men.