A genome-wide association study on obesity and obesity-related traits

Large-scale genome-wide association studies (GWAS) have identified many loci associated with body mass index (BMI), but few studies focused on obesity as a binary trait. Here we report the results of a GWAS and candidate SNP genotyping study of obesity, including extremely obese cases and never overweight controls as well as families segregating extreme obesity and thinness. We first performed a GWAS on 520 cases (BMI>35 kg/m(2)) and 540 control subjects (BMI<25 kg/m(2)), on measures of obesity and obesity-related traits. We subsequently followed up obesity-associated signals by genotyping the top ～500 SNPs from GWAS in the combined sample of cases, controls and family members totaling 2,256 individuals. For the binary trait of obesity, we found 16 genome-wide significant signals within the FTO gene (strongest signal at rs17817449, P = 2.5 × 10(-12)). We next examined obesity-related quantitative traits (such as total body weight, waist circumference and waist to hip ratio), and detected genome-wide significant signals between waist to hip ratio and NRXN3 (rs11624704, P = 2.67 × 10(-9)), previously associated with body weight and fat distribution. Our study demonstrated how a relatively small sample ascertained through extreme phenotypes can detect genuine associations in a GWAS.     

Analyzes of genome-wide association follow-up study for calving traits in dairy cattle

ABSTRACT: BACKGROUND: There is often a pronounced disagreement between results obtained from different genome-wide association studies in cattle. There are multiple reasons for this disagreement. Particularly the presence of false positives leads to a need to validate detected QTL before they are optimally incorporated or weighted in selection decisions or further studied for causal gene. In dairy cattle progeny testing scheme new data is routinely accumulated which can be used to validate previously discovered associations. However, the data is not an independent sample and the sample size may not be sufficient to have enough power to validate previous discoveries. Here we compared two strategies to validate previously detected QTL when new data is added from the same study population. We compare analyzing a combined dataset (COMB) including all data presently available to only analyzing a validation dataset (VAL) i.e. a new dataset not previously analyzed as an independent replication. Secondly, we confirm SNP detected in the Reference population (REF) (i.e. previously analyzed dataset consists of older bulls) in the VAL dataset. RESULTS: Clearly the results from the combined (COMB) dataset which had nearly twice the sample size of other two subsets allowed the detection of far more significant associations than the two smaller subsets. The number of significant SNPs in REF (older bulls) was about four times higher compare to VAL (younger bulls) though both had similar sample sizes, 2,219 and 2,039 respectively. A total of 424 SNP-trait combinations on 22 chromosomes showed genome-wide significant association involving 284 unique SNPs in the COMB dataset. In the REF data set 101 associations (73 unique SNPs) and in the VAL 24 associations (18 unique SNPs) were found genome-wide significant.Sixty-eight percent of the SNPs in the REF dataset could be confirmed in the VAL dataset. Out of 469 unique SNPs showing chromosome-wide significant association with calving traits in the REF dataset 321 could be confirmed in the VAL dataset at P < 0.05. CONCLUSIONS: The follow-up study for GWAS in cattle will depend on the aim of the study. If the aim is to discover novel QTL, analyses of the COMB dataset is recommended, while in case of identification of the causal mutation underlying a QTL, confirmation of the discovered SNPs are necessary to avoid following a false positive.     

Analyzes of genome-wide association follow-up study for calving traits in dairy cattle

Background

Genomic islands play an important role in medical, methylation and biological studies. To explore the region, we propose a CpG islands prediction analysis platform for genome sequence exploration (CpGPAP).

Results

CpGPAP is a web-based application that provides a user-friendly interface for predicting CpG islands in genome sequences or in user input sequences. The prediction algorithms supported in CpGPAP include complementary particle swarm optimization (CPSO), a complementary genetic algorithm (CGA) and other methods (CpGPlot, CpGProD and CpGIS) found in the literature. The CpGPAP platform is easy to use and has three main features (1) selection of the prediction algorithm; (2) graphic visualization of results; and (3) application of related tools and dataset downloads. These features allow the user to easily view CpG island results and download the relevant island data. CpGPAP is freely available at http://bio.kuas.edu.tw/CpGPAP/.

Conclusions

The platform's supported algorithms (CPSO and CGA) provide a higher sensitivity and a higher correlation coefficient when compared to CpGPlot, CpGProD, CpGIS, and CpGcluster over an entire chromosome.     

An imputation-based genome-wide association study for growth and fatness traits in Sujiang pigs

Sujiang pigs are a synthetic breed derived from Jiangquhai, Fengjing, and Duroc pigs. In this study, we sequenced the genome of 62 pigs with a coverage depth of 10× to 20×, including 27 Sujiang and 35 founder breed pigs, and we collected 360 global pigs’ genome sequence data from public databases including 39 Duroc pigs. We obtained a high-quality variant dataset of 365 Sujiang pigs by imputing the porcine 80 K single nucleotide polymorphism (SNP) Beadchip to the whole-genome scale with a total of 422 pigs as a reference panel. A dataset of 365 imputated Sujiang pigs was used to perform single-trait genome-wide association study (GWAS) and meta-analyses for growth and fatness traits. Single-trait GWAS identified 1 907, 18, and 14 SNPs surpassing the suggestively significant threshold for backfat thickness, chest circumference, and chest width, respectively. Meta-analyses identified 2 400 genome-wide significant SNPs and 520 suggestively significant SNPs for backfat thickness and chest circumference, and 719 genome-wide significant SNPs and 1 225 suggestively significant SNPs for all seven traits. According to the meta-analysis of backfat thickness and chest circumference, a remarkable region of 2.69 Mb on Sus scrofa chromosome 4 containing FAM110B, IMPAD1, LYN, MOS, PENK, PLAG1, SDR16C5 and XKR4 was identified as a candidate region. The haplotype heat map of the 2.69 Mb region verified that Sujiang pigs were derived from Duroc and Chinese indigenous pigs, especially Jiangquhai pigs. The Kruskal-Wallis test showed that haplotypes of the 2.69 Mb region significantly affected backfat thickness and chest circumference traits. We then focused on PLAG1, an important growth-related gene, and identified two synonymous SNPs with obvious differences among different breeds in the PLAG1 gene. We then performed genotyping of 365 Sujiang, 150 Duroc, 95 Jiangquhai, and 100 Fengjing pigs to confirm the above result and verified that the two variants significantly affected phenotypes of growth and fatness traits. Our findings not only provide insights into the genetic architecture of porcine growth and fatness traits but also provide potential markers for selective breeding of these traits in Sujiang pigs.     

Genetic dissection of wheat panicle traits using linkage analysis and a genome-wide association study

Key message

Coincident regions on chromosome 4B for GW, on 5A for SD and TSS, and on 3A for SL and GNS were detected through an integration of a linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A and 6A were identified with high PVE% on a composite map.

Abstract

The panicle traits of wheat, such as grain number per spike and 1000-grain weight, are closely correlated with grain yield. Superior and effective alleles at loci related to panicles developments play a crucial role in the progress of molecular improvement in wheat yield breeding. Here, we revealed several notable allelic variations of seven panicle-related traits through an integration of genome-wide association mapping and a linkage analysis. The linkage analysis was performed using a recombinant inbred line (RIL) population (173 lines of F_8:9) with a high-density genetic map constructed with 90K SNP arrays, Diversity Arrays Technology (DArT) and simple sequence repeat (SSR) markers in five environments. Thirty-five additive quantitative trait loci (QTL) were discovered, including eleven stable QTLs on chromosomes 1A, 2D, 4B, 5B, 6B, and 6D. The marker interval between EX_C101685 and RAC875_C27536 on chromosome 4B exhibited pleiotropic effects for GW, SL, GNS, FSN, SSN, and TSS, with the phenotypic variation explained (PVE) ranging from 5.40 to 37.70%. In addition, an association analysis was conducted using a diverse panel of 205 elite wheat lines with a composite map (24,355 SNPs) based on the Illumina Infinium assay in four environments. A total of 73 significant marker-trait associations (MTAs) were detected for panicle traits, which were distributed across all wheat chromosomes except for 4D, 5D, and 6D. Consensus regions between RAC875_C27536_611 and Tdurum_contig4974_355 on chromosome 4B for GW in multiple environments, between QTSS5A.7-43 and BS00021805_51 on 5A for SD and TSS, and between QSD3A.2-164 and RAC875_c17479_359 on 3A for SL and GNS in multiple environments were detected through linkage analysis and a genome-wide association study (GWAS). In addition, six stable QTL clusters on chromosomes 2D, 3A, 4B, 5A, and 6A were identified with high PVE% on a composite map. This study provides potentially valuable information on the dissection of yield-component traits and valuable genetic alleles for molecular-design breeding or functional gene exploration.

     

The mean and variability of a floral trait have opposing effects on fitness traits

Background and Aims Floral traits are essential for ensuring successful pollination and reproduction in flowering plants. In particular, style and anther positions are key for pollination accuracy and efficiency. Variation in these traits among individuals has been well studied, but less is known about variation within flowers and plants and its effect on pollination and reproductive success.Methods Style deflexion is responsible for herkogamy and important for pollen deposition in Passiflora incarnata. The degree of deflexion may vary among stigmas within flowers as well as among flowers. We measured the variability of style deflexion at both the flower and the plant level. The fitness consequences of the mean and variation of style deflexion were then evaluated under natural pollination by determining their relationship to pollen deposition, seed production and average seed weight using structural equation modelling. In addition, the relationship between style deflexion and self-pollen deposition was estimated in a greenhouse experiment.Key Results We found greater variation in style deflexion within flowers and plants than among plants. Variation of style deflexion at the flower and plant level was positively correlated, suggesting that variability in style deflexion may be a distinct trait in P. incarnata. Lower deflexion and reduced variation in that deflexion increased pollen deposition, which in turn increased seed number. However, lower styles also increased self-pollen deposition. In contrast, higher deflexion and greater variability of that deflexion increased variation in pollen deposition, which resulted in heavier seeds.Conclusions Variability of style deflexion and therefore stigma placement, independent from the mean, appears to be a property of individual P. incarnata plants. The mean and variability of style deflexion in P. incarnata affected seed number and seed weight in contrasting ways, through the quantity and potentially quality of pollen deposition. This antagonistic selection via different fitness components may maintain diverse style phenotypes.     

Properties of ethylmethane sulfonate-induced mutations affecting life-history traits in Caenorhabditis elegans and inferences about bivariate distributions of mutation effects

The homozygous effects of ethylmethane sulfonate (EMS)-induced mutations in Caenorhabditis elegans are compared across life-history traits. Mutagenesis has a greater effect on early than late reproductive output, since EMS-induced mutations tend to cause delayed reproduction. Mutagenesis changes the mean and variance of longevity much less than reproductive output traits. Mutations that increase total or early productivity are not detected, but the net effect of mutations is to increase and decrease late productivity to approximately equal extents. Although most mutations decrease longevity, a mutant line with increased longevity was found. A flattening of mortality curves with age is noted, particularly in EMS lines. We infer that less than one-tenth of mutations that have fitness effects in natural conditions are detected in the laboratory, and such mutations have moderately large effects ( approximately 20% of the mean). Mutational correlations for life-history traits are strong and positive. Correlations between early or late productivity and longevity are of similar magnitude. We develop a maximum-likelihood procedure to infer bivariate distributions of mutation effects. We show that strong mutation-induced genetic correlations do not necessarily imply strong directional correlations between mutational effects, since correlation is also generated by lines carrying different numbers of mutations.     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状,随着平衡育种理念的提出和发展,繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展,全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

Bayesian genome-wide association study of nut traits in Japanese chestnut

Japanese chestnut (Castanea crenata Sieb. et Zucc.) has a long juvenile phase, so breeders have to wait many years to evaluate nut traits. Molecular markers associated with genes of interest would accelerate selection in chestnut breeding programs. We evaluated five nut traits (nut harvest date, nut weight, pericarp splitting, insect infestation, and specific gravity) in 99 Japanese chestnut cultivars and selections. A wide range of phenotypic variation was observed for each of the traits, suggesting that the collection harbored sufficient genetic diversity for breeding. A Bayesian genome-wide association study was conducted using 162 simple sequence repeat markers and 741 single nucleotide polymorphism markers. To evaluate the potential of marker-assisted selection, we examined 12 molecular markers found to be associated with nut traits: 4 for nut harvest date, 4 for nut weight, 1 for pericarp splitting, and 3 for insect infestation. The percentages of phenotypic variance explained ranged from 4.8 to 37.1%. Although insect infestation showed only medium heritability (0.672), we obtained two quantitative trait loci (QTLs) with extremely high posterior probabilities (0.93 and 1.00). Out of the 12 molecular markers, 3 of the 4 markers for nut harvest time could be validated in a breeding population. Accuracies of genomic selection were extremely high for nut harvest date (0.841) and moderate for insect infestation (0.604), suggesting that genomic selection based on Bayesian regression would reduce the cost of phenotypic evaluation of these traits and possibly others.     

Testing the extreme value domain of attraction for distributions of beneficial fitness effects

In modeling evolutionary genetics, it is often assumed that mutational effects are assigned according to a continuous probability distribution, and multiple distributions have been used with varying degrees of justification. For mutations with beneficial effects, the distribution currently favored is the exponential distribution, in part because it can be justified in terms of extreme value theory, since beneficial mutations should have fitnesses in the extreme right tail of the fitness distribution. While the appeal to extreme value theory seems justified, the exponential distribution is but one of three possible limiting forms for tail distributions, with the other two loosely corresponding to distributions with right-truncated tails and those with heavy tails. We describe a likelihood-ratio framework for analyzing the fitness effects of beneficial mutations, focusing on testing the null hypothesis that the distribution is exponential. We also describe how to account for missing the smallest-effect mutations, which are often difficult to identify experimentally. This technique makes it possible to apply the test to gain-of-function mutations, where the ancestral genotype is unable to grow under the selective conditions. We also describe how to pool data across experiments, since we expect few possible beneficial mutations in any particular experiment.     

Natural selection fails to optimize mutation rates for long-term adaptation on rugged fitness landscapes

The rate of mutation is central to evolution. Mutations are required for adaptation, yet most mutations with phenotypic effects are deleterious. As a consequence, the mutation rate that maximizes adaptation will be some intermediate value. Here, we used digital organisms to investigate the ability of natural selection to adjust and optimize mutation rates. We assessed the optimal mutation rate by empirically determining what mutation rate produced the highest rate of adaptation. Then, we allowed mutation rates to evolve, and we evaluated the proximity to the optimum. Although we chose conditions favorable for mutation rate optimization, the evolved rates were invariably far below the optimum across a wide range of experimental parameter settings. We hypothesized that the reason that mutation rates evolved to be suboptimal was the ruggedness of fitness landscapes. To test this hypothesis, we created a simplified landscape without any fitness valleys and found that, in such conditions, populations evolved near-optimal mutation rates. In contrast, when fitness valleys were added to this simple landscape, the ability of evolving populations to find the optimal mutation rate was lost. We conclude that rugged fitness landscapes can prevent the evolution of mutation rates that are optimal for long-term adaptation. This finding has important implications for applied evolutionary research in both biological and computational realms.     

GWAPower: a statistical power calculation software for genome-wide association studies with quantitative traits

Background

In designing genome-wide association (GWA) studies it is important to calculate statistical power. General statistical power calculation procedures for quantitative measures often require information concerning summary statistics of distributions such as mean and variance. However, with genetic studies, the effect size of quantitative traits is traditionally expressed as heritability, a quantity defined as the amount of phenotypic variation in the population that can be ascribed to the genetic variants among individuals. Heritability is hard to transform into summary statistics. Therefore, general power calculation procedures cannot be used directly in GWA studies. The development of appropriate statistical methods and a user-friendly software package to address this problem would be welcomed.

Results

This paper presents GWAPower, a statistical software package of power calculation designed for GWA studies with quantitative traits, where genetic effect is defined as heritability. Based on several popular one-degree-of-freedom genetic models, this method avoids the need to specify the non-centrality parameter of the F-distribution under the alternative hypothesis. Therefore, it can use heritability information directly without approximation. In GWAPower, the power calculation can be easily adjusted for adding covariates and linkage disequilibrium information. An example is provided to illustrate GWAPower, followed by discussions.

Conclusions

GWAPower is a user-friendly free software package for calculating statistical power based on heritability in GWA studies with quantitative traits. The software is freely available at: http://dl.dropbox.com/u/10502931/GWAPower.zip     

Physiological effects on demography: a long-term experimental study of testosterone's effects on fitness

Understanding physiological and behavioral mechanisms underlying the diversity of observed life-history strategies is challenging because of difficulties in obtaining long-term measures of fitness and in relating fitness to these mechanisms. We evaluated effects of experimentally elevated testosterone on male fitness in a population of dark-eyed juncos studied over nine breeding seasons using a demographic modeling approach. Elevated levels of testosterone decreased survival rates but increased success of producing extra-pair offspring. Higher overall fitness for testosterone-treated males was unexpected and led us to consider indirect effects of testosterone on offspring and females. Nest success was similar for testosterone-treated and control males, but testosterone-treated males produced smaller offspring, and smaller offspring had lower postfledging survival. Older, more experienced females preferred to mate with older males and realized higher reproductive success when they did so. Treatment of young males increased their ability to attract older females yet resulted in poor reproductive performance. The higher fitness of testosterone-treated males in the absence of a comparable natural phenotype suggests that the natural phenotype may be constrained. If this phenotype were to arise, the negative social effects on offspring and mates suggest that these effects might prevent high-testosterone phenotypes from spreading in the population.     

Equilibrium between genetic drift and migration at various mutation rates: simulation analysis

Rates of approach to equilibrium values of F(ST)/R(ST) at various mutation rates and using different mutation models (K-allele model KAM and stepwise model SMM) were analyzed numerically for the finite island model and the one-dimensional stepping stone models of migration, using simulation. In the island model of migration and the KAM mutation model, the rate of approach to the equilibrium F(ST) value was appreciably higher and the equilibrium value was almost twofold lower at micro (mutation rate) = m (migration rate) than at micro < m. In the one-dimensional stepping stone model of migration and the KAM model of mutation, the mutation rate significantly affected both the rate of approaching F(ST) equilibrium and the equilibrium value. In both island and one-dimensional stepping stone models and SMM, R(ST) was not influenced by various mutation rates. The rate of approach to the equilibrium values of both F(ST) and R(ST) was lower for the stepping stone model than to the island model. RST was rather resistant to deviations from the SMM mutation model.     

A genome-wide association study of wheat yield and quality-related traits in southwest China

Wheat (Triticum aestivum L.) is one of the most productive and important crops and its yield potential and quality characteristics are tightly linked with the global food security. In this study, genome-wide association study (GWAS) was conducted for yield and quality-related traits. Based on the high-density wheat 90K Illumina iSelect SNP Array, 192 bread wheat lines from southwest China, including 25 synthetic hexaploid wheat lines, 80 landraces, and 87 cultivars were analyzed. Association analysis results indicated that there were 57, 27, 30, and 34 SNPs associated with plant height (PH), grain protein content (GPC), thousand kernel weight (TKW), and SDS sedimentation volume (SSV) have been detected, respectively. Then, integrating RNA-Seq with bioinformatics analysis, 246 candidate genes (102 for GPC, 52 for TKW, and 92 for SSV) were found. Further analysis indicated that one up-regulated and two down-regulated expression genes affect GPC. Additionally, two haplotypes significantly affecting PH were detected in a 2.2-Mb genome region encompassing a gene which encoded an ubiquitin-specific protease, TaUBP24. The functional markers of TaUBP24 have been developed, which could be used for marker-assisted selection to improve wheat quality and yield.     

A genome-wide association study of immune response traits in Canadian Holstein cattle

Background

Breeding for enhanced immune response (IR) has been suggested as a tool to improve inherent animal health. Dairy cows with superior antibody-mediated (AMIR) and cell-mediated immune responses (CMIR) have been demonstrated to have a lower occurrence of many diseases including mastitis. Adaptive immune response traits are heritable, and it is, therefore, possible to breed for improved IR, decreasing the occurrence of disease. The objective of this study was to perform genome-wide association studies to determine differences in genetic profiles among Holstein cows classified as High or Low for AMIR and CMIR. From a total of 680 cows with immune response phenotypes, 163 cows for AMIR (81 High and 82 Low) and 140 for CMIR (75 High and 65 Low) were selectively genotyped using the Illumina Bovine SNP50 BeadChip. Results were validated using an unrelated population of 164 Holstein bulls IR phenotyped for AMIR and 146 for CMIR.

Results

A generalized quasi likelihood score method was used to determine single nucleotide polymorphisms (SNP) and chromosomal regions associated with immune response. After applying a 5% chromosomal false discovery rate, 186 SNPs were significantly associated with AMIR. The majority (93%) of significant markers were on chromosome 23, with a similar peak found in the bull population. For CMIR, 21 SNP markers remained significant. Candidate genes within 250,000 base pairs of significant SNPs were identified to determine biological pathways associated with AMIR and CMIR. Various pathways were identified, including the antigen processing and presentation pathway, important in host defense. Candidate genes included those within the bovine Major Histocompatability Complex such as BoLA-DQ, BoLA-DR and the non-classical BoLA-NC1 for AMIR and BoLA-DQ for CMIR, the complement system including C2 and C4 for AMIR and C1q for CMIR, and cytokines including IL-17A, IL17F for AMIR and IL-17RA for CMIR and tumor necrosis factor for both AMIR and CMIR. Additional genes associated with CMIR included galectins 1, 2 and 3, BCL2 and β-defensin.

Conclusions

The significant genetic variation associated with AMIR and CMIR in this study may imply feasibility to include immune response in genomic breeding indices as an approach to improve inherent animal health.     

Evolutionarily unstable fitness maxima and stable fitness minima of continuous traits

Summary We present models of adaptive change in continuous traits for the following situations: (1) adaptation of a single trait within a single population in which the fitness of a given individual depends on the population's mean trait value as well as its own trait value; (2) adaptation of two (or more) traits within a single population; (3) adaptation in two or more interacting species. We analyse a dynamic model of these adaptive scenarios in which the rate of change of the mean trait value is an increasing function of the fitness gradient (i.e. the rate of increase of individual fitness with the individual's trait value). Such models have been employed in evolutionary game theory and are often appropriate both for the evolution of quantitative genetic traits and for the behavioural adjustment of phenotypically plastic traits. The dynamics of the adaptation of several different ecologically important traits can result in characters that minimize individual fitness and can preclude evolution towards characters that maximize individual fitness. We discuss biological circumstances that are likely to produce such adaptive failures for situations involving foraging, predator avoidance, competition and coevolution. The results argue for greater attention to dynamical stability in models of the evolution of continuous traits.     

Adiponectin concentrations: a genome-wide association study

Adiponectin is associated with obesity and insulin resistance. To date, there has been no genome-wide association study (GWAS) of adiponectin levels in Asians. Here we present a GWAS of a cohort of Korean volunteers. A total of 4,001 subjects were genotyped by using a genome-wide marker panel in a two-stage design (979 subjects initially and 3,022 in a second stage). Another 2,304 subjects were used for follow-up replication studies with selected markers. In the discovery phase, the top SNP associated with mean log adiponectin was rs3865188 in CDH13 on chromosome 16 (p = 1.69 × 10(-15) in the initial sample, p = 6.58 × 10(-39) in the second genome-wide sample, and p = 2.12 × 10(-32) in the replication sample). The meta-analysis p value for rs3865188 in all 6,305 individuals was 2.82 × 10(-83). The association of rs3865188 with high-molecular-weight adiponectin (p = 7.36 × 10(-58)) was even stronger in the third sample. A reporter assay that evaluated the effects of a CDH13 promoter SNP in complete linkage disequilibrium with rs3865188 revealed that the major allele increased expression 2.2-fold. This study clearly shows that genetic variants in CDH13 influence adiponectin levels in Korean adults.