Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10−9 to P = 1.8 × 10−40) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10−3 to P = 1.2 × 10−13). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

Genetic and Malt Quality Diversity of East Asian Two-Rowed Barley Accessions in Relationship to North American Malting Barley.

Because of epidemics of Fusarium head blight (FHB; caused by Fusarium graminearum Schwabe [teleomorph Gibberella zeae (Schwein.) Petch]) in the northern Great Plains of the United States and Canada in the past two decades, malting barley breeders have been forced to use nonadapted barley (Hordeum vulgare L.) accessions as sources of FHB resistance. Many of the resistant accessions are from East Asia, and limited information is available on their genetic diversity and malt quality. The objectives of this study were to determine the genetic diversity among 30 East Asian accessions and two North American cultivars. Genetic diversity was based on 49 simple-sequence repeat markers. All accessions were tested for barley grain brightness; protein content; 1,000-kernel weight; malting loss; fine-grind malt extract; content of plump kernels, free amino nitrogen, soluble protein, and wort beta-glucan; the Kolbach index (i.e., the ratio of malt soluble protein to malt total protein); a-amylase activity; diastatic power; won color; and wort viscosity. A few accessions had equal quality compared with Harrington and Conlon barley for individual traits but not for all. Qing 2, Mokkei 93-78, and Nitakia 48 could be excellent sources for increased malt extract; Nitakia 48 is a possible source for low wort viscosity; and Mokkei 93-78 and Nitakia 48 are putative sources of low beta-glucan content. The cluster analyses also implied that the malt quality of an accession cannot be predicted based on the country where it was developed.

Genetic effects and genetic relationships among shrunken (sh2) sweet corn lines and F1 hybrids

The objectives of this study were to quantify the components of genetic variance and the genetic effects, and to examine the genetic relationship of inbred lines extracted from various shrunken2(sh2) breeding populations. Ten diverse inbred lines developed from genetic background, were crossed in half diallel. Parents and their F1 hybrids were evaluated at three environments. The parents were genotyped using 20 polymorphic simple sequence repeats (SSR). Agronomic and quality traits were analysed by a mixed linear model according to additive-dominance genetic model. Genetic effects were estimated using an adjusted unbiased prediction method. Additive variance was more important than dominance variance in the expression of traits related to ear aspects (husk ratio and percentage of ear filled) and eating quality (flavour and total soluble solids). For agronomic traits, however, dominance variance was more important than additive variance. The additive genetic correlation between flavour and tenderness was strong (r = 0.84, P <0.01). Flavour, tenderness and kernel colour additive genetic effects were not correlated with yield related traits. Genetic distance (GD), estimated from SSR profiles on the basis of Jaccard's similarity coefficient varied from 0.10 to 0.77 with an average of 0.56. Cluster analysis classified parents according to their pedigree relationships. In most studied traits, F1 performance was not associated with GD.

A genome-wide linkage scan provides evidence for both new and previously reported loci influencing common migraine

Latent class analysis was performed on migraine symptom data collected in a Dutch population sample (N = 12,210, 59% female) in order to obtain empirical groupings of individuals suffering from symptoms of migraine headache. Based on these heritable groupings (h(2) = 0.49, 95% CI: 0.41-0.57) individuals were classified as affected (migrainous headache) or unaffected. Genome-wide linkage analysis was performed using genotype data from 105 families with at least 2 affected siblings. In addition to this primary phenotype, linkage analyses were performed for the individual migraine symptoms. Significance levels, corrected for the analysis of multiple traits, were determined empirically via a novel simulation approach. Suggestive linkage for migrainous headache was found on chromosomes 1 (LOD = 1.63; pointwise P = 0.0031), 13 (LOD = 1.63; P = 0.0031), and 20 (LOD = 1.85; P = 0.0018). Interestingly, the chromosome 1 peak was located close to the ATP1A2 gene, associated with familial hemiplegic migraine type 2 (FHM2). Individual symptom analysis produced a LOD score of 1.97 (P = 0.0013) on chromosome 5 (photo/phonophobia), a LOD score of 2.13 (P = 0.0009) on chromosome 10 (moderate/severe pain intensity) and a near significant LOD score of 3.31 (P = 0.00005) on chromosome 13 (pulsating headache). These peaks were all located near regions previously reported in migraine linkage studies. Our results provide important replication and support for the presence of migraine susceptibility genes within these regions, and further support the utility of an LCA-based phenotyping approach and analysis of individual symptoms in migraine genetic research. Additionally, our novel "2-step" analysis and simulation approach provides a powerful means to investigate linkage to individual trait components.

Genetic variation in tumour necrosis factor and lymphotoxin is not associated with endometriosis in an Australian sample

BACKGROUND Tumour necrosis factor (TNF) is a pleiotropic cytokine with a wide range of immunoregulatory effects. Variation in the promoter region of TNF and the neighbouring lymphotoxin alpha (LTA) gene might be associated with endometriosis. METHODS We examined the association between endometriosis and common single-nucleotide polymorphisms (SNPs) or haplotypes in the TNF/LTA region in an Australian sample by analysing 26 SNPs in 958 endometriosis cases and 959 unrelated controls. We selected functional SNPs in the coding and the promoter region of the TNF gene and HapMap tagging SNPs and typed them on a Sequenom MassARRAY platform. A key SNP (rs1800630) in the promoter region typed in previous studies did not give reliable results. Therefore, we also examined a statistically identical (r(2) = 1) SNP (siSNP) (rs2844482), identified using the web based program ssSNPer. RESULTS Genotype completion rate was 99.5% for SNPs spanning a region of 15.5 kb across the TNF/LTA locus. There was no evidence for association between endometriosis and TNF/LTA SNPs or SNP haplotypes in our case-control study. CONCLUSIONS Our data suggest both TNF and LTA genes are not major susceptibility genes for endometriosis.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.

KRAS variation and risk of endometriosis

Endometriosis is a common gynaecological disease with symptoms of pelvic pain and infertility which affects 7-10% of women in their reproductive years. Activation of an oncogenic allele of Kirsten rat sarcoma viral oncogene homologue (KRAS) in the reproductive tract of mice resulted in the development of endometriosis. We hypothesized that variation in KRAS may influence risk of endometriosis in humans. Thirty tagSNPs spanning a region of 60.7 kb across the KRAS locus were genotyped using iPLEX chemistry on a MALDI-TOF MassARRAY platform in 959 endometriosis cases and 959 unrelated controls, and data were analysed for association with endometriosis. Genotypes were obtained for most individuals with a mean completion rate of 99.1%. We identified six haplotype blocks across the KRAS locus in our sample. There were no significant differences between cases and controls in the frequencies of individual single-nucleotide polymorphisms (SNPs) or haplotypes. We also developed a rapid method to screen for 11 common KRAS and BRAF mutations on the Sequenom MassARRAY system. The assay detected all mutations previously identified by direct sequencing in a panel of positive controls. No germline variants for KRAS or BRAF were detected. Our results demonstrate that any risk of endometriosis in women because of common variation in KRAS must be very small.

ssSNPer: identifying statistically similar SNPs to aid interpretation of genetic association studies

ssSNPer is a novel user-friendly web interface that provides easy determination of the number and location of untested HapMap SNPs, in the region surrounding a tested HapMap SNP, which are statistically similar and would thus produce comparable and perhaps more significant association results. Identification of ssSNPs can have crucial implications for the interpretation of the initial association results and the design of follow-up studies. AVAILABILITY: http://fraser.qimr.edu.au/general/daleN/ssSNPer/

Genetic diversity of Atlantic salmon (Salmo salar) and European whitefish (Coregonus lavaretus) in the Baltic Sea

We described the patterns and extent of microsatellite DNA variation in historical and present-day Atlantic salmon (Salmo salar L.) stocks in the Baltic Sea and neighbouring areas, and in European whitefish (Coregonus lavaretus) ecotypes, populations and run-timing types in Finland. Moreover, the amount and pattern of genetic diversity in historical salmon populations before human impact were described, and the proportion of diversity maintained in the present hatchery stocks evaluated. Salmon populations in the Baltic Sea were, on average, significantly less variable than eastern Atlantic populations, and the diversity of landlocked populations (Lakes Vänern, Saimaa, Onega and Ladoga) was in turn significantly lower than that of anadromous salmon populations in the Baltic Sea populations. Within the Baltic Sea, the anadromous populations of Atlantic salmon formed three clear groups, corresponding to the northern (Gulf of Bothnia), eastern (Gulf of Finland and eastern Baltic Main Basin) and southern (western Baltic Main Basin) regions. Based on microsatellite data, three salmon population groups in the Baltic Sea were considered potentially different colonization lineages. In short- and long-term breeding programmes of Atlantic salmon, the average observed rate of loss of alleles was 4.9% and 2.0% per generation and the average rate of loss of heterozygosity was 1.4% and 1% per generation, respectively. When comparing the genetic parameters of stocks before and after hatchery breeding of several successive generations (Rivers Iijoki and Oulujoki), statistically significant changes in allele frequencies were common, while large wild stock in the Teno River has remained temporally very stable over 56 years. Despite the observed losses of genetic diversity in broodstock breeding, a large proportion of the genetic resources of the extirpated stocks are still conserved in the broodstocks. Genetic differentiation among European whitefish ecotypes was generally low, thus giving support to the hypothesis of one native European whitefish species in Fennoscandia. Among the ecotypes, the northern, large sparsely rakered, bottom-dwelling whitefish was the most unique. The known genetic differences in quantitative traits have thus either developed independently of potential phylogenetic lineages, or the lineages have mixed and the quantitative traits of the ecotypes, like gill-raker number, have later changed according to environment and selection pressures. Overall, genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Wild whitefish populations studied had slightly higher allelic diversity than hatchery-reared populations in corresponding rivers.

Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits

A total of 4063 young bulls of two tropical genotypes (1639 Brahman and 2424 Tropical Composite) raised in northern Australia were evaluated for a comprehensive range of production and reproduction traits up to 24 months of age. Prior to weaning, peripheral blood concentrations of luteinising hormone (LH) and inhibin were measured at 4 months of age. At weaning (6 months) blood insulin-like growth factor-1 (IGF-I) and flight time were recorded. Body composition traits of fat depth and eye-muscle area were determined by ultrasonography at 15 months of age when additional measurements of liveweight, hip height and body condition score were recorded. Bull breeding soundness was evaluated at similar to 12, 18 and 24 months of age when measurements of scrotal circumference, sheath score, semen mass activity, progressive motility of individual sperm and percent morphologically normal sperm were recorded. Magnitude of heritability and genetic correlations changed across time for some traits. Heritability of LH, inhibin, IGF-I and of 18-month scrotal circumference, mass activity, progressive motility and percent normal sperm was 0.31, 0.74, 0.44, 0.75, 0.24, 0.15 and 0.25, respectively, for Brahmans and 0.48, 0.72, 0.36, 0.43, 0.13, 0.15 and 0.20, respectively, for Tropical Composites. Inhibin and IGF-I had moderate genetic association with percent normal sperm at 24 months in Brahmans but low to negligible associations in Tropical Composites. Body condition score in Brahmans and sperm motility (mass and individual) traits in both genotypes had moderate to strong genetic correlation with percent normal sperm and may prove useful candidates for indirect selection. There is scope to increase scrotal circumference by selection and this will be associated with favourable correlated responses of improved semen quality in both genotypes. The lack of genetic antagonism among bull traits indicates that selection for improved semen quality will not adversely affect other production traits.

Predicting Additive and Non-additive Genetic Effects from Trials Where Traits Are Affected by Interplot Competition

There are two key types of selection in a plant breeding program, namely selection of hybrids for potential commercial use and the selection of parents for use in future breeding. Oakey et al. (in Theoretical and Applied Genetics 113, 809-819, 2006) showed how both of these aims could be achieved using pedigree information in a mixed model analysis in order to partition genetic effects into additive and non-additive effects. Their approach was developed for field trial data subject to spatial variation. In this paper we extend the approach for data from trials subject to interplot competition. We show how the approach may be used to obtain predictions of pure stand additive and non-additive effects. We develop the methodology in the context of a single field trial using an example from an Australian sorghum breeding program.

Migraine with aura and migraine without aura are not distinct entities: further evidence from a large Dutch population study

It is often debated whether migraine with aura (MA) and migraine without aura (MO) are etiologically distinct disorders. A previous study using latent class analysis (LCA) in Australian twins showed no evidence for separate subtypes of MO and MA. The aim of the present study was to replicate these results in a population of Dutch twins and their parents, siblings and partners (N = 10,144). Latent class analysis of International Headache Society (IHS)-based migraine symptoms resulted in the identification of 4 classes: a class of unaffected subjects (class 0), a mild form of nonmigrainous headache (class 1), a moderately severe type of migraine (class 2), typically without neurological symptoms or aura (8% reporting aura symptoms), and a severe type of migraine (class 3), typically with neurological symptoms, and aura symptoms in approximately half of the cases. Given the overlap of neurological symptoms and nonmutual exclusivity of aura symptoms, these results do not support the MO and MA subtypes as being etiologically distinct. The heritability in female twins of migraine based on LCA classification was estimated at .50 (95% confidence intervals [CI] .27 - .59), similar to IHS-based migraine diagnosis (h2 = .49, 95% CI .19-.57). However, using a dichotomous classification (affected-unaffected) decreased heritability for the IHS-based classification (h2 = .33, 95% CI .00-.60), but not the LCA-based classification (h2 = .51, 95% CI .23-.61). Importantly, use of the LCA-based classification increased the number of subjects classified as affected. The heritability of the screening question was similar to more detailed LCA and IHS classifications, suggesting that the screening procedure is an important determining factor in genetic studies of migraine.

Trait components provide tools to dissect the genetic susceptibility of migraine

The commonly used "end diagnosis" phenotype that is adopted in linkage and association studies of complex traits is likely to represent an oversimplified model of the genetic background of a disease. This is also likely to be the case for common types of migraine, for which no convincingly associated genetic variants have been reported. In headache disorders, most genetic studies have used end diagnoses of the International Headache Society (IHS) classification as phenotypes. Here, we introduce an alternative strategy; we use trait components--individual clinical symptoms of migraine--to determine affection status in genomewide linkage analyses of migraine-affected families. We identified linkage between several traits and markers on chromosome 4q24 (highest LOD score under locus heterogeneity [HLOD] 4.52), a locus we previously reported to be linked to the end diagnosis migraine with aura. The pulsation trait identified a novel locus on 17p13 (HLOD 4.65). Additionally, a trait combination phenotype (IHS full criteria) revealed a locus on 18q12 (HLOD 3.29), and the age at onset trait revealed a locus on 4q28 (HLOD 2.99). Furthermore, suggestive or nearly suggestive evidence of linkage to four additional loci was observed with the traits phonophobia (10q22) and aggravation by physical exercise (12q21, 15q14, and Xp21), and, interestingly, these loci have been linked to migraine in previous studies. Our findings suggest that the use of symptom components of migraine instead of the end diagnosis provides a useful tool in stratifying the sample for genetic studies.