A Genome-Wide Search for Linkage of Estimated Glomerular Filtration Rate (eGFR) in the Family Investigation of Nephropathy and Diabetes (FIND)

Objective

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.

Methods

Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.

Results

The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10⁻⁵) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10⁻⁴) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10⁻⁴) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.

Conclusion

The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.     

Genome Scan for Locus Involved in Mandibular Prognathism in Pedigrees from China

Background

It is well known that genetic components play an important role in the etiology of mandibular prognathism, but few susceptibility loci have been mapped.

Methodology

In order to identify linkage regions for mandibular prognathism, we analyzed two Chinese pedigrees with 6,090 genome-wide single-nucleotide polymorphism (SNP) markers from Illumina Linkage-12 DNA Analysis Kit (average spacing 0.58 cM). Multipoint parametric and non-parametric (model-free) linkage analyses were used for the pedigrees.

Principal Finding

The most statistically significant linkage results were with markers on chromosome 4 (LOD  = 3.166 and NPL = 3.65 with rs 875864, 4p16.1, 8.38 cM). Candidate genes within the 4p16.1 include EVC, EVC2.

Conclusion

We detected a novel suggestive linkage locus for mandibular prognathism in two Chinese pedigrees, and this linkage region provides target for susceptibility gene identification, a process that will provide important insights into the molecular and cellular basis of mandibular prognathism.     

Joint QTL linkage mapping for multiple-cross mating design sharing one common parent

Background

Nested association mapping (NAM) is a novel genetic mating design that combines the advantages of linkage analysis and association mapping. This design provides opportunities to study the inheritance of complex traits, but also requires more advanced statistical methods. In this paper, we present the detailed algorithm of a QTL linkage mapping method suitable for genetic populations derived from NAM designs. This method is called joint inclusive composite interval mapping (JICIM). Simulations were designed on the detected QTL in a maize NAM population and an Arabidopsis NAM population so as to evaluate the efficiency of the NAM design and the JICIM method.

Principal Findings

Fifty-two QTL were identified in the maize population, explaining 89% of the phenotypic variance of days to silking, and nine QTL were identified in the Arabidopsis population, explaining 83% of the phenotypic variance of flowering time. Simulations indicated that the detection power of these identified QTL was consistently high, especially for large-effect QTL. For rare QTL having significant effects in only one family, the power of correct detection within the 5 cM support interval was around 80% for 1-day effect QTL in the maize population, and for 3-day effect QTL in the Arabidopsis population. For smaller-effect QTL, the power diminished, e.g., it was around 50% for maize QTL with an effect of 0.5 day. When QTL were linked at a distance of 5 cM, the likelihood of mapping them as two distinct QTL was about 70% in the maize population. When the linkage distance was 1 cM, they were more likely mapped as one single QTL at an intermediary position.

Conclusions

Because it takes advantage of the large genetic variation among parental lines and the large population size, NAM is a powerful multiple-cross design for complex trait dissection. JICIM is an efficient and specialty method for the joint QTL linkage mapping of genetic populations derived from the NAM design.     

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population

Background

A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses.

Results

Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results.

Conclusions

The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1068) contains supplementary material, which is available to authorized users.     

Quantitative Trait Loci for Interhemispheric Commissure Development and Social Behaviors in the BTBR T+ tf/J Mouse Model of Autism

Background

Autism and Agenesis of the Corpus Callosum (AgCC) are interrelated behavioral and anatomic phenotypes whose genetic etiologies are incompletely understood. We used the BTBR T⁺ tf/J (BTBR) strain, exhibiting fully penetrant AgCC, a diminished hippocampal commissure, and abnormal behaviors that may have face validity to autism, to study the genetic basis of these disorders.

Methods

We generated 410 progeny from an F2 intercross between the BTBR and C57BL/6J strains. The progeny were phenotyped for social behaviors (as juveniles and adults) and commisural morphology, and genotyped using 458 markers. Quantitative trait loci (QTL) were identified using genome scans; significant loci were fine-mapped, and the BTBR genome was sequenced and analyzed to identify candidate genes.

Results

Six QTL meeting genome-wide significance for three autism-relevant behaviors in BTBR were identified on chromosomes 1, 3, 9, 10, 12, and X. Four novel QTL for commissural morphology on chromosomes 4, 6, and 12 were also identified. We identified a highly significant QTL (LOD score = 20.2) for callosal morphology on the distal end of chromosome 4.

Conclusions

We identified several QTL and candidate genes for both autism-relevant traits and commissural morphology in the BTBR mouse. Twenty-nine candidate genes were associated with synaptic activity, axon guidance, and neural development. This is consistent with a role for these processes in modulating white matter tract development and aspects of autism-relevant behaviors in the BTBR mouse. Our findings reveal candidate genes in a mouse model that will inform future human and preclinical studies of autism and AgCC.     

Concordance analysis for QTL detection in dairy cattle: a case study of leg morphology

Background

The present availability of sequence data gives new opportunities to narrow down from QTL (quantitative trait locus) regions to causative mutations. Our objective was to decrease the number of candidate causative mutations in a QTL region. For this, a concordance analysis was applied for a leg conformation trait in dairy cattle. Several QTL were detected for which the QTL status (homozygous or heterozygous for the QTL) was inferred for each individual. Subsequently, the inferred QTL status was used in a concordance analysis to reduce the number of candidate mutations.

Methods

Twenty QTL for rear leg set side view were mapped using Bayes C. Marker effects estimated during QTL mapping were used to infer the QTL status for each individual. Subsequently, polymorphisms present in the QTL regions were extracted from the whole-genome sequences of 71 Holstein bulls. Only polymorphisms for which the status was concordant with the QTL status were kept as candidate causative mutations.

Results

QTL status could be inferred for 15 of the 20 QTL. The number of concordant polymorphisms differed between QTL and depended on the number of QTL statuses that could be inferred and the linkage disequilibrium in the QTL region. For some QTL, the concordance analysis was efficient and narrowed down to a limited number of candidate mutations located in one or two genes, while for other QTL a large number of genes contained concordant polymorphisms.

Conclusions

For regions for which the concordance analysis could be performed, we were able to reduce the number of candidate mutations. For part of the QTL, the concordant analyses narrowed QTL regions down to a limited number of genes, of which some are known for their role in limb or skeletal development in humans and mice. Mutations in these genes are good candidates for QTN (quantitative trait nucleotides) influencing rear leg set side view.     

A dominant X-linked QTL regulating pubertal timing in mice found by whole genome scanning and modified interval-specific congenic strain analysis

Background

Pubertal timing in mammals is triggered by reactivation of the hypothalamic-pituitary-gonadal (HPG) axis and modulated by both genetic and environmental factors. Strain-dependent differences in vaginal opening among inbred mouse strains suggest that genetic background contribute significantly to the puberty timing, although the exact mechanism remains unknown.

Methodology/Principal Findings

We performed a genome-wide scanning for linkage in reciprocal crosses between two strains, C3H/HeJ (C3H) and C57BL6/J (B6), which differed significantly in the pubertal timing. Vaginal opening (VO) was used to characterize pubertal timing in female mice, and the age at VO of all female mice (two parental strains, F1 and F2 progeny) was recorded. A genome-wide search was performed in 260 phenotypically extreme F2 mice out of 464 female progeny of the F1 intercrosses to identify quantitative trait loci (QTLs) controlling this trait. A QTL significantly associated was mapped to the DXMit166 marker (15.5 cM, LOD = 3.86, p<0.01) in the reciprocal cross population (C3HB6F2). This QTL contributed 2.1 days to the timing of VO, which accounted for 32.31% of the difference between the original strains. Further study showed that the QTL was B6-dominant and explained 10.5% of variation to this trait with a power of 99.4% at an alpha level of 0.05.The location of the significant ChrX QTL found by genome scanning was then fine-mapped to a region of ∼2.5 cM between marker DXMit68 and rs29053133 by generating and phenotyping a panel of 10 modified interval-specific congenic strains (mISCSs).

Conclusions/Significance

Such findings in our study lay a foundation for positional cloning of genes regulating the timing of puberty, and also reveal the fact that chromosome X (the sex chromosome) does carry gene(s) which take part in the regulative pathway of the pubertal timing in mice.     

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

Background

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

Methods

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

Results

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

Conclusions

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.     

Genetic Susceptibility to Non-Necrotizing Erysipelas/Cellulitis

Background

Bacterial non-necrotizing erysipelas and cellulitis are often recurring, diffusely spreading infections of the skin and subcutaneous tissues caused most commonly by streptococci. Host genetic factors influence infection susceptibility but no extensive studies on the genetic determinants of human erysipelas exist.

Methods

We performed genome-wide linkage with the 10,000 variant Human Mapping Array (HMA10K) array on 52 Finnish families with multiple erysipelas cases followed by microsatellite fine mapping of suggestive linkage peaks. A scan with the HMA250K array was subsequently performed with a subset of cases and controls.

Results

Significant linkage was found at 9q34 (nonparametric multipoint linkage score (NPL_all) 3.84, p = 0.026), which is syntenic to a quantitative trait locus for susceptibility to group A streptococci infections on chromosome 2 in mouse. Sequencing of candidate genes in the 9q34 region did not conclusively associate any to erysipelas/cellulitis susceptibility. Suggestive linkage (NPL_all>3.0) was found at three loci: 3q22-24, 21q22, and 22q13. A subsequent denser genome scan with the HMA250K array supported the 3q22 locus, in which several SNPs in the promoter of AGTR1 (Angiotensin II receptor type I) suggestively associated with erysipelas/cellulitis susceptibility.

Conclusions

Specific host genetic factors may cause erysipelas/cellulitis susceptibility in humans.     

Genetic dynamics underlying phenotypic development of biomass yield in triticale

Background

The nature of dynamic traits with their phenotypic plasticity suggests that they are under the control of a dynamic genetic regulation. We employed a precision phenotyping platform to non-invasively assess biomass yield in a large mapping population of triticale at three developmental stages.

Results

Using multiple-line cross QTL mapping we identified QTL for each of these developmental stages which explained a considerable proportion of the genotypic variance. Some QTL were identified at each developmental stage and thus contribute to biomass yield throughout the studied developmental phases. Interestingly, we also observed QTL that were only identified for one or two of the developmental stages illustrating a temporal contribution of these QTL to the trait. In addition, epistatic QTL were detected and the epistatic interaction landscape was shown to dynamically change with developmental progression.

Conclusions

In summary, our results reveal the temporal dynamics of the genetic architecture underlying biomass accumulation in triticale and emphasize the need for a temporal assessment of dynamic traits.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-458) contains supplementary material, which is available to authorized users.     

Accounting for trait architecture in genomic predictions of US Holstein cattle using a weighted realized relationship matrix

Background

Genomic BLUP (GBLUP) can predict breeding values for non-phenotyped individuals based on the identity-by-state genomic relationship matrix (G). The G matrix can be constructed from thousands of markers spread across the genome. The strongest assumption of G and consequently of GBLUP is that all markers contribute equally to the genetic variance of a trait. This assumption is violated for traits that are controlled by a small number of quantitative trait loci (QTL) or individual QTL with large effects. In this paper, we investigate the performance of using a weighted genomic relationship matrix (wG) that takes into consideration the genetic architecture of the trait in order to improve predictive ability for a wide range of traits. Multiple methods were used to calculate weights for several economically relevant traits in US Holstein dairy cattle. Predictive performance was tested by k-means cross-validation.

Results

Relaxing the GBLUP assumption of equal marker contribution by increasing the weight that is given to a specific marker in the construction of the trait-specific G resulted in increased predictive performance. The increase was strongest for traits that are controlled by a small number of QTL (e.g. fat and protein percentage). Furthermore, bias in prediction estimates was reduced compared to that resulting from the use of regular G. Even for traits with low heritability and lower general predictive performance (e.g. calving ease traits), weighted G still yielded a gain in accuracy.

Conclusions

Genomic relationship matrices weighted by marker realized variance yielded more accurate and less biased predictions for traits regulated by few QTL. Genome-wide association analyses were used to derive marker weights for creating weighted genomic relationship matrices. However, this can be cumbersome and prone to low stability over generations because of erosion of linkage disequilibrium between markers and QTL. Future studies may include other sources of information, such as functional annotation and gene networks, to better exploit the genetic architecture of traits and produce more stable predictions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0100-1) contains supplementary material, which is available to authorized users.     

Detection and validation of stay-green QTL in post-rainy sorghum involving widely adapted cultivar,M35-1 and a popular stay-green genotype B35

Background

Sorghum [Sorghum bicolor (L.) Moench] is an important dry-land cereal of the world providing food, fodder, feed and fuel. Stay-green (delayed-leaf senescence) is a key attribute in sorghum determining its adaptation to terminal drought stress. The objective of this study was to validate sorghum stay-green quantitative trait loci (QTL) identified in the past, and to identify new QTL in the genetic background of a post-rainy adapted genotype M35-1.

Results

A genetic linkage map based on 245 F₉ Recombinant Inbred Lines (RILs) derived from a cross between M35-1 (more senescent) and B35 (less senescent) with 237 markers consisting of 174 genomic, 60 genic and 3 morphological markers was used. The phenotypic data collected for three consecutive post-rainy crop seasons on the RIL population (M35-1 × B35) was used for QTL analysis. Sixty-one QTL were identified for various measures of stay-green trait and each trait was controlled by one to ten QTL. The phenotypic variation explained by each QTL ranged from 3.8 to 18.7%. Co-localization of QTL for more than five traits was observed on two linkage groups i.e. on SBI-09-3 flanked by S18 and Xgap206 markers and, on SBI-03 flanked by XnhsbSFCILP67 and Xtxp31. QTL identified in this study were stable across environments and corresponded to sorghum stay-green and grain yield QTL reported previously. Of the 60 genic SSRs mapped, 14 were closely linked with QTL for ten traits. A genic marker, XnhsbSFCILP67 (Sb03g028240) encoding Indole-3-acetic acid-amido synthetase GH3.5, was co-located with QTL for GLB, GLM, PGLM and GLAM on SBI-03. Genes underlying key enzymes of chlorophyll metabolism were also found in the stay-green QTL regions.

Conclusions

We validated important stay-green QTL reported in the past in sorghum and detected new QTL influencing the stay-green related traits consistently. Stg2, Stg3 and StgB were prominent in their expression. Collectively, the QTL/markers identified are likely candidates for subsequent verification for their involvement in stay-green phenotype using NILs and to develop drought tolerant sorghum varieties through marker-assisted breeding for terminal drought tolerance in sorghum.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-909) contains supplementary material, which is available to authorized users.     

Evaluation of Five Candidate Genes from GWAS for Association with Oligozoospermia in a Han Chinese Population

Background

Oligozoospermia is one of the severe forms of idiopathic male infertility. However, its pathology is largely unknown, and few genetic factors have been defined. Our previous genome-wide association study (GWAS) has identified four risk loci for non-obstructive azoospermia (NOA).

Objective

To investigate the potentially functional genetic variants (including not only common variants, but also less-common and rare variants) of these loci on spermatogenic impairment, especially oligozoospermia.

Design, Setting, and Participants

A total of 784 individuals with oligozoospermia and 592 healthy controls were recruited to this study from March 2004 and January 2011.

Measurements

We conducted a two-stage study to explore the association between oligozoospermia and new makers near NOA risk loci. In the first stage, we used next generation sequencing (NGS) in 96 oligozoospermia cases and 96 healthy controls to screen oligozoospermia-susceptible genetic variants. Next, we validated these variants in a large cohort containing 688 cases and 496 controls by SNPscan for high-throughput Single Nucleotide Polymorphism (SNP) genotyping.

Results and Limitations

Totally, we observed seven oligozoospermia associated variants (rs3791185 and rs2232015 in PRMT6, rs146039840 and rs11046992 in Sox5, rs1129332 in PEX10, rs3197744 in SIRPA, rs1048055 in SIRPG) in the first stage. In the validation stage, rs3197744 in SIRPA and rs11046992 in Sox5 were associated with increased risk of oligozoospermia with an odds ratio (OR) of 4.62 (P  =  0.005, 95%CI 1.58-13.4) and 1.82 (P  =  0.005, 95%CI 1.01-1.64), respectively. Further investigation in larger populations and functional characterizations are needed to validate our findings.

Conclusions

Our study provides evidence of independent oligozoospermia risk alleles driven by variants in the potentially functional regions of genes discovered by GWAS. Our findings suggest that integrating sequence data with large-scale genotyping will serve as an effective strategy for discovering risk alleles in the future.     

Different models of genetic variation and their effect on genomic evaluation

Background

The theory of genomic selection is based on the prediction of the effects of quantitative trait loci (QTL) in linkage disequilibrium (LD) with markers. However, there is increasing evidence that genomic selection also relies on "relationships" between individuals to accurately predict genetic values. Therefore, a better understanding of what genomic selection actually predicts is relevant so that appropriate methods of analysis are used in genomic evaluations.

Methods

Simulation was used to compare the performance of estimates of breeding values based on pedigree relationships (Best Linear Unbiased Prediction, BLUP), genomic relationships (gBLUP), and based on a Bayesian variable selection model (Bayes B) to estimate breeding values under a range of different underlying models of genetic variation. The effects of different marker densities and varying animal relationships were also examined.

Results

This study shows that genomic selection methods can predict a proportion of the additive genetic value when genetic variation is controlled by common quantitative trait loci (QTL model), rare loci (rare variant model), all loci (infinitesimal model) and a random association (a polygenic model). The Bayes B method was able to estimate breeding values more accurately than gBLUP under the QTL and rare variant models, for the alternative marker densities and reference populations. The Bayes B and gBLUP methods had similar accuracies under the infinitesimal model.

Conclusions

Our results suggest that Bayes B is superior to gBLUP to estimate breeding values from genomic data. The underlying model of genetic variation greatly affects the predictive ability of genomic selection methods, and the superiority of Bayes B over gBLUP is highly dependent on the presence of large QTL effects. The use of SNP sequence data will outperform the less dense marker panels. However, the size and distribution of QTL effects and the size of reference populations still greatly influence the effectiveness of using sequence data for genomic prediction.     

Empirical and deterministic accuracies of across-population genomic prediction

Background

Differences in linkage disequilibrium and in allele substitution effects of QTL (quantitative trait loci) may hinder genomic prediction across populations. Our objective was to develop a deterministic formula to estimate the accuracy of across-population genomic prediction, for which reference individuals and selection candidates are from different populations, and to investigate the impact of differences in allele substitution effects across populations and of the number of QTL underlying a trait on the accuracy.

Methods

A deterministic formula to estimate the accuracy of across-population genomic prediction was derived based on selection index theory. Moreover, accuracies were deterministically predicted using a formula based on population parameters and empirically calculated using simulated phenotypes and a GBLUP (genomic best linear unbiased prediction) model. Phenotypes of 1033 Holstein-Friesian, 105 Groninger White Headed and 147 Meuse-Rhine-Yssel cows were simulated by sampling 3000, 300, 30 or 3 QTL from the available high-density SNP (single nucleotide polymorphism) information of three chromosomes, assuming a correlation of 1.0, 0.8, 0.6, 0.4, or 0.2 between allele substitution effects across breeds. The simulated heritability was set to 0.95 to resemble the heritability of deregressed proofs of bulls.

Results

Accuracies estimated with the deterministic formula based on selection index theory were similar to empirical accuracies for all scenarios, while accuracies predicted with the formula based on population parameters overestimated empirical accuracies by ~25 to 30%. When the between-breed genetic correlation differed from 1, i.e. allele substitution effects differed across breeds, empirical and deterministic accuracies decreased in proportion to the genetic correlation. Using a multi-trait model, it was possible to accurately estimate the genetic correlation between the breeds based on phenotypes and high-density genotypes. The number of QTL underlying the simulated trait did not affect the accuracy.

Conclusions

The deterministic formula based on selection index theory estimated the accuracy of across-population genomic predictions well. The deterministic formula using population parameters overestimated the across-population genomic accuracy, but may still be useful because of its simplicity. Both formulas could accommodate for genetic correlations between populations lower than 1. The number of QTL underlying a trait did not affect the accuracy of across-population genomic prediction using a GBLUP method.     

Genome-Wide Association Study Identified CNP12587 Region Underlying Height Variation in Chinese Females

Introduction

Human height is a highly heritable trait considered as an important factor for health. There has been limited success in identifying the genetic factors underlying height variation. We aim to identify sequence variants associated with adult height by a genome-wide association study of copy number variants (CNVs) in Chinese.

Methods

Genome-wide CNV association analyses were conducted in 1,625 unrelated Chinese adults and sex specific subgroup for height variation, respectively. Height was measured with a stadiometer. Affymetrix SNP6.0 genotyping platform was used to identify copy number polymorphisms (CNPs). We constructed a genomic map containing 1,009 CNPs in Chinese individuals and performed a genome-wide association study of CNPs with height.

Results

We detected 10 significant association signals for height (p<0.05) in the whole population, 9 and 11 association signals for Chinese female and male population, respectively. A copy number polymorphism ({"type":"entrez-protein","attrs":{"text":"CNP12587","term_id":"893597589","term_text":"CNP12587"}}CNP12587, chr18:54081842-54086942, p = 2.41×10⁻⁴) was found to be significantly associated with height variation in Chinese females even after strict Bonferroni correction (p = 0.048). Confirmatory real time PCR experiments lent further support for CNV validation. Compared to female subjects with two copies of the CNP, carriers of three copies had an average of 8.1% decrease in height. An important candidate gene, ubiquitin-protein ligase NEDD4-like (NEDD4L), was detected at this region, which plays important roles in bone metabolism by binding to bone formation regulators.

Conclusions

Our findings suggest the important genetic variants underlying height variation in Chinese.     

Linkage disequilibrium and genome-wide association analysis for anthocyanin pigmentation and fruit color in eggplant

Background

The genome-wide association (GWA) approach represents an alternative to biparental linkage mapping for determining the genetic basis of trait variation. Both approaches rely on recombination to re-arrange the genome, and seek to establish correlations between phenotype and genotype. The major advantages of GWA lie in being able to sample a much wider range of the phenotypic and genotypic variation present, in being able to exploit multiple rounds of historical recombination in many different lineages and to include multiple accessions of direct relevance to crop improvement.

Results

A 191 accessions eggplant (Solanum melongena L.) association panel, comprising a mixture of breeding lines, old varieties and landrace selections originating from Asia and the Mediterranean Basin, was SNP genotyped and scored for anthocyanin pigmentation and fruit color at two locations over two years. The panel formed two major clusters, reflecting geographical provenance and fruit type. The global level of linkage disequilibrium was 3.4 cM. A mixed linear model appeared to be the most appropriate for GWA. A set of 56 SNP locus/phenotype associations was identified and the genomic regions harboring these loci were distributed over nine of the 12 eggplant chromosomes. The associations were compared with the location of known QTL for the same traits.

Conclusion

The GWA mapping approach was effective in validating a number of established QTL and, thanks to the wide diversity captured by the panel, was able to detect a series of novel marker/trait associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-896) contains supplementary material, which is available to authorized users.     

Bayesian neural networks for detecting epistasis in genetic association studies

Background

Discovering causal genetic variants from large genetic association studies poses many difficult challenges. Assessing which genetic markers are involved in determining trait status is a computationally demanding task, especially in the presence of gene-gene interactions.

Results

A non-parametric Bayesian approach in the form of a Bayesian neural network is proposed for use in analyzing genetic association studies. Demonstrations on synthetic and real data reveal they are able to efficiently and accurately determine which variants are involved in determining case-control status. By using graphics processing units (GPUs) the time needed to build these models is decreased by several orders of magnitude. In comparison with commonly used approaches for detecting interactions, Bayesian neural networks perform very well across a broad spectrum of possible genetic relationships.

Conclusions

The proposed framework is shown to be a powerful method for detecting causal SNPs while being computationally efficient enough to handle large datasets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0368-0) contains supplementary material, which is available to authorized users.     

Genetic architecture of hybrid male sterility in Drosophila: analysis of intraspecies variation for interspecies isolation

Background

The genetic basis of postzygotic isolation is a central puzzle in evolutionary biology. Evolutionary forces causing hybrid sterility or inviability act on the responsible genes while they still are polymorphic, thus we have to study these traits as they arise, before isolation is complete.

Methodology/Principal Findings

Isofemale strains of D. mojavensis vary significantly in their production of sterile F₁ sons when females are crossed to D. arizonae males. We took advantage of the intraspecific polymorphism, in a novel design, to perform quantitative trait locus (QTL) mapping analyses directly on F₁ hybrid male sterility itself. We found that the genetic architecture of the polymorphism for hybrid male sterility (HMS) in the F₁ is complex, involving multiple QTL, epistasis, and cytoplasmic effects.

Conclusions/Significance

The role of extensive intraspecific polymorphism, multiple QTL, and epistatic interactions in HMS in this young species pair shows that HMS is arising as a complex trait in this system. Directional selection alone would be unlikely to maintain polymorphism at multiple loci, thus we hypothesize that directional selection is unlikely to be the only evolutionary force influencing postzygotic isolation.