An atlas of protein homo-oligomerization across domains of life.

Protein structures are essential to understanding cellular processes in molecular detail. While advances in artificial intelligence revealed the tertiary structure of proteins at scale, their quaternary structure remains mostly unknown. We devise a scalable strategy based on AlphaFold2 to predict homo-oligomeric assemblies across four proteomes spanning the tree of life. Our results suggest that approximately 45% of an archaeal proteome and a bacterial proteome and 20% of two eukaryotic proteomes form homomers. Our predictions accurately capture protein homo-oligomerization, recapitulate megadalton complexes, and unveil hundreds of homo-oligomer types, including three confirmed experimentally by structure determination. Integrating these datasets with omics information suggests that a majority of known protein complexes are symmetric. Finally, these datasets provide a structural context for interpreting disease mutations and reveal coiled-coil regions as major enablers of quaternary structure evolution in human. Our strategy is applicable to any organism and provides a comprehensive view of homo-oligomerization in proteomes.

Boom-bust population dynamics drive rapid genetic change.

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.

Critical amino acid residues regulating TRPA1 Zn2+ response: A comparative study across species.

Cellular zinc ions (Zn2+) are crucial for signal transduction in various cell types. The transient receptor potential (TRP) ankyrin 1 (TRPA1) channel, known for its sensitivity to intracellular Zn2+ ([Zn2+]i), has been a subject of limited understanding regarding its molecular mechanism. Here, we used metal ion-affinity prediction, three-dimensional structural modeling, and mutagenesis, utilizing data from the Protein Data Bank and AlphaFold database, to elucidate the [Zn2+]i binding domain (IZD) structure composed by specific AAs residues in human (hTRPA1) and chicken TRPA1 (gTRPA1). External Zn2+ induced activation in hTRPA1, while not in gTRPA1. Moreover, external Zn2+ elevated [Zn2+]i specifically in hTRPA1. Notably, both hTRPA1 and gTRPA1 exhibited inherent sensitivity to [Zn2+]i, as evidenced by their activation upon internal Zn2+ application. The critical AAs within IZDs, specifically histidine at 983/984, lysine at 711/717, tyrosine at 714/720, and glutamate at 987/988 in IZD1, and H983/H984, tryptophan at 710/716, E854/E855, and glutamine at 979/980 in IZD2, were identified in hTRPA1/gTRPA1. Furthermore, mutations, such as the substitution of arginine at 919 (R919) to H919, abrogated the response to external Zn2+ in hTRPA1. Among single-nucleotide polymorphisms (SNPs) at Y714 and a triple SNP at R919 in hTRPA1, we revealed that the Zn2+ responses were attenuated in mutants carrying the Y714 and R919 substitution to asparagine and proline, respectively. Overall, this study unveils the intrinsic sensitivity of hTRPA1 and gTRPA1 to [Zn2+]i mediated through IZDs. Furthermore, our findings suggest that specific SNP mutations can alter the responsiveness of hTRPA1 to extracellular and intracellular Zn2+.

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements.

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable. RESULTS: Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq). CONCLUSIONS: This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

Improved two-step testing of genome-wide gene-environment interactions.

Two-step tests for gene-environment ( G × E $G\times E$ ) interactions exploit marginal single-nucleotide polymorphism (SNP) effects to improve the power of a genome-wide interaction scan. They combine a screening step based on marginal effects used to "bin" SNPs for weighted hypothesis testing in the second step to deliver greater power over single-step tests while preserving the genome-wide Type I error. However, the presence of many SNPs with detectable marginal effects on the trait of interest can reduce power by "displacing" true interactions with weaker marginal effects and by adding to the number of tests that need to be corrected for multiple testing. We introduce a new significance-based allocation into bins for Step-2 G × E $G\times E$ testing that overcomes the displacement issue and propose a computationally efficient approach to account for multiple testing within bins. Simulation results demonstrate that these simple improvements can provide substantially greater power than current methods under several scenarios. An application to a multistudy collaboration for understanding colorectal cancer reveals a G × Sex interaction located near the SMAD7 gene.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects.

BACKGROUND: Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database. RESULTS: We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327. CONCLUSIONS: Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

Development of a coding SNP panel for tracking the origin of whole-exome sequencing samples.

Whole-exome sequencing (WES) is widely used to diagnose complex genetic diseases and rare conditions. The implementation of a robust and effective quality control system for sample identification and tracking throughout the WES process is essential. We established a multiplex panel that included 22 coding single-nucleotide polymorphism (cSNP) loci. The personal identification and paternity identification abilities of the panel were evaluated, and a preliminary validation of the practical feasibility of the panel was conducted in a clinical WES case. These results indicate that the cSNP panel could be a useful tool for sample tracking in WES.

Genetic variants affecting mitochondrial function provide further insights for kidney disease.

BACKGROUND: Chronic kidney disease (CKD) is a complex disorder that has become a high prevalence global health problem, with diabetes being its predominant pathophysiologic driver. Autosomal genetic variation only explains some of the predisposition to kidney disease. Variations in the mitochondrial genome (mtDNA) and nuclear-encoded mitochondrial genes (NEMG) are implicated in susceptibility to kidney disease and CKD progression, but they have not been thoroughly explored. Our aim was to investigate the association of variation in both mtDNA and NEMG with CKD (and related traits), with a particular focus on diabetes. METHODS: We used the UK Biobank (UKB) and UK-ROI, an independent collection of individuals with type 1 diabetes mellitus (T1DM) patients. RESULTS: Fourteen mitochondrial variants were associated with estimated glomerular filtration rate (eGFR) in UKB. Mitochondrial variants and haplogroups U, H and J were associated with eGFR and serum variables. Mitochondrial haplogroup H was associated with all the serum variables regardless of the presence of diabetes. Mitochondrial haplogroup X was associated with end-stage kidney disease (ESKD) in UKB. We confirmed the influence of several known NEMG on kidney disease and function and found novel associations for SLC39A13, CFL1, ACP2 or ATP5G1 with serum variables and kidney damage, and for SLC4A1, NUP210 and MYH14 with ESKD. The G allele of TBC1D32-rs113987180 was associated with higher risk of ESKD in patients with diabetes (OR:9.879; CI95%:4.440-21.980; P = 2.0E-08). In UK-ROI, AGXT2-rs71615838 and SURF1-rs183853102 were associated with diabetic nephropathies, and TFB1M-rs869120 with eGFR. CONCLUSIONS: We identified novel variants both in mtDNA and NEMG which may explain some of the missing heritability for CKD and kidney phenotypes. We confirmed the role of MT-ND5 and mitochondrial haplogroup H on renal disease (serum variables), and identified the MT-ND5-rs41535848G variant, along with mitochondrial haplogroup X, associated with higher risk of ESKD. Despite most of the associations were independent of diabetes, we also showed potential roles for NEMG in T1DM.

IKZF3 polymorphisms contribute to the increased risk of acute lymphoblastic leukemia in children.

BACKGROUND: Acute lymphoblastic leukemia (ALL) is the most common cancer in children. IKZF3 (IKAROS family zinc finger 3) is a hematopoietic-specific transcription factor, and it has been validated that it is involved in leukemia. However, the role of IKZF3 single-nucleotide polymorphisms (SNPs) remains unclear. In this case-control study, the authors investigated the association of IKZF3 SNPs with ALL in children. METHODS: Six IKZF3 reference SNPs (rs9635726, rs2060941, rs907092, rs12946510, rs1453559, and rs62066988) were genotyped in 692 patients who had ALL (cases) and in 926 controls. The associations between IKZF3 polymorphisms and ALL risk were determined using odds ratios (ORs) and 95% confidence intervals (CIs). The associations of rs9635726 and rs2060941 with the risk of ALL were further estimated by using false-positive report probability (FPRP) analysis. Functional analysis in silico was performed to evaluate the probability that rs9635726 and rs2060941 might influence the regulation of IKZF3. RESULTS: The authors observed that rs9635726C>T (adjusted OR, 1.49; 95% CI, 1.06-2.11; p = .023) and rs2060941G>T (adjusted OR, 1.51; 95% CI, 1.24-1.84; p = .001) were related to and increased risk of ALL in the recessive and dominant models, respectively. Furthermore, the associations of both rs9635726 (FPRP = .177) and rs2060941 (FPRP < .001) with ALL were noteworthy in the FPRP analysis. Functional analysis indicated that rs9635726 and rs2060941 might repress the transcription of IKZF3 by disrupting its binding to MLLT1, TAF1, POLR2A, and/or RAD21. CONCLUSIONS: This study revealed that IKZF3 polymorphisms were associated with increased ALL susceptibility in children and might influence the expression of IKZF3 by disrupting its binding to MLLT1, TAF1, POLR2A, and/or RAD21. IKZF3 polymorphisms were suggested as a biomarker for childhood ALL.

Causal effect of interleukin (IL)-6 on blood pressure and hypertension: A mendelian randomization study.

To examine whether circulating interleukin-6 (IL-6) levels (CirIL6) have a causal effect on blood pressure using Mendelian randomization (MR) methods. We used data from genome-wide association studies (GWAS) of European ancestry to obtain genetic instruments for circulating IL-6 levels and blood pressure measurements. We applied several robust MR methods to estimate the causal effects and to test for heterogeneity and pleiotropy. We found that circulating IL-6 had a significant positive causal effect on systolic blood pressure (SBP) and pulmonary arterial hypertension (PAH), but not on diastolic blood pressure (DBP) or hypertension. We found that as CirIL6 genetically increased, SBP increased using Inverse Variance Weighted (IVW) method (for ukb-b-20175, ß = 0.082 with SE = 0.032, P = 0.011; for ukb-a-360, ß = 0.075 with SE = 0.031, P = 0.014) and weighted median (WM) method (for ukb-b-20175, ß = 0.061 with SE = 0.022, P = 0.006; for ukb-a-360, ß = 0.065 with SE = 0.027, P = 0.014). Moreover, CirIL6 may be associated with an increased risk of PAH using WM method (odds ratio (OR) = 15.503, 95% CI, 1.025-234.525, P = 0.048), but not with IVW method. Our study provides novel evidence that circulating IL-6 has a causal role in the development of SBP and PAH, but not DBP or hypertension. These findings suggest that IL-6 may be a potential therapeutic target for preventing or treating cardiovascular diseases and metabolic disorders. However, more studies are needed to confirm the causal effects of IL-6 on blood pressure and to elucidate the underlying mechanisms and pathways.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study.

BACKGROUND: Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype. RESULTS: Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers. CONCLUSIONS: Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Potentially functional variants of CHMP4A and PANX1 in the pyroptosis-related pathway predict survival of patients with non-oropharyngeal head and neck squamous cell carcinoma.

BACKGROUND: Pyroptosis has been implicated in the advancement of various cancers. Triggering pyroptosis within tumors amplifies the immune response, thereby fostering an antitumor immune environment. Nonetheless, few published studies have evaluated associations between functional variants in the pyroptosis-related genes and clinical outcomes of patients with non-oropharyngeal head and neck squamous cell carcinoma (NON-ORO HNSCC). METHODS: We conducted an association study of 985 NON-ORO HNSCC patients who were randomly divided into two groups: the discovery group of 492 patients and the replication group of 493 patients. We used Cox proportional hazards regression analysis to examine associations between genetic variants of the pyroptosis-related genes and survival of patients with NON-ORO HNSCC. Bayesian false discovery probability (BFDP) was used for multiple testing correction. Functional annotation was applied to the identified survival-associated genetic variants. RESULTS: There are 8254 single-nucleotide polymorphisms (SNPs) located in 82 pyroptosis-related genes, of which 202 SNPs passed multiple testing correction with BFDP < 0.8 in the discovery and six SNPs retained statistically significant in the replication. In subsequent stepwise multivariable Cox regression analysis, two independent SNPs (CHMP4A rs1997996 G > A and PANX1 rs56175344 C > G) remained significant with an adjusted hazard ratios (HR) of 1.31 (95% confidence interval [CI] = 1.09-1.57, p = 0.004) and 0.65 (95% CI = 0.51-0.83, p = 0.0005) for overall survival (OS), respectively. Further analysis of the combined genotypes revealed progressively worse OS associated with the number of unfavorable genotypes (ptrend < 0.0001 and 0.021 for OS and disease-specific survival, respectively). Moreover, both PANX1 rs56175344G and CHMP4A rs1997996A alleles were correlated with reduced mRNA expression levels. CONCLUSIONS: Genetic variants in the pyroptosis pathway genes may predict the survival of NON-ORO HNSCC patients, likely by reducing the gene expression, but our findings need to be replicated by larger studies.

Functional genetic variants of GEN1 predict overall survival of Chinese epithelial ovarian cancer patients.

BACKGROUND: Inherited variations in DNA double-strand break (DSB) repair pathway are known to influence ovarian cancer occurrence, progression and treatment response. Despite its significance, survival-associated genetic variants within the DSB pathway remain underexplored. METHODS: In the present study, we performed a two-phase analysis of 19,290 single-nucleotide polymorphisms (SNPs) in 199 genes in the DSB repair pathway from a genome-wide association study (GWAS) dataset and explored their associations with overall survival (OS) in 1039 Han Chinese epithelial ovarian carcinoma (EOC) patients. After utilizing multivariate Cox regression analysis with bayesian false-discovery probability for multiple test correction, significant genetic variations were identified and subsequently underwent functional prediction and validation. RESULTS: We discovered a significant association between poor overall survival and the functional variant GEN1 rs56070363 C > T (CT + TT vs. TT, adjusted hazard ratio (HR) = 2.50, P < 0.001). And the impact of GEN1 rs56070363 C > T on survival was attributed to its reduced binding affinity to hsa-miR-1287-5p and the resultant upregulation of GEN1 mRNA expression. Overexpression of GEN1 aggregated EOC cell proliferation, invasion and migration presumably by influencing the expression of immune inhibitory factors, thereby elevating the proportion of polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) and then constructing an immunosuppressive tumor microenvironment. CONCLUSIONS: In conclusion, GEN1 rs56070363 variant could serve as a potential predictive biomarker and chemotherapeutic target for improving the survival of EOC patients.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Phylogenomics of American pika (Ochotona princeps) lineage diversification.

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Association between CBL gene polymorphism and susceptibility of microscopic polyangiitis in a Chinese population: A case-control analysis.

OBJECTIVE: To assess whether Casitas B-lineage lymphoma (CBL) gene polymorphism influences the risk of microscopic polyangiitis (MPA) in Chinese populations. METHODS: In total, 266 MPA patients and 297 healthy controls were recruited for a case-control study. Five CBL SNPs were genotyped using multiplex polymerase chain reaction and high-throughput sequencing. The relationship between SNPs and the risk of MPA under different genetic models was evaluated by SNPstats. SNP-SNP interaction was analyzed by generalized multifactor dimensionality reduction (GMDR). Finally, the association between CBL SNPs and treatment effects were assessed. RESULTS: The results showed that CBL rs2276083 was associated with decreasing MPA risk under dominant (OR: 0.53; p = 0.014) and recessive models (OR: 0.52; p = 0.0034). Stratification analysis indicated that rs2276083 and rs2509671 in age < 60 years, rs2276083 in female or in Han population were protective factors for MPA. The CBL haplotype (A-A-G-C-T) was associated with an increased risk of MPA. GMDR suggested that CBL rs2276083, phosphatidylinositol-4, 5-bisphosphate 3-kinase catalytic subunit alpha (PI3KCA) rs1607237, and autophagy-related gene 7 (ATG7) rs7549008 might interact with each other in MPA development (p = 0.0107). CBL rs1047417 with AG genotype and rs11217234 with AG genotype had better clinical treatment effects than other two genotypes (p = 0.048 and p = 0.025, respectively). CONCLUSION: The genetic polymorphism of CBL had a potential association with the risk of MPA and clinical treatment effects in Guangxi population in China.

Circulating cytokines levels and osteoarthritis: A Mendelian randomization study.

BACKGROUND: Previous traditional observational studies have suggested the contribution of several cytokines and growth factors to the development of osteoarthritis (OA). This study aimed to determine the association of circulating cytokine and growth factor levels with OA. METHODS: We used two-sample Mendelian randomization (MR) to explore the causality between circulating cytokine and growth factor levels and OA [including knee or hip OA (K/HOA), knee OA (KOA), and hip OA (HOA)]. Summary level data for circulating cytokine and growth factor levels were sourced from a genome-wide association study (GWAS) involving 8,293 participants of Finnish ancestry. Single-nucleotide polymorphisms related to K/HOA (39,427 cases and 378,169 controls), KOA (24,955 cases and 378,169 controls), and HOA (15,704 cases and 378,169 controls) were obtained from a previous GWAS. The inverse variance weighted (IVW) method was primarily used for our MR analysis. For exposures to only one relevant SNP as IV, we used the Wald ratio as the major method to assess causal effects. We also conducted a series of sensitivity analyses to improve the robustness of the results. RESULTS: Circulating vascular endothelial growth factor levels were suggestively associated with an increased risk of K/HOA (odds ratio (OR) = 1.034; 95 % confidence interval (CI) = 1.013-1.055; P = 0.001), KOA (OR = 1.034; 95 % CI = 1.014-1.065; P = 0.002), and HOA (OR = 1.039; 95 % CI = 1.003-1.067; P = 0.034). Circulating interleukin (IL)-12p70 levels was suggestively associated with K/HOA (OR = 1.047; 95 % CI = 1.018-1.077; P = 0.001), KOA (OR = 1.058; 95 % CI = 1.022-1.095; P = 0.001), and HOA (OR = 1.044; 95 % CI = 1.000-1.091; P = 0.048). Circulating IL-18 levels were suggestively associated with HOA (OR = 1.068; 95 % CI = 1.014-1.125; P = 0.012). However, limited evidence exists to support causal genetic relationships between other circulating cytokines, growth factor levels and K/HOA, KOA, and HOA. CONCLUSIONS: Our MR analysis provides suggestive evidence of causal relationships between circulating cytokines and growth factors levels and OA, providing new insights into the etiology of OA.

Effects of CASZ1, WNT2B and PTPRG SNPs on stroke susceptibility in the Chinese Han population.

BACKGROUND: Stroke is an important cause of death and disability worldwide, ranking second in the cause of death, and it is thought to be related to genetic factors. The purpose of our study is to investigate the association between CASZ1, WNT2B and PTPRG single nucleotide polymorphisms (SNPs) and stroke risk in the Chinese population. METHODS: We recruited 1418 volunteers, comprised of 710 stroke cases and 708 controls in this study. We used MassARRAY iPLEX GOLD method to genotype the three SNPs on CASZ1, WNT2B and PTPRG. Logistic regression was used to analyse the association between these SNPs and stroke, and odds ratios (ORs) and 95% confidence intervals (CIs) were then calculated. What's more, the interactions among SNPs were predicted by multi-factor dimensionality reduction (MDR) analysis. RESULTS: This research demonstrated that CASZ1 rs880315 and PTPRG rs704341 were associated with reduced stroke susceptibility. More precisely, CASZ1 rs880315 was associated with reduced stroke susceptibility in people aged ≤64 years and women. PTPRG rs704341 was associated with reduced stroke susceptibility in people aged >64 years, women, non-smokers and non-drinkers. Conversely, WNT2B rs12037987 was related to elevated stroke susceptibility in people aged >64 years, women and non-smokers. In addition, CASZ1 rs880315, WNT2B rs12037987 and PTPRG rs704341 had a strong redundancy relationship. CONCLUSION: Our study concludes that CASZ1 rs880315, WNT2B rs12037987 and PTPRG rs704341 are associated with stroke, and the study provides a basis for assessing genetic variants associated with stroke risk in the Han Chinese population.

Association of transketolase polymorphisms with diabetic polyneuropathy in the general population: The KORA F4 study.

AIMS: We recently reported that genetic variability in the TKT gene encoding transketolase, a key enzyme in the pentose phosphate pathway, is associated with measures of diabetic sensorimotor polyneuropathy (DSPN) in recent-onset diabetes. Here, we aimed to substantiate these findings in a population-based KORA F4 study. MATERIALS AND METHODS: In this cross-sectional study, we assessed seven single nucleotide polymorphisms (SNPs) in the transketolase gene in 952 participants from the KORA F4 study with normal glucose tolerance (NGT; n = 394), prediabetes (n = 411), and type 2 diabetes (n = 147). DSPN was defined by the examination part of the Michigan Neuropathy Screening Instrument (MNSI) using the original MNSI > 2 cut-off and two alternative versions extended by touch/pressure perception (TPP) (MNSI > 3) and by TPP plus cold perception (MNSI > 4). RESULTS: After adjustment for sex, age, BMI, and HbA1c, in type 2 diabetes participants, four out of seven transketolase SNPs were associated with DSPN for all three MNSI versions (all p ≤ 0.004). The odds ratios of these associations increased with extending the MNSI score, for example, OR (95% CI) for SNP rs62255988 with MNSI > 2: 1.99 (1.16-3.41), MNSI > 3: 2.27 (1.26-4.09), and MNSI > 4: 4.78 (2.22-10.26); SNP rs9284890 with MNSI > 2: 2.43 (1.42-4.16), MNSI > 3: 3.46 (1.82-6.59), and MNSI > 4: 4.75 (2.15-10.51). In contrast, no associations were found between transketolase SNPs and the three MNSI versions in the NGT and prediabetes groups. CONCLUSIONS: The link of genetic variation in transketolase enzyme to diabetic polyneuropathy corroborated at the population level strengthens the concept suggesting an important role of pathways metabolising glycolytic intermediates in the evolution of diabetic polyneuropathy.