[Educational attainment increases the risk of developing allergic rhinitis and (or) eczema: a Mendelian randomization study based on genome-wide association studies in European patients].

Objective: To determine the causal relationship between educational attainment and the risk of allergic rhinitis and (or) eczema using Mendelian randomization (MR) analyses. Methods: This study was a secondary data analysis based on the summary data of genome-wide association studies (GWAS), which involved 293 723 participants (educational attainment) from the Social Science Genetics Association Consortium and 462 013 participants [allergic rhinitis and (or) eczema] from the UK Biobank. Genetic variants that were closely related to educational attainment were identified as instrumental variables. Two-sample MR analyses, including inverse-variance weighted (IVW), MR-Egger regression, weighted median method and weighted model-based estimation, were performed to investigate the causal relationship between educational attainment and the risk of allergic rhinitis and (or) eczema, in which the odds ratio (OR) values were used as indicators. Results: A total of 70 single-nucleotide polymorphisms (SNPs) were chosen as instrumental variables. The MR-Egger regression results suggested that the genetic pleiotropy was unlikely to bias our results (P=0.107). In the univariable MR analyses, IVW regression showed that the risk of allergic rhinitis and (or) eczema was OR=1.044 (95%CI: 1.020-1.069, P<0.001) and OR=1.170 (95%CI: 1.074-1.256, P<0.001), respectively, for the increase in the duration of education by one year or one standard deviation (SD) (3.71 years). In the reverse MR analysis, IVW regression showed little evidence that allergic rhinitis and (or) eczema affected educational attainment (OR=1.020, 95%CI: 0.927-1.023, P=0.683). The results of the weighted median method and weighted mode-based estimation were consistent with the results of IVW. Conclusion: This study suggests that there is a positive causal relationship between educational attainment and the risk of allergic rhinitis and (or) eczema, which means that educational attainment can increase the occurrence of allergic rhinitis and (or) eczema.

Is Type 2 Diabetes Mellitus Associated with Spinal Degenerative Disorders?: Evidence from Observational and 2-Sample Mendelian Randomization Analyses.

BACKGROUND: Type 2 diabetes mellitus (T2DM) and spinal degenerative disorders (SDD) are common diseases that frequently coexist. However, both traditional observational studies and recent Mendelian randomization (MR) studies have demonstrated conflicting evidence on the association between T2DM and SDD. This comparative study explored and compared the association between T2DM and SDD using observational and MR analyses. METHODS: For observational analyses, cross-sectional studies (44,972 participants with T2DM and 403,095 participants without T2DM), case-control studies (38,234 participants with SDD and 409,833 participants without SDD), and prospective studies (35,550 participants with T2DM and 392,046 participants without T2DM with follow-up information until 2022) were performed to test the relationship between T2DM and SDD using individual-level data from the U.K. Biobank from 2006 to 2022. For MR analyses, the associations between single-nucleotide polymorphisms with SDD susceptibility obtained using participant data from the U.K. Biobank, which had 407,938 participants from 2006 to 2022, and the FinnGen Consortium, which had 227,388 participants from 2017 to 2022, and genetic predisposition to T2DM obtained using summary statistics from a pooled genome-wide association study involving 1,407,282 individuals were examined. The onset and severity of T2DM are not available in the databases being used. RESULTS: Participants with T2DM were more likely to have SDD than their counterparts. Logistic regression analysis identified T2DM as an independent risk factor for SDD, which was confirmed by the Cox proportional hazard model results. However, using single-nucleotide polymorphisms as instruments, the MR analyses demonstrated no causal relationship between T2DM and SDD. The lack of such an association was robust in the sensitivity analysis, and no pleiotropy was seen. CONCLUSIONS: Our results suggest that the association between T2DM and SDD may be method-dependent. Researchers and clinicians should be cautious in interpreting the association, especially the causal association, between T2DM and SDD. Our findings provide fresh insights into the association between T2DM and SDD by various analysis methods and guide future research and clinical efforts in the effective prevention and management of T2DM and SDD. LEVEL OF EVIDENCE: Prognostic Level III. See Instructions for Authors for a complete description of levels of evidence.

Exploring the relationship between life course adiposity and sepsis: insights from a two-sample Mendelian randomization analysis.

Objectives: The relationship between adiposity and sepsis has received increasing attention. This study aims to explore the causal relationship between life course adiposity and the sepsis incidence. Methods: Mendelian randomization (MR) method was employed in this study. Instrumental variants were obtained from genome-wide association studies for life course adiposity, including birth weight, childhood body mass index (BMI), childhood obesity, adult BMI, waist circumference, visceral adiposity, and body fat percentage. A meta-analysis of genome-wide association studies for sepsis including 10,154 cases and 454,764 controls was used in this study. MR analyses were performed using inverse variance weighted, MR Egger regression, weighted median, weighted mode, and simple mode. Instrumental variables were identified as significant single nucleotide polymorphisms at the genome-wide significance level (P < 5×10-8). The sensitivity analysis was conducted to assess the reliability of the MR estimates. Results: Analysis using the MR analysis of inverse variance weighted method revealed that genetic predisposition to increased childhood BMI (OR = 1.29, P = 0.003), childhood obesity (OR = 1.07, P = 0.034), adult BMI (OR = 1.38, P < 0.001), adult waist circumference (OR = 1.01, P = 0.028), and adult visceral adiposity (OR = 1.53, P < 0.001) predicted a higher risk of sepsis. Sensitivity analysis did not identify any bias in the MR results. Conclusion: The results demonstrated that adiposity in childhood and adults had causal effects on sepsis incidence. However, more well-designed studies are still needed to validate their association.

Mendelian randomization study and mediation analysis about the relation of inflammatory bowel disease and diabetic retinopathy: the further exploration of gut-retina axis.

Background: The concept of the gut-retinal axis proposed by previous scholars primarily focused on the relationship between intestinal microbiota and retinal diseases, and few further expanded the relationship between intestinal diseases and retinal diseases. To further substantiate the concept of the gut-retinal axis, we analyzed inflammatory bowel disease (IBD) and diabetic retinopathy (DR) using Mendelian randomization (MR), and use mediation analysis to further explore the potential substances that influence this causal relationship. Methods: The genome-wide association study's (GWAS) summary statistics for genetic variations were utilized in a Mendelian randomization (MR) investigation. GWAS data on IBD (including ulcerative colitis (UC), Crohn's disease (CD), and IBD) for non-Finnish Europeans (NFE) were sourced from published articles. In contrast, data on DR (including DR and diabetic maculopathy (DMP)) were obtained from FinnGen R9. The causal relationship has been investigated using inverse variance weighted (IVW), MR-Egger, and weighted median and sensitivity analysis was applied to verify the stability of the results. In addition, we applied mediation analysis to investigate whether circulating inflammatory proteins and plasma lipids played a mediating role, and calculated its effect ratio. Results: The causal relationship between IBD and DR was discovered by employing the inverse variance weighted (IVW) method and weighted median method. In forward MR, UC was significantly associated with lower risk of DR (IVW: OR=0.874; 95%CI= 0.835-0.916; P value= 1.28E-08) (Weighted median: OR=0.893; 95%CI= 0.837-0.954; P value= 7.40E-04). In reverse MR, it was shown that DR (IVW: OR=0.870; 95%CI= 0.828-0.914; P value= 2.79E-08)(Weighted median: OR=0.857; 95%CI= 0.801-0.916; P value= 6.40E-06) and DMP (IVW: OR=0.900; 95%CI= 0.865-0.937; P value= 3.34E-07)(Weighted median: OR=0.882; 95%CI= 0.841-0.924; P value= 1.82E-07) could reduce the risk of CD. What's more, DR is associated with a lower risk of IBD according to genetic prediction (IVW: OR=0.922; 95%CI= 0.873-0.972; P value= 0.002) (Weighted median: OR=0.924; 95%CI= 0.861-0.992; P value= 0.029). Fibroblast growth factor 21 (FGF21), phosphatidylcholine (PC), and triacylglycerol (TG) serve as mediators in these relationships. Conclusions: Our research offers novel insights and sources for investigating the gut-retina axis in the genetic relationship between IBD and DR. We discover four mediators and more about the association between the intestine and retinal disorders and provide more evidence for the gut-retinal axis theory.

Plasma proteins and inflammatory dermatoses: proteome-wide Mendelian randomization and colocalization analyses.

Current genome-wide association studies (GWAS) of plasma proteomes provide additional possibilities for finding new drug targets for inflammatory dermatoses. We performed proteome-wide Mendelian randomization (MR) and colocalization analyses to identify novel potential drug targets for inflammatory dermatoses. We performed MR and colocalization analysis using genetic variation as instrumental variables to determine the causal relationship between circulating plasma proteins and inflammatory dermatoses. 5 plasma proteins were found to be causally associated with dermatitis eczematosa, SLE, urticaria and psoriasis using cis-pQTLs as instrumental variables, but not found in AD and LP. 19 candidate genes with high colocalization evidence were identified. These potential drug targets still require more research and rigorous validation in future trials.

Assessing causal associations of blood counts and biochemical indicators with pulmonary arterial hypertension: a Mendelian randomization study and results from national health and nutrition examination survey 2003-2018.

Background: Blood counts and biochemical markers are among the most common tests performed in hospitals and most readily accepted by patients, and are widely regarded as reliable biomarkers in the literature. The aim of this study was to assess the causal relationship between blood counts, biochemical indicators and pulmonary arterial hypertension (PAH). Methods: A two-sample Mendelian randomization (MR) analysis was performed to assess the causal relationship between blood counts and biochemical indicators with PAH. The genome-wide association study (GWAS) for blood counts and biochemical indicators were obtained from the UK Biobank (UKBB), while the GWAS for PAH were sourced from the FinnGen Biobank. Inverse variance weighting (IVW) was used as the primary analysis method, supplemented by three sensitivity analyses to assess the robustness of the results. And we conducted an observational study using data from National Health and Nutrition Examination Survey (NHANES) 2003-2018 to verify the relationship. Results: The MR analysis primarily using the IVW method revealed genetic variants of platelet count (OR=2.51, 95% CI 1.56-4.22, P<0.001), platelet crit(OR=1.87, 95% CI1.17-7.65, P=0.022), direct bilirubin (DBIL)(OR=1.71, 95%CI 1.18-2.47,P=0.004), insulin-like growth factor (IGF-1)(OR=0.51, 95% CI 0.27-0.96, P=0.038), Lipoprotein A (Lp(a))(OR=0.66, 95% CI 0.45-0.98, P=0.037) and total bilirubin (TBIL)(OR=0.51, 95% CI 0.27-0.96, P=0.038) were significantly associated with PAH. In NHANES, multivariate logistic regression analyses revealed a significant positive correlation between platelet count and volume and the risk of PAH, and a significant negative correlation between total bilirubin and PAH. Conclusion: Our study reveals a causal relationship between blood counts, biochemical indicators and pulmonary arterial hypertension. These findings offer novel insights into the etiology and pathological mechanisms of PAH, and emphasizes the important value of these markers as potential targets for the prevention and treatment of PAH.

Identification of functional enhancer variants associated with type I diabetes in CD4+ T cells.

Type I diabetes is an autoimmune disease mediated by T-cell destruction of ß cells in pancreatic islets. Currently, there is no known cure, and treatment consists of daily insulin injections. Genome-wide association studies and twin studies have indicated a strong genetic heritability for type I diabetes and implicated several genes. As most strongly associated variants are noncoding, there is still a lack of identification of functional and, therefore, likely causal variants. Given that many of these genetic variants reside in enhancer elements, we have tested 121 CD4+ T-cell enhancer variants associated with T1D. We found four to be functional through massively parallel reporter assays. Three of the enhancer variants weaken activity, while the fourth strengthens activity. We link these to their cognate genes using 3D genome architecture or eQTL data and validate them using CRISPR editing. Validated target genes include CLEC16A and SOCS1. While these genes have been previously implicated in type 1 diabetes and other autoimmune diseases, we show that enhancers controlling their expression harbor functional variants. These variants, therefore, may act as causal type 1 diabetic variants.

A causal relationship between particulate matter 2.5 and obesity and its related indicators: a Mendelian randomization study of European ancestry.

Background: In recent years, the prevalence of obesity has continued to increase as a global health concern. Numerous epidemiological studies have confirmed the long-term effects of exposure to ambient air pollutant particulate matter 2.5 (PM2.5) on obesity, but their relationship remains ambiguous. Methods: Utilizing large-scale publicly available genome-wide association studies (GWAS), we conducted univariate and multivariate Mendelian randomization (MR) analyses to assess the causal effect of PM2.5 exposure on obesity and its related indicators. The primary outcome given for both univariate MR (UVMR) and multivariate MR (MVMR) is the estimation utilizing the inverse variance weighted (IVW) method. The weighted median, MR-Egger, and maximum likelihood techniques were employed for UVMR, while the MVMR-Lasso method was applied for MVMR in the supplementary analyses. In addition, we conducted a series of thorough sensitivity studies to determine the accuracy of our MR findings. Results: The UVMR analysis demonstrated a significant association between PM2.5 exposure and an increased risk of obesity, as indicated by the IVW model (odds ratio [OR]: 6.427; 95% confidence interval [CI]: 1.881-21.968; P FDR = 0.005). Additionally, PM2.5 concentrations were positively associated with fat distribution metrics, including visceral adipose tissue (VAT) (OR: 1.861; 95% CI: 1.244-2.776; P FDR = 0.004), particularly pancreatic fat (OR: 3.499; 95% CI: 2.092-5.855; PFDR =1.28E-05), and abdominal subcutaneous adipose tissue (ASAT) volume (OR: 1.773; 95% CI: 1.106-2.841; P FDR = 0.019). Furthermore, PM2.5 exposure correlated positively with markers of glucose and lipid metabolism, specifically triglycerides (TG) (OR: 19.959; 95% CI: 1.269-3.022; P FDR = 0.004) and glycated hemoglobin (HbA1c) (OR: 2.462; 95% CI: 1.34-4.649; P FDR = 0.007). Finally, a significant negative association was observed between PM2.5 concentrations and levels of the novel obesity-related biomarker fibroblast growth factor 21 (FGF-21) (OR: 0.148; 95% CI: 0.025-0.89; P FDR = 0.037). After adjusting for confounding factors, including external smoke exposure, physical activity, educational attainment (EA), participation in sports clubs or gym leisure activities, and Townsend deprivation index at recruitment (TDI), the MVMR analysis revealed that PM2.5 levels maintained significant associations with pancreatic fat, HbA1c, and FGF-21. Conclusion: Our MR study demonstrates conclusively that higher PM2.5 concentrations are associated with an increased risk of obesity-related indicators such as pancreatic fat content, HbA1c, and FGF-21. The potential mechanisms require additional investigation.

Role of plasma fatty acid in age-related macular degeneration: insights from a mendelian randomization analysis.

BACKGROUND: An imbalance in lipid metabolism has been linked to the development of AMD, but the causal relationship between AMD and plasma fatty acids (FAs) remains controversial. Using a two-sample Mendelian randomization (MR) approach, we sought to evaluate the impact of specific FA plasma levels on the risk of different AMD subtypes. METHODS: We analysed genome-wide association data of circulating FAs from 115,006 European-descended individuals in the UK Biobank. These data were used in a two-sample MR framework to assess the potential role of circulating FAs in developing wet and dry AMD. Sensitivity analyses were conducted to ensure the robustness of our findings. Additional multivariable and locus-specific MR analyses were conducted to evaluate direct effects of FA on AMD subtypes, minimizing biases from lipoprotein-related traits and triglycerides. RESULTS: Mendelian randomization revealed associations of omega-3 was associated with decreased wet (OR 0.78, 95%CI 0.66-0.92) and dry AMD (0.85, 0.74-0.97) risk, showed a protective effect on AMD. Notably, the omega-6 to omega-3 ratio showed potential causal effects on both wet (1.27, 1.03-1.56) and dry AMD (1.18, 1.02-1.37). Multivariable MR suggested that the causal relationship of omega-3, omega-6 to omega-3 ratio on wet AMD persists after conditioning on HDL, LDL and triglycerides, albeit with slightly diminished evidence strength. Locus-specific MR linked to omega-3(FADS1, 0.89, 0.82-0.98; FADS2, 0.88, 0.81-0.96) and omega-6 to omega-3 ratio (FADS1, 1.10, 1.02-1.20; FADS2, 1.11, 1.03-1.20) suggests causal effects of these factors on wet AMD. CONCLUSIONS: The associations between plasma FA concentrations and AMD, suggest potential causal role of omega-3, and the omega-6 to omega-3 ratio in wet AMD. These results underscore the impact of an imbalanced circulating omega-3 and omega-6 FA ratio on AMD pathophysiology from MR perspective.

The association of rs17713054 with Neanderthal origin at 3p21.31 locus with the severity of COVID-19 in Iranian patients.

Since the COVID-19 pandemic, the diversity of clinical manifestations in patients has been a tremendous challenge. It seems that genetic variations, as one of the players, contribute to the variety of symptoms. Genome-wide association studies have demonstrated the influence of certain genomic regions on the disease prognosis. Particularly, a haplotype at 3p21.31 locus, inherited from Neanderthals, showed an association with COVID-19 severity. Despite several studies regarding this haplotype, some key variants are not sufficiently addressed. In the present study, we investigated the association of rs17713054 at 3p21.31 with COVID-19 severity. We analyzed the genotype of 251 Iranian COVID-19 patients (151 patients with asymptomatic to mild form as control and 100 patients with severe to critical symptoms without any comorbidities as case group) using the ARMS-PCR method. Results demonstrated that the A allele confers an almost twofold increased risk for COVID-19 severity (P value = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (P value = 0.013). In conclusion, the A allele in rs17713054 was a risk allele in Iranian patients and was independently associated with COVID-19 severity. More studies are beneficial to confirm these findings in other populations and to develop strategies for risk assessment, prevention, and personalized medicine.

Breed of origin analysis in genome-wide association studies: enhancing SNP-based insights into production traits in a commercial Brangus population.

BACKGROUND: Carcass weight (HCW) and marbling (MARB) are critical for meat quality and market value in beef cattle. In composite breeds like Brangus, which meld the genetics of Angus and Brahman, SNP-based analyses have illuminated some genetic influences on these traits, but they fall short in fully capturing the nuanced effects of breed of origin alleles (BOA) on these traits. Focus on the impacts of BOA on phenotypic features within Brangus populations can result in a more profound understanding of the specific influences of Angus and Brahman genetics. Moreover, the consideration of BOA becomes particularly significant when evaluating dominance effects contributing to heterosis in crossbred populations. BOA provides a more comprehensive measure of heterosis due to its ability to differentiate the distinct genetic contributions originating from each parent breed. This detailed understanding of genetic effects is essential for making informed breeding decisions to optimize the benefits of heterosis in composite breeds like Brangus. OBJECTIVE: This study aims to identify quantitative trait loci (QTL) influencing HCW and MARB by utilizing SNP and BOA information, incorporating additive, dominance, and overdominance effects within a multi-generational Brangus commercial herd. METHODS: We analyzed phenotypic data from 1,066 genotyped Brangus steers. BOA inference was performed using LAMP-LD software using Angus and Brahman reference sets. SNP-based and BOA-based GWAS were then conducted considering additive, dominance, and overdominance models. RESULTS: The study identified numerous QTLs for HCW and MARB. A notable QTL for HCW was associated to the SGCB gene, pivotal for muscle growth, and was identified solely in the BOA GWAS. Several BOA GWAS QTLs exhibited a dominance effect underscoring their importance in estimating heterosis. CONCLUSIONS: Our findings demonstrate that SNP-based methods may not detect all genetic variation affecting economically important traits in composite breeds. BOA inclusion in genomic evaluations is crucial for identifying genetic regions contributing to trait variation and for understanding the dominance value underpinning heterosis. By considering BOA, we gain a deeper understanding of genetic interactions and heterosis, which is integral to advancing breeding programs. The incorporation of BOA is recommended for comprehensive genomic evaluations to optimize trait improvements in crossbred cattle populations.

Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes.

BACKGROUND: The cashmere goat industry is one of the main pillars of animal husbandry in Inner Mongolia Autonomous Region, and plays an irreplaceable role in local economic development. With the change in feeding methods and environment, the cashmere produced by Inner Mongolia cashmere goats shows a tendency of coarser, and the cashmere yield can not meet the consumption demand of people. However, the genetic basis behind these changes is not fully understood. We measured cashmere traits, including cashmere yield (CY), cashmere diameter (CD), cashmere thickness (CT), and fleece length (FL) traits for four consecutive years, and utilized Genome-wide association study of four cashmere traits in Inner Mongolia cashmere goats was carried out using new genomics tools to infer genomic regions and functional loci associated with cashmere traits and to construct haplotypes that significantly affect cashmere traits. RESULTS: We estimated the genetic parameters of cashmere traits in Inner Mongolia cashmere goats. The heritability of cashmere yield, cashmere diameter, and fleece length traits of Inner Mongolia cashmere goats were 0.229, 0.359, and 0.250, which belonged to the medium heritability traits (0.2 ~ 0.4). The cashmere thickness trait has a low heritability of 0.053. We detected 151 genome-wide significantly associated SNPs with four cashmere traits on different chromosomes, which were very close to the chromosomes of 392 genes (located within the gene or within ± 500 kb). Notch3, BMPR1B, and CCNA2 have direct functional associations with fibroblasts and follicle stem cells, which play important roles in hair follicle growth and development. Based on GO functional annotation and KEGG enrichment analysis, potential candidate genes were associated with pathways of hair follicle genesis and development (Notch, P13K-Akt, TGF-beta, Cell cycle, Wnt, MAPK). We calculated the effective allele number of the Inner Mongolia cashmere goat population to be 1.109-1.998, the dominant genotypes of most SNPs were wild-type, the polymorphic information content of 57 SNPs were low polymorphism (0 < PIC < 0.25), and the polymorphic information content of 79 SNPs were moderate polymorphism (0.25 < PIC < 0.50). We analyzed the association of SNPs with phenotypes and found that the homozygous mutant type of SNP1 and SNP3 was associated with the highest cashmere yield, the heterozygous mutant type of SNP30 was associated with the lowest cashmere thickness, the wild type of SNP76, SNP77, SNP78, SNP80, and SNP81 was associated with the highest cashmere thickness, and the wild type type of SNP137 was associated with the highest fleece length. 21 haplotype blocks and 68 haplotype combinations were constructed. Haplotypes A2A2, B2B2, C2C2, and D4D4 were associated with increased cashmere yield, haplotypes E2E2, F1F1, G5G5, and G1G5 were associated with decreased cashmere fineness, haplotypes H2H2 was associated with increased cashmere thickness, haplotypes I1I1, I1I2, J1J4, L5L3, N3N2, N3N3, O2O1, P2P2, and Q3Q3 were associated with increased cashmere length. We verified the polymorphism of 8 SNPs by KASP, and found that chr7_g.102631194A > G, chr10_g.82715068 T > C, chr1_g.124483769C > T, chr24_g.12811352C > T, chr6_g.114111249A > G, and chr6_g.115606026 T > C were significantly genotyped in verified populations (P < 0.05). CONCLUSIONS: In conclusion, the genetic effect of single SNP on phenotypes is small, and SNPs are more inclined to be inherited as a whole. By constructing haplotypes from SNPs that are significantly associated with cashmere traits, it will help to reveal the complex and potential causal variations in cashmere traits of Inner Mongolia cashmere goats. This will be a valuable resource for genomics and breeding of the cashmere goat.

Genome-wide association study provided insights into the polled phenotype and polled intersex syndrome (PIS) in goats.

BACKGROUND: Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals. METHODS: The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats. RESULTS: The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats. CONCLUSION: We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

Unveiling the molecular landscape of cognitive aging: insights from polygenic risk scores, DNA methylation, and gene expression.

BACKGROUND: Aging represents a significant risk factor for the occurrence of cerebral small vessel disease, associated with white matter (WM) lesions, and to age-related cognitive alterations, though the precise mechanisms remain largely unknown. This study aimed to investigate the impact of polygenic risk scores (PRS) for WM integrity, together with age-related DNA methylation, and gene expression alterations, on cognitive aging in a cross-sectional healthy aging cohort. The PRSs were calculated using genome-wide association study (GWAS) summary statistics for magnetic resonance imaging (MRI) markers of WM integrity, including WM hyperintensities, fractional anisotropy (FA), and mean diffusivity (MD). These scores were utilized to predict age-related cognitive changes and evaluate their correlation with structural brain changes, which distinguish individuals with higher and lower cognitive scores. To reduce the dimensionality of the data and identify age-related DNA methylation and transcriptomic alterations, Sparse Partial Least Squares-Discriminant Analysis (sPLS-DA) was used. Subsequently, a canonical correlation algorithm was used to integrate the three types of omics data (PRS, DNA methylation, and gene expression data) and identify an individual "omics" signature that distinguishes subjects with varying cognitive profiles. RESULTS: We found a positive association between MD-PRS and long-term memory, as well as a correlation between MD-PRS and structural brain changes, effectively discriminating between individuals with lower and higher memory scores. Furthermore, we observed an enrichment of polygenic signals in genes related to both vascular and non-vascular factors. Age-related alterations in DNA methylation and gene expression indicated dysregulation of critical molecular features and signaling pathways involved in aging and lifespan regulation. The integration of multi-omics data underscored the involvement of synaptic dysfunction, axonal degeneration, microtubule organization, and glycosylation in the process of cognitive aging. CONCLUSIONS: These findings provide valuable insights into the biological mechanisms underlying the association between WM coherence and cognitive aging. Additionally, they highlight how age-associated DNA methylation and gene expression changes contribute to cognitive aging.

Comparative analysis of genetic risk scores for predicting biochemical recurrence in prostate cancer patients after radical prostatectomy.

BACKGROUND: In recent years, Genome-Wide Association Studies (GWAS) has identified risk variants related to complex diseases, but most genetic variants have less impact on phenotypes. To solve the above problems, methods that can use variants with low genetic effects, such as genetic risk score (GRS), have been developed to predict disease risk. METHODS: As the GRS model with the most incredible prediction power for complex diseases has not been determined, our study used simulation data and prostate cancer data to explore the disease prediction power of three GRS models, including the simple count genetic risk score (SC-GRS), the direct logistic regression genetic risk score (DL-GRS), and the explained variance weighted GRS based on directed logistic regression (EVDL-GRS). RESULTS AND CONCLUSIONS: We used 26 SNPs to establish GRS models to predict the risk of biochemical recurrence (BCR) after radical prostatectomy. Combining clinical variables such as age at diagnosis, body mass index, prostate-specific antigen, Gleason score, pathologic T stage, and surgical margin and GRS models has better predictive power for BCR. The results of simulation data (statistical power = 0.707) and prostate cancer data (area under curve = 0.8462) show that DL-GRS has the best prediction performance. The rs455192 was the most relevant locus for BCR (p = 2.496 × 10-6) in our study.

Single-nucleus multi-omics of Parkinson's disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks.

The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

Causal association between immune cells and lung cancer risk: a two-sample bidirectional Mendelian randomization analysis.

Background: Previous studies have highlighted the crucial role of immune cells in lung cancer development; however, the direct link between immunophenotypes and lung cancer remains underexplored. Methods: We applied two-sample Mendelian randomization (MR) analysis, using genetic variants as instruments to determine the causal influence of exposures on outcomes. This method, unlike traditional randomized controlled trials (RCTs), leverages genetic variants inherited randomly at conception, thus reducing confounding and preventing reverse causation. Our analysis involved three genome-wide association studies to assess the causal impact of 731 immune cell signatures on lung cancer using genetic instrumental variables (IVs). We initially used the standard inverse variance weighted (IVW) method and further validated our findings with three supplementary MR techniques (MR-Egger, weighted median, and MR-PRESSO) to ensure robustness. We also conducted MR-Egger intercept and Cochran's Q tests to assess heterogeneity and pleiotropy. Additionally, reverse MR analysis was performed to explore potential causality between lung cancer subtypes and identified immunophenotypes, using R software for all statistical calculations. Results: Our MR analysis identified 106 immune signatures significantly associated with lung cancer. Notably, we found five suggestive associations across all sensitivity tests (P<0.05): CD25 on IgD- CD24- cells in small cell lung carcinoma (ORIVW =0.885; 95% CI: 0.798-0.983; P IVW =0.022); CD27 on IgD+ CD24+ cells in lung squamous cell carcinoma (ORIVW =1.054; 95% CI: 1.010-1.100; P IVW =0.015); CCR2 on monocyte cells in lung squamous cell carcinoma (ORIVW =0.941; 95% CI: 0.898-0.987; P IVW =0.012); CD123 on CD62L+ plasmacytoid dendritic cells (ORIVW =0.958; 95% CI: 0.924-0.992; P IVW =0.017) as well as on plasmacytoid dendritic cells (ORIVW =0.958; 95% CI: 0.924-0.992; P IVW =0.017) in lung squamous cell carcinoma. Conclusion: This study establishes a significant genomic link between immune cells and lung cancer, providing a robust basis for future clinical research aimed at lung cancer management.

Investigating the causal association between gut microbiota and type 2 diabetes: a meta-analysis and Mendelian randomization.

Background: Studies have shown that gut dysbiosis contributes to the pathophysiology of type 2 diabetes mellitus (T2DM). Identifying specific gut microbiota dysbiosis may provide insight into the pathogenesis of T2DM. Purpose: This study investigated the causal relationship between gut microbiota and T2DM using meta-analysis and Mendelian randomization (MR). Methods: In the first part, we searched for literature on gut microbiota and T2DM, and conducted a meta-analysis. We observed differences in glycosylated hemoglobin and fasting blood glucose levels in both groups. Second, we obtained GWAS data from genome-wide association study database 19 (GWAS). We used two-sample MR analysis to verify the forward and reverse causal associations between gut microbiota and T2DM. Additionally, we selected the European GWAS data from the European Bioinformatics Institute (EBI) as a validation set for external validation of the MR analysis. In the third part, we aimed to clarify which gut microbiota contribute to the degree of causal association between group disorders and T2DM through multivariate MR analysis and Bayesian model averaging (MR-BMA). Results: 1. According to the meta-analysis results, the glycated hemoglobin concentration in the gut probiotic intervention group was significantly lower than in the control group. Following treatment, fasting blood glucose levels in the intervention group were significantly lower than those in the control group. 2. The results of two samples MR analysis revealed that there were causal relationships between six gut microbiota and T2DM. Genus Haemophilus and order Pasteurellaceae were negatively correlated with T2DM. Genus Actinomycetes, class Melanobacteria and genus Lactobacillus were positively correlated. Reverse MR analysis demonstrated that T2DM and gut microbiota did not have any reverse causal relationship. The external validation data set showed a causal relationship between gut microbiota and T2DM. 3. Multivariate MR analysis and MR-BMA results showed that the independent genus Haemophilus collection had the largest PP. Conclusion: Our research results suggest that gut microbiota is closely related to T2DM pathogenesis. The results of further MR research and an analysis of the prediction model indicate that a variety of gut microbiota disorders, including genus Haemophilus, are causally related to the development of T2DM. The findings of this study may provide some insight into the diagnosis and treatment of T2DM. Systematic review registration: https://www.crd.york.ac.uk/PROSPERO.

Causal relationships between height, screen time, physical activity, sleep and myopia: univariable and multivariable Mendelian randomization.

Background: This study aims to investigate the independent causal relation between height, screen time, physical activity, sleep and myopia. Methods: Instrumental variables (IVs) for exposures and outcome were obtained from the largest publicly available genome-wide association studies (GWAS) databases. First, we performed a bidirectional univariate MR analysis using primarily the inverse variance weighted method (IVW) with height, screen time, physical activity and sleep as the exposure and myopia as the outcome to investigate the causal relationship between exposures and myopia. Sensitivity analysis was used to demonstrate its robustness. Then the multivariable MR (MVMR) and MR-based mediation approach was further used to estimate the mediating effect of potential confounders (education and time outdoors) on causality. Results: The results of univariate MR analysis showed that taller height (OR = 1.009, 95% CI = 1.005-1.012, p = 3.71 × 10-7), longer time on computer (OR = 1.048, 95% CI = 1.029-1.047, p = 3.87 × 10-7) and less moderate physical activity (OR = 0.976, 95% CI = 0.96-0.991 p = 2.37 × 10-3) had a total effect on the increased risk of developing myopia. Meanwhile our results did not have sufficient evidence to support the causal relationship between chronotype (p = 0.637), sleep duration (p = 0.952) and myopia. After adjusting for education, only taller height remains an independent risk factor for myopia. After adjusting for education, the causal relationship between height, screen and myopia still had statistical significance. A reverse causal relationship was not found in our study. Most of the sensitivity analyses showed consistent results with those of the IVW method. Conclusion: Our MR study revealed that genetically predicted taller height, longer time on computer, less moderate physical activity increased the risk of myopia. After full adjustment for confounders, only height remained independently associated with myopia. As a complement to observational studies, the results of our analysis provide strong evidence for the improvement of myopia risk factors and provide a theoretical basis for future measures to prevent and control myopia in adolescents.

Relationship between Sjogren's syndrome and gastroesophageal reflux: A bidirectional Mendelian randomization study.

The clinical incidence of sjogren's syndrome combined with gastroesophageal reflux disease is high. Existing observational studies have shown inconsistent results in the association between gastroesophageal reflux disease (GERD) and Sjogren's syndrome (SS).We observed that the symptoms of SS patients also improved after receiving GERD-related treatment. Therefore, we aimed to investigate the relationship between GERD and SS through a bidirectional two-sample Mendelian randomization (MR) study. Independent SNPs associated with GERD and SS were selected from a genome-wide association study (GWAS) as instrumental variables to conduct a bidirectional two-sample Mendelian analysis of GERD and SS. Genetic data were obtained from two databases for the following two outcomes: Gastroesophageal reflux (IEU Open GWAS) [sample size = 602,604 (patients = 129,080; nonpatients = 473,524)] and SS (FinnGen) [sample size = 392,423 (patients = 2,495; nonpatients = 389,928)]. Statistical methods for the MR analysis included the inverse-variance weighting method, weighted median, simple mode and weighted mode, as well as heterogeneity and sensitivity analyses using the Cochran Q statistic, MR‒Egger regression, outlier detection methods (MR-PRESSO). In addition, Steiger Test was conducted to test the direction of causality. MR analysis showed a positive correlation between GERD and SS risk [odds ratio (OR) = 1.3279 (95% confidence interval 1.0312-1.7099, P = 0.0280)]. However, in contrast, no significant causal effect of SS on GERD was observed [OR = 1.0024 (95% CI 0.9651-1.0412; P = 0.8995)]. This bidirectional two-sample Mendelian randomization study confirmed a causal relationship between SS and GERD, and suggested that GERD is a risk factor for SS, while SS does not affect GERD.