Association Analysis Research Articles

The joint impact of PM2.5 constituents on the risk of cerebrovascular diseases hospitalization: A large community-based cohort study

Air pollution poses significant health risks to urban areas, with limited focus on the chronic association of PM2.5 and its constituents on cerebrovascular diseases (CERs), especially regarding the joint associations. This study explores the individual and joint associations between PM2.5 constituents and CER hospitalization risks through a cohort analysis of 36,271 adults in the Pearl River Delta, South China, from 2015 to 2020. Cox proportional hazards regression and quantile-based g-computation models were used to quantify the individual and joint associations of annual mean concentrations of PM2.5 constituents with hospitalization for CERs. 1,151 participants were hospitalized due to CERs during the five-year follow-up period. Joint associations analyses identified that one quartile increase in co-exposure may result in hazard ratios of 1.530 (1.441–1.623), 1.840 (1.710–1.980), and 1.609 (1.491–1.737) for CERs, total, and ischemic stroke hospitalization, respectively. The adverse effect was primarily driven by organic matter and chlorine. Men, those with a history of tobacco or alcohol use or with low residential greenness, were more susceptible to CERs hospitalization following PM2.5 constituents co-exposure. Upcoming strategies should focus on monitoring and regulating PM2.5 constituents, encouraging healthy lifestyles, and enhancing urban greenery.

Insurance Type and Withdrawal of Life-Sustaining Therapy in Critically Injured Trauma Patients

Withdrawal of life-sustaining therapy (WLST) decisions for critically injured trauma patients are complicated and multifactorial, with potential for patients' insurance status to affect decision-making. To determine if patient insurance type (private insurance, Medicaid, and uninsured) is associated with time to WLST in critically injured adults cared for at US trauma centers. This retrospective registry-based cohort study included reported data from level I and level II trauma centers in the US that participated in the American College of Surgeons Trauma Quality Improvement Program (TQIP) registry. Participants included adult trauma patients who were injured between January 1, 2017, and December 31, 2020, and required an intensive care unit stay. Patients were excluded if they died on arrival or in the emergency department or had a preexisting do not resuscitate directive. Analyses were performed on December 12, 2023. Insurance type (private insurance, Medicaid, uninsured). An adjusted time-to-event analysis for association between insurance status and time to WLST was performed, with analyses accounting for clustering by hospital. This study included 307 731 patients, of whom 160 809 (52.3%) had private insurance, 88 233 (28.6%) had Medicaid, and 58 689 (19.1%) were uninsured. The mean (SD) age was 40.2 (14.1) years, 232 994 (75.7%) were male, 59 551 (19.4%) were African American or Black patients, and 201 012 (65.3%) were White patients. In total, 12 962 patients (4.2%) underwent WLST during their admission. Patients who are uninsured were significantly more likely to undergo earlier WLST compared with those with private insurance (HR, 1.54; 95% CI, 1.46-1.62) and Medicaid (HR, 1.47; 95% CI, 1.39-1.55). This finding was robust to sensitivity analysis excluding patients who died within 48 hours of presentation and after accounting for nonwithdrawal death as a competing risk. In this cohort study of US adult trauma patients who were critically injured, patients who were uninsured underwent earlier WLST compared with those with private or Medicaid insurance. Based on our findings, patient's ability to pay was may be associated with a shift in decision-making for WLST, suggesting the influence of socioeconomics on patient outcomes.

Bacteriophage-driven emergence and expansion of Staphylococcus aureus in rodent populations.

Human activities such as agriculturalization and domestication have led to the emergence of many new pathogens via host-switching events between humans, domesticated and wild animals. Staphylococcus aureus is a multi-host opportunistic pathogen with a global healthcare and economic burden. Recently, it was discovered that laboratory and wild rodents can be colonised and infected with S. aureus, but the origins and zoonotic potential of rodent S. aureus is unknown. In order to trace their evolutionary history, we employed a dataset of 1249 S. aureus genome sequences including 393 of isolates from rodents and other small mammals (including newly determined sequences for 305 isolates from 7 countries). Among laboratory mouse populations, we identified multiple widespread rodent-specific S. aureus clones that likely originated in humans. Phylogeographic analysis of the most common murine lineage CC88 suggests that it emerged in the 1980s in laboratory mouse facilities most likely in North America, from where it spread to institutions around the world, via the distribution of mice for research. In contrast, wild rodents (mice, voles, squirrels) were colonized with a unique complement of S. aureus lineages that are widely disseminated across Europe. In order to investigate the molecular basis for S. aureus adaptation to rodent hosts, genome-wide association analysis was carried out revealing a unique complement of bacteriophages associated with a rodent host ecology. Of note, we identified novel prophages and pathogenicity islands in rodent-derived S. aureus that conferred the potential for coagulation of rodent plasma, a key phenotype of abscess formation and persistence. Our findings highlight the remarkable capacity of S. aureus to expand into new host populations, driven by the acquisition of genes promoting survival in new host-species.

Genetic characterization of Tibetan pigs adapted to high altitude under natural selection based on a large whole-genome dataset

The Qinghai-Tibet Plateau is a valuable genetic resource pool, and the high-altitude adaptation of Tibetan pigs is a classic example of the adaptive evolution of domestic animals. Here, we report the presence of Darwinian positive selection signatures in Tibetan pigs (TBPs) using 348 genome-wide datasets (127 whole-genome sequence datasets (WGSs) and 221 whole-genome single-nucleotide polymorphism (SNP) chip datasets). We characterized a high-confidence list of genetic signatures related response to high-altitude adaptation in Tibetan pigs, including 4,598 candidate SNPs and 131 candidate genes. Functional annotation and enrichment analysis revealed that 131 candidate genes are related to multiple systems and organs in Tibetan pigs. Notably, eight of the top ten novel genes, RALB, NBEA, LIFR, CLEC17A, PRIM2, CDH7, GK5 and FAM83B, were highlighted and associated with improved adaptive heart functions in Tibetan pigs high-altitude adaptation. Moreover, genome-wide association analysis revealed that 29 SNPs were involved in 13 candidate genes associated with at least one adaptive trait. In particular, among the top ten candidate genes, CLEC17A is related to a reduction in hemoglobin (HGB) in Tibetan pigs. Overall, our study provides a robust SNP/gene list involving genetic adaptation for Tibetan pig high-altitude adaptation, and it will be a valuable resource for future Tibetan pig studies.

Association of the IceCube neutrinos with blazars in the CGRaBS sample

The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically motivated candidate populations, such as blazars -- a subclass of active galactic nuclei with jets pointed toward our line of sight -- they have not been convincingly linked to HE neutrino production yet. Here, we perform a spatio-temporal association analysis between a sample of blazars (from the CGRaBS catalog) in the radio and optical bands and the most up-to-date IceCube HE neutrino catalog. We find that if the IceCube error regions are enlarged by 1$^ in quadrature, to account for unknown systematic errors at a maximal level, a spatio-temporal correlation between the multiwavelength light curves of the CGRaBS blazars and the IceCube HE neutrinos is hinted at, least at a 2.17sigma significance level. On the other hand, when the IceCube error regions are taken as their published values, we do not find any significant correlations. A discrepancy in the blazar-neutrino correlation strengths, when using such minimal and enlarged error region scenarios, was also obtained in a recent study by the IceCube collaboration. In our study, this difference arises because several flaring blazars -- coinciding with a neutrino arrival time -- happen to narrowly miss the published 90<!PCT!>-likelihood error region of the nearest neutrino event. For all of the associations driving our most significant correlations, the flaring blazar is much less than 1$^ away from the published error regions. Therefore, our results indicate that the question of the blazar-neutrino connection is highly sensitive to the reconstruction of the neutrino error regions, whose reliability is expected to improve with the next generation of neutrino observatories.

SLC45A4 encodes a mitochondrial putrescine transporter that promotes GABA de novo synthesis.

Solute carriers (SLC) are membrane proteins that facilitate the transportation of ions and metabolites across either the plasma membrane or the membrane of intracellular organelles. With more than 450 human genes annotated as SLCs, many of them are still orphan transporters without known biochemical functions. We developed a metabolomic-transcriptomic association analysis, and we found that the expression of SLC45A4 has a strong positive correlation with the cellular level of γ-aminobutyric acid (GABA). Using mass spectrometry and the stable isotope tracing approach, we demonstrated that SLC45A4 promotes GABA de novo synthesis through the Arginine/Ornithine/Putrescine (AOP) pathway. SLC45A4 functions as a putrescine transporter localized to the mitochondrial membrane to facilitate GABA production. Taken together, our results revealed a new biochemical mechanism where SLC45A4 controls GABA production.

Assessing the carotenoid profiles and allelic diversity of yellow maize inbred lines adapted to mid-altitude subhumid maize agroecology in Ethiopia.

Biofortification of provitamin A in maize is an attractive and sustainable remedy to the problem of vitamin A deficiency in developing countries. The utilization of molecular markers represents a promising avenue to facilitate the development of provitamin A (PVA)-enriched maize varieties. We screened 752 diverse tropical yellow/orange maize lines using kompetitive allele-specific PCR (KASP) makers to validate the use of KASP markers in PVA maize breeding. To this end, a total of 161 yellow/orange inbred lines, selected from among the 752 lines, were evaluated for their endosperm PVA and other carotenoid compounds levels in two separate trials composed of 63 and 98 inbred lines in 2020 and 2021, respectively. Significant differences (p < 0.001) were observed among the yellow maize inbred lines studied for all carotenoid profiles. An inbred line TZMI1017, introduced by the International Institute of Tropical Agriculture (IITA) showed the highest level of PVA (12.99 µg/g) and β-carotene (12.08 µg/g). The molecular screening showed 43 yellow maize inbred lines carrying at least three of the favorable alleles of the KASP markers. TZMI1017 inbred line also carried the favorable alleles of almost all markers. In addition, nine locally developed inbred lines had medium to high PVA concentrations varying from 5.11 µg/g to 10.76 µg/g and harbored the favorable alleles of all the KASP PVA markers. Association analysis between molecular markers and PVA content variation in the yellow/orange maize inbred lines did not reveal a significant, predictable correlation. Further investigation is warranted to elucidate the underlying genetic architecture of the PVA content in this germplasm. However, we recommend strategic utilization of the maize-inbred lines with higher PVA content to enhance the PVA profile of the breeding program's germplasm.

Association analysis of leaf aromatic substances in cultivated and weedy types of Perilla crop using SSR markers

In East Asia, particularly South Korea, the two cultivated varieties of Perilla are commonly grown. They are clearly distinguished by their aromatic substances and have different uses as leafy vegetables or oil crop. This study was performed for the development of simple sequence repeat (SSR) markers linked to volatile compounds in Perilla leaves that show differences between cultivated var. frutescens (CF), weedy var. frutescens (WF), and weedy var. crispa (WC) of Perilla. Fifty Perilla SSR primer sets were used to analyze genetic diversity for the 80 Perilla accessions of the three types. A total of 276 alleles were detected, with an average of 5.5 alleles per locus. The average genetic diversity values for CF, WF, and WC accessions were 0.402, 0.583, and 0.437, respectively. WF accessions exhibited the highest genetic diversity among the three types of the Perilla crop. Phylogenetic tree analysis classified 80 Perilla accessions of the three types into four groups, showing 37.2 % genetic similarity. Three types of the Perilla crop were clearly distinguished except for outstanding accessions. Through the application of an association analysis involving 50 Perilla SSR primer sets and five volatile compounds (perilla aldehyde, perilla ketone, myristicin, dill apiol, (Z,E)-α-farnesene) in the three types of the Perilla accessions, we detected 11 significant marker-trait associations duplicated in both Q GLM and Q + K MLM methods. These findings serve as valuable insights for identifying the aromatic substances in Perilla plants originating from various regions of South Korea.

Metabolism of quercitrin in the colon and its beneficial regulatory effects on gut microbiota.

Quercitrin is a dietary flavonoid widely found in plants with various physiological activities. However, whether quercitrin alters gut microbiota in vivo is not well understood. The aim of this study was to investigate metabolism of quercitrin in the colon and its regulation on gut microbiota in mice. Herein, 22 flavonoids related to quercitrin metabolism were identified based on ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS). Gas chromatography and 16S rDNA gene sequencing were used to investigate short-chain fatty acid (SCFA) content and diversity of composition of gut microbiota, respectively. The results showed that quercitrin significantly alters the beta-diversity of the gut microbiota, probiotics such as Akkermansia and Lactococcus were significantly increased, and the production of propanoate, isovalerate and hexanoate of the quercitrin group were enhanced significantly. The Spearman's association analysis provided evidence that Gardnerella and Akkermansia have obvious correlations with most of quercitrin metabolites and SCFAs. Quercitrin and its metabolites in the colon altered the structure of the mice gut microbiota and increased the content of SCFAs. Our experiments provide valuable insights into quercitrin research and application. © 2024 Society of Chemical Industry.

Association analysis of indel variants and gene expression identifies MDM4 as a novel locus for skeletal muscle hypertrophy and power athlete status.

Insertions and deletions (indels) are the second most common type of variation in the human genome. However, limited data on their associations with exercise-related phenotypes have been documented. The aim of the present study was to examine the association between 18,370 indel variants and power athlete status, followed by additional studies in 357,246 individuals. In the discovery phase, the D allele of the MDM4 gene rs35493922 I/D polymorphism was over-represented in power athletes compared with control subjects (P=7.8×10-9) and endurance athletes (P=0.0012). These findings were replicated in independent cohorts, showing a higher D allele frequency in power athletes compared with control subjects (P=0.016) and endurance athletes (P=0.031). Furthermore, the D allele was positively associated (P=0.0013) with greater fat-free mass in the UK Biobank. MDM4 encodes a protein that inhibits the activity of p53, which induces muscle fibre atrophy. Accordingly, we found that MDM4 expression was significantly higher in the vastus lateralis of power athletes compared with endurance athletes (P=0.0009) and was positively correlated with the percentage of fast-twitch muscle fibres (P=0.0062) and the relative area occupied by fast-twitch muscle fibres (P=0.0086). The association between MDM4 gene expression and an increased proportion of fast-twitch muscle fibres was confirmed in two additional cohorts. Finally, we found that the MDM4 DD genotype was associated with increased MDM4 gene expression in vastus lateralis and greater cross-sectional area of fast-twitch muscle fibres. In conclusion, MDM4 is suggested to be a potential regulator of muscle fibre specification and size, with its indel variant being associated with power athlete status. HIGHLIGHTS: What is the central question of this study? Which indel variants are functional and associated with sport- and exercise-related traits? What is the main finding and its importance? Out of 18,370 tested indels, the MDM4 gene rs35493922 I/D polymorphism was found to be the functional variant (affecting gene expression) and the most significant, with the deletion allele showing associations with power athlete status, fat-free mass and cross-sectional area of fast-twitch muscle fibres. Furthermore, the expression of MDM4 was positively correlated with the percentage of fast-twitch muscle fibres and the relative area occupied by fast-twitch muscle fibres.

High-depth whole-genome sequencing identifies structure variants, copy number variants and short tandem repeats associated with Parkinson’s disease

While numerous single nucleotide variants and small indels have been identified in Parkinson’s disease (PD), the contribution of structural variants (SVs), copy number variants (CNVs), and short tandem repeats (STRs) remains poorly understood. Here we investigated the association using the high-depth whole-genome sequencing data from 466 Chinese PD patients and 513 controls. Totally, we identified 29,561 SVs, 32,153 CNVs, and 174,905 STRs, and found that CNV deletions were significantly enriched in the end-proportion of autosomal chromosomes in PD. After genome-wide association analysis and replication in an external cohort of 352 cases and 547 controls, we validated that the 1.6 kb-deletion neighboring MUC19, 12.4kb-deletion near RXFP1 and GGGAAA repeats in SLC2A13 were significantly associated with PD. Moreover, the MUC19 deletion and the SLC2A13 5-copy repeat reduced the penetrance of the LRRK2 G2385R variant. Moreover, genes with these variants were dosage-sensitive. These data provided novel insights into the genetic architecture of PD.

Cortical Morphological Networks Differ Between Gyri and Sulci.

This study explored how the human cortical folding pattern composed of convex gyri and concave sulci affected single-subject morphological brain networks, which are becoming an important method for studying the human brain connectome. We found that gyri-gyri networks exhibited higher morphological similarity, lower small-world parameters, and lower long-term test-retest reliability than sulci-sulci networks for cortical thickness- and gyrification index-based networks, while opposite patterns were observed for fractal dimension-based networks. Further behavioral association analysis revealed that gyri-gyri networks and connections between gyral and sulcal regions significantly explained inter-individual variance in Cognition and Motor domains for fractal dimension- and sulcal depth-based networks. Finally, the clinical application showed that only sulci-sulci networks exhibited morphological similarity reductions in major depressive disorder for cortical thickness-, fractal dimension-, and gyrification index-based networks. Taken together, these findings provide novel insights into the constraint of the cortical folding pattern to the network organization of the human brain.

Deciphering the genomic insights into the coexistence of congenital scoliosis and congenital anomalies of the kidney and urinary tract.

Congenital scoliosis and congenital anomalies of the kidney and urinary tract are distinct genetic disorders with differing clinical manifestations. Clinically, their coexistence is not rare, but the etiologies of these complex diseases remain largely unknown, especially their shared genetic basis. We sequenced the genomes of 40 individuals diagnosed with both CS and CAKUT, alongside 2,764 controls from a Chinese Han population cohort. Our analyses encompassed gene-based and pathway-based weighted rare variant association tests, complemented by copy number variant association analyses, aiming to unravel the shared genomic etiology underlying these congenital conditions. Gene-based analysis identified PTPN11 as a pivotal gene influencing both skeletal and urinary system development (P = 1.95E-21), participating in metabolic pathways, especially the MAPK/ERK pathway known to regulate skeletal and urinary system development. Pathway-based enrichment showed a significant signal in the MAPK/ERK pathway (P = 3E-04), reinforcing the potential role of PTPN11 and MAPK/ERK pathway in both conditions. Additionally, CNV analysis pinpointed IGFLR1 haploinsufficiency as a potential influential factor in the combined CS-CAKUT phenotypic spectrum. This study enriches our understanding of the intricate genomic interplay underlying congenital scoliosis and kidney and urinary tract anomalies, emphasizing the shared genetic foundations between these two disorders.

Sex-specific epigenetic signatures of circulating urate and its increase after BCG vaccination.

Background: Urate concentration and the physiological regulation of urate homeostasis exhibit clear sex differences. DNA methylation has been shown to explain a substantial proportion of serum urate variance, mediate the genetic effect on urate concentration, and co-regulate with cardiometabolic traits. However, whether urate concentration is associated with DNA methylation in a sex-dependent manner is unknown. Additionally, it is worth investigating if urate changes after perturbations, such as vaccination, are associated with DNA methylation in a sex-specific manner. Methods: We investigated the association between DNA methylation and serum urate concentrations in a Dutch cohort of 325 healthy individuals. Urate concentration and DNA methylation were measured before and after Bacillus Calmette-Guérin (BCG) vaccination, used as a perturbation associated with increased gout flares. The association analysis included united, interaction, and sex-stratified analysis. Validation of the identified CpG sites was conducted using three independent cohorts. Results: 215 CpG sites were associated with serum urate in males, while 5 CpG sites were associated with serum urate in females, indicating sex-specific associations. Circulating urate concentrations significantly increased after BCG vaccination, and baseline DNA methylation was associated with differences in urate concentration before and after vaccination in a sex-specific manner. The CpG sites associated with urate concentration in males were enriched in neuro-protection pathways, whereas in females, the urate change-associated CpG sites were related to lipid and glucose metabolism. Conclusion: Our study enhances the understanding of how epigenetic factors contribute to regulating serum urate levels in a sex-specific manner. These insights have significant implications for the diagnosis, prevention, and treatment of various urate-related diseases and highlight the importance of personalized and sex-specific approaches in medicine.

Estimating risk of consequences following hypoglycaemia exposure using the Hypo-RESOLVE cohort: a secondary analysis of pooled data from insulin clinical trials.

Whether hypoglycaemia increases the risk of other adverse outcomes in diabetes remains controversial, especially for hypoglycaemia episodes not requiring assistance from another person. An objective of the Hypoglycaemia REdefining SOLutions for better liVEs (Hypo-RESOLVE) project was to create and use a dataset of pooled clinical trials in people with type 1 or type 2 diabetes to examine the association of exposure to all hypoglycaemia episodes across the range of severity with incident event outcomes: death, CVD, neuropathy, kidney disease, retinal disorders and depression. We also examined the change in continuous outcomes that occurred following a hypoglycaemia episode: change in eGFR, HbA1c, blood glucose, blood glucose variability and weight. Data from 84 trials with 39,373 participants were pooled. For event outcomes, time-updated Cox regression models adjusted for age, sex, diabetes duration and HbA1c were fitted to assess association between: (1) outcome and cumulative exposure to hypoglycaemia episodes; and (2) outcomes where an acute effect might be expected (i.e. death, acute CVD, retinal disorders) and any hypoglycaemia exposure within the last 10 days. Exposures to any hypoglycaemia episode and to episodes of given severity (levels 1, 2 and 3) were examined. Further adjustment was then made for a wider set of potential confounders. The within-person change in continuous outcomes was also summarised (median of 40.4 weeks for type 1 diabetes and 26 weeks for type 2 diabetes). Analyses were conducted separately by type of diabetes. The maximally adjusted association analysis for type 1 diabetes found that cumulative exposure to hypoglycaemia episodes of any level was associated with higher risks of neuropathy, kidney disease, retinal disorders and depression, with risk ratios ranging from 1.55 (p=0.002) to 2.81 (p=0.002). Associations of a similar direction were found when level 1 episodes were examined separately but were significant for depression only. For type 2 diabetes cumulative exposure to hypoglycaemia episodes of any level was associated with higher risks of death, acute CVD, kidney disease, retinal disorders and depression, with risk ratios ranging from 2.35 (p<0.0001) to 3.00 (p<0.0001). These associations remained significant when level 1 episodes were examined separately. There was evidence of an association between hypoglycaemia episodes of any kind in the previous 10 days and death, acute CVD and retinal disorders in both type 1 and type 2 diabetes, with rate ratios ranging from 1.32 (p=0.017) to 2.68 (p<0.0001). These associations varied in magnitude and significance when examined separately by hypoglycaemia level. Within the range of hypoglycaemia defined by levels 1, 2 and 3, we could not find any evidence of a threshold at which risk of these consequences suddenly became pronounced. These data are consistent with hypoglycaemia being associated with an increased risk of adverse events across several body systems in diabetes. These associations are not confined to severe hypoglycaemia requiring assistance.

Cognitive rejuvenation in old rats by hippocampal OSKM gene therapy.

Several studies have indicated that interrupted epigenetic reprogramming using Yamanaka transcription factors (OSKM) can rejuvenate cells from old laboratory animals and humans. However, the potential of OSKM-induced rejuvenation in brain tissue has been less explored. Here, we aimed to restore cognitive performance in 25.3-month-old female Sprague-Dawley rats using OSKM gene therapy for 39days. Their progress was then compared with the cognitive performance of untreated 3.5-month-old rats as well as old control rats treated with a placebo adenovector. The Barnes maze test, used to assess cognitive performance, demonstrated enhanced cognitive abilities in old rats treated with OSKM compared to old control animals. In the treated old rats, there was a noticeable trend towards improved spatial memory relative to the old controls. Further, OSKM gene expression did not lead to any pathological alterations within the 39days. Analysis of DNA methylation following OSKM treatment yielded three insights. First, epigenetic clocks for rats suggested a marginally significant epigenetic rejuvenation. Second, chromatin state analysis revealed that OSKM treatment rejuvenated the methylome of the hippocampus. Third, an epigenome-wide association analysis indicated that OSKM expression in the hippocampus of old rats partially reversed the age-related increase in methylation. In summary, the administration of Yamanaka genes via viral vectors rejuvenates the functional capabilities and the epigenetic landscape of the rat hippocampus.

Shared genetic etiology of vessel diseases: A genome-wide multi-traits association analysis

BackgroundThe comorbidity among vascular diseases has been widely reported, however, the contribution of shared genetic components remains ambiguous. MethodsBased on genome-wide association study summary statistics, we employed statistical genetics methodologies to explore the shared genetic basis of eight vascular diseases: coronary artery disease, abdominal aortic aneurysm, ischemic stroke, peripheral artery disease, thoracic aortic aneurysm, phlebitis, varicose veins, and venous thromboembolism. We assessed global and local genetic correlations among these disorders by linkage disequilibrium score regression, high-definition likelihood, and local analysis of variant association. Cross-trait analyses conducted with CPASSOC identified pleiotropic variants and loci. Further, biological pathways at the multi-omics level were explored using multimarker analysis of genomic annotation, transcriptome-wide and proteome-wide association studies. Causal associations among the vascular diseases were evaluated by mendelian randomization and latent causal variable to assess vertical pleiotropic effects. ResultsWe found significant global genetic associations in 18 pairs of vascular diseases. Additionally, we discovered 317 unique genomic regions where at least one pair of traits demonstrated significant correlation. Multi-trait association analysis identified 19,361 significant potential pleiotropic variants in 274 independent pleiotropic loci. Multi-trait colocalization analysis revealed 56 colocalized loci in specific disease sets. Gene-based analysis identified 700 potential pleiotropic genes, which were subsequently validated at both transcriptome and protein levels. Gene-set enrichment analysis supports the role of biological pathways such as vessel wall structure, coagulation and lipid transport in vascular disease. Additionally, 7 pairs of vascular diseases have a causal relationship. ConclusionsOur study indicates a shared genetic basis and the presence of common risk genes among vascular diseases. These findings offer novel insights into potential mechanisms underlying the association between vascular diseases, as well as provide guidance for interventions and treatments of multi-vascular conditions.

Association analysis of miRNA-146a and miRNA-196a polymorphism with breast cancer risk in Khyber Pakhtunkhwa, Pakistan: a preliminary study.

microRNAs (miRNAs) are small, noncoding RNA molecules, functioning either as oncogenes or tumor suppressor genes. SNPs in miRNAs might modify the expression of genes associated with miRNAs, enhancing susceptibility to breast cancer. Both miRNA-146a (rs2910164) and miRNA-196a (rs11614913) are identified and significantly associated with breast cancer risk in several ethnicities, but remains unexplored in Khyber Pakhtunkhwa population of Pakistan. This study was aimed to check the relation of selected SNPs with breast cancer risk. The research cohort included 100 breast cancer patients and 100 healthy controls. All the participants were subjected for DNA extraction followed by T-ARMS PCR and gel electrophoresis. The results revealed a strong association between risk allele (G) of miRNA-146a and increased risk of breast cancer (OR = 2.04, P = 0.0006). Similarly, heterozygous and mutant genotypes also indicated high risk and significant association with breast cancer risk (CG; OR = 0.51, 9 P = 0.0001) (GG; OR = 3.76, P = 0.04). However, risk allele (T) of miRNA-196a (rs11614913) failed to exhibit significant association with breast cancer risk (OR = 0.92 P = 0.68). Similarly, the heterozygous and mutant genotype did not show significant association with breast cancer risk (CT; OR = 0.52, P = 0.125 (TT; OR = 0.88, P = 0.84). Furthermore, miRNA-146a (rs2910164) and miRNA-196a (rs11614913) polymorphisms exhibited non-significant associations with family history (P = 0.34, P = 0.77), PR status (P = 0.310, P = 0.397), ER status (P = 0.992, P = 0.981), nodal status (P = 0.86, P = 0.90), and menstrual status (P = 0.97, P = 0.09). Notably, miRNA-196a showed a significant association with the metastasis group (P = 0.010) and cancer stages (P = 0.047). In conclusion, this study highlights the association of miRNA-146a (rs2910164) polymorphism with breast cancer risk but suggested non-significant association of miRNA-196a (rs11614913) with breast cancer risk. However, these findings need to be confirmed through larger data set for more accurate result.

Characterization of a Salmonella enterica serovar Typhimurium lineage with rough colony morphology and multidrug resistance

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a major cause of salmonellosis, and the emergence of multidrug-resistant pathovariants has become a growing concern. Here, we investigate a distinct rough colony variant exhibiting a strong biofilm-forming ability isolated in China. Whole-genome sequencing on 2,212 Chinese isolates and 1,739 publicly available genomes reveals the population structure and evolutionary history of the rough colony variants. Characterized by macro, red, dry, and rough (mrdar) colonies, these variants demonstrate enhanced biofilm formation at 28 °C and 37 °C compared to typical rdar colonies. The mrdar variants exhibit extensive multidrug resistance, with significantly higher resistance to at least five classes of antimicrobial agents compared to non-mrdar variants. This resistance is primarily conferred by an IncHI2 plasmid harboring 19 antimicrobial resistance genes. Phylogenomic analysis divides the global collections into six lineages. The majority of mrdar variants belong to sublineage L6.5, which originated from Chinese smooth colony strains and possibly emerged circa 1977. Among the mrdar variants, upregulation of the csgDEFG operons is observed, probably due to a distinct point mutation (−44G > T) in the csgD gene promoter. Pangenome and genome-wide association analyses identify 87 specific accessory genes and 72 distinct single nucleotide polymorphisms associated with the mrdar morphotype.

Genome-wide identification of the PP1 and PPP1R3 gene family in oyster Crassostrea gigas: Unraveling their roles in glycogen metabolism

High glycogen in oysters contributes to their flavor, quality, and hardiness, making it a crucial target trait for breeding purposes. Until now, the upstream regulatory mechanism for glycogen synthesis has remained unclear. Protein phosphatase-1 regulatory subunit 3 (PPP1R3) is a member of a family that directs the protein phosphatase-1 (PP1) to glycogen, resulting in the activation of glycogen synthase (GS) and the inactivation of glycogen phosphorylase (GP). However, their functions in molluscan species have remained unclear until now. In our study, we verified the function of the PP1α/PPP1R3D pathways in the glycogen metabolism process in Crassostrea gigas. Initially, we conducted a phylogenetic analysis and identified two PP1 and PPP1R3 members in molluscan species, emphasizing copy number variations in eukaryotic animals. Additionally, we identified distinct substitutions of conserved motifs specific to bivalves, Cephalopoda, or Gastropods, revealing functional diversity in PP1 and PPP1R3. Our primary focus was on PP1α/PPP1R3D in the oyster C.gigas, detected in the muscle glycogen fraction of bivalve species, with mRNA expression showing a positive correlation with glycogen content. By manipulating CgPPP1R3D expression in oysters through adenovirus-mediated overexpression or siRNA-mediated downregulation, we established a direct positive correlation between CgPPP1R3D expression levels and glycogen content. Colocation analysis confirmed protein interactions between CgPP1α and CgPPP1R3D. Furthermore, to elucidate the correlation between gene expression and glycogen accumulation, we conducted an association analysis and identified a cluster of significant SNP loci upstream of CgPPP1R3D significantly related to glycogen content. The high polymorphism observed may be linked to differential expressions of CgPPP1R3D in oyster individuals with high and low glycogen levels. These findings serve as a crucial reference for enhancing glycogen traits in molluscan species.