Efficiency and power in genetic association studies

HASCH - A high-throughput amplicon-based SNP-platform for medicinal cannabis and industrial hemp genotyping applications

BACKGROUND

Cannabis sativa is seeing a global resurgence as a food, fiber and medicinal crop for industrial hemp and medicinal Cannabis industries respectively. However, a widespread moratorium on the use and research of C. sativa throughout most of the 20th century has seen the development of improved cultivars for specific end uses lag behind that of conventional crops. While C. sativa research and development has seen significant investments in the recent past, resulting in a suite of publicly available genomic resources and tools, a versatile and cost-effective mid-density genotyping platform for applied purposes in breeding and pre-breeding is lacking. Here we report on a first mid-density fixed-target SNP platform for C. sativa.

RESULTS

The High-throughput Amplicon-based SNP-platform for medicinal Cannabis and industrial Hemp (HASCH) was designed using a combination of filtering and Integer Linear Programming on publicly available whole-genome sequencing and RNA sequencing data, supplemented with in-house generated genotyping-by-sequencing (GBS) data. HASCH contains 1,504 genome-wide targets of high call rate (97% mean) and even distribution across the genome, designed to be highly informative (> 0.3 minor allele frequency) across both medicinal cannabis and industrial hemp gene pools. Average numbers of mismatch SNP between any two accessions were 251 for medicinal cannabis (N = 116) and 272 for industrial hemp (N = 87). Comparing HASCH data with corresponding GBS data on a collection of diverse C. sativa accessions demonstrated high concordance and resulted in comparable phylogenies and genetic distance matrices. Using HASCH on a segregating F2 population derived from a cross between a tetrahydrocannabinol (THC)-dominant and a cannabidiol (CBD)-dominant accession resulted in a genetic map consisting of 310 markers, comprising 10 linkage groups and a total size of 582.7 cM. Quantitative Trait Locus (QTL) mapping identified a major QTL for CBD content on chromosome 7, consistent with previous findings.

CONCLUSION

HASCH constitutes a versatile, easy to use and cost-effective genotyping solution for the rapidly growing Cannabis research community. It provides consistent genetic fingerprints of 1504 SNPs with wide applicability genetic resource management, quantitative genetics and breeding.

TMEM132C rs7296262 Single-Nucleotide Polymorphism Is Significantly Associated with Nausea Induced by Opioids Administered for Cancer Pain and Postoperative Pain

Opioids are almost mandatorily used for analgesia for cancer pain and postoperative pain. Opioid analgesics commonly induce nausea as a side effect. However, the genetic factors involved are still mostly unknown. To clarify the genetic background of individual differences in the occurrence of nausea during opioid administration, the incidence of nausea was investigated in 331 patients (Higashi-Sapporo Hospital [HS] group) who received morphine chronically for cancer pain treatment and in 2021 patients (Cancer Institute Hospital [CIH] group) who underwent elective surgery under general anesthesia. We conducted a genome-wide association study of nausea in HS samples. Among the top 20 candidate single-nucleotide polymorphisms (SNPs), we focused on the TMEM132C rs7296262 SNP, which has been reportedly associated with psychiatric disorders. The rs7296262 SNP was significantly associated with nausea in both the HS and CIH groups (TT+TC vs. CC; HS group, p = 0.0001; CIH group, p = 0.0064). The distribution of nausea-prone genotypes for the rs7296262 SNP was reversed between HS and CIH groups. These results suggest that the TMEM132C rs7296262 SNP is significantly associated with nausea during opioid use, and the effect of the SNP genotype on nausea is reversed between chronic and acute phases of opioid use.

Disparities in ABO blood type determination across diverse ancestries: a systematic review and validation in the All of Us Research Program

OBJECTIVES

ABO blood types have widespread clinical use and robust associations with disease. The purpose of this study is to evaluate the portability and suitability of tag single-nucleotide polymorphisms (tSNPs) used to determine ABO alleles and blood types across diverse populations in published literature.

MATERIALS AND METHODS

Bibliographic databases were searched for studies using tSNPs to determine ABO alleles. We calculated linkage between tSNPs and functional variants across inferred continental ancestry groups from 1000 Genomes. We compared r2 across ancestry and assessed real-world consequences by comparing tSNP-derived blood types to serology in a diverse population from the All of Us Research Program.

RESULTS

Linkage between functional variants and O allele tSNPs was significantly lower in African (median r2 = 0.443) compared to East Asian (r2 = 0.946, P = 1.1 × 10-5) and European (r2 = 0.869, P = .023) populations. In All of Us, discordance between tSNP-derived blood types and serology was high across all SNPs in African ancestry individuals and linkage was strongly correlated with discordance across all ancestries (ρ = -0.90, P = 3.08 × 10-23).

DISCUSSION

Many studies determine ABO blood types using tSNPs. However, tSNPs with low linkage disequilibrium promote misinference of ABO blood types, particularly in diverse populations. We observe common use of inappropriate tSNPs to determine ABO blood type, particularly for O alleles and with some tSNPs mistyping up to 58% of individuals.

CONCLUSION

Our results highlight the lack of transferability of tSNPs across ancestries and potential exacerbation of disparities in genomic research for underrepresented populations. This is especially relevant as more diverse cohorts are made publicly available.

rs12411980 single-nucleotide polymorphism related to PRTFDC1 expression is significantly associated with phantom tooth pain

Phantom tooth pain (PTP) is one type of non-odontogenic neuropathic toothache, which rarely occurs after appropriate pulpectomy or tooth extraction. The cause of PTP is unknown. We investigated pain-related genetic factors that are associated with PTP. Four pain-associated genes, including G protein-coupled receptor 158 (GPR158) and phosphoribosyl transferase domain containing 1 (PRTFDC1), are adjacent to each other on the human genome. Some of these four genes or their genomic region may be related to PTP. We statistically analyzed associations between single-nucleotide polymorphisms (SNPs) in the genomic region and PTP in patients with PTP (PTP group), other orofacial pain (OFP group), and healthy control subjects. We then performed a database search of expression quantitative trait loci (eQTLs). For the seven SNPs that were significantly associated with PTP even after Bonferroni correction, we focused on the rs12411980 tag SNP (p = 9.42 × 10-4). Statistical analyses of the PTP group and healthy subject groups (group labels: NOC and TD) revealed that the rate of the GG genotype of the rs12411980 SNP was significantly higher in the PTP group than in the healthy subject groups (PTP group vs. NOC group: p = 2.92 × 10-4, PTP group vs. TD group: p = 5.46 × 10-4; percentage of GG: 30% in PTP group, 12% in NOC group, 11% in TD group). These results suggest that the GG genotype of the rs12411980 SNP is more susceptible to PTP. The rs2765697 SNP that is in strong linkage disequilibrium with the rs12411980 SNP is an eQTL that is associated with higher PRTFDC1 expression in the minor allele homozygotes in the healthy subject groups of the rs2765697 SNP. Thus, PRTFDC1 expression similarly increases in the minor allele homozygotes (GG genotype) in the healthy subject groups of the rs12411980 SNP, which would lead to greater susceptibility to PTP.

Single-Nucleotide Polymorphisms of the PAR2 and IL-17A Genes Are Significantly Associated with Chronic Pain

Patients with chronic pain are affected psychologically and socially. There are also individual differences in treatment efficacy. Insufficient research has been conducted on genetic polymorphisms that are related to individual differences in the susceptibility to chronic pain. Autoimmune disorders can lead to inflammation and chronic pain; therefore, we focused on the autoimmune-related protease-activated receptor 2 (PAR2/F2RL1) and interleukin 17A (IL-17A/IL17A) genes. PAR2 and IL-17A are associated with autoimmune diseases that lead to chronic pain, and PAR2 regulates T-helper (Th) cell activation and differentiation. We hypothesized that the PAR2 and IL-17A genes are associated with chronic pain. The present study used a case-control design to statistically examine associations between genetic polymorphisms and the vulnerability to chronic pain. The rs2243057 polymorphism of the PAR2 gene and rs3819025 polymorphism of the IL-17A gene were previously reported to be associated with pain- or autoimmune-related phenotypes. Thus, these polymorphisms were investigated in the present study. We found that both rs2243057 and rs3819025 were significantly associated with a susceptibility to chronic pain. The present findings revealed autoimmune-related genetic factors that are involved in individual differences in chronic pain, further aiding understanding of the pathomechanism that underlies chronic pain and possibly contributing to future personalized medicine.

Haplotype blocks for genomic prediction: a comparative evaluation in multiple crop datasets

In modern plant breeding, genomic selection is becoming the gold standard for selection of superior genotypes. The basis for genomic prediction models is a set of phenotyped lines along with their genotypic profile. With high marker density and linkage disequilibrium (LD) between markers, genotype data in breeding populations tends to exhibit considerable redundancy. Therefore, interest is growing in the use of haplotype blocks to overcome redundancy by summarizing co-inherited features. Moreover, haplotype blocks can help to capture local epistasis caused by interacting loci. Here, we compared genomic prediction methods that either used single SNPs or haplotype blocks with regards to their prediction accuracy for important traits in crop datasets. We used four published datasets from canola, maize, wheat and soybean. Different approaches to construct haplotype blocks were compared, including blocks based on LD, physical distance, number of adjacent markers and the algorithms implemented in the software "Haploview" and "HaploBlocker". The tested prediction methods included Genomic Best Linear Unbiased Prediction (GBLUP), Extended GBLUP to account for additive by additive epistasis (EGBLUP), Bayesian LASSO and Reproducing Kernel Hilbert Space (RKHS) regression. We found improved prediction accuracy in some traits when using haplotype blocks compared to SNP-based predictions, however the magnitude of improvement was very trait- and model-specific. Especially in settings with low marker density, haplotype blocks can improve genomic prediction accuracy. In most cases, physically large haplotype blocks yielded a strong decrease in prediction accuracy. Especially when prediction accuracy varies greatly across different prediction models, prediction based on haplotype blocks can improve prediction accuracy of underperforming models. However, there is no "best" method to build haplotype blocks, since prediction accuracy varied considerably across methods and traits. Hence, criteria used to define haplotype blocks should not be viewed as fixed biological parameters, but rather as hyperparameters that need to be adjusted for every dataset.

The peroxisome proliferator-activated receptor γ coactivator-1 alpha rs8192678 (Gly482Ser) variant and hepatitis B virus clearance

BACKGROUND

Chronic hepatitis B virus (CHB) infection is still incurable a major public health problem. It is yet unclear how host genetic factors influence the development of HBV infection. The peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PPARGC1A) has been shown to regulate hepatitis B virus (HBV). Several reports found that PPARGC1A variants are involved in a number of distinct liver diseases. Here we investigate whether the PPARGC1A rs8192678 (Gly482Ser) variant is involved in the spontaneous clearance of acute HBV infection and if it participates in chronic disease progression in Moroccan patients.

METHODS

Our study included 292 chronic hepatitis B (CHB) patients and 181 individuals who spontaneously cleared-HBV infection. We genotyped the rs8192678 SNP using a TaqMan allelic discrimination assay and then explored its association with spontaneous HBV clearance and CHB progression.

RESULTS

Our data showed that individuals carrying CT and TT genotypes were more likely to achieve spontaneous clearance (OR = 0.48, 95% CI (0.32-0.73), p = 0.00047; OR = 0.28, 95% CI (0.15-0.53), p = 0.00005, respectively). Subjects carrying the mutant allele T were more likely to achieve spontaneous clearance (OR = 0.51, 95% CI (0.38-0.67), P = 2.68E-06). However, when we investigated the impact of rs8192678 on the progression of liver diseases, we neither observe any influence (p > 0.05) nor found any significant association between ALT, AST, HBV viral loads, and the PPARGC1A rs8192678 genotypes in patients with CHB (p > 0.05).

CONCLUSION

Our result suggests that PPARGC1A rs8192678 may modulate acute HBV infection, and could therefore represent a potential predictive marker in the Moroccan population.

Circadian Gene Variants in Diseases

The circadian rhythm is a self-sustaining 24 h cycle that regulates physiological processes within the body, including cycles of alertness and sleepiness. Cells have their own intrinsic clock, which consists of several proteins that regulate the circadian rhythm of each individual cell. The core of the molecular clock in human cells consists of four main circadian proteins that work in pairs. The CLOCK-BMAL1 heterodimer and the PER-CRY heterodimer each regulate the other pair's expression, forming a negative feedback loop. Several other proteins are involved in regulating the expression of the main circadian genes, and can therefore also influence the circadian rhythm of cells. This review focuses on the existing knowledge regarding circadian gene variants in both the main and secondary circadian genes, and their association with various diseases, such as tumors, metabolic diseases, cardiovascular diseases, and sleep disorders.

A diverse ancestrally-matched reference panel increases genotype imputation accuracy in a underrepresented population

Variant imputation, a common practice in genome-wide association studies, relies on reference panels to infer unobserved genotypes. Multiple public reference panels are currently available with variations in size, sequencing depth, and represented populations. Currently, limited data exist regarding the performance of public reference panels when used in an imputation of populations underrepresented in the reference panel. Here, we compare the performance of various public reference panels: 1000 Genomes Project, Haplotype Reference Consortium, GenomeAsia 100 K, and the recent Trans-Omics for Precision Medicine (TOPMed) program, when used in an imputation of samples from the Thai population. Genotype yields were assessed, and imputation accuracies were examined by comparison with high-depth whole genome sequencing data of the same sample. We found that imputation using the TOPMed panel yielded the largest number of variants (~ 271 million). Despite being the smallest in size, GenomeAsia 100 K achieved the best imputation accuracy with a median genotype concordance rate of 0.97. For rare variants, GenomeAsia 100 K also offered the best accuracy, although rare variants were less accurately imputable than common variants (30.3% reduction in concordance rates). The high accuracy observed when using GenomeAsia 100 K is likely attributable to the diverse representation of populations genetically similar to the study cohort emphasizing the benefits of sequencing populations classically underrepresented in human genomics.

Identification of PROK2 gene polymorphisms as predictors of methamphetamine use disorder risk and indicators of craving scale in the Chinese Han population

Background: Methamphetamine use disorder (MUD) has become a global problem due to the highly addictive nature of methamphetamine. Earlier research have demonstrated that PROK2 functions as a compensatory and protective response against neurotoxic stress by stimulating astrocyte reactivity. The aim of our study was to evaluate the correlation between the PROK2 gene and both MUD risk susceptibility and craving scale in the Chinese Han population. Methods: A total of 5,282 participants (1,796 MUD patients and 3,486 controls) were recruited. Seven tag SNPs of the PROK2 gene were chosen and genotyped in the samples. Genetic association analyses were performed to capture the significant SNPs. To investigate the relationship between PROK2 levels and craving scores with the associated-SNP genotypes, we conducted a linear model. Results: SNP rs75433452 was significantly linked with MUD risk (p-value = 1.54 × 10^-8), with the A allele being positively correlated with an increased risk of MUD. Moreover, the average serum level of PROK2 decreased when more copies of the A allele were presented in both MUD patients (p-value = 4.57 × 10^-6) and controls (p-value = 1.13 × 10^-5). Furthermore, the genotypes of SNP rs75433452 were strongly correlated with the craving scores in MUD patients (p-value = 4.05 × 10^-4). Conclusion: Our study identified a significant association signal of the PROK2 gene with MUD risk susceptibility and methamphetamine craving scores in the Chinese Han population, providing potential valuable insights into the underlying mechanisms of METH dependence.

Haplotype-tagged SNPs improve genomic prediction accuracy for Fusarium head blight resistance and yield-related traits in wheat

Linkage disequilibrium (LD)-based haplotyping with subsequent SNP tagging improved the genomic prediction accuracy up to 0.07 and 0.092 for Fusarium head blight resistance and spike width, respectively, across six different models. Genomic prediction is a powerful tool to enhance genetic gain in plant breeding. However, the method is accompanied by various complications leading to low prediction accuracy. One of the major challenges arises from the complex dimensionality of marker data. To overcome this issue, we applied two pre-selection methods for SNP markers viz. LD-based haplotype-tagging and GWAS-based trait-linked marker identification. Six different models were tested with preselected SNPs to predict the genomic estimated breeding values (GEBVs) of four traits measured in 419 winter wheat genotypes. Ten different sets of haplotype-tagged SNPs were selected by adjusting the level of LD thresholds. In addition, various sets of trait-linked SNPs were identified with different scenarios from the training-test combined and only from the training populations. The BRR and RR-BLUP models developed from haplotype-tagged SNPs had a higher prediction accuracy for FHB and SPW by 0.07 and 0.092, respectively, compared to the corresponding models developed without marker pre-selection. The highest prediction accuracy for SPW and FHB was achieved with tagged SNPs pruned at weak LD thresholds (r² < 0.5), while stringent LD was required for spike length (SPL) and flag leaf area (FLA). Trait-linked SNPs identified only from training populations failed to improve the prediction accuracy of the four studied traits. Pre-selection of SNPs via LD-based haplotype-tagging could play a vital role in optimizing genomic selection and reducing genotyping costs. Furthermore, the method could pave the way for developing low-cost genotyping methods through customized genotyping platforms targeting key SNP markers tagged to essential haplotype blocks.

P2X3-P2X7 SNPs and gene-gene and gene-environment interactions on pediatric asthma

BACKGROUND

To investigate the relationship between polymorphisms of P2X3, P2X7 genes and environment interaction with susceptibility of childhood asthma.

METHODS

We conducted a matched case-control study with 170 cases and 175 healthy controls. The rs10896611, rs2276038, rs3781899 in P2X3 and rs1718119, rs3751143 in P2X7 polymorphisms were genotyped using the technique of an improved multiplex ligation detection reaction. Gene-gene, gene-environment and haplotype-environment interactions were tested using the generalized multi-factor dimensionality reduction method.

RESULTS

There were no differences between cases and controls in allele or genotype frequencies of P2X3 and P2X7. The C/C, G/C genotypes of rs10896611, and C/C, C/T genotypes of rs2276038 and G/G, G/A genotypes of rs3781899 were associated with asthmatic cough (p > 0.05). The haplotype GCT of P2X3 reduced the risk of asthma (OR = 0.48, p = 0.048), and the haplotypes AGT (OR = 0.45, p = 0.001) and GCC (OR = 2.16, p = 0.002) were associated with asthmatic cough. The haplotype AA of P2X7 increased risk of asthma severity (p < 0.05). The three-locus model indicated a potential haplotype-environment interaction in GCT, ETS, and pet (p = 0.001).

CONCLUSIONS

The rs10896611, rs2276038 and rs3781899 of P2X3 minor alleles increased the risk of asthmatic cough. Haplotype GCT of P2X3 was a protective factor for asthma, the haplotype AGT was a protective factor and GCC was a risk factor for asthma with cough. In addition, the interactions of haplotype GCT of P2X3, ETS and pet may increase an individual's susceptibility to asthma.

The rs216009 single-nucleotide polymorphism of the CACNA1C gene is associated with phantom tooth pain

Phantom tooth pain (PTP) is a rare and specific neuropathic pain that occurs after pulpectomy and tooth extraction, but its cause is not understood. We hypothesized that there is a genetic contribution to PTP. The present study focused on the CACNA1C gene, which encodes the α1C subunit of the Cav1.2 L-type Ca2+ channel (LTCC) that has been reported to be associated with neuropathic pain in previous studies. We investigated genetic polymorphisms that contribute to PTP. We statistically examined the association between genetic polymorphisms and PTP vulnerability in 33 patients with PTP and 118 patients without PTP but with pain or dysesthesia in the orofacial region. From within and around the CACNA1C gene, 155 polymorphisms were selected and analyzed for associations with clinical data. We found that the rs216009 single-nucleotide polymorphism (SNP) of the CACNA1C gene in the recessive model was significantly associated with the vulnerability to PTP. Homozygote carriers of the minor C allele of rs216009 had a higher rate of PTP. Nociceptive transmission in neuropathic pain has been reported to involve Ca2+ influx from LTCCs, and the rs216009 polymorphism may be involved in CACNA1C expression, which regulates intracellular Ca2+ levels, leading to the vulnerability to PTP. Furthermore, psychological factors may lead to the development of PTP by modulating the descending pain inhibitory system. Altogether, homozygous C-allele carriers of the rs216009 SNP were more likely to be vulnerable to PTP, possibly through the regulation of intracellular Ca2+ levels and affective pain systems, such as those that mediate fear memory recall.

Variants Tagging LGALS-3 Haplotype Block in Association with First Myocardial Infarction and Plasma Galectin-3 Six Months after the Acute Event

Galectin-3 is encoded by LGALS-3, located in a unique haplotype block in Caucasians. According to the Tagger server, rs4040064, rs11628437, and rs7159490 cover 82% (r² > 0.8) of the genetic variance of this HapBlock. Our aims were to examine the association of their haplotypes with first myocardial infarction (MI), changes in left ventricular echocardiographic parameters over time, and impact on plasma galectin-3 and LGALS-3 mRNA in peripheral blood mononuclear cells, both 6 months post-MI. The study group consisted of 546 MI patients and 323 controls. Gene expression was assessed in 92 patients and plasma galectin-3 in 189 patients. Rs4040064, rs11628437, rs7159490, and LGALS-3 mRNA expression were detected using TaqMan^® technology. Plasma galectin-3 concentrations were determined by the ELISA method. We found that the TGC haplotype could have a protective effect against MI (adjusted OR 0.19 [0.05-0.72], p = 0.015) and that the GAC haplotype had significantly higher galectin-3 concentrations (48.3 [37.3-59.4] ng/mL vs. 18.9 [14.5-23.4] ng/mL, p < 0.0001), both in males and compared to the referent haplotype GGC. Higher plasma Gal-3 was also associated with higher NYHA class and systolic dysfunction. Our results suggest that variants tagging LGALS-3 HapBlock could reflect plasma Gal-3 levels 6 months post-MI and may have a potential protective effect against MI in men. Further replication, validation, and functional studies are needed.

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ∼5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

The Roles of MTRR and MTHFR Gene Polymorphisms in Colorectal Cancer Survival

Background: Paradoxically epidemiological data illustrate a negative relationship between dietary folate intake and colorectal cancer (CRC) risk. The occurrence and progression of CRC may be influenced by variants in some key enzyme coding genes in the folate metabolic pathway. We investigated the correlation between genetic variants in methionine synthase reductase (MTRR) and methylenetetrahydrofolate reductase (MTHFR) and CRC survival. Methods: This study used data collected from the Newfoundland Familial Colorectal Cancer Study. A total of 532 patients diagnosed with CRC for the first time from 1999 to 2003 were enrolled, and their mortality were tracked until April 2010. DNA samples were genotyped by Illumina’s integrated quantum 1 million chip. Cox models were established to assess 33 tag single-nucleotide polymorphisms in MTRR and MTHFR in relation to overall survival (OS), disease-free survival (DFS) and CRC-specific survival. Results: The MTRR and MTHFR genes were associated with DFS and CRC-specific survival in CRC patients at the gene level. After multiple comparison adjustment, MTRR rs1801394 A (vs. G) allele was associated with increased DFS (p = 0.024), while MTHRT rs3737966 (G vs. A), rs4846049 (T vs. G), rs1476413 (A vs. G), rs1801131 (C vs. A), rs12121543 (A vs. C), rs1801133 (C vs. T), rs4846052 (T vs. C), rs2066471 (A vs. G) and rs7533315 (T vs. C) were related to worse CRC-specific survival. Additionally, significant interactions were seen among pre-diagnostic alcohol consumption with MTRR rs1801394, rs3776467, rs326124, rs162040, and rs3776455, with superior OS associated with those protective variant alleles limited to patients with alcohol consumption under the median. The MTHFR rs3737966 (G vs. A) allele seemed to be detrimental to CRC survival only among subjects with fruit intake below the median. Conclusions: Polymorphic variants in MTRR and MTHFR genes that code for key enzymes for folate metabolism may be associated with survival in patients with CRC. The gene-CRC outcome association seems modulated by alcohol drinking and fruit intake.

Polymorphic variants of IGF2BP3 and SENCR have an impact on predisposition and/or progression of Ewing sarcoma

Ewing sarcoma (EWS), the second most common malignant bone tumor in children and adolescents, occurs abruptly without clear evidence of tumor history or progression. Previous association studies have identified some inherited variants associated with the risk of developing EWS but a common picture of the germline susceptibility to this tumor remains largely unclear. Here, we examine the association between thirty single nucleotide polymorphisms (SNPs) of the IGF2BP3, a gene that codes for an oncofetal RNA-binding protein demonstrated to be important for EWS patient’s risk stratification, and five SNPs of SENCR, a long non-coding RNA shown to regulate IGF2BP3. An association between polymorphisms and EWS susceptibility was observed for three IGF2BP3 SNPs - rs112316332, rs13242065, rs12700421 - and for four SENCR SNPs - rs10893909, rs11221437, rs12420823, rs4526784 -. In addition, IGF2BP3 rs34033684 and SENCR rs10893909 variants increased the risk for female respect to male subgroup when carried together, while IGF2BP3 rs13242065 or rs76983703 variants reduced the probability of a disease later onset (> 14 years). Moreover, the absence of IGF2BP3 rs10488282 variant and the presence of rs199653 or rs35875486 variant were significantly associated with a worse survival in EWS patients with localized disease at diagnosis. Overall, our data provide the first evidence linking genetic variants of IGF2BP3 and its modulator SENCR to the risk of EWS development and to disease progression, thus supporting the concept that heritable factors can influence susceptibility to EWS and may help to predict patient prognosis.

[Association study between haplotypes of WNT signaling pathway genes and nonsyndromic oral clefts among Chinese Han populations]

OBJECTIVE

To explore whether WNT signaling pathway genes were associated with non-syndromic oral clefts (NSOC) based on haplotypes analyses among 1 008 Chinese NSOC case-parent trios.

METHODS

The genome-wide association study (GWAS) data of 806 Chinese non-syndromic cleft lip with or without cleft palate (NSCL/P) trios and 202 Chinese non-syndromic cleft palate (NSCP) case-parent trios were drawn from the International Consortium to Identify Genes and Interactions Controlling Oral Clefts (ICOCs) study GWAS data set, whose Chinese study population were recruited from four provinces in China, namely Taiwan, Shandong, Hubei, and Sichuan provinces. The process of DNA genotyping was conducted by the Center for Inherited Disease Research in the Johns Hopkins University, using Illumina Human610-Quad v.1_B Bead Chip. The method of sliding windows was used to determine the haplotypes for analyses, including 2 SNPs haplotypes and 3 SNPs haplotypes. Haplotypes with a frequency lower than 1% were excluded for further analyses. To further assess the association between haplotypes and NSOC risks, and the transmission disequilibrium test (TDT) was performed. The Bonferroni method was adopted to correct multiple tests in the study, with which the threshold of statistical significance level was set as P < 0.05 divided by the number of tests, e.g P < 3.47×10^-4 in the current stu-dy. All the statistical analyses were performed by using plink (v1.07).

RESULTS

After quality control, a total of 144 single nucleotide polymorphisms (SNPs) mapped in seven genes in WNT signaling pathway were included for the analyses among the 806 Chinese NSCL/P trios and 202 Chinese NSCP trios. A total of 1 042 haplotypes with frequency higher than 1% were included for NSCL/P analyses and another 1 057 haplotypes with frequency higher than 1% were included for NSCP analyses. Results from the TDT analyses showed that a total of 69 haplotypes were nominally associated with the NSCL/P risk among Chinese (P < 0.05). Another 34 haplotypes showed nominal significant association with the NSCP risk among Chinese (P < 0.05). However, none of these haplotypes reached pre-defined statistical significance level after Bonferroni correction (P>3.47×10^-4).

CONCLUSION

This study failed to observe any statistically significant associations between haplotypes of seven WNT signaling pathway genes and the risk of NSOC among Chinese. Further studies are warranted to replicate the findings here.

Single-nucleotide polymorphisms of the SLC17A9 and P2RY12 genes are significantly associated with phantom tooth pain

Phantom tooth pain (PTP) is a rare and specific neuropathic pain that occurs after pulpectomy and tooth extraction, but its cause is not understood. We hypothesized that there is a genetic contribution to PTP. We focused on solute carrier family 17 member 9 (SLC17A9)/vesicular nucleotide transporter (VNUT) and purinergic receptor P2Y12 (P2RY12), both of which have been associated with neuropathic pain and pain transduction signaling in the trigeminal ganglion in rodents. We sought to corroborate these associations in humans. We investigated gene polymorphisms that contribute to PTP. We statistically examined the association between genetic polymorphisms and PTP vulnerability in 150 patients with orofacial pain, including PTP, and 500 healthy subjects. We found that the rs735055 polymorphism of the SLC17A9 gene and rs3732759 polymorphism of the P2RY12 gene were associated with the development of PTP. Carriers of the minor allele of rs735055 and individuals who were homozygous for the major allele of rs3732759 had a higher rate of PTP. Carriers of the minor allele of rs735055 reportedly had high SLC17A9 mRNA expression in the spinal cord, which may increase the storage and release of adenosine triphosphate. Individuals who were homozygous for the major allele of rs3732759 may have higher P2RY12 expression that is more active in microglia. Therefore, these carriers may be more susceptible to PTP. These results suggest that specific genetic polymorphisms of the SLC17A9 and P2RY12 genes are involved in PTP. This is the first report on genes that are associated with PTP in humans.

Gene-environment interactions between CREB1 and childhood maltreatment on aggression among male Chinese adolescents

Both the genetic and environmental factors may affect aggression susceptibility. However, the conclusions of these associations remain discrepant. In addition, studies that explored the association between CREB1 and aggression were meager. The aim of our present study was to assess whether CREB1 polymorphisms were related to aggression and also to explore the interactive effects of CREB1 variants and childhood maltreatment on aggression. A total of 488 individuals with aggressive behavior and 488 controls were recruited. Aggression and childhood maltreatment were surveyed by standardized self-administered questionnaires. Buccal cells were also obtained and genotyping was conducted using SNPscan. Logistic regressions were applied to investigate both individual effects of CREB1 polymorphisms and the interactive influences with childhood maltreatment on aggression. We found that adolescents who carried the rs4675690 T allele in CREB1 showed a higher level of aggression compared with those who carried wildtype genotypes (CC) under the dominant model (OR = 1.67, 95% CI, 1.16–2.40) after controlling for age and childhood maltreatment. Moreover, we also found that rs4675690 T allele had a synergic additive interaction with childhood sexual abuse and emotional neglect on aggression. The significant interactive effects of CREB1 polymorphisms and childhood maltreatment on aggression were reported for the first time.

Using population-specific add-on polymorphisms to improve genotype imputation in underrepresented populations

Genome-wide association studies rely on the statistical inference of untyped variants, called imputation, to increase the coverage of genotyping arrays. However, the results are often suboptimal in populations underrepresented in existing reference panels and array designs, since the selected single nucleotide polymorphisms (SNPs) may fail to capture population-specific haplotype structures, hence the full extent of common genetic variation. Here, we propose to sequence the full genomes of a small subset of an underrepresented study cohort to inform the selection of population-specific add-on tag SNPs and to generate an internal population-specific imputation reference panel, such that the remaining array-genotyped cohort could be more accurately imputed. Using a Tanzania-based cohort as a proof-of-concept, we demonstrate the validity of our approach by showing improvements in imputation accuracy after the addition of our designed add-on tags to the base H3Africa array.

Genome-wide association studies, which study the association between genetic variants and various phenotypes, typically rely on genotyping arrays. Only a small proportion of genetic variants within the genome are typed on genotyping arrays. Untyped variants are statistically inferred through a process known as genotype imputation, where correlations between variants (haplotypes) observed in external reference panels are leveraged to infer untyped variants in the study population. However, for study populations that are underrepresented in existing reference panels, the quality of imputation is often sub-optimal. This is because typed variants incorporated on existing genotyping arrays can be unsuitable for the study population, and haplotype structures can be different between the reference and the study population. Here, we illustrate an approach to select a custom set of population-specific typed variants to improve genotype imputation in such underrepresented populations.

Polymorphisms in the mTOR-PI3K-Akt pathway, energy balance-related exposures and colorectal cancer risk in the Netherlands Cohort Study

BACKGROUND

The mTOR-PI3K-Akt pathway influences cell metabolism and (malignant) cell growth. We generated sex-specific polygenic risk scores capturing natural variation in 7 out of 10 top-ranked genes in this pathway. We studied the scores directly and in interaction with energy balance-related factors (body mass index (BMI), trouser/skirt size, height, physical activity, and early life energy restriction) in relation to colorectal cancer (CRC) risk in the Netherlands Cohort Study (NLCS) (n=120,852). The NLCS has a case-cohort design and 20.3 years of follow-up. Participants completed a baseline questionnaire on diet and cancer in 1986 when 55-69 years old. ~75% of the cohort returned toenail clippings used for DNA isolation and genotyping (n subcohort=3,793, n cases=3,464). To generate the scores, the dataset was split in two and risk alleles were defined and weighted based on sex-specific associations with CRC risk in the other dataset half, because there were no SNPs in the top-ranked genes associated with CRC risk in previous genome-wide association studies at a significance level p<1*10^-5.

RESULTS

Cox regression analyses showed positive associations between the sex-specific polygenic risk scores and colon but not rectal cancer risk in men and women, with hazard ratios for continuously modeled scores close to 1.10. There was no modifying effect observed of the scores on associations between the energy balance-related factors and CRC risk. However, BMI (in men), non-occupational physical activity (in women), and height (in men and women) were associated with the risk of CRC, in particular (proximal and distal) colon cancer, in the direction as expected in the lower tertiles of the sex-specific polygenic risk scores.

CONCLUSIONS

Current data suggest that the mTOR-PI3K-Akt pathway may be involved in colon cancer development. This study thereby sheds more light on colon cancer etiology through use of genetic variation in the mTOR-PI3K-Akt pathway.

Genetic Analysis Reveals the Prognostic Significance of the DNA Mismatch Repair Gene MSH2 in Advanced Prostate Cancer

Simple Summary

Androgen deprivation therapy is the most effective and widely used treatment for advanced prostate cancer, but its efficacy is highly variable among patients. Therefore, the identification of potent prognostic biomarkers is needed to determine patients at risk. We demonstrated that MSH2 rs1400633 was notably associated with patient survival during androgen deprivation therapy even after adjustment for clinical predictors and false discovery rate correction. Furthermore, our meta-analyses demonstrated that the MSH2 gene is highly expressed in prostate cancer and correlates positively with poor prognosis for this disease.

Abstract

DNA damage repair is frequently dysregulated in advanced prostate cancer and has been linked to cancer susceptibility and survival outcomes. The aim of this study is to assess the influence of genetic variants in DNA damage repair pathways on the prognosis of prostate cancer. Specifically, 167 single nucleotide polymorphisms (SNPs) in 18 DNA damage repair pathway genes were assessed for association with cancer-specific survival (CSS), overall survival (OS), and progression-free survival (PFS) in a cohort of 630 patients with advanced prostate cancer receiving androgen deprivation therapy. Univariate analysis identified four SNPs associated with CSS, four with OS, and two with PFS. However, only MSH2 rs1400633 C > G showed a significant association upon multivariate analysis and multiple testing adjustments (hazard ratio = 0.75, 95% confidence interval = 0.63–0.90, p = 0.002). Furthermore, rs1400633 risk allele C increased MSH2 expression in the prostate and other tissues, which correlated with more aggressive prostate cancer characteristics. A meta-analysis of 31 gene expression datasets revealed significantly higher MSH2 expression in prostate cancer than in normal tissues (p < 0.001), and this high expression was associated with a poor prognosis of prostate cancer (p = 0.002). In summary, we identified MSH2 rs1400633 as an independent prognostic biomarker for prostate cancer survival, and the association of MSH2 with cancer progression lends relevance to our findings.

Interpretation of Manhattan Plots and Other Outputs of Genome-Wide Association Studies

With increasing marker density, estimation of recombination rate between a marker and a causal mutation using linkage analysis becomes less important. Instead, linkage disequilibrium (LD) becomes the major indicator for gene mapping through genome-wide association studies (GWAS). In addition to the linkage between the marker and the causal mutation, many other factors may contribute to the LD, including population structure and cryptic relationships among individuals. As statistical methods and software evolve to improve statistical power and computing speed in GWAS, the corresponding outputs must also evolve to facilitate the interpretation of input data, the analytical process, and final association results. In this chapter, our descriptions focus on (1) considerations in creating a Manhattan plot displaying the strength of LD and locations of markers across a genome; (2) criteria for genome-wide significance threshold and the different appearance of Manhattan plots in single-locus and multiple-locus models; (3) exploration of population structure and kinship among individuals; (4) quantile-quantile (QQ) plot; (5) LD decay across the genome and LD between the associated markers and their neighbors; (6) exploration of individual and marker information on Manhattan and QQ plots via interactive visualization using HTML. The ultimate objective of this chapter is to help users to connect input data to GWAS outputs to balance power and false positives, and connect GWAS outputs to the selection of candidate genes using LD extent.

Tag Variants of LGALS-3 Containing Haplotype Block in Advanced Carotid Atherosclerosis

OBJECTIVES

Galectin-3 affects a variety of biological processes. It is encoded by LGALS-3, located in unique haplotype block in Caucasians. Most of the studies regarding the gal-3 role in atherosclerosis are focused exclusively on protein/mRNA levels. Genetic analyses of LGALS-3 are scarce. We sought to thoroughly examine the genetic background of gal-3 and to analyze tag variants that cover more than 80% variability of the LGALS-3 containing hap-block in association with carotid plaque presence (CPP). According to Tagger server, rs4040064 G/T, rs11628437 G/A and rs7159490 C/T cover 82% (r² > 0.8) of the genetic variance of this hap-block. Our aims were to investigate possible association of rs4040064, rs11628437 and rs7159490 haplotypes with CPP in patients with advanced carotid atherosclerosis (CA) and to analyze their possible effect on LGALS-3 mRNA expression in carotid plaques.

MATERIALS AND METHODS

Study group consisted of 468 patients and 296 controls. Rs4040064, rs11628437, rs7159490 and LGALS-3 mRNA expression were detected by TaqMan® technology.

RESULTS

We have found that haplotype TAC was associated with the cerebrovascular insult (CVI) occurrence (OR = 1.68, 95% CI = 1.09-2.58, p = 0.02), compared to the referent haplotype. OR was adjusted for hypertension, age and BMI. TAC also showed higher, but not statistically significant, LGALS-3 expression in carotid plaques.

CONCLUSIONS

Our results suggest that rs4040064, rs11628437 and rs7159490 bear no association with CPP, neither they affect LGALS-3 mRNA in carotid plaques. However, we showed a significant association of haplotype TAC with the CVI occurrence in CA patients from Serbia. Replication and validation of our results are required.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P &lt; 10−9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P &lt; 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Genome-wide association study identifies new loci associated with noise-induced tinnitus in Chinese populations

Background

Tinnitus is an auditory phantom sensation in the absence of an acoustic stimulus, which affects nearly 15% of the population. Excessive noise exposure is one of the main causes of tinnitus. To now, the knowledge of the genetic determinants of susceptibility to tinnitus remains limited.

Results

We performed a two-stage genome-wide association study (GWAS) and identified that two single nucleotide polymorphisms (SNPs), rs2846071 located in the intergenic region at 11q13.5 (odds ratio [OR] = 2.14, 95% confidence interval [CI] = 1.96–3.40, combined P = 4.89 × 10^− 6) and rs4149577 located in the intron of TNFRSF1A gene at 12p13.31 (OR = 2.05, 95% CI = 1.89–2.51, combined P = 6.88 × 10^− 6), are significantly associated with the susceptibility to noise-induced tinnitus. Furthermore, the expression quantitative trait loci (eQTL) analyses revealed that rs2846071 is significantly correlated with the expression of WNT11 gene, and rs4149577 with the expression of TNFRSF1A gene in multiple brain tissues (all P < 0.05). The newly identified candidate gene WNT11 is involved in Wnt pathway, and TNFRSF1A in the tumor necrosis factor pathway, respectively. Pathway enrichment analyses also showed that these two pathways are closely relevant to tinnitus.

Conclusions

Our findings highlight two novel loci at 11q13.5 and 12p13.31 conferring susceptibility to noise-induced tinnitus. and suggest that the WNT11 and TNFRSF1A genes might be the candidate causal targets of 11q13.5 and 12p13.31 loci, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00987-y.

SNX29, a new susceptibility gene shared with major mental disorders in Han Chinese population

OBJECTIVES

Environmental and genetic factors play important roles in the development of schizophrenia (SCZ), bipolar disorder (BPD) or major depressive disorder (MDD). Some risk loci are identified with shared genetic effects on major psychiatric disorders. To investigate whether SNX29 gene played a significant role in these psychiatric disorders in the Han Chinese population.

METHODS

We focussed on 11 single-nucleotide polymorphisms (SNPs) harbouring SNX29 gene and carried out case-control studies in patients with SCZ (n = 1248), BPD (n = 1344), or MDD (n = 1056), and 1248 healthy controls (HC) recruited from the Han Chinese population. We constructed weighted gene co-expression network analysis (WGCNA) and extracted significant modules by R package.

RESULTS

We found that rs3743592 was significantly associated with MDD and rs6498263 with BPD in both allele and genotype distributions. Before correction, rs3743592 showed allelic and genotypic significance with SCZ, rs6498263 showed allelic significance with SCZ. WGCNA identified top 10 modules of co-expressed genes. Gene Ontology (GO) and pathway analysis were used to examine the functions of SNX29, which revealed that SNX29 was involved in the regulation of a number of biological processes, such as TGF-beta, ErbB, and Wnt signalling pathway, etc.

CONCLUSIONS

Our results supported common risk factors in SNX29 might share among these three mental disorders in the Han Chinese population.

Genetic susceptibility to multiple sclerosis in African Americans

Objective

To explore the nature of genetic-susceptibility to multiple sclerosis (MS) in African-Americans.

Background

Recently, the number of genetic-associations with MS has exploded although the MS-associations of specific haplotypes within the major histocompatibility complex (MHC) have been known for decades. For example, the haplotypes HLA-DRB1*15:01~HLA-DQB1*06:02, and HLA-DRB1*03:01~ HLA-DQB1*02:01 have odds ratios (ORs) for an MS-association orders of magnitude stronger than many of these newly-discovered associations. Nevertheless, all these haplotypes are part of much larger conserved extended haplotypes (CEHs), which span both the Class I and Class II MHC regions. African-Americans are at greater risk of developing MS compared to a native Africans but at lesser risk compared to Europeans. It is the purpose of this manuscript to explore the relationship between MS-susceptibility and the CEH make-up of our African-American cohort.

Design/methods

The African-American (AA) cohort consisted of 1,305 patients with MS and 1,155 controls, who self-identified as being African-American. For comparison, we used the 18,492 controls and 11,144 MS-cases from the predominantly European Wellcome Trust Case Control Consortium (WTCCC) and the 28,557 phased native Africans from the multinational “Be the Match” registry. The WTCCC and the African-Americans were phased at each of five HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRB1 and HLA-DQB1) and the at 11 SNPs (10 of which were in non-coding regions) surrounding the Class II region of the DRB1 gene using previously-published probabilistic phasing algorithms.

Results

Of the 32 most frequent CEHs, 18 (56%) occurred either more frequently or exclusively in Africans) whereas 9 (28%) occurred more frequently or exclusively in Europeans. The remaining 5 CEHs occurred in neither control group although, likely, these were African in origin. Eight of these CEHs carried the DRB1*15:03~DQB1*06:02~a36 haplotype and three carried the DRB1*15:01~DQB1*06:02~a1 haplotype. In African Americans, a single-copy of the European CEH (03:01_07:02_07:02_15:01_06:02_a1) was associated with considerable MS-risk (OR = 3.30; p = 0.0001)–similar to that observed in the WTCCC (OR = 3.25; p<10⁻¹⁶⁸). By contrast, the MS-risk for the European CEH (02:01_07:02_07:02_15:01_06:02_a1) was less (OR = 1.49; ns)–again, similar to the WTCCC (OR = 2.2; p<10⁻³⁸). Moreover, four African haplotypes were “protective” relative to a neutral reference, to three European CEHs, and also to the five other African CEHs.

Conclusions

The common CEHs in African Americans are divisible into those that are either African or European in origin, which are derived without modification from their source population. European CEHs, linked to MS-risk, in general, had similar impacts in African-Americans as they did in Europeans. By contrast, African CEHs had mixed MS-risks. For a few, the MS-risk exceeded that in a neutral-reference group whereas, for many others, these CEHs were “protective”–perhaps providing a partial rationale for the lower MS-risk in African-Americans compared to European-Americans.

Interactions between exposure to polycyclic aromatic hydrocarbons and xenobiotic metabolism genes, and risk of breast cancer

PURPOSE

Polycyclic aromatic hydrocarbons (PAHs) are a group of environmental pollutants associated with multiple cancers, including female breast cancer. Several xenobiotic metabolism genes (XMGs), including the CYP450 family, play an important role in activating and detoxifying PAHs, and variations in the activity of the enzymes they encode can impact this process. This study aims to examine the association between XMGs and breast cancer, and to assess whether these variants modify the effects of PAH exposure on breast cancer risk.

METHODS

In a case-control study in Vancouver, British Columbia, and Kingston, Ontario, 1037 breast cancer cases and 1046 controls had DNA extracted from blood or saliva and genotyped for 138 single nucleotide polymorphisms (SNPs) and tagSNPs in 27 candidate XMGs. Occupational PAH exposure was assessed using a measurement-based job-exposure matrix.

RESULTS

An association between genetic variants and breast cancer was observed among six XMGs, including increased risk among the minor allele carriers of AKR1C3 variant rs12387 (OR 2.71, 95% CI 1.42-5.19) and AKR1C4 variant rs381267 (OR 2.50, 95% CI 1.23-5.07). Heterogeneous effects of occupational PAH exposure were observed among carriers of AKR1C3/4 variants, as well as the PTGS2 variant rs5275.

CONCLUSION

Our findings support an association between SNPs of XMGs and female breast cancer, including novel genetic variants that modify the toxicity of PAH exposure. These results highlight the interplay between genetic and environmental factors, which can be helpful in understanding the modifiable risks of breast cancer and its complex etiology.

Genetic susceptibility to multiple sclerosis: interactions between conserved extended haplotypes of the MHC and other susceptibility regions

BACKGROUND

To study the accumulation of MS-risk resulting from different combinations of MS-associated conserved-extended-haplotypes (CEHs) of the MHC and three non-MHC "risk-haplotypes" nearby genes EOMES, ZFP36L1, and CLEC16A. Many haplotypes are MS-associated despite having population-frequencies exceeding the percentage of genetically-susceptible individuals. The basis of this frequency-disparity requires explanation.

METHODS

The SNP-data from the WTCCC was phased at the MHC and three non-MHC susceptibility-regions. CEHs at the MHC were classified into five haplotype-groups: (HLA-DRB1*15:01 ~ DQB1*06:02 ~ a1)-containing (H +); extended-risk (ER); all-protective (AP); neutral (0); and the single-CEH (c1). MS-associations for different "risk-combinations" at the MHC and other non-MHC "risk-loci" and the appropriateness of additive and multiplicative risk-accumulation models were assessed.

RESULTS

Different combinations of "risk-haplotypes" produce a final MS-risk closer to additive rather than multiplicative risk-models but neither model was consistent. Thus, (H +)-haplotypes had greater impact when combined with (0)-haplotypes than with (H +)-haplotypes, whereas, (H +)-haplotypes had greater impact when combined with a (c1)-haplotypes than with (0)-haplotypes. Similarly, risk-genotypes (0,H +), (c1,H +), (H + ,H +) and (0,c1) were additive with risks from non-MHC risk-loci, whereas risk-genotypes (ER,H +) and (AP,c1) were unaffected.

CONCLUSIONS

Genetic-susceptibility to MS is essential for MS to develop but actually developing MS depends heavily upon both an individual's particular combination of "risk-haplotypes" and how these loci interact.

Diet Modulates the Effects of Genetic Variants on the Vitamin D Metabolic Pathway and Bone Mineral Density in Mexican Postmenopausal Women

BACKGROUND

Macro- and micronutrients, such as proteins, vitamin D, and calcium (Ca), are important dietary factors that can modify bone mineral density (BMD). Genetic factors can interact with diet, affecting an individual's predisposition to osteoporosis.

OBJECTIVES

This study aimed to evaluate the associations between macro- and micronutrient intakes and BMD in Mexican postmenopausal women, and their interactions with genetic polymorphisms involved in the vitamin D metabolic pathway.

METHODS

We analyzed data from 317 postmenopausal women from the Health Workers Cohort Study, a longitudinal cohort studied in Cuernavaca, Mexico. Postmenopausal women participated in 2 data collection waves (2004-2006 and 2010-2011), with a mean time of 6.4 years. Dietary intake was assessed with a semi-quantitative FFQ. BMD (femoral neck, hip, and lumbar spine) was measured by DXA. Hybrid mixed-effects regression models were used to assess the associations of dietary macro- and micronutrients on BMD, after adjusting for confounding factors and for diet and single nucleotide polymorphism interactions.

RESULTS

At baseline, the median age was 57 years (IQR, 50-64). Mean femoral neck, hip, and lumbar spine BMDs decreased over time. We observed statistically significant longitudinal associations for diet (Ca, vitamin D, magnesium, phosphorus, and protein intake) and BMD. Increases of vitamin D, Ca, and protein intakes by 1 SD were associated with mean increases in the femoral neck BMD (0.083 SD, 0.064 SD, and 0.130 SD, respectively). Multiple significant interactions were identified between several loci (CYP2R1, CYP24A1, CYP27B1, VDR, and DHCR7/NADSYN1) and diet for BMDs (femoral neck, hip, and lumbar spine), mainly for protein intake.

CONCLUSIONS

Our data support associations of vitamin D, Ca, protein, phosphorous, and magnesium consumption with BMD in Mexican postmenopausal women and suggest possible gene-diet interactions. These results could facilitate future personalized nutrition recommendations to help prevent low BMD.

Polymorphisms in eicosanoid-related biosynthesis enzymes associated with acute urticaria/angioedema induced by nonsteroidal anti-inflammatory drug hypersensitivity

BACKGROUND

Nonsteroidal anti-inflammatory drugs (NSAIDs) are the main triggers of drug hypersensitivity, with NSAID-induced acute urticaria/angioedema (NIUA) the most frequent phenotype. NSAID hypersensitivity is caused by cyclooxygenase 1 inhibition, which leads to an imbalance in prostaglandin (PG) and cysteinyl leukotriene (CysLT) synthesis. As only susceptible individuals develop NSAID hypersensitivity, genetic factors are believed to be involved; however, no study has assessed the overall genetic variability of key enzymes in PG and CysLT synthesis in NSAID hypersensitivity.

OBJECTIVES

To evaluate simultaneously variants in the main genes involved in PG and CysLT biosynthesis in NIUA.

METHODS

Two independent cohorts of patients were recruited in Spain, alongside NSAID-tolerant controls. The discovery cohort included only patients with NIUA; the replication cohort included patients with NSAID-exacerbated respiratory disease (NERD). A set of tagging single-nucleotide polymorphisms (tagSNPs) in PTGS1, PTGS2, ALOX5 and LTC4S was genotyped using mass spectrometry coupled with endpoint polymerase chain reaction.

RESULTS

The study included 1272 individuals. Thirty-five tagSNPs were successfully genotyped in the discovery cohort, with three being significantly associated after Bonferroni correction (rs10306194 and rs1330344 in PTGS1; rs28395868 in ALOX5). These polymorphisms were genotyped in the replication cohort: rs10306194 and rs28395868 remained associated with NIUA, and rs28395868 was marginally associated with NERD. Odds ratios (ORs) in the combined analysis (discovery and replication NIUA populations) were 1·7 for rs10306194 [95% confidence interval (CI) 1·34-2·14; P_corrected = 2·83 × 10^-4 ) and 2·19 for rs28395868 (95% CI 1·43-3·36; P_corrected = 0·002).

CONCLUSIONS

Variants of PTGS1 and ALOX5 may play a role in NIUA and NERD, supporting the proposed mechanisms of NSAID-hypersensitivity and shedding light on their genetic basis.

Rare variants discovery by extensive whole-genome sequencing of the Han Chinese population in Taiwan: Applications to cardiovascular medicine

Introduction

A population-specific genomic reference is important for research and clinical practice, yet it remains unavailable for Han Chinese (HC) in Taiwan.

Objectives

We report the first whole genome sequencing (WGS) database of HC (1000 Taiwanese genome (1KTW-WGS)) and demonstrate several applications to cardiovascular medicine.

Methods

Whole genomes of 997 HC were sequenced to at least 30X depth. A total of 20,117 relatively healthy HC individuals were genotyped using a customized Axiom GWAS array. We performed a genome-wide genotype imputation technique using IMPUTE2.

Results

We identified 26.7 million single-nucleotide variants (SNVs) and 4.2 million insertions-deletions. Of the SNVs, 16.1% were novel relative to dbSNP (build 152), and 34.2% were novel relative to gnomAD. A total of 18,450 healthy HC individuals were genotyped using a customized Genome-Wide Association Study (GWAS) array. We identified hypertension-associated variants and developed a hypertension prediction model based on the correlation between the WGS data and GWAS data (combined clinical and genetic models, AUC 0.887), and also identified 3 novel hyperlipidemia-associated variants. Each individual carried an average of 16.42 (SD = 3.72) disease-causing variants. Additionally, we established an online SCN5A (an important cardiac gene) database that can be used to explore racial differences. Finally, pharmacogenetics studies identified HC population-specific SNVs in genes (CYP2C9 and VKORC1) involved in drug metabolism and blood clotting.

Conclusion

This research demonstrates the benefits of constructing a population-specific genomic reference database for precision medicine.

Genetic Variants in Cytosolic Phospholipase A2 Associated With Nonsteroidal Anti-Inflammatory Drug-Induced Acute Urticaria/Angioedema

Nonsteroidal anti-inflammatory drugs (NSAIDs) are among the main triggers of drug hypersensitivity reactions, probably due to their high consumption worldwide. The most frequent type of NSAID hypersensitivity is NSAID cross-hypersensitivity, in which patients react to NSAIDs from different chemical groups in the absence of a specific immunological response. The underlying mechanism of NSAID cross-hypersensitivity has been linked to cyclooxygenase (COX)-1 inhibition causing an imbalance in the arachidonic acid pathway. Despite NSAID-induced acute urticaria/angioedema (NIUA) being the most frequent clinical phenotype, most studies have focused on NSAID-exacerbated respiratory disease. As NSAID cross-hypersensitivity reactions are idiosyncratic, only appearing in some subjects, it is believed that individual susceptibility is under the influence of genetic factors. Although associations with polymorphisms in genes from the AA pathway have been described, no previous study has evaluated the potential role of cytosolic phospholipase A2 (cPLA2) variants. This enzyme catalyzes the initial hydrolysis of membrane phospholipids to release AA, which can be subsequently metabolized into eicosanoids. Here, we analyzed for the first time the overall genetic variation in the cPLA2 gene (PLA2G4A) in NIUA patients. For this purpose, a set of tagging single nucleotide polymorphisms (tagSNPs) in PLA2G4A were selected using data from Europeans subjects in the 1,000 Genomes Project, and genotyped with the iPlex Sequenom MassArray technology. Two independent populations, each comprising NIUA patients and NSAID-tolerant controls, were recruited in Spain, for the purposes of discovery and replication, comprising a total of 1,128 individuals. Fifty-eight tagSNPs were successfully genotyped in the discovery cohort, of which four were significantly associated with NIUA after Bonferroni correction (rs2049963, rs2064471, rs12088010, and rs12746200). These polymorphisms were then genotyped in the replication cohort: rs2049963 was associated with increased risk for NIUA after Bonferroni correction under the dominant and additive models, whereas rs12088010 and rs12746200 were protective under these two inheritance models. Our results suggest a role for PLA2G4A polymorphisms in NIUA. However, further studies are required to replicate our findings, elucidate the mechanistic role, and evaluate the participation of PLA2G4A variants in other phenotypes induced by NSAID cross-hypersensitivity.

How Genetics Might Explain the Unusual Link Between Malaria and COVID-19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-associated coronavirus disease 2019 (COVID-19) pandemic has been the subject of a large number of studies in recent times. Here, starting from the evidence that in Italy, the areas with the lowest number of COVID-19 cases were those with the highest incidence of malaria in the early 1900's, we explore possible inverse relationships between malaria and COVID-19. Indeed, some genetic variants, which have been demonstrated to give an advantage against malaria, can also play a role in the incidence and severity of SARS-CoV-2 infections (e.g., the ACE2 receptor). To verify this scientific hypothesis, we here use public data from whole-genome sequencing (WGS) experiments to extrapolate the genetic information of 46 world populations with matched COVID-19 data. In particular, we focus on 47 genes, including ACE2 and genes which have previously been reported to play a role in malaria. Only common variants (>5%) in at least 30% of the selected populations were considered, and, for this subset, we correlate the intra-population allele frequency with the COVID-19 data (cases/million inhabitants), eventually pinpointing meaningful variants in 6 genes. This study allows us to distinguish between positive and negative correlations, i.e., variants whose frequency significantly increases with increasing or decreasing COVID-19 cases. Finally, we discuss the possible molecular mechanisms associated with these variants and advance potential therapeutic options, which may help fight and/or prevent COVID-19.

Association between MIF gene promoter rs755622 and susceptibility to coronary artery disease and inflammatory cytokines in the Chinese Han population

Macrophage migration inhibitory factor (MIF) is an essential mediator of atherosclerotic plaque progression and instability leading to intracoronary thrombosis, therefore contributing to coronary artery disease (CAD). In this study, we investigated the relationship between MIF gene polymorphism and CAD in Chinese Han population. Three single nucleotide polymorphisms (SNP, rs755622, rs1007888 and rs2096525) of MIF gene were genotyped by TaqMan genotyping assay in 1120 control participants and 1176 CAD patients. Coronary angiography was performed in all CAD patients and Gensini score was used to assess the severity of coronary artery lesions. The plasma levels of MIF and other inflammatory mediators were measured by ELISA. The CAD patients had a higher frequency of CC genotype and C allele of rs755622 compared with that in control subjects (CC genotype: 6.5% vs. 3.9%, P = 0.008, C allele: 24.0% vs. 20.6%, P = 0.005). The rs755622 CC genotype was associated with an increased risk of CAD (OR: 1.804, 95%CI: 1.221-2.664, P = 0.003). CAD patients with a variation of rs755622 CC genotype had significantly higher Gensini score compared with patients with GG or CG genotype (all P < 0.05). In addition, the circulating MIF level was highest in CAD patients carrying rs755622 CC genotype (40.7 ± 4.2 ng/mL) and then followed by GC (37.9 ± 3.4 ng/mL) or GG genotype (36.9 ± 3.7 ng/mL, all P < 0.01). Our study showed an essential relationship between the MIF gene rs755622 variation and CAD in Chinese Han population. Individuals who carrying MIF gene rs755622 CC genotype were more susceptible to CAD and had more severe coronary artery lesion. This variation also had a potential influence in circulating MIF levels.

The nature of genetic and environmental susceptibility to multiple sclerosis

OBJECTIVE

To understand the nature of genetic and environmental susceptibility to multiple sclerosis (MS) and, by extension, susceptibility to other complex genetic diseases.

BACKGROUND

Certain basic epidemiological parameters of MS (e.g., population-prevalence of MS, recurrence-risks for MS in siblings and twins, proportion of women among MS patients, and the time-dependent changes in the sex-ratio) are well-established. In addition, more than 233 genetic-loci have now been identified as being unequivocally MS-associated, including 32 loci within the major histocompatibility complex (MHC), and one locus on the X chromosome. Despite this recent explosion in genetic associations, however, the association of MS with the HLA-DRB1*15:01~HLA-DQB1*06:02~a1 (H+) haplotype has been known for decades.

DESIGN/METHODS

We define the "genetically-susceptible" subset (G) to include everyone with any non-zero life-time chance of developing MS. Individuals who have no chance of developing MS, regardless of their environmental experiences, belong to the mutually exclusive "non-susceptible" subset (G-). Using these well-established epidemiological parameters, we analyze, mathematically, the implications that these observations have regarding the genetic-susceptibility to MS. In addition, we use the sex-ratio change (observed over a 35-year interval in Canada), to derive the relationship between MS-probability and an increasing likelihood of a sufficient environmental exposure.

RESULTS

We demonstrate that genetic-susceptibitly is confined to less than 7.3% of populations throughout Europe and North America. Consequently, more than 92.7% of individuals in these populations have no chance whatsoever of developing MS, regardless of their environmental experiences. Even among carriers of the HLA-DRB1*15:01~HLA-DQB1*06:02~a1 haplotype, far fewer than 32% can possibly be members the (G) subset. Also, despite the current preponderance of women among MS patients, women are less likely to be in the susceptible (G) subset and have a higher environmental threshold for developing MS compared to men. Nevertheless, the penetrance of MS in susceptible women is considerably greater than it is in men. Moreover, the response-curves for MS-probability in susceptible individuals increases with an increasing likelihood of a sufficient environmental exposure, especially among women. However, these environmental response-curves plateau at under 50% for women and at a significantly lower level for men.

CONCLUSIONS

The pathogenesis of MS requires both a genetic predisposition and a suitable environmental exposure. Nevertheless, genetic-susceptibility is rare in the population (< 7.3%) and requires specific combinations of non-additive genetic risk-factors. For example, only a minority of carriers of the HLA-DRB1*15:01~HLA-DQB1*06:02~a1 haplotype are even in the (G) subset and, thus, genetic-susceptibility to MS in these carriers must result from the combined effect this haplotype together with the effects of certain other (as yet, unidentified) genetic factors. By itself, this haplotype poses no MS-risk. By contrast, a sufficient environmental exposure (however many events are involved, whenever these events need to act, and whatever these events might be) is common, currently occurring in, at least, 76% of susceptible individuals. In addition, the fact that environmental response-curves plateau well below 50% (especially in men), indicates that disease pathogenesis is partly stochastic. By extension, other diseases, for which monozygotic-twin recurrence-risks greatly exceed the disease-prevalence (e.g., rheumatoid arthritis, diabetes, and celiac disease), must have a similar genetic basis.

A Workflow for Selection of Single Nucleotide Polymorphic Markers for Studying of Genetics of Ischemic Stroke Outcomes

In this paper we propose a workflow for studying the genetic architecture of ischemic stroke outcomes. It develops further the candidate gene approach. The workflow is based on the animal model of brain ischemia, comparative genomics, human genomic variations, and algorithms of selection of tagging single nucleotide polymorphisms (tagSNPs) in genes which expression was changed after ischemic stroke. The workflow starts from a set of rat genes that changed their expression in response to brain ischemia and results in a set of tagSNPs, which represent other SNPs in the human genes analyzed and influenced on their expression as well.

Common variants in FAN1, located in 15q13.3, confer risk for schizophrenia and bipolar disorder in Han Chinese

Multiple genetic risk factors have been associated with psychiatric disorders which provides the genetic insight to these disorders; however, the etiology of these disorders is still elusive. 15q13.3 was previously associated with schizophrenia, bipolar and other neurodevelopmental disorders. Whereas, the FAN1 which encodes the Fanconi anemia associated nuclease 1 was suggested to be causal gene for 15q13.3 related psychiatric disorders. This study aimed to investigate the association of FAN1 with three major psychiatric disorders. Herein, we conducted a case-control study with the Chinese Han population. Three single nucleotide polymorphisms (SNPs) of FAN1 were genotyped in 1248 schizophrenia cases, 1344 bipolar disorder cases, 1056 major depressive disorder cases and 1248 normal controls. We found that SNPs rs7171212 was associated with bipolar (p_allele = 0.023, p_genotype = 0.022, OR = 0.658) and schizophrenia (p_allele = 0.021, p_genotype = 0.019, OR = 0.645). Whereas, rs4779796 was associated with schizophrenia (p_genotype = 0.001, adjusted p_genotype = 0.003, OR = 1.089). In addition, rs7171212 (adjusted p_allele = 0.018, adjusted p_genotype = 0.018, OR = 0.652) and rs4779796 (adjusted p_genotype = 0.024, OR = 1.12) showed significantly associated with combined cases of schizophrenia and bipolar disorder. Further, meta-analysis was performed with the case-control data and dataset extracted from previously reported genome-wide association study to validate the promising SNPs. Our results provide the new evidence that FAN1 may be a common susceptibility gene for schizophrenia and bipolar disorder in Han Chinese. These novel findings need further validation with larger sample size and functional characterization to understand the underlying pathogenic mechanism behind FAN1 in the prevalence of schizophrenia and bipolar disorders.

A Variant in the Nicotinic Acetylcholine Receptor Alpha 3 Subunit Gene Is Associated With Hypertension Risks in Hypogonadic Patients

Ephb6 gene knockout causes hypertension in castrated mice. EPHB6 controls catecholamine secretion by adrenal gland chromaffin cells (AGCCs) in a testosterone-dependent way. Nicotinic acetylcholine receptor (nAChR) is a ligand-gated Ca²⁺/Na⁺ channel, and its opening is the first signaling event leading to catecholamine secretion by AGCCs. There is a possibility that nAChR might be involved in EPHB6 signaling, and thus sequence variants of its subunit genes are associated with hypertension risks. CHRNA3 is the major subunit of nAChR used in human and mouse AGCCs. We conducted a human genetic study to assess the association of CHRNA3 variants with hypertension risks in hypogonadic males. The study cohort included 1,500 hypogonadic Chinese males with (750 patients) or without (750 patients) hypertension. The result revealed that SNV rs3743076 in the fourth intron of CHRNA3 was significantly associated with hypertension risks in the hypogonadic males. We further showed that EPHB6 physically interacted with CHRNA3 in AGCCs, providing a molecular basis for nAChR being in the EPHB6 signaling pathway.

Muscle-Liver Trafficking of BCAA-Derived Nitrogen Underlies Obesity-Related Glycine Depletion

Glycine levels are inversely associated with branched-chain amino acids (BCAAs) and cardiometabolic disease phenotypes, but biochemical mechanisms that explain these relationships remain uncharted. Metabolites and genes related to BCAA metabolism and nitrogen handling were strongly associated with glycine in correlation analyses. Stable isotope labeling in Zucker fatty rats (ZFRs) shows that glycine acts as a carbon donor for the pyruvate-alanine cycle in a BCAA-regulated manner. Inhibition of the BCAA transaminase (BCAT) enzymes depletes plasma pools of alanine and raises glycine levels. In high-fat-fed ZFRs, dietary glycine supplementation raises urinary acyl-glycine content and lowers circulating triglycerides but also results in accumulation of long-chain acyl-coenzyme As (acyl-CoAs), lower 5' adenosine monophosphate-activated protein kinase (AMPK) phosphorylation in muscle, and no improvement in glucose tolerance. Collectively, these studies frame a mechanism for explaining obesity-related glycine depletion and also provide insight into the impact of glycine supplementation on systemic glucose, lipid, and amino acid metabolism.

The Genetic Polymorphism of CYP3A4 rs2242480 is Associated with Sirolimus Trough Concentrations Among Adult Renal Transplant Recipients

BACKGROUND

The large interindividual variability in the genetic polymorphisms of sirolimus (SIR)- metabolizing enzymes, transporters, and receptors can lead to qualitatively and quantitatively distinct therapeutic responses.

OBJECTIVE

We examined the impact of numerous candidate single-nucleotide polymorphisms (SNPs) involved in the trough concentration of SIR-based immunosuppressant regimen.

METHODS

This is a retrospective, long-term cohort study involving 69 renal allograft recipients. Total DNA was isolated from recipient blood samples and trough SIR concentrations were measured by microparticle enzyme immunoassay. Genome sequence reading was targeted based on next-generation sequencing. The association of tagger SNPs to SIR trough concentrations with non-genetic covariate adjusting was analyzed using logistic regression.

RESULTS

A total of 300 SNPs were genotyped in the recipient DNA samples using target sequencing analysis. Only the SNP of CYP3A4 (Ch7: 99361466 C>T, rs2242480) had a significantly higher association with SIR trough concentration as compared to the other 36 tagger SNPs. The mean trough SIR concentration of patients in the CYP3A4 rs2242480-CC group was more significant compared to that of the CYP3A4 rs2242480-TC and TT group, respectively 533.3; 157.4 and 142.5 (ng/ml)/mg/kg, P<0.0001. After adjusting the SNPs, there was no significant association between clinical factors such as age, follow-up period, the incidence of delayed graft function, immunosuppression protocol, and sex with SIR trough concentration.

CONCLUSION

These findings indicated a significant association of polymorphism in the CYP3A4 (Ch7: 99361466 C>T, rs2242480) with SIR trough concentration after 1-year administration in patients who have undergone kidney transplantation.

Association Analysis Between Common Variants of the TRPM1 Gene and Three Mental Disorders in the Han Chinese Population

Objective: Our study was designed to determine if the TRPM1 gene is associated with any of three mental disorders. The project included a cross disorder meta-analysis and association analysis including 141701 cases and 175248 controls. Materials and Methods: We genotyped eight tag single nucleotide polymorphisms (SNPs) in 1248 unrelated schizophrenia (SCZ) patients, 1056 major depressive disorder patients, 1344 bipolar disorder patients, and 1248 normal controls. We then performed a meta-analysis of 10 GWASs to identify common genetic factors among these three mental disorders. Finally, we performed a meta-analysis of six GWASs to explore the role of rs10162727 in SCZ. Result: Although two haplotypes of the TRPM1 gene, rs1035706-rs10162727 and rs10162727-rs3784599, were identified in the control group, as well as all three disease groups, none of the eight tag SNP associations remained significant after correction for multiple tests. In this cross-disorder meta-analysis of the three diseases, none of the tag SNPs were confirmed to be common among the diseases. In addition, in the meta-analysis conducted for the rs10162727 locus in SCZ, there was no significant association (p-value = 0.84, odds ratio = 1.02 [95% CI = 0.87-1.19]). Conclusion: In the Han Chinese population, no significant evidence was found linking variants of the TRPM1 gene with any of the mental disorders examined.

Characterizing the impact of an exotic soybean line on elite cultivar development

The genetic diversity of North American soybean cultivars has been largely influenced by a small number of ancestors. High yielding breeding lines that possess exotic pedigrees have been developed, but identifying beneficial exotic alleles has been difficult as a result of complex interactions of yield alleles with genetic backgrounds and environments as well as the highly quantitative nature of yield. PI 416937 has been utilized in the development of many high yielding lines that have been entered into the USDA Southern States Uniform Tests over the past ~20 years. The primary goal of this research was to identify genomic regions under breeding selection from PI 416937 and introduce a methodology for identifying and potentially utilizing beneficial diversity from lines prevalent in the ancestry of elite cultivars. Utilizing SoySNP50K Infinium BeadChips, 52 high yielding PI 416937-derived lines as well as their parents were genotyped to identify PI 416937 alleles under breeding selection. Nine genomic regions across three chromosomes and 17 genomic regions across seven chromosomes were identified where PI 416937 alleles were under positive or negative selection. Minimal significant associations between PI 416937 alleles and yield were observed in replicated yield trials of five RIL populations, highlighting the difficulty of consistently detecting yield associations.

Clinical relevance of single nucleotide polymorphisms in the CXCL1 and CXCL12 genes in patients with major trauma

BACKGROUND

Genetic backgrounds have been recognized as significant determinants of susceptibility to sepsis. CXC chemokines play a significant role in innate immunity against infectious diseases. Genetic polymorphisms of CXC chemokine genes have been widely studied in inflammatory and infectious diseases but not in sepsis. Thus, we aimed to investigate the clinical relevance of CXC chemokine gene polymorphisms and susceptibility to sepsis in a traumatically injured population.

METHODS

Thirteen tag single nucleotide polymorphisms were selected from CXC chemokine genes using a multimarker tagging algorithm in the Tagger software. Three independent cohorts of injured patients (n = 1700) were prospectively recruited. Selected single nucleotide polymorphisms were genotyped using an improved multiplex ligation detection reaction method. Cytokine production in lipopolysaccharide-stimulated whole blood was measured using an enzyme-linked immunosorbent assay.

RESULTS

Among the 13 tag single nucleotide polymorphisms, four single nucleotide polymorphisms (rs1429638, rs266087, rs2297630, and rs2839693) were significantly associated with the susceptibility to sepsis, and three (rs3117604, rs1429638, and rs4074) were significantly associated with an increased multiple organ dysfunction score in the derivation cohort. However, only the clinical relevance of rs1429638 and rs266087 was confirmed in the validation cohorts. In addition, rs2297630 was significantly associated with interleukin 6 production.

CONCLUSION

The rs1429638 polymorphism in the CXCL1 gene and the rs2297630 polymorphism in the CXCL12 gene were associated with altered susceptibility to sepsis and might be used as important genetic markers to assess the risks of sepsis in trauma patients.

LEVEL OF EVIDENCE

Prognostic and epidemiologic study, level II.

Association of DEAR1 Tagging Single Nucleotide Polymorphisms With Breast Cancer in a Sample of Colombian Population: A Case Control Study

Purpose:

Ubiquitin ligase genes can act as oncogenes or tumor suppressor genes. They play a role in various diseases, including development and progression of breast cancer; the objective of this study was to evaluate the association of common variants in the ductal-epithelium-associated RING chromosome 1 (DEAR1) gene with breast cancer risk in a sample of Colombian population.

Methods:

We carried out a case-control study to investigate associations of variants in DEAR1 with breast cancer in women from Colombia. Single nucleotide polymorphisms (SNPs) rs584298, rs2927970, rs59983645, and rs599167 were genotyped in 1022 breast cancer cases and 1023 healthy controls using the iPLEX® and Kompetitive Allele Specific PCR (polymerase chain reaction) (KASP) method. The associations between SNPs and breast cancer were examined by conditional logistic regression. The associations between SNPs and epidemiological/histopathological variables were examined by multinomial logistic regression.

Results:

Associations were found between tag SNPs and breast cancer adjusted for the epidemiological risk factors rs584298 genotypes AG and GG (P = .048 and P = .004, respectively). The analysis of the disease characteristics showed that SNP rs584298 (genotype AG) (P = .015) shows association with progesterone receptor (PR) status and (genotype AA) (P = .048) shows association with human epidermal growth factor receptor 2 (HER2) status.

Conclusions:

The SNP rs584298 in DEAR1 showed associations with breast cancer and the expression of HER2 receptor; when this receptor is amplified, the result is aggressive tumoral subtype and expression of PR receptor that is associated with high-proliferative tumor grade. Validation of this SNP is important to establish whether this variant or the tagged variant is the cause for the risk association showed.