Haplotype-association analysis

Molecular Characteristics, Functional Definitions, and Regulatory Mechanisms for Cross-Presentation Mediated by the Major Histocompatibility Complex: A Comprehensive Review

The major histocompatibility complexes of vertebrates play a key role in the immune response. Antigen-presenting cells are loaded on MHC I molecules, which mainly present endogenous antigens; when MHC I presents exogenous antigens, this is called cross-presentation. The discovery of cross-presentation provides an important theoretical basis for the study of exogenous antigens. Cross-presentation is a complex process in which MHC I molecules present antigens to the cell surface to activate CD8+ T lymphocytes. The process of cross-representation includes many components, and this article briefly outlines the origins and development of MHC molecules, gene structures, functions, and their classical presentation pathways. The cross-presentation pathways of MHC I molecules, the cell lines that support cross-presentation, and the mechanisms of MHC I molecular transporting are all reviewed. After more than 40 years of research, the specific mechanism of cross-presentation is still unclear. In this paper, we summarize cross-presentation and anticipate the research and development prospects for cross-presentation.

Polymorphisms in the promoter regions of RHD and RHCE genes in the Chinese Han population

BACKGROUND AND OBJECTIVES

The Rh blood group system is the most polymorphic human blood group system. Previous studies have investigated variants in the RHD and RHCE promoter. The relevance of these variants to the Chinese Han population is further clarified in this study.

MATERIALS AND METHODS

In total, 317 donors (223 Rh D-positive [D+], including 20 Del and 94 Rh D-negative [D-]) were randomly selected. The promoter regions and exon 1 of RHD and RHCE were amplified through polymerase chain reaction (PCR) whose products were directly sequenced using forward and reverse primers.

RESULTS

Expected PCR products of the RHD promoter and exon 1 were amplified in 223 D+ individuals, including 20 Del individuals, and were absent in 81 of 94 D- individuals. Expected PCR products of RHCE were observed in all donors. Two single nucleotide variants (SNVs) were observed in the RHD promoter region. Moreover, 11 SNVs were observed in the promoter and exon 1 of RHCE. rs4649082, rs2375313, rs2281179, rs2072933, rs2072932, rs2072931 and rs586178 with strong linkage disequilibria were significantly different between the D+ and D- groups. [A;C] was the most common haplotype in the RHD promoter (NC_000001.11:g.[-1033A>G;-831C>T]). [G;T;T;A;T;A;C;G;A;C;G] was the most predominant haplotype in both total and D- groups. In D+ individuals, [A;C;T;G;C;G;C;G;C;C;C] was the most frequent haplotype in the RHCE promoter (NC_000001.11:g.[-1080A>G;-958C>T;-390T>C;-378G>A;-369C>T;-296G>A;-144C>G;-132G>A;-122C>A;28C>T;48C>G]).

CONCLUSION

We speculate that the SNVs/haplotypes found in this article cannot significantly affect gene expression. The present study findings should help elucidate the molecular basis of the polymorphic expression of RHD and RHCE promoter regions.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains

BACKGROUND

The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD.

METHODS

Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs).

RESULTS

Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes.

CONCLUSIONS

To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Identification of common genetic polymorphisms associated with down-regulated gonadotropin levels in an exome-wide association study

OBJECTIVE

To investigate whether common genetic polymorphisms are associated with gonadotropin levels after down-regulation with daily gonadotropin-releasing hormone agonist and whether the polymorphisms of candidate variants influence the ovarian response to exogenous gonadotropins.

DESIGN

Genetic association study.

SETTING

University-affiliated in vitro fertilization center.

PATIENTS

Subjects enrolled in an exploratory exome-wide association study (n = 862), a replication exome-wide association study (n = 86), and a classifier validation study (n = 148) were recruited from September 2016 to October 2018, September 2019 to September 2020, and January 2021 to December 2021, respectively. The included patients were aged ≤40 years and had a basal follicle-stimulating hormone (FSH) ≤12 IU/L.

INTERVENTIONS

All participants received a luteal phase down-regulation long protocol. Genome DNA was extracted from the peripheral blood leukocytes. For the exploratory and replication cohorts, exome sequencing was conducted on a HiSeq 2500 sequencing platform. The multiplex polymerase chain reaction amplification technique and next-generation sequencing also were performed in the exploratory and replication cohorts. For the samples of the validation cohort, Sanger sequencing was performed.

MAIN OUTCOME MEASURES

The primary endpoint was the gonadotropin levels after down-regulation, and the secondary endpoints were hormone levels and follicle diameters during stimulation, the total dose of FSH, duration of FSH stimulation, number of oocytes retrieved, and clinical pregnancy rate.

RESULTS

In the exploratory cohort, we identified that FSHB rs6169 (P=2.71 × 10-24) and its single-nucleotide polymorphisms in high linkage disequilibrium were associated with the down-regulated FSH level. The same locus was confirmed in the replication cohort. Women carrying the C allele of FSHB rs6169 exhibited higher average estradiol level during stimulation (P=6.82 × 10-5), shorter duration of stimulation, and less amount of exogenous FSH (Pduration=0.0002; Pdose=0.0024). In the independent validation set, adding rs6169 genotypes into the prediction model for FSH level after down-regulation enhanced the area under the curve from 0.560 to 0.712 in a logistic regression model, and increased prediction accuracy by 41.05% when a support vector machine classifier was applied.

CONCLUSION

The C allele of FSHB rs6169 is a susceptibility site for the relatively high level of FSH after down-regulation, which may be associated with increased ovarian FSH sensitivity.

TP53 and p21 (CDKN1A) polymorphisms and the risk of systemic lupus erythematosus

BACKGROUND

The p53 and p21 proteins are important regulators of cell cycle and apoptosis and may contribute to autoimmune diseases, such as systemic lupus erythematosus (SLE). As genetic polymorphisms may cause changes in protein levels and functions, we investigated associations of TP53 and p21 (CDKN1A) polymorphisms (p53 72 G > C-rs1042522; p53 PIN3-rs17878362; p21 31 C > A-rs1801270; p21 70 C > T-rs1059234) with the development of systemic lupus erythematosus (SLE) in a Southeastern Brazilian population.

METHODS

Genotyping of 353 female volunteers (cases, n = 145; controls, n = 208) was performed by polymerase chain reaction, restriction fragment length polymorphism and/or DNA sequencing. Associations between TP53 and p21 polymorphisms and SLE susceptibility and clinical manifestations of SLE patients were assessed by logistic regression analysis.

RESULTS

Protective effect was observed for the genotype combinations p53 PIN3 A1/A1-p21 31 C/A, in the total study population (OR 0.45), and p53 PIN3 A1/A2-p21 31 C/C, in non-white women (OR 0.28). In Whites, p53 72 C-containing (OR 3.06) and p53 PIN3 A2-containing (OR 6.93) genotypes were associated with SLE risk, and higher OR value was observed for the combined genotype p53 72 G/C-p53 PIN3 A1/A2 (OR 9.00). Further, p53 PIN3 A1/A2 genotype was associated with serositis (OR 2.82), while p53 PIN3 A2/A2 and p53 72 C/C genotypes were associated with neurological disorders (OR 4.69 and OR 3.34, respectively).

CONCLUSIONS

Our findings showed that the TP53 and p21 polymorphisms included in this study may have potential to emerge as SLE susceptibility markers for specific groups of patients. Significant interactions of the TP53 polymorphisms with serositis and neurological disorders were also observed in SLE patients.

Single-nucleotide polymorphisms in ABC drug transporters alter expression and circulating tenofovir in healthy South African women exposed to pre-exposure prophylaxis

Aim: We investigated if single-nucleotide polymorphisms (SNPs) in ATP-binding cassette (ABC) drug transporters alter gene expression and tenofovir disposition in South African women taking Truvada® for HIV prevention. Materials & methods: In 393 women, real-time PCR was used to determine the associations between six SNPs in ABC transporter genes, mRNA expression and circulating-tenofovir. Results: Univariable and multivariable analyses showed that CT and TT relative to CC genotypes for the ABCC4(3463C/T) SNP had significantly higher tenofovir levels. In contrast, the AA genotype for the ABCC4(4976A/G) SNP showed significantly less tenofovir, while mRNA expression was increased. Conclusion: SNPs in the ABCC4 gene may differentially affect gene expression and circulating tenofovir. Their impact may inform on low pre-exposure prophylaxis efficacy and discern effective drugs in clinical trials of African women enriched for certain genotypes.

Exploring the association between fat-related traits in chickens and the RGS16 gene: insights from polymorphism and functional validation analysis

Introduction

Excessive fat deposition in chickens can lead to reduced feed utilization and meat quality, resulting in significant economic losses for the broiler industry. Therefore, reducing fat deposition has become an important breeding objective in addition to achieving high broiler weight, growth rate, and feed conversion efficiency. In our previous studies, we observed high expression of Regulators of G Protein Signaling 16 Gene (RGS16) in high-fat individuals. This led us to speculate that RGS16 might be involved in the process of fat deposition in chickens.

Methods

Thus, we conducted a polymorphism and functional analysis of the RGS16 gene to investigate its association with fat-related phenotypic traits in chickens. Using a mixed linear model (MLM), this study explored the relationship between RGS16 gene polymorphisms and fat-related traits for the first time. We identified 30 SNPs of RGS16 in a population of Wens Sanhuang chickens, among which 8 SNPs were significantly associated with fat-related traits, including sebum thickness (ST), abdominal fat weight (AFW), and abdominal fat weight (AFR). Furthermore, our findings demonstrated that AFW, AFR, and ST showed significant associations with at least two or more out of the eight identified SNPs of RGS16. We also validated the role of RGS16 in ICP-1 cells through various experimental methods, including RT-qPCR, CCK- 8, EdU assays, and oil red O staining.

Results

Our functional validation experiments showed that RGS16 was highly expressed in the abdominal adipose tissue of high-fat chickens and played a critical role in the regulation of fat deposition by promoting preadipocyte differentiation and inhibiting their proliferation. Taken together, our findings suggest that RGS16 polymorphisms are associated with fat-related traits in chickens. Moreover, the ectopic expression of RGS16 could inhibit preadipocyte proliferation but promote preadipocyte differentiation.

Discussion

Based on our current findings, we propose that the RGS16 gene could serve as a powerful genetic marker for marker-assisted breeding of chicken fat-related traits.

Comprehensive investigating of mismatch repair genes (MMR) polymorphisms in participants with chronic hepatitis B virus infection

Background and aim: In this study, we focused on the relationship between single nucleotide polymorphisms in MMR genes and the occurrence and development of HBV infection. Materials and methods: A total of 3,128 participants were divided into five groups: negative control group (NeC), spontaneous clearance group (SC), chronic hepatitis B group (CHB), liver cirrhosis group (LC) and hepatocellular carcinoma group (HCC), CHB, liver cirrhosis and hepatocellular carcinoma constitute HLD. We conducted three case-control studies: NeC (840 cases) vs. HLD (1792 cases), SC (486 cases) vs. HLD (1792 cases) and CHB + LC (1,371 cases) vs. HCC (421 cases). 11 polymorphic loci in MLH1, MLH3, MSH5, PMS1 and PMS2 were involved in genotyping by Sequenom MassArray. The SNPStats performed Hardy-Weinberg equilibrium test. Linkage disequilibrium patterns were visualized using Haploview4.2. The GMDR (v0.9) was conducted to generalized multifactor dimension reduction analysis. The correlation, multiplicative interaction and additive interaction analyses were calculated by Logistic Regression through SPSS21.0. Matrix and programmed excel were also involved in the calculation of additive interaction. Results: In NeC vs. HLD group, MSH5-rs1150793(G) was a risk base to HBV susceptibility (nominal p = 0.002, OR = 1.346). We found multiplicative interaction between MLH1-rs1540354 (AA + AT) and PMS1-rs1233255 (AA) (nominal p = 0.024, OR = 1.240). There was additive interaction between PMS1-rs1233255 (AA) and PMS1-rs256554(CA + CC). In SC vs. HLD group, MLH1-rs1540354 (TT) was a risk genotype (nominal p < 0.05, OR>1). Through haplotype analysis, we found the linkage disequilibrium of three loci in MLH1. The results of GMDR showed the optimal five-locus model about the spontaneous clearance of HBV. In CHB + LC vs. HCC group, PMS2-rs12112229(A) was related to the cancerization of liver. Conclusion: We found rs1150793(G), rs1540354(T) and rs12112229(A) were significantly related to HBV susceptibility, spontaneous clearance of HBV and cancerization after infection, respectively.

Pharmacogenomic and Statistical Analysis

Genetic variants can alter response to drugs and other therapeutic interventions. The study of this phenomenon, called pharmacogenomics, is similar in many ways to other types of genetic studies but has distinct methodological and statistical considerations. Genetic variants involved in the processing of exogenous compounds exhibit great diversity and complexity, and the phenotypes studied in pharmacogenomics are also more complex than typical genetic studies. In this chapter, we review basic concepts in pharmacogenomic study designs, data generation techniques, statistical analysis approaches, and commonly used methods and briefly discuss the ultimate translation of findings to clinical care.

Developing best practices for genotyping-by-sequencing analysis in the construction of linkage maps

BACKGROUND

Genotyping-by-sequencing (GBS) provides affordable methods for genotyping hundreds of individuals using millions of markers. However, this challenges bioinformatic procedures that must overcome possible artifacts such as the bias generated by polymerase chain reaction duplicates and sequencing errors. Genotyping errors lead to data that deviate from what is expected from regular meiosis. This, in turn, leads to difficulties in grouping and ordering markers, resulting in inflated and incorrect linkage maps. Therefore, genotyping errors can be easily detected by linkage map quality evaluations.

RESULTS

We developed and used the Reads2Map workflow to build linkage maps with simulated and empirical GBS data of diploid outcrossing populations. The workflows run GATK, Stacks, TASSEL, and Freebayes for single-nucleotide polymorphism calling and updog, polyRAD, and SuperMASSA for genotype calling, as well as OneMap and GUSMap to build linkage maps. Using simulated data, we observed which genotype call software fails in identifying common errors in GBS sequencing data and proposed specific filters to better handle them. We tested whether it is possible to overcome errors in a linkage map using genotype probabilities from each software or global error rates to estimate genetic distances with an updated version of OneMap. We also evaluated the impact of segregation distortion, contaminant samples, and haplotype-based multiallelic markers in the final linkage maps. Through our evaluations, we observed that some of the approaches produce different results depending on the dataset (dataset dependent) and others produce consistent advantageous results among them (dataset independent).

CONCLUSIONS

We set as default in the Reads2Map workflows the approaches that showed to be dataset independent for GBS datasets according to our results. This reduces the number of required tests to identify optimal pipelines and parameters for other empirical datasets. Using Reads2Map, users can select the pipeline and parameters that best fit their data context. The Reads2MapApp shiny app provides a graphical representation of the results to facilitate their interpretation.

Interleukin-38 promoter variants and risk of COVID-19 among Iraqis

Coronavirus disease-19 (COVID-19) has recently emerged as a respiratory infection with a significant impact on health and society. The pathogenesis is primarily attributed to a dysregulation of cytokines, especially those with pro-inflammatory and anti-inflammatory effects. Interleukin-38 (IL-38) is a recently identified anti-inflammatory cytokine with a proposed involvement in mediating COVID-19 pathogenesis, while the association between IL38 gene variants and disease susceptibility has not been explored. Therefore, a pilot study was designed to evaluate the association of three gene variants in the promoter region of IL38 gene (rs7599662 T/A/C/G, rs28992497 T/C and rs28992498 C/A/T) with COVID-19 risk. DNA sequencing was performed to identify these variants. The study included 148 Iraqi patients with COVID-19 and 113 healthy controls (HC). Only rs7599662 showed a significant negative association with susceptibility to COVID-19. The mutant T allele was presented at a significantly lower frequency in patients compared to HC. Analysis of recessive, dominant and codominant models demonstrated that rs7599662 TT genotype frequency was significantly lower in patients than in HC. In terms of haplotypes (in order: rs7599662, rs28992497 and rs28992498), frequency of CTC haplotype was significantly increased in patients compared to HC, while TTC haplotype showed significantly lower frequency in patients. The three SNPs influenced serum IL-38 levels and homozygous genotypes of mutant alleles were associated with elevated levels. In conclusion, this study indicated that IL38 gene in terms of promoter variants and haplotypes may have important implications for COVID-19 risk.

The effect of polymorphisms in FSHR gene on late follicular phase progesterone and estradiol serum levels in predicted normoresponders

STUDY QUESTION

Does the presence of FSHR single-nucleotide polymorphisms (SNPs) affect late follicular phase progesterone and estradiol serum levels in predicted normoresponders treated with rFSH?

SUMMARY ANSWER

The presence of FSHR SNPs (rs6165, rs6166, rs1394205) had no clinically significant impact on late follicular phase serum progesterone and estradiol levels in predicted normoresponders undergoing a GnRH antagonist protocol with a fixed daily dose of 150 IU rFSH.

WHAT IS KNOWN ALREADY

Previous studies have shown that late follicular phase serum progesterone and estradiol levels are significantly correlated with the magnitude of ovarian response. Several authors have proposed that individual variability in the response to ovarian stimulation (OS) could be explained by variants in FSHR. However, so far, the literature is scarce on the influence of this genetic variability on late follicular phase steroidogenic response. Our aim is to determine whether genetic variants in the FSHR gene could modulate late follicular phase serum progesterone and estradiol levels.

STUDY DESIGN, SIZE, DURATION

In this multicenter multinational prospective study conducted from November 2016 to June 2019, 366 patients from Vietnam, Belgium and Spain (166 from Europe and 200 from Asia) underwent OS followed by oocyte retrieval in a GnRH antagonist protocol with a fixed daily dose of 150 IU rFSH. All patients were genotyped for 3 FSHR SNPs (rs6165, rs6166, rs1394205) and had a serum progesterone and estradiol measurement on the day of trigger.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Included patients were predicted normal responder women <38 years old undergoing their first or second OS cycle. The prevalence of late follicular phase progesterone elevation (PE), as well as mean serum progesterone and estradiol levels on the day of trigger were compared between the different FSHR SNPs genotypes. PE was defined as >1.50 ng/ml.

MAIN RESULTS AND THE ROLE OF CHANCE

The overall prevalence of PE was 15.8% (n = 58). No significant difference was found in the prevalence of PE in Caucasian and Asian patients (17.5% versus 14.5%). Estradiol levels on the day of trigger and the number of retrieved oocytes were significantly higher in patients with PE (4779 ± 6236.2 versus 3261 ± 3974.5 pg/ml, P = 0.003, and 16.1 ± 8.02 versus 13.5 ± 6.66, P = 0.011, respectively). Genetic model analysis, adjusted for patient age, body mass index, number of retrieved oocytes and continent (Asia versus Europe), revealed a similar prevalence of PE in co-dominant, dominant and recessive models for variants FSHR rs6166, rs6165 and rs1394205. No statistically significant difference was observed in the mean late follicular phase progesterone serum levels according to the genotypes of FSHR rs6166 (P = 0.941), rs6165 (P = 0.637) and rs1394205 (P = 0.114) in the bivariate analysis. Also, no difference was found in the genetic model analysis regarding mean late follicular phase progesterone levels across the different genotypes. Genetic model analysis has also revealed no statistically significant difference regarding mean estradiol levels on the day of trigger in co-dominant, dominant and recessive models for variants FSHR rs6166, rs6165 and rs1394205. Haplotype analysis revealed a statistically significant lower estradiol level on the day of trigger for rs6166/rs6165 haplotypes GA, AA and GG when compared to AG (respectively, estimated mean difference (EMD) -441.46 pg/ml (95% CI -442.47; -440.45), EMD -673.46 pg/ml (95% CI -674.26; -672.67) and EMD -582.10 pg/ml (95% CI -584.92; -579.28)). No statistically significant differences were found regarding the prevalence of PE nor late follicular phase progesterone levels according to rs6166/rs6165 haplotypes.

LIMITATIONS, REASONS FOR CAUTION

Results refer to a population of predicted normal responders treated with a normal/low fixed dose of 150 IU rFSH throughout the whole OS. Consequently, caution is needed before generalizing our results to all patient categories.

WIDER IMPLICATIONS OF THE FINDINGS

Based on our results, FSHR SNPs rs6165, rs6166 and rs1394205 do not have any clinically significant impact neither on late follicular phase serum progesterone nor on estradiol levels in predicted normal responders. These findings add to the controversy in the literature regarding the impact of individual genetic susceptibility in response to OS in this population.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by an unrestricted grant by Merck Sharp & Dohme (MSD, IISP56222). N.P.P. reports grants and/or personal fees from MSD, Merck Serono, Roche Diagnostics, Ferring International, Besins Healthcare, Gedeon Richter, Organon, Theramex and Institut Biochimique SA (IBSA). C.A. reports conference fees from Merck Serono, Medea and Event Planet. A.R.N., C.B., C.S., P.Q.M.M., H.T., C.B., N.L.V., M.T.H. and S.G. report no conflict of interests related to the content of this article.

TRIAL REGISTRATION NUMBER

NCT03007043.

The Influence of FTO, FABP2, LEP, LEPR, and MC4R Genes on Obesity Parameters in Physically Active Caucasian Men

Obesity is a complex multifactorial abnormality that has a well-confirmed genetic basis. However, the problem still lies in identifying the polymorphisms linked to body mass and composition. Therefore, this study aimed to analyze associations between FTO (rs9939609), FABP2 (rs1799883), and LEP (rs2167270), LEPR (rs1137101), and MC4R (rs17782313) polymorphisms and obesity-related parameters. Unrelated Caucasian males (n = 165) were recruited. All participants had similar physical activity levels. The participants were divided into two groups depending on their body mass index (BMI) and fat mass index (FMI). All samples were genotyped using real-time polymerase chain reaction (real-time PCR). When tested individually, only one statistically significant result was found. The FTO A/T polymorphism was significantly associated with FMI (p = 0.01). The chance of having increased FMI was >2-fold higher for the FTO A allele carriers (p < 0.01). Gene−gene interaction analyses showed the additional influence of all investigated genes on BMI and FMI. In summary, it was demonstrated that harboring the FTO A allele might be a risk factor for elevated fat mass. Additionally, this study confirmed that all five polymorphisms are involved in the development of common obesity in the studied population and the genetic risk of obesity is linked to the accumulation of numerous variants.

Impact of the DRD2 Polymorphisms on the Effectiveness of the Training Program

Dopamine receptor D2 gene (DRD2) polymorphisms have been associated with cognitive abilities, obesity, addictions, and physical-activity-related behaviors, which may underlie differences in the effectiveness of training programs. What is not yet clear is the impact of DRD2 polymorphisms on the effectiveness of exercise programs. Thus, the aim of this study was to investigate the association between the DRD2 polymorphic sites (rs1076560, rs12364283, rs1799732, rs1800497, and rs1800498) and the body's response to regular physical activity. We studied genotypes and haplotypes distribution in a group of 165 females measured for body mass and body composition measurements, lipid profile, and glucose levels before and after realization of a 12-week training program. When tested individually, statistical analyses revealed one significant genotype by training interaction under the general model (for the basal metabolic rate, BMR, p = 0.033). Carriers of the rs1076560 CC genotype exhibited a decrease in BMR in response to training (p = 0.006). Haplotype analyses also showed that (i) the CACCC and CACTT haplotypes were associated with a post-training decrease in glucose level (β = -4.11, p = 0.032; β = -6.86, p = 0.020, respectively); (ii) the CGCCT with an increase in BMR (β = 0.65, p = 0.003) and fat free mass (FFM, β = 1.20, p = 0.009); (iii) the CA-CT with a decrease in low-density lipoprotein cholesterol (LDL, β = -17.26, p = 0.046). These results provide some evidence that the DRD2 polymorphisms may play a role in post-training changes in lipid and carbohydrate metabolism, and, as a consequence, in the effectiveness of training programs.

circTAF8 Regulates Myoblast Development and Associated Carcass Traits in Chicken

Recent studies have shown that circular RNAs (circRNAs) play important roles in skeletal muscle development. CircRNA biogenesis is dependent on the genetic context. Single-nucleotide polymorphisms in the introns flanking circRNAs may be intermediate-inducible factors between circRNA expression and phenotypic traits. Our previous study showed that circTAF8 is an abundantly and differentially expressed circRNA in leg muscle during chicken embryonic development. Here, we aimed to investigate circTAF8 function in muscle development and the association of the SNPs in the circTAF8 flanking introns with carcass traits. In this study, we observed that overexpression of circTAF8 could promote the proliferation of chicken primary myoblasts and inhibit their differentiation. In addition, the SNPs in the introns flanking the circTAF8 locus and those associated with chicken carcass traits were analyzed in 335 partridge chickens. A total of eight SNPs were found associated with carcass traits such as leg muscle weight, live weight, and half and full-bore weight. The association analysis results of haplotype combinations were consistent with the association analysis of a single SNP. These results suggest that circTAF8 plays a regulatory role in muscle development. These identified SNPs were found correlated with traits to muscle development and carcass muscle weight in chickens.

Association of SMAD7 genetic markers and haplotypes with colorectal cancer risk

Purpose

Colorectal cancer (CRC) is one of the common cancers with a high mortality rate worldwide. In Iran, there has been a trend of increased incidence of colorectal cancer in the last three decades that necessitates the early diagnosis. Genetic factors have an influential role in its etiology along with the conventional risk factors such as age, diet, and lifestyle. Results from GWAS have shown significant associations between SMAD7 gene variants and risk of CRC. This study aimed to assess the association of certain polymorphisms as well as haplotypes of this gene and risk of colorectal cancer.

Methods and materials

This study was designed as a case–control association study. After obtaining ethical approval and informed consent, blood samples from 209 patients with colorectal cancer were collected and DNA was extracted. Four variants: rs4939827, rs34007497, rs8085824 and rs8088297 were genotyped using ARMS-PCR method.

Results

SMAD7 rs4939827 in the recessive and co-dominant models was associated with colorectal cancer risk [TT/CT + CC: OR = 2.90, 95%CI (1.38–6.09), p = 0.005; CC + TT/CT: OR = 1.66, 95%CI (1.00–2.75), p = 0.01]. Haplotype analysis indicated that some SNP combinations including two for-SNPs haplotypes of T-T-C-C and T-C-C-A were significantly associated with CRC risk.

Conclusion

Based on the identified association of SMAD7 gene variations and haplotypes with colorectal cancer risk in our population, genetic variations in this gene region may have a role in CRC development. This data may shed light on the genetic predisposition of CRC which involves different pathways including TGF-β.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01150-3.

Is BRD7 associated with spermatogenesis impairment and male infertility in humans? A case-control study in a Han Chinese population

Background

Bromodomain-containing protein 7 (BRD7), a member of the bromodomain-containing protein family, plays important roles in chromatin modification and transcriptional regulation. A recent model of Brd7-knockout mice presented azoospermia and male infertility, implying the potential role of BRD7 in spermatogenic failure in humans. This case-control study aimed to explore the association of the BRD7 gene with spermatogenic efficiency and the risk of spermatogenic defects in humans.

Results

A total of six heterozygous variants were detected in the coding and splicing regions of the BRD7 gene in patients with azoospermia. For each of four rare variants predicted to potentially damage BRD7 function, we further identified these four variants in oligozoospermia and normozoospermia as well. However, no difference in the allele and genotype frequencies of rare variants were observed between cases with spermatogenic failure and controls with normozoospermia; the sperm products of variant carriers were similar to those of noncarriers. Moreover, similar distribution of the alleles, genotypes and haplotypes of seven tag single nucleotide polymorphisms (tagSNPs) was observed between the cases with azoospermia and oligozoospermia and controls with normozoospermia; associations of tagSNP-distinguished BRD7 alleles with sperm products were not identified.

Conclusions

The lack of an association of BRD7-linked rare and common variants with spermatogenic failure implied a limited contribution of the BRD7 gene to spermatogenic efficiency and susceptibility to male infertility in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12610-021-00139-3.

Association of eNOS and MCP-1 Genetic Variants with Type 2 Diabetes and Diabetic Nephropathy Susceptibility: A Case-Control and Meta-Analysis Study

Type 2 diabetes (T2D) and its secondary complications result from the complex interplay of genetic and environmental factors. To understand the role of these factors on disease susceptibility, the present study was conducted to assess the association of eNOS and MCP-1 variants with T2D and diabetic nephropathy (DN) in two ethnically and geographically different cohorts from North India. A total of 1313 subjects from two cohorts were genotyped for eNOS (rs2070744, rs869109213 and rs1799983) and MCP-1 (rs1024611 and rs3917887) variants. Cohort-I (Punjab) comprised 461 T2D cases (204 T2D with DN and 257 T2D without DN) and 315 healthy controls. Cohort-II (Jammu and Kashmir) included 337 T2D (150 T2D with DN and 187 T2D without DN) and 200 controls. Allele, genotype and haplotype frequencies were compared among the studied participants, and phenotype-genotype interactions were determined. Meta-analysis was performed to investigate the association between the selected variants and disease susceptibility. All three eNOS variants were associated with 1.5-4.0-fold risk of DN in both cohorts. MCP-1 rs1024611 conferred twofold risk towards DN progression in cohort-II, while rs3917887 provided twofold risk for both T2D and DN in both cohorts. eNOS and MCP-1 haplotypes conferred risk for T2D and DN susceptibility. Phenotype-genotype interactions showed significant associations between the studied variants and anthropometric and biochemical parameters. In meta-analysis, all eNOS variants conferred risk towards DN progression, whereas no significant association was observed for MCP-1 rs1024611. We show evidences for an association of eNOS and MCP-1 variants with T2D and DN susceptibility.

Genetic variants in S-adenosyl-methionine synthesis pathway and nonsyndromic cleft lip with or without cleft palate in Chile

BACKGROUND

The S-adenosyl-methionine (SAM) availability is crucial for DNA methylation, an epigenetic mechanism involved in nonsyndromic cleft lip with or without cleft palate (NSCL/P) expression. The aim of this study was to assess the association between single-nucleotide polymorphisms (SNPs) of genes involved in SAM synthesis and NSCL/P in a Chilean population.

METHODS

In 234 cases and 309 controls, 18 SNPs in AHCY, MTR, MTRR, and MAT2A were genotyped, and the association between them and the phenotype was evaluated based on additive (allele), dominant, recessive and haplotype models, by odds ratio (OR) computing.

RESULTS

Three deep intronic SNPs of MTR showed a protective effect on NSCL/P expression: rs10925239 (OR 0.68; p = 0.0032; q = 0.0192), rs10925254 (OR 0.66; p = 0.0018; q = 0.0162), and rs3768142 (OR 0.66; p = 0.0015; q = 0.0162). Annotations in expression database demonstrate that the protective allele of the three SNPs is associated with a reduction of MTR expression summed to the prediction by bioinformatic tools of its potentiality to modify splicing sites.

CONCLUSIONS

The protective effect against NSCL/P of these intronic MTR SNPs seems to be related to a decrease in MTR enzyme expression, modulating the SAM availability for proper substrate methylation. However, functional analyses are necessary to confirm our findings.

IMPACT

SAM synthesis pathway genetic variants are factors associated to NSCL/P. This article adds new evidence for folate related genes in NSCL/P in Chile. Its impact is to contribute with potential new markers for genetic counseling.

Polymorphisms and haplotypes in the Interleukin 17 Alfa gene: potential effect of smoking habits in the association with periodontitis and type 2 diabetes mellitus

Few studies evaluate interrelationships between periodontitis (P) and Type 2 Diabetes Mellitus (T2DM). The aim of this study is to investigate the genetic susceptibility to periodontitis alone, or concomitant with T2DM (as comorbidities), analyzing single nucleotide polymorphisms (SNPs) in the Interleukin 17 alpha (IL17A) gene, considering the biochemical profile and smoking habits on the subjects' periodontal status. We investigated 879 individuals divided into: T2DM subjects also affected by severe or moderate periodontitis (T2DM-P, n = 199); non-diabetics with severe or moderate periodontitis (PERIODONTITIS, n = 342); and healthy subjects (HEALTHY, n = 338). Subjects underwent complete periodontal examination, history of smoking habits, glycemic and lipid biochemical evaluation. DNA from buccal cells was utilized to genotype the SNPs rs2275913, rs3819024 and rs10484879. The impact of the subjects' biochemical profile was analyzed in their periodontal status. Each SNP was analyzed independently, and as haplotypes, by multiple logistic regressions, adjusted for covariates, and also stratifying the groups by age, sex and smoking habits. Independently of the periodontitis degree, poorly-controlled T2DM subjects showed worse glycemic and lipid profile. Multiple logistic regressions demonstrated that smokers and former-smokers carrying the GG genotype of rs3819024 seemed to have higher risk for T2DM-Periodontitis (OR = 6.33; 95% CI = 1.26-31.77, p = 0.02), and mainly for T2DM alone (OR = 5.11; 95% CI = 1.37-19.06, p = 0.01), than never smokers. We found the potential effect of smoking habits in the association of IL17A-rs3819024-GG with diseased phenotypes. Because the observed wide confidence intervals, further studies enrolling larger populations, and SNPs' functional evaluations are needed to better understand our findings.

APOE: The New Frontier in the Development of a Therapeutic Target towards Precision Medicine in Late-Onset Alzheimer's

Alzheimer's disease (AD) has a critical unmet medical need. The consensus around the amyloid cascade hypothesis has been guiding pre-clinical and clinical research to focus mainly on targeting beta-amyloid for treating AD. Nevertheless, the vast majority of the clinical trials have repeatedly failed, prompting the urgent need to refocus on other targets and shifting the paradigm of AD drug development towards precision medicine. One such emerging target is apolipoprotein E (APOE), identified nearly 30 years ago as one of the strongest and most reproduceable genetic risk factor for late-onset Alzheimer's disease (LOAD). An exploration of APOE as a new therapeutic culprit has produced some very encouraging results, proving that the protein holds promise in the context of LOAD therapies. Here, we review the strategies to target APOE based on state-of-the-art technologies such as antisense oligonucleotides, monoclonal antibodies, and gene/base editing. We discuss the potential of these initiatives in advancing the development of novel precision medicine therapies to LOAD.

Association of Vitamin D Metabolism Gene Polymorphisms with Autoimmunity: Evidence in Population Genetic Studies

A high prevalence of vitamin D (calcidiol) serum deficiency has been described in several autoimmune diseases, including multiple sclerosis (MS), rheumatoid arthritis (AR), and systemic lupus erythematosus (SLE). Vitamin D is a potent immunonutrient that through its main metabolite calcitriol, regulates the immunomodulation of macrophages, dendritic cells, T and B lymphocytes, which express the vitamin D receptor (VDR), and they produce and respond to calcitriol. Genetic association studies have shown that up to 65% of vitamin D serum variance may be explained due to genetic background. The 90% of genetic variability takes place in the form of single nucleotide polymorphisms (SNPs), and SNPs in genes related to vitamin D metabolism have been linked to influence the calcidiol serum levels, such as in the vitamin D binding protein (VDBP; rs2282679 GC), 25-hydroxylase (rs10751657 CYP2R1), 1α-hydroxylase (rs10877012, CYP27B1) and the vitamin D receptor (FokI (rs2228570), BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) VDR). Therefore, the aim of this comprehensive literature review was to discuss the current findings of functional SNPs in GC, CYP2R1, CYP27B1, and VDR associated to genetic risk, and the most common clinical features of MS, RA, and SLE.

A chaotic viewpoint-based approach to solve haplotype assembly using hypergraph model

Decreasing the cost of high-throughput DNA sequencing technologies, provides a huge amount of data that enables researchers to determine haplotypes for diploid and polyploid organisms. Although various methods have been developed to reconstruct haplotypes in diploid form, their accuracy is still a challenging task. Also, most of the current methods cannot be applied to polyploid form. In this paper, an iterative method is proposed, which employs hypergraph to reconstruct haplotype. The proposed method by utilizing chaotic viewpoint can enhance the obtained haplotypes. For this purpose, a haplotype set was randomly generated as an initial estimate, and its consistency with the input fragments was described by constructing a weighted hypergraph. Partitioning the hypergraph specifies those positions in the haplotype set that need to be corrected. This procedure is repeated until no further improvement could be achieved. Each element of the finalized haplotype set is mapped to a line by chaos game representation, and a coordinate series is defined based on the position of mapped points. Then, some positions with low qualities can be assessed by applying a local projection. Experimental results on both simulated and real datasets demonstrate that this method outperforms most other approaches, and is promising to perform the haplotype assembly.

Haplotype of RNASE 3 polymorphisms is associated with severe malaria in an Indian population

Severe malaria (SM) caused by Plasmodium falciparum (Pf) infection has been associated with life-threatening anemia, metabolic acidosis, cerebral malaria and multiorgan dysfunction. It may lead to death if not treated promptly. RNASE 3 has been linked to Pf growth inhibition and its polymorphisms found associated with SM and cerebral malaria in African populations. This study aimed to assess the association of RNASE 3 polymorphisms with SM in an Indian population. RNASE 3 gene and flanking regions were amplified followed by direct DNA sequencing in 151 Indian patients who visited Wenlock District Government Hospital, Mangalore, Karnataka, India. Allele, genotype and haplotype frequencies were compared between patients with SM (n = 47) and uncomplicated malaria (UM; n = 104). Homozygous mutant genotype was only found for rs2233860 (+ 499G > C) polymorphism (< 1% frequency). No significant genetic associations were found for RNASE 3 polymorphism genotypes and alleles in Indian SM patients using the Fisher's exact test. C-G-G haplotype of rs2233859 (− 38C > A), rs2073342 (+ 371C > G) and rs2233860 (+ 499G > C) polymorphisms was correlated significantly with SM patients (OR = 3.03; p = 0.008) after Bonferroni correction. A haplotype of RNASE 3 gene was found associated with an increased risk of SM and confirming that RNASE 3 gene plays a role in susceptibility to SM.

NCMHap: a novel method for haplotype reconstruction based on Neutrosophic c-means clustering

BACKGROUND

Single individual haplotype problem refers to reconstructing haplotypes of an individual based on several input fragments sequenced from a specified chromosome. Solving this problem is an important task in computational biology and has many applications in the pharmaceutical industry, clinical decision-making, and genetic diseases. It is known that solving the problem is NP-hard. Although several methods have been proposed to solve the problem, it is found that most of them have low performances in dealing with noisy input fragments. Therefore, proposing a method which is accurate and scalable, is a challenging task.

RESULTS

In this paper, we introduced a method, named NCMHap, which utilizes the Neutrosophic c-means (NCM) clustering algorithm. The NCM algorithm can effectively detect the noise and outliers in the input data. In addition, it can reduce their effects in the clustering process. The proposed method has been evaluated by several benchmark datasets. Comparing with existing methods indicates when NCM is tuned by suitable parameters, the results are encouraging. In particular, when the amount of noise increases, it outperforms the comparing methods.

CONCLUSION

The proposed method is validated using simulated and real datasets. The achieved results recommend the application of NCMHap on the datasets which involve the fragments with a huge amount of gaps and noise.

Beyond Single Nucleotide Polymorphisms: CYP3A5∗3∗6∗7 Composite and ABCB1 Haplotype Associations to Tacrolimus Pharmacokinetics in Black and White Renal Transplant Recipients

Interpatient variability in tacrolimus pharmacokinetics is attributed to metabolism by cytochrome P-450 3A5 (CYP3A5) isoenzymes and membrane transport by P-glycoprotein. Interpatient pharmacokinetic variability has been associated with genotypic variants for both CYP3A5 or ABCB1. Tacrolimus pharmacokinetics was investigated in 65 stable Black and Caucasian post-renal transplant patients by assessing the effects of multiple alleles in both CYP3A5 and ABCB1. A metabolic composite based upon the CYP3A5 polymorphisms: ^∗3(rs776746), ^∗6(10264272), and ^∗7(41303343), each independently responsible for loss of protein expression was used to classify patients as extensive, intermediate and poor metabolizers. In addition, the role of ABCB1 on tacrolimus pharmacokinetics was assessed using haplotype analysis encompassing the single nucleotide polymorphisms: 1236C > T (rs1128503), 2677G > T/A(rs2032582), and 3435C > T(rs1045642). Finally, a combined analysis using both CYP3A5 and ABCB1 polymorphisms was developed to assess their inter-related influence on tacrolimus pharmacokinetics. Extensive metabolizers identified as homozygous wild type at all three CYP3A5 loci were found in 7 Blacks and required twice the tacrolimus dose (5.6 ± 1.6 mg) compared to Poor metabolizers [2.5 ± 1.1 mg (P < 0.001)]; who were primarily Whites. These extensive metabolizers had 2-fold faster clearance (P < 0.001) with 50% lower AUC^∗ (P < 0.001) than Poor metabolizers. No differences in C_{12 h} were found due to therapeutic drug monitoring. The majority of blacks (81%) were classified as either Extensive or Intermediate Metabolizers requiring higher tacrolimus doses to accommodate the more rapid clearance. Blacks who were homozygous for one or more loss of function SNPS were associated with lower tacrolimus doses and slower clearance. These values are comparable to Whites, 82% of who were in the Poor metabolic composite group. The ABCB1 haplotype analysis detected significant associations of the wildtype 1236T-2677T-3435T haplotype to tacrolimus dose (P = 0.03), CL (P = 0.023), CL/LBW (P = 0.022), and AUC^∗ (P = 0.078). Finally, analysis combining CYP3A5 and ABCB1 genotypes indicated that the presence of the ABCB1 3435 T allele significantly reduced tacrolimus clearance for all three CPY3A5 metabolic composite groups. Genotypic associations of tacrolimus pharmacokinetics can be improved by using the novel composite CYP3A5^∗3^∗4^∗5 and ABCB1 haplotypes. Consideration of multiple alleles using CYP3A5 metabolic composites and drug transporter ABCB1 haplotypes provides a more comprehensive appraisal of genetic factors contributing to interpatient variability in tacrolimus pharmacokinetics among Whites and Blacks.

Association analysis in a Latin American population revealed ethnic differences in rheumatoid arthritis-associated SNPs in Caucasian and Asian populations

Large genome-wide association studies (GWAS) have increased our knowledge of the genetic risk factors of rheumatoid arthritis (RA). However, little is known about genetic susceptibility in populations with a large admixture of Amerindian ancestry. The aim of the present study was to test the generalizability of previously reported RA loci in a Latin American (LA) population with admixed ancestry. We selected 128 single nucleotide polymorphisms (SNPs) in linkage equilibrium, with high association to RA in multiple populations of non-Amerindian origin. Genotyping of 118 SNPs was performed in 313 RA patients/487 healthy control subjects by mid-density arrays of polymerase chain reaction (PCR). Some of the identified associations were validated in an additional cohort (250 cases/290 controls). One marker, the SNP rs2451258, located upstream of T Cell Activation RhoGTPase Activating Protein (TAGAP) gene, showed significant association with RA (p = 5 × 10-3), whereas 18 markers exhibited suggestive associations (p < 0.05). Haplotype testing showed association of some groups of adjacent SNPs around the signal transducer and activator of transcription 4 (STAT4) gene (p = 9.82 × 10-3 to 2.04 × 10-3) with RA. Our major finding was little replication of previously reported genetic associations with RA. These results suggest that performing GWAS and admixture mapping in LA populations has the potential to reveal novel loci associated with RA. This in turn might help to gain insight into the 'pathogenomics' of this disease and to explore trans-population differences for RA in general.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

Comparing Single-SNP, Multi-SNP, and Haplotype-Based Approaches in Association Studies for Major Traits in Barley

The multiple single nucleotide polymorphism (multi-SNP) and haplotype-based approaches that jointly consider multiple markers unveiled a larger number of associations, some of which were shared with the single-SNP approach. A larger overlap of quantitative trait loci (QTLs) between the single-SNP and haplotype-based approaches was obtained than with the multi-SNP approach. Despite a limited overlap between the QTLs detected by these approaches, each uncovered QTLs reported previously, suggesting that each approach is capable of uncovering a different subset of QTLs. We demonstrated the efficiency of an integrated genome-wide association study (GWAS) procedure, combining single-locus and multilocus approaches to improve the capacity and reliability of association analysis to detect key QTLs. The efficiency of barley breeding programs may be improved by the practical use of QTLs identified in this study. Genome-wide association studies (GWAS) have been widely used to identify quantitative trait loci (QTLs) underlying complex agronomic traits. The conventional GWAS model is based on a single-locus model, which may prove inaccurate if a trait is controlled by multiple loci, which is the case for most agronomic traits in barley (Hordeum vulgare L.). Additionally, an individual single nucleotide polymorphism (SNP) will prove incapable of capturing underlying allelic diversity. A multilocus model could potentially represent a better alternative for QTL identification. This study aimed to explore different GWAS approaches (single-SNP, multi-SNP, and haplotype-based) to establish SNP-trait associations and to potentially describe the complex genetic architecture of seven key traits in spring barley. The multi-SNP and haplotype-based approaches unveiled a larger number of significant associations, some of which were shared with the single-SNP approach. Globally, the multi-SNP approach explained more of the phenotypic variance (cumulative R² ) and provided the best fit with the genetic model [Bayesian information criterion (BIC)]. Compared with the multi-SNP approach, the single-SNP and haplotype-based approaches were relatively similar in terms of cumulative R² and BIC, with an improvement with the haplotype-based approach. Despite limited overlap between detected QTLs, each approach discovered QTLs that had been validated previously, suggesting that each approach can uncover a different subset of QTLs. An integrated GWAS procedure, considering single-locus and multilocus GWAS approaches jointly, may improve the capacity of association studies to detect key QTLs and to provide a more complete picture of the genetic architecture of complex traits in barley.

Genome-wide association studies using binned genotypes

Linear mixed models (LMM) that tests trait association one marker at a time have been the most popular methods for genome-wide association studies. However, this approach has potential pitfalls: over conservativeness after Bonferroni correction, ignorance of linkage disequilibrium (LD) between neighboring markers, and power reduction due to overfitting SNP effects. So, multiple locus models that can simultaneously estimate and test all markers in the genome are more appropriate. Based on the multiple locus models, we proposed a bin model that combines markers into bins based on their LD relationships. A bin is treated as a new synthetic marker and we detect the associations between bins and traits. Since the number of bins can be substantially smaller than the number of markers, a penalized multiple regression method can be adopted by fitting all bins to a single model. We developed an innovative method to bin the neighboring markers and used the least absolute shrinkage and selection operator (LASSO) method. We compared BIN-Lasso with SNP-Lasso and Q + K-LMM in a simulation experiment, and showed that the new method is more powerful with less Type I error than the other two methods. We also applied the bin model to a Chinese Simmental beef cattle population for bone weight association study. The new method identified more significant associations than the classical LMM. The bin model is a new dimension reduction technique that takes advantage of biological information (i.e., LD). The new method will be a significant breakthrough in associative genomics in the big data era.

CSHAP: efficient haplotype frequency estimation based on sparse representation

MOTIVATION

Estimating haplotype frequencies from genotype data plays an important role in genetic analysis. In silico methods are usually computationally involved since phase information is not available. Due to tight linkage disequilibrium and low recombination rates, the number of haplotypes observed in human populations is far less than all the possibilities. This motivates us to solve the estimation problem by maximizing the sparsity of existing haplotypes. Here, we propose a new algorithm by applying the compressive sensing (CS) theory in the field of signal processing, compressive sensing haplotype inference (CSHAP), to solve the sparse representation of haplotype frequencies based on allele frequencies and between-allele co-variances.

RESULTS

Our proposed approach can handle both individual genotype data and pooled DNA data with hundreds of loci. The CSHAP exhibits the same accuracy compared with the state-of-the-art methods, but runs several orders of magnitude faster. CSHAP can also handle with missing genotype data imputations efficiently.

AVAILABILITY AND IMPLEMENTATION

The CSHAP is implemented in R, the source code and the testing datasets are available at http://home.ustc.edu.cn/∼zhouys/CSHAP/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Association of HLA-DPA1 and HLA-DPB1 polymorphisms with spontaneous HBsAg seroclearance in Caucasians

BACKGROUND & AIMS

Acute hepatitis B virus (HBV) infections may clear spontaneously or become chronic and run through different phases. The single nucleotide polymorphisms (SNPs) rs3077, rs9277535 and rs9277534 within the human leucocyte antigen (HLA)-DP gene have been found to be associated with HBV susceptibility and persistence in Asians. However, evidence for the influence of these variants in Caucasians has been limited so far. The aim of our study was to investigate the impact of these polymorphisms on the outcome of HBV infections in a large Caucasian population.

METHODS

In this case-control study, we retrospectively analysed 1111 Caucasian individuals, including 618 with chronic HBV infections (CHB), 239 with spontaneous HBsAg seroclearance (SC) and 254 healthy controls (HC). The rs3077, rs9277535 and rs9277534 SNPs were genotyped by a polymerase chain reaction from blood samples and melting curve analysis.

RESULTS

A significant difference in the allele distributions was observed only for the rs3077 SNP between the HC and the CHB group as well as between the SC and CHB groups. The rs3077-C allele was associated with a lower probability for spontaneous HBsAg seroclearance in comparison with the rs3077-T allele (OR 0.704, 95% CI 0.509-0.974; P = 0.033). No association of the three SNPs with the stages of chronic HBV infection was found.

CONCLUSION

This is the first study demonstrating an association of the rs3077-T allele with spontaneous HBsAg seroclearance in Caucasians. Further studies are needed to elucidate the role of HLA-DP variants in disease pathogenesis and their potential role for individualized disease management.

Differential Methylation in APOE (Chr19; Exon Four; from 44,909,188 to 44,909,373/hg38) and Increased Apolipoprotein E Plasma Levels in Subjects with Mild Cognitive Impairment

BACKGROUND

Biomarkers are essential for identification of individuals at high risk of mild cognitive impairment (MCI) for potential prevention of dementia. We investigated DNA methylation in the APOE gene and apolipoprotein E (ApoE) plasma levels as MCI biomarkers in Colombian subjects with MCI and controls.

METHODS

In total, 100 participants were included (71% women; average age, 70 years; range, 43⁻91 years). MCI was diagnosed by neuropsychological testing, medical and social history, activities of daily living, cognitive symptoms and neuroimaging. Using multivariate logistic regression models adjusted by age and gender, we examined the risk association of MCI with plasma ApoE and APOE methylation.

RESULTS

MCI was diagnosed in 41 subjects (average age, 66.5 ± 9.6 years) and compared with 59 controls. Elevated plasma ApoE and APOE methylation of CpGs 165, 190, and 198 were risk factors for MCI (p < 0.05). Higher CpG-227 methylation correlated with lower risk for MCI (p = 0.002). Only CpG-227 was significantly correlated with plasma ApoE levels (correlation coefficient = -0.665; p = 0.008).

CONCLUSION

Differential APOE methylation and increased plasma ApoE levels were correlated with MCI. These epigenetic patterns require confirmation in larger samples but could potentially be used as biomarkers to identify early stages of MCI.

Haplotype and Haplotype-Environment Interaction Analysis Revealed Roles of SPRY2 for NSCL/P among Chinese Populations

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is one of common birth defects in China, with genetic and environmental components contributing to the etiology. Genome wide association studies (GWASs) have identified SPRY1 and SPRY2 to be associated with NSCL/P among Chinese populations. This study aimed to further explore potential genetic effect and gene—environment interaction among SPRY genes based on haplotype analysis, using 806 Chinese case—parent NSCL/P trios drawn from an international consortium which conducted a genome-wide association study. After the process of quality control, 190 single nucleotide polymorphisms (SNPs) of SPRY genes were included for analyses. Haplotype and haplotype—environment interaction analyses were conducted in Population-Based Association Test (PBAT) software. A 2-SNP haplotype and three 3-SNP haplotypes showed a significant association with the risk of NSCL/P after Bonferroni correction (corrected significance level = 2.6 × 10⁻⁴). Moreover, haplotype—environment interaction analysis identified these haplotypes respectively showing statistically significant interactions with maternal multivitamin supplementation or maternal environmental tobacco smoke. This study showed SPRY2 to be associated with NSCL/P among the Chinese population through not only gene effects, but also a gene—environment interaction, highlighting the importance of considering environmental exposures in the genetic etiological study of NSCL/P.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Polymorphisms in the Toll-like receptor 3 (TLR3) gene are associated with the natural course of hepatitis B virus infection in Caucasian population

Innate immunity can induce spontaneous hepatitis B surface antigen (HBsAg) seroclearance (SC) of hepatitis B virus (HBV) infection or transition towards an inactive carrier state. Toll-like receptor (TLR) 3 signalling has been linked to these processes. Alterations in the TLR3 gene might impair immune responses against HBV. In our study, we analysed the impact of the TLR3 polymorphisms rs3775291 and rs5743305 on the natural course of HBV infection. In this retrospective study, a Caucasian cohort of 621 patients with chronic HBV infection (CHB), 239 individuals with spontaneous HBsAg SC, and 254 healthy controls were enrolled. In the CHB group, 49% of patients were inactive carriers, and 17% were HBeAg-positive. The TLR3 rs3775291 A allele was associated with a reduced likelihood of spontaneous HBsAg SC and HBeAg SC, and an increased risk of developing chronic hepatitis B. In haplotype analysis, the haplotype including both risk variants rs3775291A and rs5743305A had the lowest likelihood of HBsAg SC. Further research in larger cohorts and functional analyses are needed to shed light on the impact of TLR3 signalling.

Genome-wide haplotype association analysis of primary biliary cholangitis risk in Japanese

Primary biliary cholangitis (PBC) susceptibility loci have largely been discovered through single SNP association testing. In this study, we report genic haplotype patterns associated with PBC risk genome-wide in two Japanese cohorts. Among the 74 genic PBC risk haplotype candidates we detected with a novel methodological approach in a discovery cohort of 1,937 Japanese, nearly two-thirds were replicated (49 haplotypes, Bonferroni-corrected P < 6.8 × 10^-4) in an independent Japanese cohort (N = 949). Along with corroborating known PBC-associated loci (TNFSF15, HLA-DRA), risk haplotypes may potentially model cis-interactions that regulate gene expression. For example, one replicated haplotype association (9q32-9q33.1, OR = 1.7, P = 3.0 × 10^-21) consists of intergenic SNPs outside of the human leukocyte antigen (HLA) region that overlap regulatory histone mark peaks in liver and blood cells, and are significantly associated with TNFSF8 expression in whole blood. We also replicated a novel haplotype association involving non-HLA SNPs mapped to UMAD1 (7p21.3; OR = 15.2, P = 3.9 × 10^-9) that overlap enhancer peaks in liver and memory T_h cells. Our analysis demonstrates the utility of haplotype association analyses in discovering and characterizing PBC susceptibility loci.

Gene-gene interaction for nonsyndromic cleft lip with or without cleft palate in Chilean case-parent trios

OBJECTIVE

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a birth defect for which several genes susceptibility genes been proposed. Consequently, it has been suggested that many of these genes belong to common inter-related pathways during craniofacial development gene-gene interaction. We evaluated the presence of gene-gene interaction for single nucleotide polymorphisms within interferon regulatory factor 6 (IRF6), muscle segment homeobox 1 (MSX1), bone morphogenetic protein 4 (BMP4) and transforming growth factor 3 (TGFB3) genes in NSCL/P risk in Chilean case-parent trios.

DESIGN

From previous studies, we retrieved genotypes for 13 polymorphic variants within these four genes in 152 case-parent trios. Using the trio package (R) we evaluate the gene-gen interaction in genetic markers pairs applying a 1°-of-freedom test (1df) and a confirmatory 4°-of-freedom (4df) test for epistasis followed by both a permutation test and a Benjamini-Hochberg test for multiple comparisons adjustment.

RESULTS

We found evidence of gene-gene interaction for rs6446693 (MSX1) and rs2268625 (TGFB3) (4df p = 0.024; permutation p = 0.015, Benjamini-Hochberg p = 0.001).

CONCLUSIONS

A significant gene-gene interaction was detected for rs6446693 (MSX1) and rs2268625 (TGFB3). This finding is concordant with research in animal models showing that MSX1 and TGFB3 are expressed in common molecular pathways acting in an epistatic manner during maxillofacial development.

A method combining a random forest-based technique with the modeling of linkage disequilibrium through latent variables, to run multilocus genome-wide association studies

BACKGROUND

Genome-wide association studies (GWASs) have been widely used to discover the genetic basis of complex phenotypes. However, standard single-SNP GWASs suffer from lack of power. In particular, they do not directly account for linkage disequilibrium, that is the dependences between SNPs (Single Nucleotide Polymorphisms).

RESULTS

We present the comparative study of two multilocus GWAS strategies, in the random forest-based framework. The first method, T-Trees, was designed by Botta and collaborators (Botta et al., PLoS ONE 9(4):e93379, 2014). We designed the other method, which is an innovative hybrid method combining T-Trees with the modeling of linkage disequilibrium. Linkage disequilibrium is modeled through a collection of tree-shaped Bayesian networks with latent variables, following our former works (Mourad et al., BMC Bioinformatics 12(1):16, 2011). We compared the two methods, both on simulated and real data. For dominant and additive genetic models, in either of the conditions simulated, the hybrid approach always slightly performs better than T-Trees. We assessed predictive powers through the standard ROC technique on 14 real datasets. For 10 of the 14 datasets analyzed, the already high predicted power observed for T-Trees (0.910-0.946) can still be increased by up to 0.030. We also assessed whether the distributions of SNPs' scores obtained from T-Trees and the hybrid approach differed. Finally, we thoroughly analyzed the intersections of top 100 SNPs output by any two or the three methods amongst T-Trees, the hybrid approach, and the single-SNP method.

CONCLUSIONS

The sophistication of T-Trees through finer linkage disequilibrium modeling is shown beneficial. The distributions of SNPs' scores generated by T-Trees and the hybrid approach are shown statistically different, which suggests complementary of the methods. In particular, for 12 of the 14 real datasets, the distribution tail of highest SNPs' scores shows larger values for the hybrid approach. Thus are pinpointed more interesting SNPs than by T-Trees, to be provided as a short list of prioritized SNPs, for a further analysis by biologists. Finally, among the 211 top 100 SNPs jointly detected by the single-SNP method, T-Trees and the hybrid approach over the 14 datasets, we identified 72 and 38 SNPs respectively present in the top25s and top10s for each method.

The effects of the TOMM40 poly-T alleles on Alzheimer's disease phenotypes

The TOMM40 poly-T is a polymorphism in intron 6 of the TOMM40 gene, which is adjacent to and in linkage disequilibrium with APOE. Roses et al. identified the association between the length of TOMM40 poly-T with the risk and age of onset of late-onset Alzheimer's disease (LOAD). Following the original discovery, additional studies found associations between the TOMM40 poly-T and LOAD-related phenotypes independent of APOE genotypes, while others did not replicate these associations. Furthermore, the identity of the TOMM40 poly-T risk allele has been controversial between different LOAD-related phenotypes. Here, we propose a framework to address the conflicting findings with respect to the TOMM40 poly-T allele associations with LOAD phenotypes and their functional effects. The framework is used to interpret previous studies as means to gain insights regarding the nature of the risk allele, very long versus short. We suggest that the identity of the TOMM40 poly-T risk allele depends on the phenotype being evaluated, the ages of the study subjects at the time of assessment, and the context of the APOE genotypes. In concluding remarks, we outline future studies that will inform the mechanistic interpretation of the genetic data.

An ancestral human genetic variant linked to an ancient disease: A novel association of FMO2 polymorphisms with tuberculosis (TB) in Ethiopian populations provides new insight into the differential ethno-geographic distribution of FMO2*1

The human FMO2 (flavin-containing monooxygenase 2) gene has been shown to be involved in innate immunity against microbial infections, including tuberculosis (TB), via the modulation of oxidative stress levels. It has also been found to possess a curious loss-of-function mutation (FMO2*1/FMO2*2) that demonstrates a distinctive differentiation in expression, function and ethno-geographic distribution. However, despite evidences of ethnic-specific genetic associations in the inflammatory profile of TB, no studies were done to investigate whether these patterns of variations correlate with evidences for the involvement of FMO2 in antimicrobial immune responses and ethnic differences in the distribution of FMO2 polymorphisms except for some pharmacogenetic data that suggest a potentially deleterious role for the functional variant (FMO2*1). This genetic epidemiological study was designed to investigate whether there is an association between FMO2 polymorphisms and TB, an ancient malady that remains a modern global health concern, in a sub-Saharan Africa setting where there is not only a relatively high co-prevalence of the disease and the ancestral FMO2*1 variant but also where both Mycobcaterium and Homo sapiens are considered to have originated and co-evolved. Blood samples and TB related clinical data were collected from ascertained TB cases and unrelated household controls (n = 292) from 3 different ethnic groups in Ethiopia. Latent Mtb infection was determined using Quantiferon to develop reliable TB progression phenotypes. We sequenced exonic regions of FMO2.We identified for the first time an association between FMO2 and TB both at the SNP and haplotype level. Two novel SNPs achieved a study-wide significance [chr1:171181877(A), p = 3.15E-07, OR = 4.644 and chr1:171165749(T), p = 3.32E-06, OR = 6.825] while multiple SNPs (22) showed nominal signals. The pattern of association suggested a protective effect of FMO2 against both active and latent TB with distinct genetic variants underlying the TB-progression pathway. The results were robust for population stratification. Haplotype-based tests confirmed the SNP-based results with a single haplotype bearing the ancestral-and-functional FMO2*1 "C" allele ("AGCTCTACAATCCCCTCGTTGCGC") explaining the overall association (haplotype-specific-p = 0.000103). Strikingly, not only was FMO2*1 nominally associated with reduced risk to "Active TB" (p = 0.0118, OR = 0.496) but it also does not co-segregate with the 5'-3' flanking top high-TB-risk alleles. The study provides an evidence for the existence of an evolutionary adaptation to an ancient disease based on an ancestral genetic variant acting in a haplotypic framework in Ethiopian populations.

HPA axis in major depression: cortisol, clinical symptomatology and genetic variation predict cognition

The hypothalamic-pituitary-adrenal (HPA) axis has been implicated in the pathophysiology of a variety of mood and cognitive disorders. Neuroendocrine studies have demonstrated HPA axis overactivity in major depression, a relationship of HPA axis activity to cognitive performance and a potential role of HPA axis genetic variation in cognition. The present study investigated the simultaneous roles HPA axis activity, clinical symptomatology and HPA genetic variation play in cognitive performance. Patients with major depression with psychotic major depression (PMD) and with nonpsychotic major depression (NPMD) and healthy controls (HC) were studied. All participants underwent a diagnostic interview and psychiatric ratings, a comprehensive neuropsychological battery, overnight hourly blood sampling for cortisol and genetic assessment. Cognitive performance differed as a function of depression subtype. Across all subjects, cognitive performance was negatively correlated with higher cortisol, and PMD patients had higher cortisol than did NPMDs and HCs. Cortisol, clinical symptoms and variation in genes, NR3C1 (glucocorticoid receptor; GR) and NR3C2 (mineralocorticoid receptor; MR) that encode for GRs and MRs, predicted cognitive performance. Beyond the effects of cortisol, demographics and clinical symptoms, NR3C1 variation predicted attention and working memory, whereas NR3C2 polymorphisms predicted memory performance. These findings parallel the distribution of GR and MR in primate brain and their putative roles in specific cognitive tasks. HPA axis genetic variation and activity were important predictors of cognition across the entire sample of depressed subjects and HR. GR and MR genetic variation predicted unique cognitive functions, beyond the influence of cortisol and clinical symptoms. GR genetic variation was implicated in attention and working memory, whereas MR was implicated in verbal memory.

Characterizing a region on BTA11 affecting β-lactoglobulin content of milk using high-density genotyping and haplotype grouping

Background

Milk β-lactoglobulin (β-LG) content is of interest as it is associated with nutritional and manufacturing properties. It is known that milk β-LG content is strongly affected by genetic factors. In cattle, most of the genetic differences are associated with a chromosomal region on BTA11, which contains the β-LG gene. The aim of this study was to characterize this region using 777 k SNP data (BovineHDbeadChip) and perform a haplotype-based association study. A statistical approach was developed to build haplotypes that capture the genetic variation associated with this genomic region.

Results

The SNP with the most significant effect on β-lactoglobulin content was one of the 2 causal mutations responsible for the β-lactoglobulin protein variants A/B. Haplotypes based on 2 to 5 selected lead SNP were clustered in groups with different effects on β-lactoglobulin content. Four different groups were identified suggesting that β-lactoglobulin variant A and B can be further refined in A₁, A₂, B₁ and B₂.

Conclusions

This study showed that β-lactoglobulin protein variants A/B do not explain all genetic variation associated with the tail part of BTA11 but this region contains more than one mutation with an effect on β-lactoglobulin content. These findings can be used for selection of cows with higher cheese yield, which is desirable for the dairy industry.

A Genome-Wide Association Study for Agronomic Traits in Soybean Using SNP Markers and SNP-Based Haplotype Analysis

Mapping quantitative trait loci through the use of linkage disequilibrium (LD) in populations of unrelated individuals provides a valuable approach for dissecting the genetic basis of complex traits in soybean (Glycine max). The haplotype-based genome-wide association study (GWAS) has now been proposed as a complementary approach to intensify benefits from LD, which enable to assess the genetic determinants of agronomic traits. In this study a GWAS was undertaken to identify genomic regions that control 100-seed weight (SW), plant height (PH) and seed yield (SY) in a soybean association mapping panel using single nucleotide polymorphism (SNP) markers and haplotype information. The soybean cultivars (N = 169) were field-evaluated across four locations of southern Brazil. The genome-wide haplotype association analysis (941 haplotypes) identified eleven, seventeen and fifty-nine SNP-based haplotypes significantly associated with SY, SW and PH, respectively. Although most marker-trait associations were environment and trait specific, stable haplotype associations were identified for SY and SW across environments (i.e., haplotypes Gm12_Hap12). The haplotype block 42 on Chr19 (Gm19_Hap42) was confirmed to be associated with PH in two environments. These findings enable us to refine the breeding strategy for tropical soybean, which confirm that haplotype-based GWAS can provide new insights on the genetic determinants that are not captured by the single-marker approach.

Towards precision medicine in Alzheimer's disease: deciphering genetic data to establish informative biomarkers

INTRODUCTION

Developing biomarker tools for identification of individuals at high-risk for late-onset Alzheimer's disease (LOAD) is important for prognosis and early treatment. This review focuses on genetic factors and their potential role for precision medicine in LOAD.

AREAS COVERED

APOEe4 is the strongest genetic risk factor for non-Mendelian LOAD, and the APOE-linkage disequilibrium (LD) region has produced the most significant association signal in multi-center genome-wide-association-studies (GWAS). Consideration of extended haplotypes in the APOE-LD region and specifically, non-coding variants in putative enhancer elements, such as the TOMM40-polyT, in-addition to the coding variants that comprise the APOE-genotypes, may be useful for predicting subjects at high-risk of developing LOAD and estimating age-of-onset of early disease-stage symptoms. A genetic-biomarker based on APOE-TOMM40-polyT haplotypes, and age is currently applied in a clinical trial for prevention/delay of LOAD onset. Additionally, we discuss LOAD-GWAS discoveries and the development of new genetic risk scores based on LOAD-GWAS findings other than the APOE-LD region.

EXPERT COMMENTARY

Deciphering the precise causal genetic-variants within LOAD-GWAS regions will advance the development of genetic-biomarkers to complement and refine the APOE-LD region based prediction model. Collectively, the genetic-biomarkers will be translational for early diagnosis and enrichment of clinical trials with subjects at high-risk.

Haplotype-based gene-gene interaction of bone morphogenetic protein 4 and interferon regulatory factor 6 in the etiology of non-syndromic cleft lip with or without cleft palate in a Chilean population

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans, the etiology of which can be dependent on the interactions of multiple genes. We previously reported haplotype associations for polymorphic variants of interferon regulatory factor 6 (IRF6), msh homeobox 1 (MSX1), bone morphogenetic protein 4 (BMP4), and transforming growth factor beta 3 (TGFB3) in Chile. Here, we analyzed the haplotype-based gene-gene interaction for markers of these genes and NSCL/P risk in the Chilean population. We genotyped 15 single nucleoptide polymorphisms (SNPs) in 152 Chilean patients and 164 controls. Linkage disequilibrium (LD) blocks were determined using the Haploview software, and phase reconstruction was performed by the Phase program. Haplotype-based interactions were evaluated using the multifactor dimensionality reduction (MDR) method. We detected two LD blocks composed of two SNPs from BMP4 (Block 1) and three SNPs from IRF6 (Block 2). Although MDR showed no statistical significance for the global interaction model involving these blocks, we found four combinations conferring a statistically significantly increased NSCL/P risk (Block 1-Block 2): T-T/T-G C-G-T/G-A-T; T-T/T-G C-G-C/C-G-C; T-T/T-G G-A-T/G-A-T; and T-T/C-G G-A-T/G-A-T. These findings may reflect the presence of a genomic region containing potential causal variants interacting in the etiology of NSCL/P and may contribute to disentangling the complex etiology of this birth defect.

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.