Two-stage Genome-wide Association Study Research Articles

Genome-Wide Association Analysis Identifies LILRB2 Gene for Pathological Myopia.

Pathological myopia (PM) is one of the leading causes of blindness, especially in Asia. To identify the genetic risk factors of PM, a two-stage genome-wide association study (GWAS) and replication analysis in East Asian populations is conducted. The analysis identified LILRB2 in 19q13.42 as a new candidate locus for PM. The increased protein expression of LILRB2/Pirb (mouse orthologous protein) in PM patients and myopia mouse models is validated. It is further revealed that the increase in LILRB2/Pirb promoted fatty acid synthesis and lipid accumulation, leading to the destruction of choroidal function and the development of PM. This study revealed the association between LILRB2 and PM, uncovering the molecular mechanism of lipid metabolism disorders leading to the pathogenesis of PM due to LILRB2 upregulation.

Genome-wide meta-analysis identifies ancestry-specific loci for Alzheimer's disease.

Alzheimer's disease (AD) is a devastating neurological disease with complex genetic etiology. Yet most known loci have only identified from the late-onset type AD in populations of European ancestry. We performed a two-stage genome-wide association study (GWAS) of AD totaling 6878 Chinese and 63,926 European individuals. In addition to the apolipoprotein E (APOE) locus, our GWAS of two independent Chinese samples uncovered three novel AD susceptibility loci (KIAA2013, SLC52A3, and TCN2) and a novel ancestry-specific variant within EGFR (rs1815157). More replicated variants were observed in the Chinese (31%) than in the European samples (15%). In combining genome-wide associations and functional annotations, EGFR and TCN2 were prioritized as two of the most biologically significant genes. Phenome-wide Mendelian randomization suggests that high mean corpuscular hemoglobin concentration might protect against AD. The current study reveals novel AD susceptibility loci, emphasizes the importance of diverse populations in AD genetic research, and advances our understanding of disease etiology. Loci KIAA2013, SLC52A3, and TCN2 were associated with Alzheimer's disease (AD) in Chinese populations. rs1815157 within the EGFR locus was associated with AD in Chinese populations. The genetic architecture of AD varied between Chinese and European populations. EGFR and TCN2 were prioritized as two of the most biologically significant genes. High mean corpuscular hemoglobin concentrations might have protective effects against AD.

Integrative Multi-omics Analysis Identifies Genetic Variants Contributing to Non-syndromic Cleft Lip with or without Cleft Palate.

To provide novel insights into the aetiology of non-syndromic cleft lip with or without cleft palate (NSCL/P) by integrating multi-omics data and exploring susceptibility genes associated with NSCL/P. A two-stage genome-wide association study (GWAS) of NSCL/P was performed, involving a total of 1,069 cases and 1,724 controls. Using promoter capture Hi-C (pCHi-C) datasets in human embryonic stem cells (hESC) and chromatin immunoprecipitation sequencing (ChIP-seq) in craniofacial tissues, we filtered out single nucleotide polymorphisms (SNPs) with active cis-regulation and their target genes. Additionally, we employed expression quantitative trait loci (eQTL) analysis to identify candidate genes. Thirteen SNPs were identified as cis-regulation units associated with the risk of NSCL/P. Five of these were proven to be active in chromatin states in early human craniofacial development (rs7218002: odds ratio [OR] 1.50, P = 8.14E-08; rs835367: OR 0.78, P = 3.48E- 05; rs77022994: OR 0.55, P = 1.05E-04; rs961470: OR 0.73, P = 1.38E-04; rs17314727: OR 0.73, P = 1.85E-04). Additionally, pCHi-C and eQTL analysis prioritised three candidate genes (rs7218002: NTN1, rs835367: FGGY, LINC01135). NTN1 and FGGY were expressed in mouse orofacial development. Deficiencies in NTN1, FGGY and LINC01135 were associated with cleft palate and cleft lip, abnormal facial shape and bifid uvula, and abnormality of the face, respectively. Our study identified five SNPs (rs7218002, rs835367, rs77022994, rs961470 and rs17314727) and three susceptibility genes (NTN1, FGGY and LINC01135) associated with NSCL/P. These findings contribute to a better understanding of the genetic factors involved.

Identification of Gene-Allele System Conferring Alkali-Tolerance at Seedling Stage in Northeast China Soybean Germplasm.

Salinization of cultivated soils may result in either high salt levels or alkaline conditions, both of which stress crops and reduce performance. We sampled genotypes included in the Northeast China soybean germplasm population (NECSGP) to identify possible genes that affect tolerance to alkaline soil conditions. In this study, 361 soybean accessions collected in Northeast China were tested under 220 mM NaHCO3:Na2CO3 = 9:1 (pH = 9.8) to evaluate the alkali-tolerance (ATI) at the seedling stage in Mudanjiang, Heilongjiang, China. The restricted two-stage multi-locus model genome-wide association study (RTM-GWAS) with gene-allele sequences as markers (6503 GASMs) based on simplified genome resequencing (RAD-sequencing) was accomplished. From this analysis, 132 main effect candidate genes with 359 alleles and 35 Gene × Environment genes with 103 alleles were identified, explaining 90.93% and 2.80% of the seedling alkali-tolerance phenotypic variation, respectively. Genetic variability of ATI in NECSGP was observed primarily within subpopulations, especially in ecoregion B, from which 80% of ATI-tolerant accessions were screened out. The biological functions of 132 candidate genes were classified into eight functional categories (defense response, substance transport, regulation, metabolism-related, substance synthesis, biological process, plant development, and unknown function). From the ATI gene-allele system, six key genes-alleles were identified as starting points for further study on understanding the ATI gene network.

Genome-wide association with transcriptomics reveals a shade-tolerance gene network in soybean

Shade tolerance is essential for soybeans in inter/relay cropping systems. A genome-wide association study (GWAS) integrated with transcriptome sequencing was performed to identify genes and construct a genetic network governing the trait in a set of recombinant inbred lines derived from two soybean parents with contrasting shade tolerance. An improved GWAS procedure, restricted two-stage multi-locus genome-wide association study based on gene/allele sequence markers (GASM-RTM-GWAS), identified 140 genes and their alleles associated with shade-tolerance index (STI), 146 with relative pith cell length (RCL), and nine with both. Annotation of these genes by biological categories allowed the construction of a protein–protein interaction network by 187 genes, of which half were differentially expressed under shading and non-shading conditions as well as at different growth stages. From the identified genes, three ones jointly identified for both traits by both GWAS and transcriptome and two genes with maximum links were chosen as beginners for entrance into the network. Altogether, both STI and RCL gene systems worked for shade-tolerance with genes interacted each other, this confirmed that shade-tolerance is regulated by more than single group of interacted genes, involving multiple biological functions as a gene network.

A genome-wide association study of 24-hour urinary excretion of endocrine disrupting chemicals

Ubiquitous exposure to environmental endocrine disrupting chemicals (EDCs) instigates a major public health problem, but much remains unknown on the inter-individual differences in metabolism and excretion of EDCs. To examine this we performed a two-stage genome-wide association study (GWAS) for 24-hour urinary excretions of four parabens, two bisphenols, and nine phthalate metabolites. Results showed five genome-wide significant (p-value < 5x10-8) and replicated single nucleotide polymorphisms (SNPs) representing four independent signals that associated with mono-(2-ethyl-5-carboxypentyl) phthalate (MECPP) and mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP). Three of the four signals were located on chromosome 10 in a locus harboring the cytochrome P450 (CYP) genes CYP2C9, CYP2C58P, and CYP2C19 (rs117529685, pMECPP = 5.38x10-25; rs117033379, pMECPP = 1.96x10-19; rs4918798, pMECPP = 4.01x10-71; rs7895726, pMEHHP = 1.37x10-15, r2 with rs4918798 = 0.93). The other signal was on chromosome 6 close to the solute carrier (SLC) genes SLC17A1, SLC17A3, SLC17A4, and SCGN (rs1359232, pMECPP = 7.6x10-16). These four SNPs explained a substantial part (8.3 % - 9.2 %) of the variance in MECPP in the replication cohort. Bioinformatics analyses supported a likely causal role of CYP2C9 and SLC17A1 in metabolism and excretion of MECPP and MEHHP. Our results provide biological insights into mechanisms of phthalate metabolism and excretion with a likely causal role for CYP2C9 and SLC17A1.

The Identification of a Quantative Trait Loci-Allele System of Antixenosis against the Common Cutworm (Spodoptera litura Fabricius) at the Seedling Stage in the Chinese Soybean Landrace Population.

Common cutworm (CCW) is an omnivorous insect causing severe yield losses in soybean crops. The seedling-stage mini-tray identification system with the damaged leaf percentage (DLP) as an indicator was used to evaluate antixenosis against CCW in the Chinese soybean landrace population (CSLRP) under three environments. Using the innovative restricted two-stage multi-locus genome-wide association study procedure (RTM-GWAS), 86 DLP QTLs with 243 alleles (2-11/QTL) were identified, including 66 main-effect loci with 203 alleles and 57 QTL-environment interaction loci with 172 alleles. Among the main-effect loci, 12 large-contribution loci (R2 ≥ 1%) explained 25.45% of the phenotypic variation (PV), and 54 small-contribution loci (R2 < 1%) explained 16.55% of the PV. This indicates that the CSLRP can be characterized with a DLP QTL-allele system complex that has not been found before, except for a few individual QTLs without alleles involved. From the DLP QTL-allele matrix, the recombination potentials expressed in the 25th percentile of the DLP of all possible crosses were predicted to be reduced by 41.5% as the maximum improvement and 14.2% as the maximum transgression, indicating great breeding potential in the antixenosis of the CSLRP. From the QTLs, 62 candidate genes were annotated, which were involved in eight biological function categories as a gene network of the DLP. Changing from susceptible to moderate plus resistant varieties in the CSLRP, 26 QTLs had 32 alleles involved, in which 19 genes were annotated from 25 QTL-alleles, including eight increased negative alleles on seven loci and 11 decreased positive alleles on 11 loci, showing the major genetic constitution changes for the antixenosis enhancement at the seedling stage in the CSLRP.

Gene-allele system of shade tolerance in southern China soybean germplasm revealed by genome-wide association study using gene-allele sequence as markers.

Fifty-three shade tolerance genes with 281 alleles in the SCSGP were identified directly using gene-allele sequence as markers in RTM GWAS, from which optimized crosses, evolutionary motivators, and gene-allele networks were explored. Shade tolerance is a key for optimal cultivation of soybean inter/relay-cropped with corn. To explore the shade tolerance gene-allele system in the southern China soybean germplasm, we proposed using gene-allele sequence markers (GASMs) in a restricted two-stage multi-locus model genome-wide association study (GASM-RTM-GWAS). A representative sample with 394 accessions was tested for their shade tolerance index (STI), in Nanning, China. Through whole-genome re-sequencing, 47,586 GASMs were assembled. From GASM-RTM-GWAS, 53 main-effect STI genes with 281 alleles (2-13 alleles/gene) (totally 63 genes with 308 alleles, including 38 G × E genes with 191 alleles) were identified and then organized into a gene-allele matrix composed of eight submatrices corresponding to geo-seasonal subpopulations. The population featured mild STI changes (1.69 → 1.56-1.82) and mild gene-allele changes (92.5% alleles inherited, 0% alleles excluded, 7.5% alleles emerged) from the primitive (SAIII) to the derived seven subpopulations, but large transgressive recombination potentials and optimal crosses were predicted. The 63 STI genes were annotated into six biological categories (metabolic process, catalytic activity, response to stresses, transcription and translation, signal transduction and transport and unknown functions), interacted as gene networks. From the STI gene-allele system, 38 important alleles of 22 genes were nominated for further in-depth study. GASM-RTM-GWAS performed powerful and efficient in germplasm population genetic study comparing to other procedures through facilitating direct and thorough identification of its gene-allele system, from which genome-wide breeding by design could be achieved, and evolutionary motivators and gene-allele networks could be explored.

An Improved Genome-Wide Association Procedure Explores Gene-Allele Constitutions and Evolutionary Drives of Growth Period Traits in the Global Soybean Germplasm Population.

In soybeans (Glycine max (L.) Merr.), their growth periods, DSF (days of sowing-to-flowering), and DFM (days of flowering-to-maturity) are determined by their required accumulative day-length (ADL) and active temperature (AAT). A sample of 354 soybean varieties from five world eco-regions was tested in four seasons in Nanjing, China. The ADL and AAT of DSF and DFM were calculated from daily day-lengths and temperatures provided by the Nanjing Meteorological Bureau. The improved restricted two-stage multi-locus genome-wide association study using gene-allele sequences as markers (coded GASM-RTM-GWAS) was performed. (i) For DSF and its related ADLDSF and AATDSF, 130-141 genes with 384-406 alleles were explored, and for DFM and its related ADLDFM and AATDFM, 124-135 genes with 362-384 alleles were explored, in a total of six gene-allele systems. DSF shared more ADL and AAT contributions than DFM. (ii) Comparisons between the eco-region gene-allele submatrices indicated that the genetic adaptation from the origin to the geographic sub-regions was characterized by allele emergence (mutation), while genetic expansion from primary maturity group (MG)-sets to early/late MG-sets featured allele exclusion (selection) without allele emergence in addition to inheritance (migration). (iii) Optimal crosses with transgressive segregations in both directions were predicted and recommended for breeding purposes, indicating that allele recombination in soybean is an important evolutionary drive. (iv) Genes of the six traits were mostly trait-specific involved in four categories of 10 groups of biological functions. GASM-RTM-GWAS showed potential in detecting directly causal genes with their alleles, identifying differential trait evolutionary drives, predicting recombination breeding potentials, and revealing population gene networks.

A two-stage genome-wide association study identified four potential early-onset nonsmall cell lung cancer risk loci based on 26,652 participants in Chinese population.

Early-onset lung cancer is rare with an increasing incidence rate. Although several genetic variants have been identified for it with candidate gene approaches, no genome-wide association study (GWAS) has been reported. In this study, a two-stage strategy was adopted: firstly we performed a GWAS to identify variants associated with early-onset nonsmall-cell lung cancer (NSCLC) risk using 2556 cases (age ≤ 50 years) and 13,327 controls by logistic regression model. To further discriminate younger cases from older ones, we took a case-case analysis for the promising variants with above early-onset cases and 10,769 cases (age > 50 years) by Cox regression model. After combining these results, we identified four early-onset NSCLC susceptibility loci at 5p15.33 (rs2853677, odds ratio[OR] = 1.48, 95% confidence interval[CI]: 1.36-1.60, Pcase-control = 3.58 × 10-21 ; hazard ratio[HR] = 1.10, 95% CI: 1.04-1.16, Pcase-case = 6.77 × 10-4 ), 5p15.1 (rs2055817, OR = 1.24, 95% CI: 1.15-1.35, Pcase-control = 1.39 × 10-7 ; HR = 1.08, 95% CI: 1.02-1.14, Pcase-case = 6.90 × 10-3 ), 6q24.2 (rs9403497, OR = 1.24, 95% CI: 1.15-1.35, Pcase-control = 1.61 × 10-7 ; HR = 1.11, 95% CI: 1.05-1.17, Pcase-case = 3.60 × 10-4 ) and 12q14.3 (rs4762093, OR = 1.31, 95% CI: 1.18-1.45, Pcase-control = 1.90 × 10-7 ; HR = 1.10, 95% CI: 1.03-1.18, Pcase-case = 7.49 × 10-3 ). Except for 5p15.33, other loci were found to be associated with NSCLC risk for the first time. All of them had stronger effects in younger patients than in older ones. These results provide a promising overview for early-onset NSCLC genetics.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Identification of a Common Variant for Coronary Heart Disease at PDE1A Contributes to Individualized Treatment Goals and Risk Stratification of Cardiovascular Complications in Chinese Patients With Type 2 Diabetes.

In this study we aim to unravel genetic determinants of coronary heart disease (CHD) in type 2 diabetes (T2D) and explore their applications. We performed a two-stage genome-wide association study for CHD in Chinese patients with T2D (3,596 case and 8,898 control subjects), followed by replications in European patients with T2D (764 case and 4,276 control subjects) and general populations (n = 51,442-547,261). Each identified variant was examined for its association with a wide range of phenotypes and its interactions with glycemic, blood pressure (BP), and lipid controls in incident cardiovascular diseases. We identified a novel variant (rs10171703) for CHD (odds ratio 1.21 [95% CI 1.13-1.30]; P = 2.4 × 10-8) and BP (β ± SE 0.130 ± 0.017; P = 4.1 × 10-14) at PDE1A in Chinese T2D patients but found only a modest association with CHD in general populations. This variant modulated the effects of BP goal attainment (130/80 mmHg) on CHD (Pinteraction = 0.0155) and myocardial infarction (MI) (Pinteraction = 5.1 × 10-4). Patients with CC genotype of rs10171703 had >40% reduction in either cardiovascular events in response to BP control (2.9 × 10-8 < P < 3.6 × 10-5), those with CT genotype had no difference (0.0726 < P < 0.2614), and those with TT genotype had a threefold increase in MI risk (P = 6.7 × 10-3). We discovered a novel CHD- and BP-related variant at PDE1A that interacted with BP goal attainment with divergent effects on CHD risk in Chinese patients with T2D. Incorporating this information may facilitate individualized treatment strategies for precision care in diabetes, only when our findings are validated.

Genome-wide association study of tuberculosis in the western Chinese Han and Tibetan population.

Tuberculosis (TB) remains a serious global public health threat. Accumulated evidence has demonstrated that human susceptibility to TB has a strong genetic basis. And different susceptibility single nucleotide polymorphisms (SNP) have been reported in different studies. To gain greater insight into the host susceptibility to TB, we perform a two-stage genome-wide association study to identify the susceptible loci of TB. In the discovery stage, 3116 (1532TB patients and 1584 healthy controls) and 439 (211TB patients and 228 healthy controls) individuals were genome-wide genotyped from a western Chinese Han and Tibetan population, respectively. Based on the additive genetic model, we discovered 14 and three independent loci that had potential associations with TB susceptibility in the Chinese Han and Tibetan populations, respectively (p<1×10-5). Furthermore, we conducted an imputation-based meta-analysis on another two East Asia cohorts to replicate our findings. We identified one independent locus harbored by the human leukocyte antigen (HLA) class II genes that was genome-wide significantly associated with TB (lead SNP rs111875628 with a p-value of 2.20×10-9). Our findings suggest a novel mechanism of the interaction with the HLA class II genes and reinforce the importance of the HLA class II alleles in response to TB.

Genome-wide association study identifies genetic variants which predict the response of bone mineral density to teriparatide therapy

ObjectivesTeriparatide (TPTD) is an effective treatment for osteoporosis but the individual response to therapy is variable for reasons that are unclear. This study aimed to determine whether the response to...

Abstract P3-05-31: PVT1, CASC16, TNFRSF13B and NAMPT may be the key genes correlated with recurrence and metastasis of breast cancer

Abstract Background: It is vital to find the risk genes and loci related to recurrence and metastasis in breast cancer (BC) patients. The Genome-wide association study (GWAS) can detect the risk of single nucleotide polymorphism (SNP) and help us understand the mechanism of disease progression by analyzing the SNP information of the whole genome. At present, GWAS have identified several genetic susceptibility SNPs for BC, but it is unclear about its significance in the recurrence and metastasis of BC patients. Therefore, in this study we aimed to identify the risk SNPs related to recurrence and metastasis of BC based on GWAS analysis. Methods: This study adopts a two-stage GWAS research strategy. In the first stage, according to propensity score matching, we chose 97 pairs of BC patients with or without recurrence and metastasis who were treated in Shandong Cancer Hospital and Institute from November 2013 to April 2014. Their peripheral blood samples were collected and DNA was extracted for genome-wide detection by Illumina ASA chip. The SNP loci related to recurrence and metastasis of BC were obtained by logistic regression analysis of GWAS. In the second stage, 854 BC patients from the same hospital were verified from May 2014 to June 2015. The SNP genotyping was detected by time-of-flight mass spectrometry to verify its correlation with recurrence and metastasis of BC. The SNP loci and their corresponding genes and pathways were analyzed using DAVID (https://david.ncifcrf.gov/) online enrichment analysis website. Results: In the first stage, 191 SNPs were obtained from GWAS analysis (P&amp;lt; 0.001). Based on these SNPs and their corresponding genes, we further found that glutamatergic synaptic transmission may be related to BC recurrence and metastasis through GO analysis. Similar results were found in calcium signaling pathway and insulin secretion pathway through KEGG enrichment analysis. And 21 SNPs were screened for their association with the tumor evolution process and whether they are located in the gene coding functional region. In the second stage, Cox regression analysis showed that the genotyping of the four SNPs was consistent with the results of GWAS analysis, and these SNPs are located in the intron region of the gene. Specifically, compared with homozygous AA, BC patients carrying TNFRSF13B rs4273077(AG/GG) and PVT1 rs10108514(AG/GG) had a 26% (HR: 0.74, 95%CI: 0.54-1.01) and 27% (HR: 0.73, 95%CI: 0.54-0.98) reduced risk of recurrence and metastasis, respectively; compared with homozygous CC, patients with BC carrying NAMPT rs4730155(CT/TT) had a 35% (HR: 0.65, 95%CI: 0.45-0.93) reduced risk of recurrence and metastasis; while compared with CASC16 rs12920540(CA/AA), BC patients with homozygous CC has a 52% (HR: 1.52, 95%CI: 1.06-2.18) increased risk of recurrence and metastasis. Analysis of TCGA and GTEX databases showed that the expression levels of PVT1 and CASC16 in breast tumors were significantly higher than those in normal tissues and matched para-cancerous tissues (P&amp;lt; 0.001), whereas the expression level of NAMPT in breast tumors was significantly lower than that in normal tissues (P&amp;lt; 0.001). No significant difference was shown about the expression level of TNFRSF13B between breast tumor and normal tissues or matched para-cancerous tissues (P≥0.05). Conclusion: Our study suggested that PVT1, CASC16, TNFRSF13B and NAMPT may be correlated with the risk of recurrence and metastasis of BC. Glutamatergic synaptic transmission, calcium signaling pathway and insulin secretion pathway may be involved in BC recurrence and metastasis. Citation Format: Huihui Li, Sha Yin, Shujuan Sun, Qiaorui Tan, Xiaochu Man, Jie Huang, Xuemei Xie. PVT1, CASC16, TNFRSF13B and NAMPT may be the key genes correlated with recurrence and metastasis of breast cancer [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P3-05-31.

Genome-Wide Association Study in Acute Tubulointerstitial Nephritis.

Polymorphisms of HLA genes may confer susceptibility to acute tubulointerstitial nephritis (ATIN), but small sample sizes and candidate gene design have hindered their investigation. The first genome-wide association study of ATIN identified two significant loci, risk haplotype DRB1*14-DQA1*0101-DQB1*0503 (DR14 serotype) and protective haplotype DRB1*1501-DQA1*0102-DQB1*0602 (DR15 serotype), with amino acid position 60 in the peptide-binding groove P10 of HLA-DR β 1 key. Risk alleles were shared among different causes of ATIN and HLA genotypes associated with kidney injury and immune therapy response. HLA alleles showed the strongest association. The findings suggest that a genetically conferred risk of immune dysregulation is part of the pathogenesis of ATIN. Acute tubulointerstitial nephritis (ATIN) is a rare immune-related disease, accounting for approximately 10% of patients with unexplained AKI. Previous elucidation of the relationship between genetic factors that contribute to its pathogenesis was hampered because of small sample sizes and candidate gene design. We undertook the first two-stage genome-wide association study and meta-analysis involving 544 kidney biopsy-defined patients with ATIN and 2346 controls of Chinese ancestry. We conducted statistical fine-mapping analysis, provided functional annotations of significant variants, estimated single nucleotide polymorphism (SNP)-based heritability, and checked genotype and subphenotype correlations. Two genome-wide significant loci, rs35087390 of HLA-DQA1 ( P =3.01×10 -39 ) on 6p21.32 and rs2417771 of PLEKHA5 on 12p12.3 ( P =2.14×10 -8 ), emerged from the analysis. HLA imputation using two reference panels suggested that HLA-DRB1*14 mainly drives the HLA risk association . HLA-DRB1 residue 60 belonging to pocket P10 was the key amino acid position. The SNP-based heritability estimates with and without the HLA locus were 20.43% and 10.35%, respectively. Different clinical subphenotypes (drug-related or tubulointerstitial nephritis and uveitis syndrome) seemed to share the same risk alleles. However, the HLA risk genotype was associated with disease severity and response rate to immunosuppressive therapy. We identified two candidate genome regions associated with susceptibility to ATIN. The findings suggest that a genetically conferred risk of immune dysregulation is involved in the pathogenesis of ATIN.

Comparison of the loci associated with HbA1c and blood glucose levels identified by a genome-wide association study in the Japanese population.

Hemoglobin A1c (HbA1c) levels are widely employed to diagnose diabetes. However, estimates of the heritability of HbA1c and glucose levels are different. Therefore, we explored HbA1c- and blood glucose-associated loci in a non-diabetic Japanese population. We conducted a two-stage genome-wide association study (GWAS) on variants associated with HbA1c and blood glucose levels in a Japanese population. In the initial stage, data of 4911 participants of the Japan Multi-Institutional Collaborative Cohort (J-MICC) were subjected to discovery analysis. In the second stage, two datasets from the Tohoku Medical Megabank project, with 8175 and 40,519 participants, were used for the replication study. Association of the imputed variants with HbA1c and blood glucose levels was determined via linear regression analyses adjusted for age, sex, body mass index (BMI), smoking, and genetic principal components (PC1-PC10). Moreover, we performed a BMI-stratified GWAS on HbA1c levels in the J-MICC. The discovery analysis and BMI-stratified GWAS results were validated with re-analyses of normalized HbA1c levels adjusted for site in addition to the above, and blood glucose adjusted for fasting time as an additional covariate. Genetic variants associated with HbA1c levels were identified in KCNQ1 and TMC6. None of the genetic variants associated with blood glucose levels in the discovery analysis were replicated. Association of rs2299620 in KCNQ1 with HbA1c levels showed heterogeneity between individuals with BMI ≥ 25kg/m2 and BMI < 25kg/m2. The variant rs2299620 in KCNQ1 might affect HbA1c levels differentially based on BMI grouping in the Japanese population. The online version contains supplementary material available at 10.1007/s13340-023-00618-0.

A two-stage genome-wide association study identifies novel germline genetic variations in CACNA2D3 associated with radiotherapy response in nasopharyngeal carcinoma

BackgroundRadiotherapy (RT) is the standard treatment for nasopharyngeal carcinoma (NPC). However, due to individual differences in radiosensitivity, biomarkers are needed to tailored radiotherapy to cancer patients. However, comprehensive genome-wide radiogenomic studies on them are still lacking. The aim of this study was to identify genetic variants associated with radiotherapy response in patients with NPC.MethodsThis was a large‑scale genome-wide association analysis (GWAS) including a total of 981 patients. 319 individuals in the discovery stage were genotyped for 688,783 SNPs using whole genome-wide screening microarray. Significant loci were further genotyped using MassARRAY system and TaqMan SNP assays in the validation stages of 847 patients. This study used logistic regression analysis and multiple bioinformatics tools such as PLINK, LocusZoom, LDBlockShow, GTEx, Pancan-meQTL and FUMA to examine genetic variants associated with radiotherapy efficacy in NPC.ResultsAfter genome-wide level analysis, 19 SNPs entered the validation stage (P < 1 × 10− 6), and rs11130424 ultimately showed statistical significance among these SNPs. The efficacy was better in minor allele carriers of rs11130424 than in major allele carriers. Further stratified analysis showed that the association existed in patients in the EBV-positive, smoking, and late-stage (III and IV) subgroups and in patients who underwent both concurrent chemoradiotherapy and induction/adjuvant chemotherapy.ConclusionOur study showed that rs11130424 in the CACNA2D3 gene was associated with sensitivity to radiotherapy in NPC patients.Trial registration number: Effect of genetic polymorphism on nasopharyngeal carcinoma chemoradiotherapy reaction, ChiCTR-OPC-14005257, Registered 18 September 2014, http://www.chictr.org.cn/showproj.aspx?proj=9546.

Genome-Wide Association Studies (GWAS).

Most of the breeding targets are quantitative traits. In exploring the quantitative trait locus (QTL) system of a trait, linkage mapping was established using sparse polymerase chain reaction (PCR) markers. With the genome-wide sequencing technology advanced, genome-wide association study (GWAS) was developed for natural (germplasm) populations using dense genomic markers, which facilitates the identification of the complete QTL system with their multiple alleles on genomic locations. GWAS makes use of the linkage disequilibrium (LD) due to historical saturate recombination and high-density genomic markers to detect QTLs through statistical test for the association between molecular markers and phenotypes. However, due to inbreeding and mixture of source populations, the germplasm population often has complex and unknown structure, which leads to false positives/negatives in GWAS. Various GWAS methods have been proposed to reduce false positives/negatives, including those of the general linear model and the mixed linear model, which focused mainly on finding a handful of major QTLs under single-locus model for major gene cloning and could not detect directly the multiple alleles using bi-allelic single-nucleotide polymorphism (SNP) marker. As a relatively thorough detection of QTLs with their multiple alleles is required for germplasm population, the restricted two-stage multi-locus multi-allele GWAS (RTM-GWAS) procedure was proposed for identifying the QTL system with varying multiple alleles. From the RTM-GWAS results, a QTL-allele matrix is constructed as a compact form of the population genetic constitution, which can be further used for crop genetic and breeding studies, including major gene mining, population evolution, and breeding by genetic design.

Gamma-glutamyl-leucine levels are causally associated with elevated cardio-metabolic risks.

Gamma-glutamyl dipeptides are bioactive peptides involved in inflammation, oxidative stress, and glucose regulation. Gamma-glutamyl-leucine (Gamma-Glu-Leu) has been extensively reported to be associated with the risk of cardio-metabolic diseases, such as obesity, metabolic syndrome, and type 2 diabetes. However, the causality remains to be uncovered. The aim of this study was to explore the causal-effect relationships between Gamma-Glu-Leu and metabolic risk. In this study, 1,289 subjects were included from a cross-sectional survey on metabolic syndrome (MetS) in eastern China. Serum Gamma-Glu-Leu levels were measured by untargeted metabolomics. Using linear regressions, a two-stage genome-wide association study (GWAS) for Gamma-Glu-Leu was conducted to seek its instrumental single nucleotide polymorphisms (SNPs). One-sample Mendelian randomization (MR) analyses were performed to evaluate the causality between Gamma-Glu-Leu and the metabolic risk. Four SNPs are associated with serum Gamma-Glu-Leu levels, including rs12476238, rs56146133, rs2479714, and rs12229654. Out of them, rs12476238 exhibits the strongest association (Beta = -0.38, S.E. = 0.07 in discovery stage, Beta = -0.29, S.E. = 0.14 in validation stage, combined P-value = 1.04 × 10-8). Each of the four SNPs has a nominal association with at least one metabolic risk factor. Both rs12229654 and rs56146133 are associated with body mass index, waist circumference (WC), the ratio of WC to hip circumference, blood pressure, and triglyceride (5 × 10-5 < P < 0.05). rs56146133 also has nominal associations with fasting insulin, glucose, and insulin resistance index (5 × 10-5 < P < 0.05). Using the four SNPs serving as the instrumental SNPs of Gamma-Glu-Leu, the MR analyses revealed that higher Gamma-Glu-Leu levels are causally associated with elevated risks of multiple cardio-metabolic factors except for high-density lipoprotein cholesterol and low-density lipoprotein cholesterol (P > 0.05). Four SNPs (rs12476238, rs56146133, rs2479714, and rs12229654) may regulate the levels of serum Gamma-Glu-Leu. Higher Gamma-Glu-Leu levels are causally linked to cardio-metabolic risks. Future prospective studies on Gamma-Glu-Leu are required to explain its role in metabolic disorders.