Evaluating associations of haplotypes with traits

Haplotype-based breeding: A new insight in crop improvement

Haplotype-based breeding (HBB) is one of the cutting-edge technologies in the realm of crop improvement due to the increasing availability of Single Nucleotide Polymorphisms identified by Next Generation Sequencing technologies. The complexity of the data can be decreased with fewer statistical tests and a lower probability of spurious associations by combining thousands of SNPs into a few hundred haplotype blocks. The presence of strong genomic regions in breeding lines of most crop species facilitates the use of haplotypes to improve the efficiency of genomic and marker-assisted selection. Haplotype-based breeding as a Genomic Assisted Breeding (GAB) approach harnesses the genome sequence data to pinpoint the allelic variation used to hasten the breeding cycle and circumvent the challenges associated with linkage drag. This review article demonstrates ways to identify candidate genes, superior haplotype identification, haplo-pheno analysis, and haplotype-based marker-assisted selection. The crop improvement strategies that utilize superior haplotypes will hasten the breeding progress to safeguard global food security.

Seq2Sat and SatAnalyzer toolkit: Towards comprehensive microsatellite genotyping from sequencing data

Accurate and efficient microsatellite loci genotyping is an essential process in population genetics that is also used in various demographic analyses. Protocols for next-generation sequencing of microsatellite loci enable high-throughput and cross-compatible allele scoring, common issues that are not addressed by conventional capillary-based approaches. To improve this process, we have developed an all-in-one software, called Seq2Sat (sequence to microsatellite), in C++ to support automated microsatellite genotyping. It directly takes raw reads of microsatellite amplicons and conducts read quality control before inferring genotypes based on depth-of-read, read ratio, sequence composition and length. We have also developed a module for sex identification based on sex chromosome-specific locus amplicons. To allow for greater user access and complement autoscoring, we developed SatAnalyzer (microsatellite analyzer), a user-friendly web-based platform that conducts reads-to-report analyses by calling Seq2Sat for genotype autoscoring and produces interactive genotype graphs for manual editing. SatAnalyzer also allows users to troubleshoot multiplex optimization by analysing read quality and distribution across loci and samples in support of high-quality library preparation. To evaluate its performance, we benchmarked our toolkit Seq2Sat/SatAnalyzer against a conventional capillary gel method and existing microsatellite genotyping software, MEGASAT, using two datasets. Results showed that SatAnalyzer can achieve >99.70% genotyping accuracy and Seq2Sat is ~5 times faster than MEGASAT despite many more informative tables and figures being generated. Seq2Sat and SatAnalyzer are freely available on github (https://github.com/ecogenomicscanada/Seq2Sat) and dockerhub (https://hub.docker.com/r/rocpengliu/satanalyzer).

Refining the genetic risk of breast cancer with rare haplotypes and pattern mining

Hundreds of common variants have been found to confer small but significant differences in breast cancer risk, supporting the widely accepted polygenic model of inherited predisposition. Using a novel closed-pattern mining algorithm, we provide evidence that rare haplotypes may refine the association of breast cancer risk with common germline alleles. Our method, called Chromosome Overlap, consists in iteratively pairing chromosomes from affected individuals and looking for noncontiguous patterns of shared alleles. We applied Chromosome Overlap to haplotypes of genotyped SNPs from female breast cancer cases from the UK Biobank at four loci containing common breast cancer-risk SNPs. We found two rare (frequency <0.1%) haplotypes bearing a GWAS hit at 11q13 (hazard ratio = 4.21 and 16.7) which replicated in an independent, European ancestry population at P < 0.05, and another at 22q12 (frequency <0.2%, hazard ratio = 2.58) which expanded the risk pool to noncarriers of a GWAS hit. These results suggest that rare haplotypes (or mutations) may underlie the "synthetic association" of breast cancer risk with at least some common variants.

Fetal hemoglobin-boosting haplotypes of BCL11A gene and HBS1L-MYB intergenic region in the prediction of clinical and hematological outcomes in a cohort of children with sickle cell anemia

Single nucleotide polymorphisms (SNPs) of BCL11A gene and HBS1L-MYB intergenic region (named HMIP-2) affect both fetal hemoglobin (HbF) concentration and clinical outcomes in patients with sickle cell anemia (SCA). However, no previous study has examined the interaction among these SNPs in the regulation of HbF. We examined whether HbF-boosting haplotypes combining alleles of functional SNPs of BCL11A and HMIP-2 were associated with clinical outcomes and hematological parameters, and whether they interact to regulate HbF in a cohort of Brazilian children with SCA. The minor haplotype of BCL11A ("TCA", an allele combination of rs1427407, rs766432, and rs4671393) was associated with higher HbF, hemoglobin and lower reticulocytes count compared to reference haplotype "GAG". The minor haplotype of HMIP-2 ("CGC", an allele combination of rs9399137, rs4895441, and rs9494145) was associated with higher HbF and hemoglobin compared to reference haplotype "TAT". Subjects carrying minor haplotypes showed reduced rate of clinical complications compared to reference haplotypes. Non-carriers of both minor haplotypes for BCL11A and HMIP-2 showed the lowest HbF concentration. Subjects carrying only the minor haplotype of BCL11A showed significantly higher HbF concentration than non-carriers of any minor haplotype, which showed no significant difference compared to subjects carrying only the minor haplotype of HMIP-2. Interestingly, subjects carrying both minor haplotypes of BCL11A ("TCA") and HMIP-2 ("CGC") showed significantly higher HbF levels than subjects carrying only the minor haplotype of BCL11A. Our novel findings suggest that HbF-boosting haplotypes of BCL11A and HMIP-2 can predict clinical outcomes and may interact to regulate HbF in patients with SCA.

FNDC5/Irisin polymorphisms are associated with osteopenia in postmenopausal Mayan-Mestizo women

OBJECTIVE

This study aimed to analyze the association between rs3480 and rs16835198 of FNDC5/Irisin and their haplotypes with variations in bone mineral density (BMD) and osteopenia/osteoporosis in postmenopausal Mayan-Mestizo women.

METHODS

We studied 547 postmenopausal women of Maya-Mestizo origin. BMD was measured in the lumbar spine and total hip by dual-energy X-ray absorptiometry. DNA was obtained from blood leukocytes. rs3480 and rs16835198 of FNDC5/Irisin were studied using real-time PCR allelic discrimination. Differences between the means of BMD according to genotype were analyzed with covariance. Allele frequency differences were assessed by χ² and logistic regression was used to test for associations. Pairwise linkage disequilibrium between polymorphisms was calculated by direct correlation r², and haplotype analysis was conducted.

RESULTS

Under a recessive model, we observed a significant association of rs3480 with the presence of osteopenia at the total hip and femoral neck (p = 0.008 and p = 0.003, respectively). For rs16835198, we found an association with osteopenia at the total hip and femoral neck in a dominant model (p = 0.043 and p = 0.009, respectively).

CONCLUSIONS

We found an association of rs3480 with risk to present osteopenia at the total hip and femoral neck, while rs16835198 was associated as a protector for presence of osteopenia only at the femoral neck.

Confirming the TMEM232 gene associated with atopic dermatitis through targeted capture sequencing

Atopic dermatitis (AD) is a common and complex skin disorder, and the 5q22.1 region had been reported to be associated with AD. To confirm the susceptibility gene for AD in the 5q22.1 region by haplotype and targeted capture sequencing. The haplotypes were reconstructed with the genotyping data of four SNPs and six deletions from 3624 Chinese Hans AD patients and 5076 controls. The targeted capture sequencing spanning 5q22.1 region was performed in the selected samples. The gene level enrichment analysis was done using loss of function variants. A total of 62 haplotypes were found, and the H15 haplotype had the strongest association with AD (P = 3.92 × 10^-10, OR 0.17, 95% CI 0.09-0.32). However, no co-segregation mutation sites were found in the sequencing analysis within the 16 selected samples, while the enrichment analysis indicated that TMEM232 was significantly associated with AD (P = 7.33 × 10^-5, OR 0.33, 95% CI 0.19-0.58). This study confirms previous findings that the TMEM232 gene is associated with AD by haplotype analysis and targeted capture sequencing.

Statistically efficient association analysis of quantitative traits with haplotypes and untyped SNPs in family studies

BACKGROUND

Associations between haplotypes and quantitative traits provide valuable information about the genetic basis of complex human diseases. Haplotypes also provide an effective way to deal with untyped SNPs. Two major challenges arise in haplotype-based association analysis of family data. First, haplotypes may not be inferred with certainty from genotype data. Second, the trait values within a family tend to be correlated because of common genetic and environmental factors.

RESULTS

To address these challenges, we present an efficient likelihood-based approach to analyzing associations of quantitative traits with haplotypes or untyped SNPs. This approach properly accounts for within-family trait correlations and can handle general pedigrees with arbitrary patterns of missing genotypes. We characterize the genetic effects on the quantitative trait by a linear regression model with random effects and develop efficient likelihood-based inference procedures. Extensive simulation studies are conducted to examine the performance of the proposed methods. An application to family data from the Childhood Asthma Management Program Ancillary Genetic Study is provided. A computer program is freely available.

CONCLUSIONS

Results from extensive simulation studies show that the proposed methods for testing the haplotype effects on quantitative traits have correct type I error rates and are more powerful than some existing methods.

Selecting Closely-Linked SNPs Based on Local Epistatic Effects for Haplotype Construction Improves Power of Association Mapping

Genome-wide association studies (GWAS) have gained central importance for the identification of candidate loci underlying complex traits. Single nucleotide polymorphism (SNP) markers are mostly used as genetic variants for the analysis of genotype-phenotype associations in populations, but closely linked SNPs that are grouped into haplotypes are also exploited. The benefit of haplotype-based GWAS approaches vs. SNP-based approaches is still under debate because SNPs in high linkage disequilibrium provide redundant information. To overcome some constraints of the commonly-used haplotype-based GWAS in which only consecutive SNPs are considered for haplotype construction, we propose a new method called functional haplotype-based GWAS (FH GWAS). FH GWAS is featured by combining SNPs into haplotypes based on the additive and epistatic effects among SNPs. Such haplotypes were termed functional haplotypes (FH). As shown by simulation studies, the FH GWAS approach clearly outperformed the SNP-based approach unless the minor allele frequency of the SNPs making up the haplotypes is low and the linkage disequilibrium between them is high. Applying FH GWAS for the trait flowering time in a large Arabidopsis thaliana population with whole-genome sequencing data revealed its potential empirically. FH GWAS identified all candidate regions which were detected in SNP-based and two other haplotype-based GWAS approaches. In addition, a novel region on chromosome 4 was solely detected by FH GWAS. Thus both the results of our simulation and empirical studies demonstrate that FH GWAS is a promising method and superior to the SNP-based approach even if almost complete genotype information is available.

Genetic Variants in TMEM39A Gene Are Associated with Autoimmune Thyroid Diseases

Transmembrane protein 39A (TMEM39A) gene polymorphisms have been related to various autoimmune diseases, but the relationship between TMEM39A polymorphisms and autoimmune thyroid disease (AITD) remains unknown. This study was aimed at investigating whether the polymorphisms of the TMEM39A were associated with AITD in the Chinese Han population. A case-control study was performed in a total of 906 AITD patients and 744 healthy controls. Four single nucleotide polymorphisms, including rs1132200, rs12492609, rs2282175, and rs7629750, in TMEM39A were examined by polymerase chain reaction-ligase detection reaction. We found that the allele T of rs12492609 in TMEM39A was associated with AITD and Hashimoto's thyroiditis (HT) (p = 0.023; p = 0.028 respectively). Compared with controls, the frequency of haplotype CCA in AITD patients was higher than that in controls (p = 0.036), but the frequency of haplotype CTA in AITD and HT patients was lower than that in controls (p = 0.046; p = 0.047 respectively). In addition, the frequency of allele A at rs7629750 in AITD patients with onset age ≤18 years old was higher than that in AITD patients with onset age ≥19 (p = 0.046). Further, there was an obvious difference in the genotype distributions of rs12492609 and rs7629750 between HT patients with hypothyroidism and those without hypothyroidism (p = 0.034 and p = 0.023, respectively). Our study first reveals that the polymorphisms of the TMEM39A gene are associated with the susceptibility to AITD, especially for early-onset AITD and HT with hypothyroidism.

Bivariate logistic Bayesian LASSO for detecting rare haplotype association with two correlated phenotypes

In genetic association studies, joint modeling of related traits/phenotypes can utilize the correlation between them and thereby provide more power and uncover additional information about genetic etiology. Moreover, detecting rare genetic variants are of current scientific interest as a key to missing heritability. Logistic Bayesian LASSO (LBL) has been proposed recently to detect rare haplotype variants using case-control data, that is, a single binary phenotype. As there is currently no haplotype association method that can handle multiple binary phenotypes, we extend LBL to fill this gap. We develop a bivariate model by using a latent variable to induce correlation between the two outcomes. We carry out extensive simulations to investigate the bivariate LBL and compare with the univariate LBL. The bivariate LBL performs better or similar to the univariate LBL in most settings. It has the highest gain in power when a haplotype is associated with both traits and it affects at least one trait in a direction opposite to the direction of the correlation between the traits. We analyze two data sets-Genetic Analysis Workshop 19 sequence data on systolic and diastolic blood pressures and a genome-wide association data set on lung cancer and smoking and detect several associated rare haplotypes.

E-selectin gene haplotypes are associated with the risk of myocardial infarction

INTRODUCTION

Endothelial dysfunction is one of the most important factors implicated in the pathogenesis of coronary artery disease (CAD). The aim of this study was to investigate the association of the E-selectin gene (SELE) with CAD and CAD-related traits using tagging polymorphisms.

MATERIAL AND METHODS

A total of 379 Polish patients who had undergone angiography were included: 261 patients with angiographically documented CAD, 202 CAD patients without myocardial infarction (CAD/MI(-) group) and 59 patients with myocardial infarction (CAD/MI(+) group) as well as 118 healthy control subjects (non-CAD). Eight tagging single nucleotide polymorphisms (SNPs) in the SELE gene were selected using genotype data from HapMap. Genotyping was performed using PCR-RFLP and PCR-DHPLC methods.

RESULTS

The most common SELE haplotype in this analysis ([C;G;T;C;G;T], 31.2%) showed a negative association with myocardial infarction (MI) (CAD/MI(+) vs. non-CAD) under the additive (p = 0.001), dominant (p = 0.006) and recessive (p = 0.012) model. Two other haplotypes ([C;G;C;C;A;C], [C;A;C;A;G;T], 5.73% and 18.1%, respectively) were also negatively associated with MI under the additive and dominant model. We also found two haplotypes ([T;G;T;C;G;T], [C;G;C;C;A;T], 1.52% and 6.71%, respectively) associated with the risk for MI (CAD/MI(+) vs. CAD/MI(-)), acting in both additive (p = 0.04, p = 0.007, respectively) and dominant (p = 0.04, p = 0.004, respectively) manner. There was no association with either CAD/MI(-) or with severity of CAD expressed as the number of vessels involved.

CONCLUSIONS

Our results suggest that SELE is one of the independent genetic factors modifying the risk of myocardial infarction.

Application of Chaotic Laws to Improve Haplotype Assembly Using Chaos Game Representation

Sequence data are deposited in the form of unphased genotypes and it is not possible to directly identify the location of a particular allele on a specific parental chromosome or haplotype. This study employed nonlinear time series modeling approaches to analyze the haplotype sequences obtained from the NGS sequencing method. To evaluate the chaotic behavior of haplotypes, we analyzed their whole sequences, as well as several subsequences from distinct haplotypes, in terms of the SNP distribution on their chromosomes. This analysis utilized chaos game representation (CGR) followed by the application of two different scaling methods. It was found that chaotic behavior clearly exists in most haplotype subsequences. For testing the applicability of the proposed model, the present research determined the alleles in gap positions and positions with low coverage by using chromosome subsequences in which 10% of each subsequence's alleles are replaced by gaps. After conversion of the subsequences' CGR into the coordinate series, a Local Projection (LP) method predicted the measure of ambiguous positions in the coordinate series. It was discovered that the average reconstruction rate for all input data is more than 97%, demonstrating that applying this knowledge can effectively improve the reconstruction rate of given haplotypes.

DHA intake interacts with ELOVL2 and ELOVL5 genetic variants to influence polyunsaturated fatty acids in human milk

Endogenous synthesis of PUFAs is mediated by genes controlling fatty acid elongases 2 and 5 (ELOVL2 and ELOVL5) and by exogenous DHA intake. Associations between elongases and PUFA levels probably involve genetic variants of ELOVL and changes in DHA intake, but data about their combined effect on PUFA levels are sparse. We hypothesized that each factor would directly affect PUFAs and that interactions between haplotypes and DHA intake would influence PUFAs. We explored four levels of DHA intake in pregnant Chinese Han women and 10 SNPs in the ELOVL genes to determine associations with PUFAs in breast milk. The SNP rs3798713 and 3-SNP haplotype (rs2281591, rs12332786, and rs3798713) in ELOVL2 were associated with linoleic acid (LA) concentrations. However, carriers of the 3-SNP haplotype with higher DHA intake (second quartile: 14.58-43.15 mg/day) had higher concentrations of LA, arachidonic acid, EPA, and DHA compared with the interaction baseline. In ELOVL5, five SNPs (rs2294867, rs9357760, rs2397142, rs209512, and rs12207094) correlated with PUFA changes. Compared with those who had the 5-SNP haplotype C-A-C-G-A and low DHA intake (<14.58 mg/day), carriers with other haplotypes (A-A-C-A-A or C-A-C-A-A) and high DHA intake (≥118.82 mg/day) had increased EPA levels after adjustments for age and BMI. This study showed that maternal genetic variants in ELOVL2 and ELOVL5 were associated with PUFA levels in breast milk and that the combination of SNP haplotypes and higher DHA intake increased PUFA concentrations.

Powerful testing via hierarchical linkage disequilibrium in haplotype association studies

Marginal tests based on individual SNPs are routinely used in genetic association studies. Studies have shown that haplotype-based methods may provide more power in disease mapping than methods based on single markers when, for example, multiple disease-susceptibility variants occur within the same gene. A limitation of haplotype-based methods is that the number of parameters increases exponentially with the number of SNPs, inducing a commensurate increase in the degrees of freedom and weakening the power to detect associations. To address this limitation, we introduce a hierarchical linkage disequilibrium model for disease mapping, based on a reparametrization of the multinomial haplotype distribution, where every parameter corresponds to the cumulant of each possible subset of a set of loci. This hierarchy present in the parameters enables us to employ flexible testing strategies over a range of parameter sets: from standard single SNP analyses through the full haplotype distribution tests, reducing degrees of freedom and increasing the power to detect associations. We show via extensive simulations that our approach maintains the type I error at nominal level and has increased power under many realistic scenarios, as compared to single SNP and standard haplotype-based studies. To evaluate the performance of our proposed methodology in real data, we analyze genome-wide data from the Wellcome Trust Case-Control Consortium.

Main and interaction effects of childhood trauma and the MAOA uVNTR polymorphism on psychopathy

Psychopathy is characterized by callous affect, interpersonal manipulation, a deviant lifestyle, and antisocial behavior. Previous research has linked psychopathic traits to childhood trauma, but also to the upstream variable number tandem repeat (uVNTR) polymorphism of the monoamine oxidase A (MAOA) gene. An interaction between childhood trauma and MAOA genotype has been associated with antisocial behavior, but so far little is known about interaction effects of childhood trauma and the MAOA uVNTR on psychopathy. In order to bridge this gap, we used data of 1531 male and 1265 female twins and their siblings from a Finnish community sample to estimate structural equation models. The psychopathy and childhood trauma constructs were conceptualized as bifactor models with one general and two orthogonal group factors. Data comprised self-reports on childhood trauma and psychopathic traits as well as MAOA uVNTR genotype. In both genders, childhood trauma was associated with the general factor that represents the overarching psychopathy construct, and with the group factor that captures social deviance, but not with the group factor capturing psychopathic core personality traits. Women with a low activity variant of the MAOA uVNTR reported slightly higher levels of psychopathy than those with a high activity allele, but only with respect to the general psychopathy factor. There was no evidence for an interaction effect between MAOA uVNTR genotype and childhood trauma on psychopathy in either gender. Our results suggest that psychopathy in general and social deviance in particular are associated with childhood trauma in men and women, and that psychopathic traits are subject to variation in the MAOA uVNTR genotype in women.

Association between 1p13 polymorphisms and peripheral arterial disease in a Chinese population with diabetes

AIMS/INTRODUCTION

Variants on chromosome 1p13 have been associated with coronary artery disease and acute myocardial infarction risk in different ethnic groups. The present study aimed to investigate the association between 1p13 polymorphisms and the development of peripheral artery disease (PAD) in a Chinese population with type 2 diabetes mellitus.

MATERIALS AND METHODS

1p13 polymorphisms, rs599839, rs646776 and rs12740374, were assessed in a cohort of 882 type 2 diabetes mellitus patients including 440 type 2 diabetes mellitus patients with PAD (DM + PAD group) and 442 patients without PAD (DM group). Genotyping was carried out using TaqMan assay.

RESULTS

Compared with the DM group, the frequencies of the minor G allele of both rs599839 and rs646776 and the minor T allele of rs12740374 decreased (P = 0.013, P = 0.019 and P = 0.005, respectively), and the frequencies of rs599839 AG + GG, rs646776 AG + GG and rs12740374 CT+TT genotypes were statistically significantly decreased as well (P = 0.017, P = 0.011 and P = 0.007, respectively) in the dominant model in the DM + PAD group than in the DM group. Multivariate unconditional logistic regression analyses adjusted for age, glycated hemoglobin, triglyceride, low-density lipoprotein cholesterol, smoking, hypertension, diabetes duration, coronary heart disease and cerebral infarction showed that the genotypic distribution of rs599839 AG + GG, rs646776 AG + GG and rs12740374 CT + TT remained statistically different between the DM and DM + PAD group (P = 0.014, P = 0.003 and P = 0.004, respectively). The frequencies of haplotype GGT were statistically significantly different between groups (P = 0.08).

CONCLUSIONS

The present study strongly supports that genotypes of rs599839, rs646776 and rs12740374 on 1p13 are protective factors for diabetic PAD in a Chinese population. Haplotype GGT generated by rs599839, rs646776 and rs12740374 might also decrease the risk of the disease.

Identifying new susceptibility genes on dopaminergic and serotonergic pathways for the framing effect in decision-making

The framing effect refers the tendency to be risk-averse when options are presented positively but be risk-seeking when the same options are presented negatively during decision-making. This effect has been found to be modulated by the serotonin transporter gene (SLC6A4) and the catechol-o-methyltransferase gene (COMT) polymorphisms, which are on the dopaminergic and serotonergic pathways and which are associated with affective processing. The current study aimed to identify new genetic variations of genes on dopaminergic and serotonergic pathways that may contribute to individual differences in the susceptibility to framing. Using genome-wide association data and the gene-based principal components regression method, we examined genetic variations of 26 genes on the pathways in 1317 Chinese Han participants. Consistent with previous studies, we found that the genetic variations of the SLC6A4 gene and the COMT gene were associated with the framing effect. More importantly, we demonstrated that the genetic variations of the aromatic-L-amino-acid decarboxylase (DDC) gene, which is involved in the synthesis of both dopamine and serotonin, contributed to individual differences in the susceptibility to framing. Our findings shed light on the understanding of the genetic basis of affective decision-making.

DHPLC Elution Patterns of VDR PCR Products Can Predict Prostate Cancer Susceptibility in African American Men

BACKGROUND/AIM

Denaturing high-performance liquid chromatography (DHPLC) is a technique that is used to detect mutations. The aim of the present study was to determine whether DHPLC elution patterns of vitamin D receptor (VDR) gene PCR products can serve as indicators of susceptibility to prostate cancer (PCa) risk.

MATERIALS AND METHODS

DNA samples of PCa cases and controls were screened for mutations and/or polymorphisms in coding exons of VDR gene using DHPLC analysis. Logistic regression, phi-coefficient (ϕ), and Backward Wald models were used to analyze the data.

RESULTS

Similar elution patterns of exons 1, 6, 7 and 9 along with higher prevalence of heteroduplex DNA were observed in PCa samples than in controls. Exons 4 and 8 had highly significant protective effects (p<0.05). Whereas, exons 5, 7, and 9 were perfectly positively correlated with PCa risk (ϕ=1), thus presenting candidate exons significantly associated with susceptibility to PCa.

CONCLUSION

DHPLC elution patterns of the selected exons could be useful to predict susceptibility to develop PCa.

Effects of genetic variants of the bovine WNT8A gene on nine important growth traits in beef cattle

WNT-β-catenin-TCF pathway is involved in carcinogenesis and foetal development. As a member of the WNT gene family, Wnt8A encodes secreted signalling proteins and responds to many biological processes.However, similar research on the effects of genetic variations of Wnt8A gene on growth traits is lacking. Therefore, in this study, polymorphisms of Wnt8A were detected in 396 animals from Chinese Qinchuan cattle using DNA pool sequencing and PCR-RFLP methods. Four novel single-nucleotide polymorphisms (SNPs) of Wnt8A gene were identified, including three mutations in introns (g.T-445C, g.G244C and g.G910A) and one in exon (g.T4922C). Additionally, we examined the associations of four SNPs with growth traits. The results revealed that SNP2 (g.G244C) was significantly associated with shoulder height, hip height, body length, hip width, and body weight (P < 0.05). SNP3 (g.G910A) also displayed notable effects on hip width (P < 0.05). Meanwhile, the haplotype combination CC-GC-GA-CC was strongly associated with heavier, taller and longer animals (P < 0.05). These results show that the Wnt8A gene may be a potential candidate gene, and the SNPs could be used as molecular markers in early marker-assisted selection in beef cattle breeding programmes.

Evaluation of novel SNPs and haplotypes within the <i>ATBF1</i> gene and their effects on economically important production traits in cattle

Abstract. AT motif binding factor 1 (ATBF1) gene can promote the expression level of the growth hormone 1 (GH1) gene by binding to the enhancers of the POU1F1 and PROP1 genes; thus, it affects the growth and development of livestock. Considering that the ATBF1 gene also has a close relationship with the Janus kinase–signal transductor and activator of transcription (JAK–STAT) pathway, the objective of this work was to identify novel single-nucleotide polymorphism (SNP) variations and their association with growth traits in native Chinese cattle breeds. Five novel SNPs within the ATBF1 gene were found in 644 Qinchuan and Jinnan cattle for first time using 25 pairs of screening and genotyping primers. The five novel SNPs were named as AC_000175:g.140344C>G (SNP1), g.146573T>C (SNP2), g.205468C>T (SNP3), g.205575A>G (SNP4) and g.297690C<T (SNP5). Among them, SNP1 and SNP2 were synonymous coding SNPs, while SNP5 was a missense coding SNP, and the other SNPs were intronic. Haplotype analysis found 18 haplotypes in the two breeds, and three and five closely linked loci were revealed in Qinchuan and Jinnan breeds, respectively. Association analysis revealed that SNP1 was significantly associated with the height across the hip in Qinchuan cattle. SNP2 was found to be significantly related to chest circumference and body side length traits in Jinnan cattle. SNP3 was found to have significant associations with four growth traits in Qinchuan cattle. Moreover, the different combined genotypes, SNP1–SNP3, SNP1–SNP4 and SNP2–SNP5 were significantly associated with the growth traits in cattle. These findings indicated that the bovine ATBF1 gene had marked effects on growth traits, and the growth-trait-related loci can be used as DNA markers for maker-assisted selection (MAS) breeding programs in cattle.

Single Marker and Haplotype-Based Association Analysis of Semolina and Pasta Colour in Elite Durum Wheat Breeding Lines Using a High-Density Consensus Map

Association mapping is usually performed by testing the correlation between a single marker and phenotypes. However, because patterns of variation within genomes are inherited as blocks, clustering markers into haplotypes for genome-wide scans could be a worthwhile approach to improve statistical power to detect associations. The availability of high-density molecular data allows the possibility to assess the potential of both approaches to identify marker-trait associations in durum wheat. In the present study, we used single marker- and haplotype-based approaches to identify loci associated with semolina and pasta colour in durum wheat, the main objective being to evaluate the potential benefits of haplotype-based analysis for identifying quantitative trait loci. One hundred sixty-nine durum lines were genotyped using the Illumina 90K Infinium iSelect assay, and 12,234 polymorphic single nucleotide polymorphism (SNP) markers were generated and used to assess the population structure and the linkage disequilibrium (LD) patterns. A total of 8,581 SNPs previously localized to a high-density consensus map were clustered into 406 haplotype blocks based on the average LD distance of 5.3 cM. Combining multiple SNPs into haplotype blocks increased the average polymorphism information content (PIC) from 0.27 per SNP to 0.50 per haplotype. The haplotype-based analysis identified 12 loci associated with grain pigment colour traits, including the five loci identified by the single marker-based analysis. Furthermore, the haplotype-based analysis resulted in an increase of the phenotypic variance explained (50.4% on average) and the allelic effect (33.7% on average) when compared to single marker analysis. The presence of multiple allelic combinations within each haplotype locus offers potential for screening the most favorable haplotype series and may facilitate marker-assisted selection of grain pigment colour in durum wheat. These results suggest a benefit of haplotype-based analysis over single marker analysis to detect loci associated with colour traits in durum wheat.

Multi-strategy genome-wide association studies identify the DCAF16-NCAPG region as a susceptibility locus for average daily gain in cattle

Average daily gain (ADG) is the most economically important trait in beef cattle industry. Using genome-wide association study (GWAS) approaches, previous studies have identified several causal variants within the PLAG1, NCAPG and LCORL genes for ADG in cattle. Multi-strategy GWASs were implemented in this study to improve detection and to explore the causal genes and regions. In this study, we conducted GWASs based on the genotypes of 1,173 Simmental cattle. In the SNP-based GWAS, the most significant SNPs (rs109303784 and rs110058857, P = 1.78 × 10⁻⁷) were identified in the NCAPG intron on BTA6 and explained 4.01% of the phenotypic variance, and the independent and significant SNP (rs110406669, P = 5.18 × 10⁻⁶) explained 3.32% of the phenotypic variance. Similarly, in the haplotype-based GWAS, the most significant haplotype block, Hap-6-N1416 (P = 2.56 × 10⁻⁸), spanned 12.7 kb on BTA6 and explained 4.85% of the phenotypic variance. Also, in the gene-based GWAS, seven significant genes were obtained which included DCAF16 and NCAPG. Moreover, analysis of the transcript levels confirmed that transcripts abundance of NCAPG (P = 0.046) and DCAF16 (P = 0.046) were significantly correlated with the ADG trait. Overall, our results from the multi-strategy GWASs revealed the DCAF16-NCAPG region to be a susceptibility locus for ADG in cattle.

ABO blood type is associated with ovarian reserve in Chinese women with subfertility

Ovarian reserve reflects both the quantity and quality of oocytes available for procreation, and is affected by many known and unknown factors. ABO blood type is related to a number of infertility processes, but it is unclear whether and how ABO blood type affects ovarian reserve. Here, we explored the relationship between ABO blood type and ovarian reserve in Chinese women with subfertility. Day-3 serum follicle-stimulating hormone (FSH) levels and blood type were examined in 14,875 women who underwent IVF or ICSI treatment. Blood type proportions in the patient population were as follows: 30.98% type A, 24.54% type B, 7.57% type AB, and 36.91% type O. A higher percentage of women with diminished ovarian reserve (DOR) were blood type O, while a lower percentage had the B antigen (B and AB). Multiple logistic regression analysis revealed that blood type O was associated with a greater risk of DOR than blood type B and B antigen-positive types. By contrast, the B antigen (B and AB) was associated with a lower incidence of DOR than blood type O. These results suggest that blood type O is a risk factor for DOR while the B antigen (blood type B or AB) is a protective factor for ovarian reserve in Chinese women with subfertility. Further studies are needed to confirm this effect and identify the underlying mechanisms.

Causal Genetic Inference Using Haplotypes as Instrumental Variables

In genomic studies with both genotypes and gene or protein expression profile available, causal effects of gene or protein on clinical outcomes can be inferred through using genetic variants as instrumental variables (IVs). The goal of introducing IV is to remove the effects of unobserved factors that may confound the relationship between the biomarkers and the outcome. A valid inference under the IV framework requires pairwise associations and pathway exclusivity. Among these assumptions, the IV expression association needs to be strong for the casual effect estimates to be unbiased. However, a small number of single nucleotide polymorphisms (SNPs) often provide limited explanation of the variability in the gene or protein expression and can only serve as weak IVs. In this study, we propose to replace SNPs with haplotypes as IVs to increase the variant-expression association and thus improve the casual effect inference of the expression. In the classical two-stage procedure, we developed a haplotype regression model combined with a model selection procedure to identify optimal instruments. The performance of the new method was evaluated through simulations and compared with the IV approaches using observed multiple SNPs. Our results showed the gain of power to detect a causal effect of gene or protein on the outcome using haplotypes compared with using only observed SNPs, under either complete or missing genotype scenarios. We applied our proposed method to a study of the effect of interleukin-1 beta (IL-1β) protein expression on the 90-day survival following sepsis and found that overly expressed IL-1β is likely to increase mortality.

Polymorphism of IL37 gene as a protective factor for autoimmune thyroid disease

Autoimmune thyroid disease (AITD) comprises Graves' disease (GD) and Hashimoto's thyroiditis (HT). IL37 has been recently proved to be a natural suppressor for innate immunity and acquired immunity. Therefore, this study was conducted to identify the association of IL37 genetic polymorphisms with AITD in Chinese Han population. Polymorphisms of rs3811046/rs3811047/rs2723176/rs272186 in the IL37 gene were assessed in a case-control study comprising 701 GD patients, 301 HT patients and 939 controls. Genetic variants were genotyped by multiplex polymerase chain reaction and ligase detection reaction. The frequencies of the minor allele A of rs2723176 and A of rs2723186 were significantly lower in the GD patients than in the controls (P=0.014, OR=0.774; P=0.014, OR=0.777). After gender stratification, the rs3811046 G allele and the rs3811047/rs2723186 A allele were both significantly associated with a decreased risk of GD in female patients (P=0.030, OR=0.777; P=0.023, OR=0.774; P=0.029, OR=0.761). However, none of the four single nucleotide polymorphisms of IL37 gene showed any significant association with HT. Moreover, haplotype analysis revealed the GCG haplotype conferred increased risk for GD as a whole and in female GD patients (OR=1.213; OR=1.320). The ACG haplotype was associated with an increased risk of HT as a whole (OR=1.567) and in male GD patients (OR=1.820). In contrast, the AAA haplotype showed a protective role for GD as a whole (OR=0.760) and in female GD patients (OR=0.765). Our study strongly supports that the IL37 gene variants are associated with the susceptibility to AITD.

ABCB1 gene variants, digoxin and risk of sudden cardiac death in a general population

OBJECTIVE

The ATP-binding cassette B1 (ABCB1) gene encodes P-glycoprotein, a transport protein, which plays an important role in the bioavailability of digoxin. We aimed to investigate the interaction between variants within the ABCB1 gene and digoxin on the risk of sudden cardiac death (SCD).

METHODS

Within the Rotterdam Study, a population-based cohort study in persons 45 years of age and older, we used Cox regression to analyse the association between three polymorphisms that have been associated with digoxin bioavailability, extracted from 1000-Genomes imputed ABCB1 genotypes and the risk of SCD, stratified by digoxin use.

RESULTS

In a total study population of 10,932 persons, 419 SCDs occurred during a median follow-up of 9.8 years. In non-users of digoxin, the risk of SCD was not different across genotypes. In digoxin users, homozygous T allele carriers of C1236T (HR 1.90; 95% CI 1.09 to 3.30; allele frequency 0.43), G2677T (HR 1.89; 95% CI 1.10 to 3.24; allele frequency 0.44) and C3435T (HR 1.72; 95% CI 1.03 to 2.87; allele frequency 0.53) had a significantly increased risk of SCD in a recessive model. Interaction between the ABCB1 polymorphisms and digoxin use was significant for C1236T and G2677T in the age-adjusted and sex-adjusted model.

CONCLUSIONS

In this study, we showed that in digoxin users variant alleles at each of the three loci in the ABCB1 gene were associated with an increased risk of SCD compared with digoxin users with none or one T allele. If replicated, the findings imply that the ABCB1 genotype modifies the risk of cardiac digoxin toxicity.

Tetra-primer ARMS-PCR identified four pivotal genetic variations in bovine PNPLA3 gene and its expression patterns

Patatin-like phospholipase domain-containing protein 3 (PNPLA3), a member of the patatin like phospholipase domain-containing (PNPLA) family, plays an important role in energy balance, fat metabolism regulation, glucose metabolism and fatty liver disease. Tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR) is a new method offering fast detection and extreme simplicity at a negligible cost for SNP genotyping. In this paper, we investigated the genetic variations at different ages of 660 Chinese indigenous cattle belonging to three breeds (QC, NY, JX) and applied T-ARMS-PCR and PCR-RFLP methods to genotype four SNPs, SNP1: g.A2980G, SNP2: g.A2996T, SNP3: g.A36718G, SNP4: g.G36850A. The statistical analyses indicated that these 4 SNPs affected growth traits markedly (P<0.05) in QC population, whereas combined haplotypes were not (P>0.05). The qPCR (quantitative PCR) indicated that bovine PNPLA3 gene was exclusively expressed in fat tissues. Besides, the analysis between SNP and mRNA expression revealed that, in SNP1, the expression of AG was much higher than AA and GG (P<0.05), which was in accordance with the results of growth traits association analysis, while the results of SNP4 was not. These results supported high potential that SNPs of bovine PNPLA3 gene might be utilized as genetic markers in marker-assisted selection (MAS) for Chinese cattle breeding programs.

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.

Analysis of the role of interleukin 6 receptor haplotypes in the regulation of circulating levels of inflammatory biomarkers and risk of coronary heart disease

Variants at the interleukin 6 receptor (IL6R) gene regulate inflammation and are associated with risk of coronary heart disease (CHD). The aim of the present study was to investigate the effects of IL6R haplotypes on circulating levels of inflammatory biomarkers and risk of CHD. We performed a discovery analysis in SHEEP, a myocardial infarction (MI) case control study (n = 2,774) and replicated our results in two large, independent European populations, PROCARDIS, a CHD case control study (n = 7,998), and IMPROVE (n = 3,711) a prospective cardiovascular cohort study. Two major haplotype blocks (rs12083537A/G and rs4075015A/T—block 1; and rs8192282G/A, rs4553185T/C, rs8192284A/C, rs4240872T/C and rs7514452T/C—block 2) were identified in the IL6R gene. IL6R haplotype associations with C-reactive protein (CRP), fibrinogen, IL6, soluble IL6R (sIL6R), IL6, IL8 and TNF-α in SHEEP, CRP and fibrinogen in PROCARDIS and CRP in IMPROVE as well as association with risk of MI and CHD, were analyzed by THESIAS. Haplotypes in block 1 were associated neither with circulating inflammatory biomarkers nor with the MI/CHD risk. Haplotypes in block 2 were associated with circulating levels of CRP, in all three study populations, with fibrinogen in SHEEP and PROCARDIS, with IL8 and sIL6Rin SHEEP and with a modest, non significant, increase (7%) in MI/CHD risk in the three populations studied. Our results indicate that IL6R haplotypes regulate the circulating levels of inflammatory biomarkers. Lack of association with the risk of CHD may be explained by the combined effect of SNPs with opposite effect on the CHD risk, the sample size as well as by structural changes affecting sIL6R stability in the circulation.

A penalized likelihood approach for investigating gene-drug interactions in pharmacogenetic studies

Pharmacogenetics investigates the relationship between heritable genetic variation and the variation in how individuals respond to drug therapies. Often, gene-drug interactions play a primary role in this response, and identifying these effects can aid in the development of individualized treatment regimes. Haplotypes can hold key information in understanding the association between genetic variation and drug response. However, the standard approach for haplotype-based association analysis does not directly address the research questions dictated by individualized medicine. A complementary post-hoc analysis is required, and this post-hoc analysis is usually under powered after adjusting for multiple comparisons and may lead to seemingly contradictory conclusions. In this work, we propose a penalized likelihood approach that is able to overcome the drawbacks of the standard approach and yield the desired personalized output. We demonstrate the utility of our method by applying it to the Scottish Randomized Trial in Ovarian Cancer. We also conducted simulation studies and showed that the proposed penalized method has comparable or more power than the standard approach and maintains low Type I error rates for both binary and quantitative drug responses. The largest performance gains are seen when the haplotype frequency is low, the difference in effect sizes are small, or the true relationship among the drugs is more complex.

Transethnic replication of association of CTG18.1 repeat expansion of TCF4 gene with Fuchs' corneal dystrophy in Chinese implies common causal variant

PURPOSE

To test the association between the CTG18.1 trinucleotide repeat expansion of TCF4 gene and Fuchs' endothelial corneal dystrophy (FECD) in a Chinese population.

METHODS

The trinucleotide repeat polymorphism CTG18.1 was genotyped using short tandem repeat and triplet repeat primed polymerase chain reaction assays in 57 Chinese subjects with FECD and 121 controls. Statistical association of the expanded CTG18.1 allele and 18 single nucleotide polymorphisms (SNPs) across TCF4 with FECD was evaluated. To investigate the linkage disequilibrium structure of the TCF4 region, haplotype analysis was performed on our study subjects and compared with genotyping data of 97 Han Chinese and 85 Caucasians in the 1000 Genomes Project.

RESULTS

The expanded CTG18.1 allele was associated with FECD (P = 4.7 × 10(-14)), with the odds ratio of each copy of the expanded allele estimated to be 66.5 (95% confidence interval: 12.6-350.1). Five TCF4 SNPs showed association with FECD at a nominal level (P < 5.0 × 10(-2)); however, conditional on the expanded CTG18.1 polymorphism, none of the SNPs showed association with FECD. The only haplotype associated with the disease was the one with the expansion at the CTG18.1 locus.

CONCLUSIONS

Transethnic replication of the association between the CTG18.1 repeat expansion in the TCF4 gene and FECD suggests it is a common, causal variant shared in Eurasian populations conferring significant risk for the development of FECD. Our data suggest that the expanded CTG18.1 allele is the main, if not sole, causal variant at this gene locus in the Chinese population.

Integrating dilution-based sequencing and population genotypes for single individual haplotyping

Background

Haplotype information is useful for many genetic analyses and haplotypes are usually inferred using computational approaches. Among such approaches, the importance of single individual haplotyping (SIH), which infers individual haplotypes from sequence fragments, has been increasing with the advent of novel sequencing techniques, such as dilution-based sequencing. These techniques could produce virtual long read fragments by separating DNA fragments into multiple low-concentration aliquots, sequencing and mapping each aliquot, and merging clustered short reads. Although these experimental techniques are sophisticated, they have the problem of producing chimeric fragments whose left and right parts match different chromosomes. In our previous research, we found that chimeric fragments significantly decrease the accuracy of SIH. Although chimeric fragments can be removed by using haplotypes which are determined from pedigree genotypes, pedigree genotypes are generally not available. The length of reads cluster and heterozygous calls were also used to detect chimeric fragments. Although some chimeric fragments will be removed with these features, considerable number of chimeric fragments will be undetected because of the dispersion of the length and the absence of SNPs in the overlapped regions. For these reasons, a general method to detect and remove chimeric fragments is needed.

Results

In this paper, we propose a general method to detect chimeric fragments. The basis of our method is that a chimeric fragment would correspond to an artificial recombinant haplotype and would differ from biological haplotypes. To detect differences from biological haplotypes, we integrated statistical phasing, which is a haplotype inference approach from population genotypes, into our method. We applied our method to two datasets and detected chimeric fragments with high AUC. AUC values of our method are higher than those of just using cluster length and heterozygous calls. We then used multiple SIH algorithm to compare the accuracy of SIH before and after removing the chimeric fragment candidates. The accuracy of assembled haplotypes increased significantly after removing chimeric fragment candidates.

Conclusions

Our method is useful for detecting chimeric fragments and improving SIH accuracy. The Ruby script is available at https://sites.google.com/site/hmatsu1226/software/csp.

Electronic supplementary material

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

Using haplotypes for the prediction of allelic identity to fine-map QTL: characterization and properties

BACKGROUND

Numerous methods have been developed over the last decade to predict allelic identity at unobserved loci between pairs of chromosome segments along the genome. These loci are often unobserved positions tested for the presence of quantitative trait loci (QTL). The main objective of this study was to understand from a theoretical standpoint the relation between linkage disequilibrium (LD) and allelic identity prediction when using haplotypes for fine mapping of QTL. In addition, six allelic identity predictors (AIP) were also compared in this study to determine which one performed best in theory and application.

RESULTS

A criterion based on a simple measure of matrix distance was used to study the relation between LD and allelic identity prediction when using haplotypes. The consistency of this criterion with the accuracy of QTL localization, another criterion commonly used to compare AIP, was evaluated on a set of real chromosomes. For this set of chromosomes, the criterion was consistent with the mapping accuracy of a simulated QTL with either low or high effect. As measured by the matrix distance, the best AIP for QTL mapping were those that best captured LD between a tested position and a QTL. Moreover the matrix distance between a tested position and a QTL was shown to decrease for some AIP when LD increased. However, the matrix distance for AIP with continuous predictions in the [0,1] interval was algebraically proven to decrease less rapidly up to a lower bound with increasing LD in the simplest situations, than the discrete predictor based on identity by state between haplotypes (IBS hap), for which there was no lower bound. The expected LD between haplotypes at a tested position and alleles at a QTL is a quantity that increases naturally when the tested position gets closer to the QTL. This behavior was demonstrated with pig and unrelated human chromosomes.

CONCLUSIONS

When the density of markers is high, and therefore LD between adjacent loci can be assumed to be high, the discrete predictor IBS hap is recommended since it predicts allele identity correctly when taking LD into account.

A new anti-lipopolysaccharide factor (ALF) gene with its SNP polymorphisms related to WSSV-resistance of Litopenaeus vannamei

Anti-lipopolysaccharide factors (ALFs) of crustacean play an important role against bacteria or virus infection. In this study, the cDNA sequence and genomic sequence of one new isoform of ALF designated as nLvALF1 were reported. The open reading frame (ORF) of nLvALF1 consisted of 369 bp encoding 123 amino acids and the genomic structure of nLvALF1 comprised four introns and three exons. The predicted pI of the deduced protein was 8.82 and the molecular weight (MW) was 13.72 KDa. The deduced amino acid sequence of nLvALF1 contained a typical functional domain of ALF: LPS-binding domain. Phylogenetic analysis showed that nLvALF1 had the closest relationship with FcALF1 from Fenneropenaeus chinensis. The nLvALF1 was specifically expressed in lymphoid organ (Oka) of shrimp. Its transcriptional level was significantly up-regulated after white spot syndrome virus (WSSV) challenge, suggesting that nLvALF1 might participate in defense against WSSV in Litopenaeus vannamei. In order to search potential genetic markers associated with WSSV-resistance, we scanned the polymorphisms of the genomic fragment with 397 bp where the LPS-binding domain encoding sequence located and 18 SNPs were found. The distribution frequency of these SNPs was analyzed in WSSV susceptible shrimp and resistant shrimp separately. Significant differences existed in allelic frequencies at loci g.1361-T > C, g.1370-T > C, g.1419-T > A between the WSSV-resistant group and the WSSV-susceptible/susceptible group (P < 0.05). The specific haplotype CT consisted of g.1415-C > A and g.1419-T > A was associated with susceptibility to WSSV (P < 0.05). These findings provide theoretical support for selection of WSSV-resistant varieties of L. vannamei.

Tetra-primer ARMS-PCR identifies the novel genetic variations of bovine HNF-4α gene associating with growth traits

Hepatocyte nuclear factor-4α (HNF-4α), a member of the hepatocyte nuclear factor family, plays an important role in regulating the expression of genes involved in the development, differentiation and normal function of liver and pancreatic β cells, as well as the maintenance of glucose homeostasis. Tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR) is a new method offering fast detection and extreme simplicity at a negligible cost for SNP genotyping. In this paper, we characterize the polymorphisms of the bovine HNF-4α gene in three Chinese indigenous cattle breeds (n=660). Six novel SNPs were identified including 1 mutation in the coding region and others in introns. The statistical analyses indicated that 4 SNPs (g.T53729C, g.A53861G, g.A65188C and g.T65444C) affected growth traits markedly (P<0.05) in Qinchuan cattle (2 years after birth). Besides, haplotypes involving these 4 SNP sites in the bovine HNF-4α gene were identified and their effects on growth traits were also analyzed. The results showed that haplotypes 2, 7, 9 and 11 were predominant and accounted for 73.2%, 59.6%, and 67.1% in Qinchuan, Nanyang and Jiaxian cattle breeds, respectively. Hap9 (TAAT) was extremely predominant in all test populations, which suggested that individuals with Hap9 were more adapted to the environment. Furthermore, 4 combined haplotypes were constructed to guarantee the reliability of analysis results in Qinchuan cattle. There were also significant differences in body length (P<0.05). These findings will benefit for the application of DNA marker related to the growth traits on marker-assisted selection (MAS), and improve the performance of beef cattle.

Detecting local haplotype sharing and haplotype association

A novel haplotype association method is presented, and its power is demonstrated. Relying on a statistical model for linkage disequilibrium (LD), the method first infers ancestral haplotypes and their loadings at each marker for each individual. The loadings are then used to quantify local haplotype sharing between individuals at each marker. A statistical model was developed to link the local haplotype sharing and phenotypes to test for association. We devised a novel method to fit the LD model, reducing the complexity from putatively quadratic to linear (in the number of ancestral haplotypes). Therefore, the LD model can be fitted to all study samples simultaneously, and, consequently, our method is applicable to big data sets. Compared to existing haplotype association methods, our method integrated out phase uncertainty, avoided arbitrariness in specifying haplotypes, and had the same number of tests as the single-SNP analysis. We applied our method to data from the Wellcome Trust Case Control Consortium and discovered eight novel associations between seven gene regions and five disease phenotypes. Among these, GRIK4, which encodes a protein that belongs to the glutamate-gated ionic channel family, is strongly associated with both coronary artery disease and rheumatoid arthritis. A software package implementing methods described in this article is freely available at http://www.haplotype.org.

Haplotype association analysis of combining unrelated case-control and triads with consideration of population stratification

Combining data when data are collected under different study designs, such as family trios and unrelated case-control samples, gains more power and is cost-effective than analyzing each data separately. However, a potential concern is population stratification (PS) among unrelated case-control samples and analyses integrating data should address this confounding effect. In this paper, we develop a simpler method, haplotype generalized linear model (HGLM), that tests and estimates haplotype effects on disease risk and allows for modification against PS for combining data. We proposed to combine information across aggregations of haplotype weighted-counts estimated from population case-control data and trio data separately, and to perform subsequent GLM analysis. Furthermore, we present a framework of analysis of variance based on haplotype weighted-counts for detecting whether it is appropriate to combine two data sources, as well as the modified HGLM with clustering methods for addressing PS. We evaluate the statistical properties in terms of the accuracy, false positive rate (FPR) and empirical power using simulated data with regard to various disease risks, sample sizes, multi-SNP haplotypes and the presence of PS. Our simulation results indicate that HGLM performs comparably well with the likelihood-based haplotype association analysis, particularly when the haplotype effects are moderate, but may not perform well when dealing with lengthy haplotypes for small sample sizes. In the presence of PS, the modified HGLM remains valid and has satisfactory nominal level and small bias. Overall, HGLM appears to be successful in combining data and is simple to implement in standard statistical software.

Identification of bovine NPC1 gene cSNPs and their effects on body size traits of Qinchuan cattle

NPC1 gene is an important gene closely related to the Niemann-Pick type C (NPC). Mutations in the NPC1 gene tend to cause Niemann-Pick type C, a lysosomal storage disorder. Previous studies have shown that NPC1 protein plays an important role in subcellular lipid transport, homeostasis, platelet function and formation, which are basic metabolic activities in the process of development. In this study, to explore the association between the NPC1 gene variation and body size traits in Qinchuan cattle, we detected four novel coding single nucleotide polymorphisms (cSNPs) in the bovine NPC1 gene, including one missense mutation (SNP1) and three synonymous mutations (SNP2, SNP3 and SNP4). Population genetic analyses of 518 individuals and association correlations between cSNPs and bovine body size traits were conducted in this research. A missense mutation at SNP1 locus was found to be significantly related to the heart girth, hip width and body weight (P<0.01 or P<0.05, 3.5-year-old). Two synonymous mutations at SNP2 and SNP3 loci also showed significant effects on hip width (P<0.05, 3.5-year-old). One synonymous mutation at SNP4 locus showed significant effect on body weight (P<0.05, 2.0-year-old). Combined haplotypes H2H6 and H6H6 showed significant effects on body size traits such as heart girth, hip width, and body weight (3.5-year-old, P<0.01 or P<0.05). This study provides evidence that the NPC1 gene might be involved in the regulation of bovine growth and body development, and may be considered as a candidate gene for marker assisted selection (MAS) in beef cattle breeding industry.

DNA mismatch repair MSH2 gene-based SNP associated with different populations

We screened for the major essential single-nucleotide polymorphism (SNP) variant that might be associated with the MSH2 gene based on the data available from three types of human tissue samples [156 lymphoblastoid cell variations (LCL), 160 epidermis, 166 fat]. An association analysis confirmed that the KCNK12 SNP variant (rs748780) was highly associated (p value 9 × 10(-4)) with the MSH2 gene for all three samples. Using SNP identification, we further found that the recognized SNP was also relevant among Hapmap populations. Techniques that display specific SNPs associated with the gene of interest or nearby genes provide more reliable genetic associations than techniques that rely on data from individual SNPs. We investigated the MSH2 gene regional linkage association with the determined SNP (rs748780), KCNK12 variant (Allele T>C) in the intronic region, in HapMap3 full dataset populations, Yoruba in Ibadan, Nigeria (YRI), Utah residents with ancestry from northern Europe (CEU), Han Chinese in Beijing, China (CHB), and a population of Mexican ancestry in Los Angeles, California (MEX). A gene-based SNP association analysis analyzes the combined impact of every variant within the gene while creating referrals to linkage disequilibrium or connections between markers. Our results indicated that among the four populations studied, this association was highest in the MEX population based on the r(2) value; a similar pattern was also observed in the other three populations. The relevant SNP rs748780 in KCNK12 is related to a superfamily of potassium channel pore-forming P-domain proteins as well as to other non-pore-forming proteins and has been shown to be relevant to neurological disorder predisposition in MEX as well as in other populations.

Dental caries and enamelin haplotype

In the literature, the enamelin gene ENAM has been repeatedly designated as a possible candidate for caries susceptibility. Here, we checked whether ENAM variants could increase caries susceptibility. To this aim, we sequenced coding exons and exon-intron boundaries of ENAM in 250 children with a severe caries phenotype and in 149 caries-free patients from 9 French hospital groups. In total, 23 single-nucleotide polymorphisms (SNPs) were found, but none appeared to be responsible for a direct change of ENAM function. Six SNPs had a high minor allele frequency (MAF) and 6 others were identified for the first time. Statistical and evolutionary analyses showed that none of these SNPs was associated with caries susceptibility or caries protection when studied separately and challenged with environmental factors. However, haplotype interaction analysis showed that the presence, in a same variant, of 2 exonic SNPs (rs7671281 and rs3796704; MAF 0.12 and 0.10, respectively), both changing an amino acid in the protein region encoded by exon 10 (p.I648T and p.R763Q, respectively), increased caries susceptibility 2.66-fold independent of the environmental risk factors. These findings support ENAM as a gene candidate for caries susceptibility in the studied population.

Significant interactions between maternal PAH exposure and haplotypes in candidate genes on B[a]P-DNA adducts in a NYC cohort of non-smoking African-American and Dominican mothers and newborns

Polycyclic aromatic hydrocarbons (PAH) are a class of chemicals common in the environment. Certain PAH are carcinogenic, although the degree to which genetic variation influences susceptibility to carcinogenic PAH remains unclear. Also unknown is the influence of genetic variation on the procarcinogenic effect of in utero exposures to PAH. Benzo[a]pyrene (B[a]P) is a well-studied PAH that is classified as a probable human carcinogen. Within our New York City-based cohort, we explored interactions between maternal exposure to airborne PAH during pregnancy and maternal and newborn haplotypes (and in one case, a single-nucleotide polymorphism) in key B[a]P metabolism genes on B[a]P-DNA adducts in paired cord blood samples. The study subjects included non-smoking African-American (n = 132) and Dominican (n = 235) women with available data on maternal PAH exposure, paired cord adducts and genetic data who resided in the Washington Heights, Central Harlem and South Bronx neighborhoods of New York City. We selected seven maternal and newborn genes related to B[a]P metabolism, detoxification and repair for our analyses: CYP1A1, CYP1A2, CYP1B1, GSTM3, GSTT2, NQO1 and XRCC1. We found significant interactions between maternal PAH exposure and haplotype on cord B[a]P-DNA adducts in the following genes: maternal CYP1B1, XRCC1 and GSTM3, and newborn CYP1A2 and XRCC1 in African-Americans; and maternal XRCC1 and newborn NQO1 in Dominicans. These novel findings highlight differences in maternal and newborn genetic contributions to B[a]P-DNA adduct formation, as well as ethnic differences in gene-environment interactions, and have the potential to identify at-risk subpopulations who are susceptible to the carcinogenic potential of B[a]P.

Nonparametric tests of associations with disease based on U-statistics

In case-control studies, association analysis was designed to test whether genetic variants were associated with human diseases. To evaluate the association, analysing one genetic marker at a time suffered from weak power, because of the correction for multiple testing and possibly small genetic effects. An alternative strategy was to test simultaneous effects of multiple markers, which was believed to be more powerful. However, when the number of markers under investigation was large, they would be subjected to weak power as well, because of the greater degrees of freedom. To conquer these limitations in case-control studies, we proposed a novel method that could test joint association of several loci (i.e. haplotype), with only a single degree of freedom. In this research, we developed a nonparametric approach, which was based on U-statistics. We also introduced a new kernel for U-statistic, which could combine the haplotype structure information, and was expected to enhance the power. Simulations indicated that our proposed approach offered merits in identifying the associations between diseases and haplotypes. Application of our method to a study of candidate genes for internalising disorder illustrated its virtue in utility and interpretation, and provided an excellent result in detecting the associations.

Novel polymorphisms of the APOA2 gene and its promoter region affect body traits in cattle

Apolipoprotein A-II (APOA2) is one of the major constituents of high-density lipoprotein and plays a critical role in lipid metabolism and obesity. However, similar research for the bovine APOA2 gene is lacking. In this study, polymorphisms of the bovine APOA2 gene and its promoter region were detected in 1021 cows from four breeds by sequencing and PCR-RFLP methods. Totally, we detected six novel mutations which included one mutation in the promoter region, two mutations in the exons and three mutations in the introns. There were four polymorphisms within APOA2 gene were analyzed. The allele A, T, T and G frequencies of the four loci were predominant in the four breeds when in separate or combinations analysis which suggested cows with those alleles to be more adapted to the steppe environment. The association analysis indicated three SVs in Nangyang cows, two SVs in Qinchun cows and the 9 haplotypes in Nangyang cows were significantly associated with body traits (P<0.05 or P<0.01). The results of this study suggested the bovine APOA2 gene may be a strong candidate gene for body traits in the cattle breeding program.

CYP19 genetic polymorphism haplotype AASA is associated with a poor prognosis in premenopausal women with lymph node-negative, hormone receptor-positive breast cancer

Given the critical role of CYP19 in estrogen synthesis, we investigated the influence of CYP19 gene polymorphisms on the clinical outcome of lymph node- (LN-) negative, hormone receptor- (HR-) positive early breast cancers. Genotyping for the CYP19 polymorphisms rs4646 (A/C), rs1065779 (A/C), CYP19 (TTTA)n (short allele/long (S/L) allele using the 7 TTTA repeat polymorphism as the cut-off), and rs1870050 (A/C) was performed on 296 patients with LN-negative, HR-positive breast cancers. All patients received adjuvant hormonal therapy. Associations were examined between these 4 genotypes and 6 common haplotypes of CYP19 and distant disease-free survival (DDFS), disease-free survival (DFS), and overall survival (OS). Patients were divided into the 6 subhaplotypes of CCLA (41.1%), AASA (17.1%), CASA (11.9%), CCLC (8.9%), CCSA (7.5%), AASC (8.9%), and others (4.6%). In premenopausal patients, haplotype AASA was significantly associated with a poor DDFS (adjusted hazard ratio (aHR), 3.3; P = 0.001), DFS (aHR, 2.5; P = 0.0008), and OS (aHR, 2.9; P = 0.0004) after adjusting for age, tumor size, tumor grade, estrogen receptor status, progesterone receptor status, chemotherapy, pathology, adjuvant hormone therapy, menopausal status, and radiotherapy. Furthermore, haplotype AASA remained a negative prognostic factor for premenopausal patients receiving adjuvant chemotherapy in terms of DDFS (aHR, 4.5; P = 0.0005), DFS (HR, 3.2; P = 0.003), and OS (HR, 6.4; P = 0.0009). However, in postmenopausal patients, haplotype AASA was not associated with a poor prognosis, whereas the AASC haplotype was significantly associated with a poor DFS (aHR, 3.1; P = 0.03) and OS (aHR, 4.4; P = 0.01). Our results indicate that, in patients with LN-negative, HR-positive breast cancers, genetic polymorphism haplotype AASA is associated with poor survival of premenopausal women but does not affect survival of postmenopausal women.

Impact of genetic variants of IL-6, IL6R, LRP5, ESR1 and SP7 genes on bone mineral density in postmenopausal Mexican-Mestizo women with obesity

BACKGROUND

Since obesity and osteoporosis present a high genetic predisposition and polymorphisms of IL-6, IL6R, LRP5, ESR1 and SP7 may influence the risk of both diseases, the aim of this study was to analyze the possible association of polymorphisms in these genes, as well as their haplotypes, with BMD variations in postmenopausal Mexican-Mestizo women with grade 2 or grade 3 obesity.

METHODS

One hundred eighty unrelated postmenopausal women with grade 2 or grade 3 obesity were included. BMD was measured in total hip and lumbar spine by dual-energy X-ray absorptiometry. DNA was obtained from blood leukocytes. Rs1800795 of IL-6, rs2228145 of IL6R, rs3736228 of LRP5, rs9340799 (XbaI) and rs2234693 (PvuII), of ESR1, rs10876432 and rs2016266, of SP7 (and their haplotypes), were studied by real-time PCR allelic discrimination. Deviations from Hardy-Weinberg equilibrium were tested. Pairwise linkage disequilibrium between single nucleotide polymorphisms was calculated by direct correlation r(2), and haplotype analysis was conducted.

RESULTS

Using WHO criteria, 54.5% had grade 2 obesity, and 45.5% had grade 3 obesity. Regarding DXA results, 11.1% women had osteoporosis, 41.7% had osteopenia, and 47.2% had normal BMD. Genotype and haplotype analysis showed no significant differences with BMD variations at the lumbar spine, total hip or femoral neck.

CONCLUSIONS

We did not find a significant association between the polymorphisms analyzed or their haplotypes and BMD variations in postmenopausal women with obesity. The higher BMD observed in women with obesity could be the result of an adaptive response to the higher loading of the skeleton.

An EM algorithm based on an internal list for estimating haplotype distributions of rare variants from pooled genotype data

Background

Pooling is a cost effective way to collect data for genetic association studies, particularly for rare genetic variants. It is of interest to estimate the haplotype frequencies, which contain more information than single locus statistics. By viewing the pooled genotype data as incomplete data, the expectation-maximization (EM) algorithm is the natural algorithm to use, but it is computationally intensive. A recent proposal to reduce the computational burden is to make use of database information to form a list of frequently occurring haplotypes, and to restrict the haplotypes to come from this list only in implementing the EM algorithm. There is, however, the danger of using an incorrect list, and there may not be enough database information to form a list externally in some applications.

Results

We investigate the possibility of creating an internal list from the data at hand. One way to form such a list is to collapse the observed total minor allele frequencies to “zero” or “at least one”, which is shown to have the desirable effect of amplifying the haplotype frequencies. To improve coverage, we propose ways to add and remove haplotypes from the list, and a benchmarking method to determine the frequency threshold for removing haplotypes. Simulation results show that the EM estimates based on a suitably augmented and trimmed collapsed data list (ATCDL) perform satisfactorily. In two scenarios involving 25 and 32 loci respectively, the EM-ATCDL estimates outperform the EM estimates based on other lists as well as the collapsed data maximum likelihood estimates.

Conclusions

The proposed augmented and trimmed CD list is a useful list for the EM algorithm to base upon in estimating the haplotype distributions of rare variants. It can handle more markers and larger pool size than existing methods, and the resulting EM-ATCDL estimates are more efficient than the EM estimates based on other lists.