Detection of quantitative trait loci in Bos indicus and Bos taurus cattle using genome-wide association studies

Background

The apparent effect of a single nucleotide polymorphism (SNP) on phenotype depends on the linkage disequilibrium (LD) between the SNP and a quantitative trait locus (QTL). However, the phase of LD between a SNP and a QTL may differ between Bos indicus and Bos taurus because they diverged at least one hundred thousand years ago. Here, we test the hypothesis that the apparent effect of a SNP on a quantitative trait depends on whether the SNP allele is inherited from a Bos taurus or Bos indicus ancestor.

Methods

Phenotype data on one or more traits and SNP genotype data for 10 181 cattle from Bos taurus, Bos indicus and composite breeds were used. All animals had genotypes for 729 068 SNPs (real or imputed). Chromosome segments were classified as originating from B. indicus or B. taurus on the basis of the haplotype of SNP alleles they contained. Consequently, SNP alleles were classified according to their sub-species origin. Three models were used for the association study: (1) conventional GWAS (genome-wide association study), fitting a single SNP effect regardless of subspecies origin, (2) interaction GWAS, fitting an interaction between SNP and subspecies-origin, and (3) best variable GWAS, fitting the most significant combination of SNP and sub-species origin.

Results

Fitting an interaction between SNP and subspecies origin resulted in more significant SNPs (i.e. more power) than a conventional GWAS. Thus, the effect of a SNP depends on the subspecies that the allele originates from. Also, most QTL segregated in only one subspecies, suggesting that many mutations that affect the traits studied occurred after divergence of the subspecies or the mutation became fixed or was lost in one of the subspecies.

Conclusions

The results imply that GWAS and genomic selection could gain power by distinguishing SNP alleles based on their subspecies origin, and that only few QTL segregate in both B. indicus and B. taurus cattle. Thus, the QTL that segregate in current populations likely resulted from mutations that occurred in one of the subspecies and can have both positive and negative effects on the traits. There was no evidence that selection has increased the frequency of alleles that increase body weight.     

A multi-SNP association test for complex diseases incorporating an optimal P-value threshold algorithm in nuclear families

Background

Genome-wide association studies (GWAS) have become a common approach to identifying single nucleotide polymorphisms (SNPs) associated with complex diseases. As complex diseases are caused by the joint effects of multiple genes, while the effect of individual gene or SNP is modest, a method considering the joint effects of multiple SNPs can be more powerful than testing individual SNPs. The multi-SNP analysis aims to test association based on a SNP set, usually defined based on biological knowledge such as gene or pathway, which may contain only a portion of SNPs with effects on the disease. Therefore, a challenge for the multi-SNP analysis is how to effectively select a subset of SNPs with promising association signals from the SNP set.

Results

We developed the Optimal P-value Threshold Pedigree Disequilibrium Test (OPTPDT). The OPTPDT uses general nuclear families. A variable p-value threshold algorithm is used to determine an optimal p-value threshold for selecting a subset of SNPs. A permutation procedure is used to assess the significance of the test. We used simulations to verify that the OPTPDT has correct type I error rates. Our power studies showed that the OPTPDT can be more powerful than the set-based test in PLINK, the multi-SNP FBAT test, and the p-value based test GATES. We applied the OPTPDT to a family-based autism GWAS dataset for gene-based association analysis and identified MACROD2-AS1 with genome-wide significance (p-value= 2.5 × 10^− 6).

Conclusions

Our simulation results suggested that the OPTPDT is a valid and powerful test. The OPTPDT will be helpful for gene-based or pathway association analysis. The method is ideal for the secondary analysis of existing GWAS datasets, which may identify a set of SNPs with joint effects on the disease.

Electronic supplementary material

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

Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 (HMGA2) as a positional candidate gene

Background

Ear size and shape are distinct conformation characteristics of pig breeds. Previously, we identified a significant quantitative trait locus (QTL) influencing ear surface on pig chromosome 5 in a White Duroc × Erhualian F₂ resource population. This QTL explained more than 17% of the phenotypic variance.

Methods

Four new markers on pig chromosome 5 were genotyped across this F₂ population. RT-PCR was performed to obtain expression profiles of different candidate genes in ear tissue. Standard association test, marker-assisted association test and F-drop test were applied to determine the effects of single nucleotide polymorphisms (SNP) on ear size. Three synthetic commercial lines were also used for the association test.

Results

We refined the QTL to an 8.7-cM interval and identified three positional candidate genes i.e. HMGA2, SOX5 and PTHLH that are expressed in ear tissue. Seven SNP within these three candidate genes were selected and genotyped in the F₂ population. Of the seven SNP, HMGA2 SNP ({"type":"entrez-nucleotide","attrs":{"text":"JF748727","term_id":"347311053","term_text":"JF748727"}}JF748727: g.2836 A > G) showed the strongest association with ear size in the standard association test and marker-assisted association test. With the F-drop test, F value decreased by more than 97% only when the genotypes of HMGA2 g.2836 A > G were included as a fixed effect. Furthermore, the significant association between g.2836 A > G and ear size was also demonstrated in the synthetic commercial Sutai pig line. The haplotype-based association test showed that the phenotypic variance explained by HMGA2 was similar to that explained by the QTL and at a much higher level than by SOX5. More interestingly, HMGA2 is also located within the dog orthologous chromosome region, which has been shown to be associated with ear type and size.

Conclusions

HMGA2 was the closest gene with a potential functional effect to the QTL or marker for ear size on chromosome 5. This study will contribute to identify the causative gene and mutation underlying this QTL.     

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.     

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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

Genome-wide association study of footrot in Texel sheep

Background

This is the first study based on a genome-wide association approach that investigates the links between ovine footrot scores and molecular polymorphisms in Texel sheep using the ovine 50 K SNP array (42 883 SNPs (single nucleotide polymorphisms) after quality control). Our aim was to identify molecular predictors of footrot resistance.

Methods

This study used data from animals selected from a footrot-phenotyped Texel sheep population of 2229 sheep with an average of 1.60 scoring records per animal. From these, a subset of 336 animals with extreme trait values for footrot was selected for genotyping based on their phenotypic records. De-regressed estimated breeding values (EBV) for footrot were used as pseudo-phenotypes in the genome-wide association analysis.

Results

Seven SNPs were significant on a chromosome-wise level but the association analysis did not reveal any genome-wise significant SNPs associated with footrot. Based on the current state of knowledge of the ovine genome, it is difficult to clearly link the function of the genes that contain these significant SNPs with a potential role in resistance/susceptibility to footrot. Linkage disequilibrium (LD) was analysed as one of the factors that influence the power of detecting QTL (quantitative trait loci). A mean LD of 0.20 (r² at a distance of 50 kb between two SNPs) in the population analysed was estimated. LD declined from 0.15 to 0.07 and to 0.04 at distances between two SNPs of 100, 1000 and 2000 kb, respectively.

Conclusions

Based on a relatively small number of genotyped animals, this study is a first step to search for genomic regions that are involved in resistance to footrot using the ovine 50 K SNP array. Seven SNPs were found to be significant on a chromosome-wise level. No major genome-wise significant QTL were identified.     

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

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

Local genealogies in a linear mixed model for genome-wide association mapping in complex pedigreed populations

Introduction

The state-of-the-art for dealing with multiple levels of relationship among the samples in genome-wide association studies (GWAS) is unified mixed model analysis (MMA). This approach is very flexible, can be applied to both family-based and population-based samples, and can be extended to incorporate other effects in a straightforward and rigorous fashion. Here, we present a complementary approach, called ‘GENMIX (genealogy based mixed model)’ which combines advantages from two powerful GWAS methods: genealogy-based haplotype grouping and MMA.

Subjects and Methods

We validated GENMIX using genotyping data of Danish Jersey cattle and simulated phenotype and compared to the MMA. We simulated scenarios for three levels of heritability (0.21, 0.34, and 0.64), seven levels of MAF (0.05, 0.10, 0.15, 0.20, 0.25, 0.35, and 0.45) and five levels of QTL effect (0.1, 0.2, 0.5, 0.7 and 1.0 in phenotypic standard deviation unit). Each of these 105 possible combinations (3 h² x 7 MAF x 5 effects) of scenarios was replicated 25 times.

Results

GENMIX provides a better ranking of markers close to the causative locus'' location. GENMIX outperformed MMA when the QTL effect was small and the MAF at the QTL was low. In scenarios where MAF was high or the QTL affecting the trait had a large effect both GENMIX and MMA performed similarly.

Conclusion

In discovery studies, where high-ranking markers are identified and later examined in validation studies, we therefore expect GENMIX to enrich candidates brought to follow-up studies with true positives over false positives more than the MMA would.     

Mapping of a milk production quantitative trait locus to a 1.056 Mb region on bovine chromosome 5 in the Fleckvieh dual purpose cattle breed

Background

In a previous study in the Fleckvieh dual purpose cattle breed, we mapped a quantitative trait locus (QTL) affecting milk yield (MY1), milk protein yield (PY1) and milk fat yield (FY1) during first lactation to the distal part of bovine chromosome 5 (BTA5), but the confidence interval was too large for positional cloning of the causal gene. Our objective here was to refine the position of this QTL and to define the candidate region for high-throughput sequencing.

Methods

In addition to those previously studied, new Fleckvieh families were genotyped, in order to increase the number of recombination events. Twelve new microsatellites and 240 SNP markers covering the most likely QTL region on BTA5 were analysed. Based on haplotype analysis performed in this complex pedigree, families segregating for the low frequency allele of this QTL (minor allele) were selected. Single- and multiple-QTL analyses using combined linkage and linkage disequilibrium methods were performed.

Results

Single nucleotide polymorphism haplotype analyses on representative family sires and their ancestors revealed that the haplotype carrying the minor QTL allele is rare and most probably originates from a unique ancestor in the mapping population. Analyses of different subsets of families, created according to the results of haplotype analysis and availability of SNP and microsatellite data, refined the previously detected QTL affecting MY1 and PY1 to a region ranging from 117.962 Mb to 119.018 Mb (1.056 Mb) on BTA5. However, the possibility of a second QTL affecting only PY1 at 122.115 Mb was not ruled out.

Conclusion

This study demonstrates that targeting families segregating for a less frequent QTL allele is a useful method. It improves the mapping resolution of the QTL, which is due to the division of the mapping population based on the results of the haplotype analysis and to the increased frequency of the minor allele in the families. Consequently, we succeeded in refining the region containing the previously detected QTL to 1 Mb on BTA5. This candidate region contains 27 genes with unknown or partially known function(s) and is small enough for high-throughput sequencing, which will allow future detailed analyses of candidate genes.     

QTL fine mapping with Bayes C(π): a simulation study

Background

Accurate QTL mapping is a prerequisite in the search for causative mutations. Bayesian genomic selection models that analyse many markers simultaneously should provide more accurate QTL detection results than single-marker models. Our objectives were to (a) evaluate by simulation the influence of heritability, number of QTL and number of records on the accuracy of QTL mapping with Bayes Cπ and Bayes C; (b) estimate the QTL status (homozygous vs. heterozygous) of the individuals analysed. This study focussed on the ten largest detected QTL, assuming they are candidates for further characterization.

Methods

Our simulations were based on a true dairy cattle population genotyped for 38 277 phased markers. Some of these markers were considered biallelic QTL and used to generate corresponding phenotypes. Different numbers of records (4387 and 1500), heritability values (0.1, 0.4 and 0.7) and numbers of QTL (10, 100 and 1000) were studied. QTL detection was based on the posterior inclusion probability for individual markers, or on the sum of the posterior inclusion probabilities for consecutive markers, estimated using Bayes C or Bayes Cπ. The QTL status of the individuals was derived from the contrast between the sums of the SNP allelic effects of their chromosomal segments.

Results

The proportion of markers with null effect (π) frequently did not reach convergence, leading to poor results for Bayes Cπ in QTL detection. Fixing π led to better results. Detection of the largest QTL was most accurate for medium to high heritability, for low to moderate numbers of QTL, and with a large number of records. The QTL status was accurately inferred when the distribution of the contrast between chromosomal segment effects was bimodal.

Conclusions

QTL detection is feasible with Bayes C. For QTL detection, it is recommended to use a large dataset and to focus on highly heritable traits and on the largest QTL. QTL statuses were inferred based on the distribution of the contrast between chromosomal segment effects.     

A gene frequency model for QTL mapping using Bayesian inference

Background

Information for mapping of quantitative trait loci (QTL) comes from two sources: linkage disequilibrium (non-random association of allele states) and cosegregation (non-random association of allele origin). Information from LD can be captured by modeling conditional means and variances at the QTL given marker information. Similarly, information from cosegregation can be captured by modeling conditional covariances. Here, we consider a Bayesian model based on gene frequency (BGF) where both conditional means and variances are modeled as a function of the conditional gene frequencies at the QTL. The parameters in this model include these gene frequencies, additive effect of the QTL, its location, and the residual variance. Bayesian methodology was used to estimate these parameters. The priors used were: logit-normal for gene frequencies, normal for the additive effect, uniform for location, and inverse chi-square for the residual variance. Computer simulation was used to compare the power to detect and accuracy to map QTL by this method with those from least squares analysis using a regression model (LSR).

Results

To simplify the analysis, data from unrelated individuals in a purebred population were simulated, where only LD information contributes to map the QTL. LD was simulated in a chromosomal segment of 1 cM with one QTL by random mating in a population of size 500 for 1000 generations and in a population of size 100 for 50 generations. The comparison was studied under a range of conditions, which included SNP density of 0.1, 0.05 or 0.02 cM, sample size of 500 or 1000, and phenotypic variance explained by QTL of 2 or 5%. Both 1 and 2-SNP models were considered. Power to detect the QTL for the BGF, ranged from 0.4 to 0.99, and close or equal to the power of the regression using least squares (LSR). Precision to map QTL position of BGF, quantified by the mean absolute error, ranged from 0.11 to 0.21 cM for BGF, and was better than the precision of LSR, which ranged from 0.12 to 0.25 cM.

Conclusions

In conclusion given a high SNP density, the gene frequency model can be used to map QTL with considerable accuracy even within a 1 cM region.     

Overlapping mouse subcongenic strains successfully separate two linked body fat QTL on distal MMU 2

Background

Mouse chromosome 2 is linked to growth and body fat phenotypes in many mouse crosses. With the goal to identify the underlying genes regulating growth and body fat on mouse chromosome 2, we developed five overlapping subcongenic strains that contained CAST/EiJ donor regions in a C57BL/6J^hg/hg background (hg is a spontaneous deletion of 500 Kb on mouse chromosome 10). To fine map QTL on distal mouse chromosome 2 a total of 1,712 F2 mice from the five subcongenic strains, plus 278 F2 mice from the HG2D founder congenic strain were phenotyped and analyzed. Interval mapping (IM) and composite IM (CIM) were performed on body weight and body fat traits on a combination of SNP and microsatellite markers, which generated a high-density genotyping panel.

Results

Phenotypic analysis and interval mapping of total fat mass identified two QTL on distal mouse chromosome 2. One QTL between 150 and 161 Mb, Fatq2a, and the second between 173.3 and 175.6 Mb, Fatq2b. The two QTL reside in different congenic strains with significant total fat differences between homozygous cast/cast and b6/b6 littermates. Both of these QTL were previously identified only as a single QTL affecting body fat, Fatq2. Furthermore, through a novel approach referred here as replicated CIM, Fatq2b was mapped to the Gnas imprinted locus.

Conclusions

The integration of subcongenic strains, high-density genotyping, and CIM succesfully partitioned two previously linked QTL 20 Mb apart, and the strongest QTL, Fatq2b, was fine mapped to a ~2.3 Mb region interval encompassing the Gnas imprinted locus.

Electronic supplementary material

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

Multiple pathway-based genetic variations associated with tobacco related multiple primary neoplasms

Background

In order to elucidate a combination of genetic alterations that drive tobacco carcinogenesis we have explored a unique model system and analytical method for an unbiased qualitative and quantitative assessment of gene-gene and gene-environment interactions. The objective of this case control study was to assess genetic predisposition in a biologically enriched clinical model system of tobacco related cancers (TRC), occurring as Multiple Primary Neoplasms (MPN).

Methods

Genotyping of 21 candidate Single Nucleotide Polymorphisms (SNP) from major metabolic pathways was performed in a cohort of 151 MPN cases and 210 cancer-free controls. Statistical analysis using logistic regression and Multifactor Dimensionality Reduction (MDR) analysis was performed for studying higher order interactions among various SNPs and tobacco habit.

Results

Increased risk association was observed for patients with at least one TRC in the upper aero digestive tract (UADT) for variations in SULT1A1 Arg²¹³His, mEH Tyr¹¹³His, hOGG1 Ser³²⁶Cys, XRCC1 Arg²⁸⁰His and BRCA2 Asn³⁷²His. Gene - environment interactions were assessed using MDR analysis. The overall best model by MDR was tobacco habit/p53(Arg/Arg)/XRCC1(Arg³⁹⁹His)/mEH(Tyr¹¹³His) that had highest Cross Validation Consistency (8.3) and test accuracy (0.69). This model also showed significant association using logistic regression analysis.

Conclusion

This is the first Indian study on a multipathway based approach to study genetic susceptibility to cancer in tobacco associated MPN. This approach could assist in planning additional studies for comprehensive understanding of tobacco carcinogenesis.     

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

Background

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

Results

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

Conclusions

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Electronic supplementary material

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

On the association of common and rare genetic variation influencing body mass index: a combined SNP and CNV analysis

Background

As the architecture of complex traits incorporates a widening spectrum of genetic variation, analyses integrating common and rare variation are needed. Body mass index (BMI) represents a model trait, since common variation shows robust association but accounts for a fraction of the heritability. A combined analysis of single nucleotide polymorphisms (SNP) and copy number variation (CNV) was performed using 1850 European and 498 African-Americans from the Study of Addiction: Genetics and Environment. Genetic risk sum scores (GRSS) were constructed using 32 BMI-validated SNPs and aggregate-risk methods were compared: count versus weighted and proxy versus imputation.

Results

The weighted SNP-GRSS constructed from imputed probabilities of risk alleles performed best and was highly associated with BMI (p = 4.3×10⁻¹⁶) accounting for 3% of the phenotypic variance. In addition to BMI-validated SNPs, common and rare BMI/obesity-associated CNVs were identified from the literature. Of the 84 CNVs previously reported, only 21-kilobase deletions on 16p12.3 showed evidence for association with BMI (p = 0.003, frequency = 16.9%), with two CNVs nominally associated with class II obesity, 1p36.1 duplications (OR = 3.1, p = 0.009, frequency 1.2%) and 5q13.2 deletions (OR = 1.5, p = 0.048, frequency 7.7%). All other CNVs, individually and in aggregate, were not associated with BMI or obesity. The combined model, including covariates, SNP-GRSS, and 16p12.3 deletion accounted for 11.5% of phenotypic variance in BMI (3.2% from genetic effects). Models significantly predicted obesity classification with maximum discriminative ability for morbid-obesity (p = 3.15×10⁻¹⁸).

Conclusion

Results show that incorporating validated effect sizes and allelic probabilities improve prediction algorithms. Although rare-CNVs did not account for significant phenotypic variation, results provide a framework for integrated analyses.

Electronic supplementary material

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

Nitric oxide compounds have different effects profiles on human articular chondrocyte metabolism

Introduction

The pathogenesis of osteoarthritis (OA) is characterized by the production of high amounts of nitric oxide (NO), as a consequence of up-regulation of chondrocyte-inducible nitric oxide synthase (iNOS) induced by inflammatory cytokines. NO donors represent a powerful tool for studying the role of NO in the cartilage in vitro. There is no consensus about NO effects on articular cartilage in part because the differences between the NO donors available. The aim of this work is to compare the metabolic profile of traditional and new generation NO donors to see which one points out the osteoarthritic process in the best way.

Methods

Human healthy and OA chondrocytes were isolated from patients undergoing joint replacement surgery, and primary cultured. Cells were stimulated with NO donors (NOC-12 or SNP). NO production was evaluated by the Griess method, and apoptosis was quantified by flow cytometry. Mitochondrial function was evaluated by analysing respiratory chain enzyme complexes, citrate synthase (CS) activities by enzymatic assay, mitochondrial membrane potential (Δψm) by JC-1 using flow cytometry, and ATP levels were measured by luminescence assays. Glucose transport was measured as the uptake of 2-deoxy-[³H]glucose (2-[³H]DG). Statistical analysis was performed using the Mann-Whitney U test.

Results

NOC-12 liberates approximately ten times more NO₂^- than SNP, but the level of cell death induced was not as profound as that produced by SNP. Normal articular chondrocytes stimulated with NOC-12 had reduced activity from complexes I, III y IV, and the mitochondrial mass was increased in these cells. Deleterious effects on ΔΨm and ATP levels were more profound with SNP, and this NO donor was able to reduce 2-[³H]DG levels. Both NO donors had opposite effects on lactate release, SNP diminished the levels and NOC-12 lead to lactate accumulation. OA chondrocytes incorporate significantly more 2-[³H]DG than healthy cells.

Conclusions

These findings suggest that the new generation donors, specifically NOC-12, mimic the OA metabolic process much better than SNP. Previous results using SNP have to be considered prudently since most of the effects observed can be induced by the interactions of secondary products of NO.     

Macrophage Phenotype Is Associated with Disease Severity in Preterm Infants with Chronic Lung Disease

Background

The etiology of persistent lung inflammation in preterm infants with chronic lung disease of prematurity (CLD) is poorly characterized, hampering efforts to stratify prognosis and treatment. Airway macrophages are important innate immune cells with roles in both the induction and resolution of tissue inflammation.

Objectives

To investigate airway innate immune cellular phenotypes in preterm infants with respiratory distress syndrome (RDS) or CLD.

Methods

Bronchoalveolar lavage (BAL) fluid was obtained from term and preterm infants requiring mechanical ventilation. BAL cells were phenotyped by flow cytometry.

Results

Preterm birth was associated with an increase in the proportion of non-classical CD14+/CD16+ monocytes on the day of delivery (58.9±5.8% of total mononuclear cells in preterm vs 33.0±6.1% in term infants, p = 0.02). Infants with RDS were born with significantly more CD36⁺ macrophages compared with the CLD group (70.3±5.3% in RDS vs 37.6±8.9% in control, p = 0.02). At day 3, infants born at a low gestational age are more likely to have greater numbers of CD14⁺ mononuclear phagocytes in the airway (p = 0.03), but fewer of these cells are functionally polarized as assessed by HLA-DR (p = 0.05) or CD36 (p = 0.05) positivity, suggesting increased recruitment of monocytes or a failure to mature these cells in the lung.

Conclusions

These findings suggest that macrophage polarization may be affected by gestational maturity, that more immature macrophage phenotypes may be associated with the progression of RDS to CLD and that phenotyping mononuclear cells in BAL could predict disease outcome.     

Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Background

In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs.

Methods

In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees.

Results

Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL.

Conclusions

Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.     

Detecting parent of origin and dominant QTL in a two-generation commercial poultry pedigree using variance component methodology

Introduction

Variance component QTL methodology was used to analyse three candidate regions on chicken chromosomes 1, 4 and 5 for dominant and parent-of-origin QTL effects. Data were available for bodyweight and conformation score measured at 40 days from a two-generation commercial broiler dam line. One hundred dams were nested in 46 sires with phenotypes and genotypes on 2708 offspring. Linear models were constructed to simultaneously estimate fixed, polygenic and QTL effects. Different genetic models were compared using likelihood ratio test statistics derived from the comparison of full with reduced or null models. Empirical thresholds were derived by permutation analysis.

Results

Dominant QTL were found for bodyweight on chicken chromosome 4 and for bodyweight and conformation score on chicken chromosome 5. Suggestive evidence for a maternally expressed QTL for bodyweight and conformation score was found on chromosome 1 in a region corresponding to orthologous imprinted regions in the human and mouse.

Conclusion

Initial results suggest that variance component analysis can be applied within commercial populations for the direct detection of segregating dominant and parent of origin effects.