Human Genome-wide Association Studies Research Articles

PandasGWAS: a Python package for easy retrieval of GWAS catalog data

BackgroundSince the NHGRI-EBI Catalog of human genome-wide association studies was established by NHGRI in 2008, research on it has attracted more and more researchers as the amount of data has grown rapidly. Easy-to-use, open-source, general-purpose programs for accessing the NHGRI-EBI Catalog of human genome-wide association studies are in great demand for current Python data analysis pipeline.ResultsIn this work we present pandasGWAS, a Python package that provides programmatic access to the NHGRI-EBI Catalog of human genome-wide association studies. Instead of downloading all data locally, pandasGWAS queries data based on input criteria and handles paginated data gracefully. The data is then transformed into multiple associated pandas.DataFrame objects according to its hierarchical relationships, which makes it easy to integrate into current Python-based data analysis toolkits.ConclusionspandasGWAS is an open-source Python package that provides the first Python client interface to the GWAS Catalog REST API. Compared with existing tools, the data structure of pandasGWAS is more consistent with the design specification of GWAS Catalog REST API, and provides many easy-to-use mathematical symbol operations.

Whole-genome sequencing study to identify candidate markers indicating susceptibility to type 2 diabetes in Bama miniature pigs.

Hundreds of single-nucleotide polymorphism (SNP) sites have been found to be potential genetic markers of type 2 diabetes mellitus (T2DM). However, SNPs related to T2DM in minipigs have been less reported. This study aimed to screen the T2DM-susceptible candidate SNP loci in Bama minipigs so as to improve the success rate of the minipig T2DM model. The genomic DNAs of three Bama minipigs with T2DM, six sibling low-susceptibility minipigs with T2DM, and three normal control minipigs were compared by whole-genome sequencing. The T2DM Bama minipig-specific loci were obtained, and their functions were annotated. Meanwhile, the Biomart software was used to perform homology alignment with T2DM-related loci obtained from the human genome-wide association study to screen candidate SNP markers for T2DM in Bama miniature pigs. Whole-genome resequencing detected 6960 specific loci in the minipigs with T2DM, and 13 loci corresponding to 9 diabetes-related genes were selected. Further, a set of 122 specific loci in 69 orthologous genes of human T2DM candidate genes were obtained in the pigs. Collectively, a batch of T2DM-susceptible candidate SNP markers in Bama minipigs, covering 16 genes and 135 loci, was established. Whole-genome sequencing and comparative genomics analysis of the orthologous genes in pigs that corresponded to the human T2DM-related variant loci successfully screened out T2DM-susceptible candidate markers in Bama miniature pigs. Using these loci to predict the susceptibility of the pigs before constructing an animal model of T2DM may help to establish an ideal animal model.

Assessment of the Response of Young and Aged Sarcospan-Deficient Mice to Nutrient Excess

This study follows up on human genome-wide association studies that have identified variants in the SSPN locus. This locus encodes the protein sarcospan that segregates with several anthropometric traits. One of these traits includes an increase in waist circumference. In this study, we compared the response of young and old SSPN-deficient (SSPN-/-) and WT control mice to high fat diet (HFD) (60% fat). The HFD was initiated in 2 month and 1 year old male and female for 4 months and 5 months respectively. This study showed that SSPN-deficiency provided partial protection from weight gain in young male mice when compared to age-matched WT mice and a significantly smaller increase in body weight (BW) compared to controls. Also, young female SSPN-/- mice gained 50% less BW than WT age matched mice compared to controls and visceral WAT fat pad weights were significantly lower. To examine whether these differences are lost in aging (e.g. due to hormonal decline), we subjected 12-month old WT and SSPN-/- mice to 5 months of HFD and assessed anthropometric traits, blood lipids, glucose tolerance, body composition, and alterations of cytokine/adipokine/chemokine levels in blood serum and adipose tissue. There were signs of inflammation/fibrosis in young WT male WAT. This was largely absent in SSPN-/- mice. Interestingly, a specific subset of upregulated inflammatory mediators was found in WT WAT including leptin. However, another distinct subset of adipokines/chemokines were upregulated including adiponectin in SSPN-/- WAT. Therefore, we are investigating whether SSPN affects obesity-associated inflammatory signaling in adipocytes that influences adipose tissue expansion and ectopic triglyceride accumulation in the liver. FDOH JEK 21K12 This is the full abstract presented at the American Physiology Summit 2023 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

A human FLII gene variant alters sarcomeric actin thin filament length and predisposes to cardiomyopathy

To better understand the genetic basis of heart disease, we identified a variant in the Flightless-I homolog (FLII) gene that generates a R1243H missense change and predisposes to cardiac remodeling across multiple previous human genome-wide association studies (GWAS). Since this gene is of unknown function in the mammalian heart we generated gain- and loss-of-function genetically altered mice, as well as knock-in mice with the syntenic R1245H amino acid substitution, which showed that Flii protein binds the sarcomeric actin thin filament and influences its length. Deletion of Flii from the heart, or mice with the R1245H amino acid substitution, show cardiomyopathy due to shortening of the actin thin filaments. Mechanistically, Flii is a known actin binding protein that we show associates with tropomodulin-1(TMOD1) to regulate sarcomere thin filament length. Indeed, overexpression of leiomodin-2 in the heart, which lengthens the actin-containing thin filaments, partially rescued disease due to heart-specific deletion of Flii. Collectively, the identified FLII human variant likely increases cardiomyopathy risk through an alteration in sarcomere structure and associated contractile dynamics, like other sarcomere gene-based familial cardiomyopathies.

Zebrafish as a model for drug induced liver injury: state of the art and beyond

Zebrafish as a preclinical drug induced liver injury (DILI) model provides multiple advantages ranging from ease of breeding and maintenance, availability of different strains and transgenic fish amenable to study liver function, and highly conserved liver structure and function with the human liver. In this review, the authors have aimed to provide an account of the metabolic enzymes that take roles in drug detoxification in both human and zebrafish in a comparative manner and exemplify several recent models in studying liver functionality. Moreover, the authors emphasize the difficulties associated with studying idiosyncratic DILI in preclinical models and propose that zebrafish could be an important complement to mice in testing functions of genes that are associated with DILI with respect to different drugs in human genome-wide association studies (GWAS) Catalog. Finally, this review highlights the state-of-the-art in the development of novel transgenic reporter strains that can be used to study degree and molecular mechanisms of hepatotoxicity caused by drugs in zebrafish. All of these will help researchers to use effectively the available resources in the zebrafish DILI models, while advocating potential leads that can be taken to provide advancements in a better understanding and treatment of DILI.

Genome-wide CRISPR screening of chondrocyte maturation newly implicates genes in skeletal growth and height-associated GWAS loci

Alterations in the growth and maturation of chondrocytes can lead to variation in human height, including monogenic disorders of skeletal growth. We aimed to identify genes and pathways relevant to human growth by pairing human height genome-wide association studies (GWASs) with genome-wide knockout (KO) screens of growth-plate chondrocyte proliferation and maturation invitro. We identified 145 genes that alter chondrocyte proliferation and maturation at early and/or late time points in culture, with 90% of genes validating in secondary screening. These genes are enriched in monogenic growth disorder genes and in KEGG pathways critical for skeletal growth and endochondral ossification. Further, common variants near these genes capture height heritability independent of genes computationally prioritized from GWASs. Our study emphasizes the value of functional studies in biologically relevant tissues as orthogonal datasets to refine likely causal genes from GWASs and implicates new genetic regulators of chondrocyte proliferation and maturation.

A gene-level test for directional selection on gene expression.

Most variants identified in human genome-wide association studies and scans for selection are non-coding. Interpretation of their effects and the way in which they contribute to phenotypic variation and adaptation in human populations is therefore limited by our understanding of gene regulation and the difficulty of confidently linking non-coding variants to genes. To overcome this, we developed a gene-wise test for population-specific selection based on combinations of regulatory variants. Specifically, we use the QX statistic to test for polygenic selection on cis-regulatory variants based on whether the variance across populations in predicted expression of a particular gene is higher than expected under neutrality. We then applied this approach to human data, testing for selection on 17,388 protein-coding genes in 26 populations from the Thousand Genomes Project. We identified 45 genes with significant evidence (FDR<0.1) for selection, including FADS1, KHK, SULT1A2, ITGAM, and several genes in the HLA region. We further confirm that these signals correspond to plausible population-level differences in predicted expression. While the small number of significant genes (0.2%) is consistent with most cis-regulatory variation evolving under genetic drift or stabilizing selection, it remains possible that there are effects not captured in this study. Our gene-level QX score is independent of standard genomic tests for selection, and may therefore be useful in combination with traditional selection scans to specifically identify selection on regulatory variation. Overall, our results demonstrate the utility of combining population-level genomic data with functional data to understand the evolution of gene expression.

Polygenic influences on the behavioral effects of alcohol withdrawal in a mixed-ancestry population from the collaborative study on the genetics of alcoholism (COGA)

Alcohol withdrawal (AW) is a feature of alcohol use disorder that may occur in up to half of individuals with chronic, heavy alcohol consumption whenever alcohol use is abruptly stopped or significantly reduced. To date, few genes have been robustly associated with AW; this may be partly due to most studies defining AW as a binary construct despite the multiple symptoms and their range in severity from mild to severe. The current study examined the effects of genome-wide loci on a factor score for AW in high risk and community family samples in the Collaborative Study for the Genetics of Alcoholism (COGA). In addition, we tested whether differentially expressed genes associated with alcohol withdrawal in model organisms are enriched in human genome-wide association study (GWAS) effects. Analyses employed roughly equal numbers of males and females (mean age 35, standard deviation = 15; total N = 8009) and included individuals from multiple ancestral backgrounds. Genomic data were imputed to the HRC reference panel and underwent strict quality control procedures using Plink2. Analyses controlled for age, sex, and population stratification effects using ancestral principal components. We found support that AW is a polygenic disease (SNP-heritability = 0.08 [95 % CI = 0.01, 0.15; pedigree-based heritability = 0.12 [0.08,0.16]. We identified five single nucleotide variants that met genomewide significance, some of which have previously been associated with alcohol phenotypes. Gene-level analyses suggest a role for COL19A1 in AW; H-MAGMA analyses implicated 12 genes associated with AW. Cross-species enrichment analyses indicated that variation within genes identified in model organism studies explained <1 % of the phenotypic variability in human AW. Notably, the surrounding regulatory regions of model organism genes explained more variance than expected by chance, indicating that these regulatory regions and gene sets may be important for human AW. Lastly, when comparing the overlap in genes identified from the human GWAS and H-MAGMA analyses with the genes identified from the animal studies, there was modest overlap, indicating some convergence between the methods and organisms.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases.

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

An exponential increase in QTL detection with an increased sample size.

Power analyses are often used to determine the number of animals required for a genome-wide association study (GWAS). These analyses are typically intended to estimate the sample size needed for at least 1 locus to exceed a genome-wide significance threshold. A related question that is less commonly considered is the number of significant loci that will be discovered with a given sample size. We used simulations based on a real data set that consisted of 3,173 male and female adult N/NIH heterogeneous stock rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in subsamples of the full data set. The subsampling analysis was conducted for 4 traits with low (0.15 ± 0.03), medium (0.31 ± 0.03 and 0.36 ± 0.03), and high (0.46 ± 0.03) SNP-based heritabilities. For each trait, we subsampled the data 100 times at different sample sizes (500, 1,000, 1,500, 2,000, and 2,500). We observed an exponential increase in the number of significant loci with larger sample sizes. Our results are consistent with similar observations in human GWAS and imply that future rodent GWAS should use sample sizes that are significantly larger than those needed to obtain a single significant result.

Genetic Networks of Alzheimer's Disease, Aging, and Longevity in Humans.

Human genomic analysis and genome-wide association studies (GWAS) have identified genes that are risk factors for early and late-onset Alzheimer's disease (AD genes). Although the genetics of aging and longevity have been extensively studied, previous studies have focused on a specific set of genes that have been shown to contribute to or are a risk factor for AD. Thus, the connections among the genes involved in AD, aging, and longevity are not well understood. Here, we identified the genetic interaction networks (referred to as pathways) of aging and longevity within the context of AD by using a gene set enrichment analysis by Reactome that cross-references more than 100 bioinformatic databases to allow interpretation of the biological functions of gene sets through a wide variety of gene networks. We validated the pathways with a threshold of p-value < 1.00 × 10-5 using the databases to extract lists of 356 AD genes, 307 aging-related (AR) genes, and 357 longevity genes. There was a broad range of biological pathways involved in AR and longevity genes shared with AD genes. AR genes identified 261 pathways within the threshold of p < 1.00 × 10-5, of which 26 pathways (10% of AR gene pathways) were further identified by overlapping genes among AD and AR genes. The overlapped pathways included gene expression (p = 4.05 × 10-11) including ApoE, SOD2, TP53, and TGFB1 (p = 2.84 × 10-10); protein metabolism and SUMOylation, including E3 ligases and target proteins (p = 1.08 × 10-7); ERBB4 signal transduction (p = 2.69 × 10-6); the immune system, including IL-3 and IL-13 (p = 3.83 × 10-6); programmed cell death (p = 4.36 × 10-6); and platelet degranulation (p = 8.16 × 10-6), among others. Longevity genes identified 49 pathways within the threshold, of which 12 pathways (24% of longevity gene pathways) were further identified by overlapping genes among AD and longevity genes. They include the immune system, including IL-3 and IL-13 (p = 7.64 × 10-8), plasma lipoprotein assembly, remodeling and clearance (p < 4.02 × 10-6), and the metabolism of fat-soluble vitamins (p = 1.96 × 10-5). Thus, this study provides shared genetic hallmarks of aging, longevity, and AD backed up by statistical significance. We discuss the significant genes involved in these pathways, including TP53, FOXO, SUMOylation, IL4, IL6, APOE, and CEPT, and suggest that mapping the gene network pathways provide a useful basis for further medical research on AD and healthy aging.

BACH1 deficiency prevents neointima formation and maintains the differentiated phenotype of vascular smooth muscle cells by regulating chromatin accessibility.

The transcription factor BTB and CNC homology 1(BACH1) has been linked to coronary artery disease risk by human genome-wide association studies, but little is known about the role of BACH1 in vascular smooth muscle cell (VSMC) phenotype switching and neointima formation following vascular injury. Therefore, this study aims to explore the role of BACH1 in vascular remodeling and its underlying mechanisms. BACH1 was highly expressed in human atherosclerotic plaques and has high transcriptional factor activity in VSMCs of human atherosclerotic arteries. VSMC-specific loss of Bach1 in mice inhibited the transformation of VSMC from contractile to synthetic phenotype and VSMC proliferation and attenuated the neointimal hyperplasia induced by wire injury. Mechanistically, BACH1 suppressed chromatin accessibility at the promoters of VSMC marker genes via recruiting histone methyltransferase G9a and cofactor YAP and maintaining the H3K9me2 state, thereby repressing VSMC marker genes expression in human aortic smooth muscle cells (HASMCs). BACH1-induced repression of VSMC marker genes was abolished by the silencing of G9a or YAP. Thus, these findings demonstrate a crucial regulatory role of BACH1 in VSMC phenotypic transition and vascular homeostasis and shed light on potential future protective vascular disease intervention via manipulation of BACH1.

Identification of candidate genes for nicotine withdrawal in C57BL/6J×DBA/2J recombinant inbred mice.

Nicotine is the reinforcing ingredient in tobacco. Following chronic exposure, sudden cessation of nicotine use produces negative symptoms of withdrawal that contribute to dependence. The molecular mechanisms underlying nicotine withdrawal behaviors, however, are poorly understood. Using recombinant inbred mice, chronic nicotine was delivered by minipump and withdrawal induced using mecamylamine. Somatic signs of withdrawal, and anxiety-like behavior using elevated plus maze, were then assessed. Interval mapping was used to identify associations between genetic variation and withdrawal behaviors, and with basal gene expression. Differential gene expression following nicotine exposure and withdrawal was also assessed in progenitor mice using microarrays. Quantitative trait loci mapping identified chromosome intervals with significant genetic associations to somatic signs of withdrawal or withdrawal-induced anxiety-like behavior. Using bioinformatics, and association with basal gene expression in nucleus accumbens, we implicated Rb1, Bnip3l, Pnma2, Itm2b, and Kif13b as candidate genes for somatic signs of withdrawal, and Galr1, which showed trans-regulation from a region of chromosome 14 that was associated with somatic signs of withdrawal. Candidate genes within the chromosome 9 region associated with anxiety-like withdrawal behavior included Dixdc1, Ncam1, and Sorl1. Bioinformatics identified six genes that were also significantly associated with nicotine or alcohol traits in recent human genome-wide association studies. Withdrawal-associated somatic signs and anxiety-like behavior had strong non-overlapping genetic associations, respectively, with regions of chromosome 14 and chromosome 9. Genetic, behavioral and gene expression correlations, and bioinformatics analysis identified several candidate genes that may represent novel molecular targets for modulating nicotine withdrawal symptoms.

Knockout of murine Lyplal1 confers sex-specific protection against diet-induced obesity.

Human genome-wide association studies found single-nucleotide polymorphisms (SNPs) near LYPLAL1 (Lysophospholipase-like protein 1) that have sex-specific effects on fat distribution and metabolic traits. To determine whether altering LYPLAL1 affects obesity and metabolic disease, we created and characterized a mouse knockout (KO) of Lyplal1. We fed the experimental group of mice a high-fat, high-sucrose (HFHS) diet for 23 weeks, and the controls were fed regular chow diet. Here, we show that CRISPR-Cas9 whole-body Lyplal1 KO mice fed an HFHS diet showed sex-specific differences in weight gain and fat accumulation as compared to chow diet. Female, not male, KO mice weighed less than WT mice, had reduced body fat percentage, had white fat mass, and had adipocyte diameter not accounted for by changes in the metabolic rate. Female, but not male, KO mice had increased serum triglycerides, decreased aspartate, and decreased alanine aminotransferase. Lyplal1 KO mice of both sexes have reduced liver triglycerides and steatosis. These diet-specific effects resemble the effects of SNPs near LYPLAL1 in humans, suggesting that LYPLAL1 has an evolutionary conserved sex-specific effect on adiposity. This murine model can be used to study this novel gene-by-sex-by-diet interaction to elucidate the metabolic effects of LYPLAL1 on human obesity.

Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10−8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats.

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

Long-Term Simulation of Microgravity Induces Changes in Gene Expression in Breast Cancer Cells.

Microgravity changes the gene expression pattern in various cell types. This study focuses on the breast cancer cell lines MCF-7 (less invasive) and MDA-MB-231 (triple-negative, highly invasive). The cells were cultured for 14 days under simulated microgravity (s-µg) conditions using a random positioning machine (RPM). We investigated cytoskeletal and extracellular matrix (ECM) factors as well as focal adhesion (FA) and the transmembrane proteins involved in different cellular signaling pathways (MAPK, PAM and VEGF). The mRNA expressions of 24 genes of interest (TUBB, ACTB, COL1A1, COL4A5, LAMA3, ITGB1, CD44, VEGF, FLK1, EGFR, SRC, FAK1, RAF1, AKT1, ERK1, MAPK14, MAP2K1, MTOR, RICTOR, VCL, PXN, CDKN1, CTNNA1 and CTNNB1) were determined by quantitative real-time PCR (qPCR) and studied using STRING interaction analysis. Histochemical staining was carried out to investigate the morphology of the adherent cells (ADs) and the multicellular spheroids (MCSs) after RPM exposure. To better understand this experimental model in the context of breast cancer patients, a weighted gene co-expression network analysis (WGCNA) was conducted to obtain the expression profiles of 35 breast cell lines from the HMS LINCS Database. The qPCR-verified genes were searched in the mammalian phenotype database and the human genome-wide association studies (GWAS) Catalog. The results demonstrated the positive association between the real metastatic microtumor environment and MCSs with respect to the extracellular matrix, cytoskeleton, morphology, different cellular signaling pathway key proteins and several other components. In summary, the microgravity-engineered three-dimensional MCS model can be utilized to study breast cancer cell behavior and to assess the therapeutic efficacies of drugs against breast cancer in the future.

Diabetes Mellitus and Pancreatic Cancer: Investigation of Causal Pathways Through Mendelian Randomization Analysis

Abstract Objective: This study was aimed at investigating the association between diabetes mellitus and susceptibility to pancreatic cancer by using Mendelian randomization (MR) methods and an extensive human genome-wide association study (GWAS) dataset. Methods: The publicly accessible MR Base database was used to obtain the complete genome, relevant research findings, and summary data pertaining to diabetes mellitus and pancreatic cancer. Genetic variables, specifically single-nucleotide polymorphisms closely associated with diabetes mellitus, were selected for analysis. Four methods—inverse variance weighted (IVW) analysis, weighted median analysis, weighted mode, and MR-Egger regression—were used. Statistical analysis was conducted to explore the potential association between diabetes mellitus and susceptibility to pancreatic cancer. Results: The results of the IVW analysis (OR = 11.56519319, 95% CI 1.275068624–104.8992116, P = 0.0296) indicated a significant causal relationship between diabetes and elevated pancreatitis risk. Furthermore, the absence of horizontal pleiotropic effects (Egger intercept = 0.29, P = 0.384) and heterogeneity (P = 0.126) suggested that the observed association was not influenced by confounding factors. Sensitivity analysis and other statistical methods also supported the conclusion that genetic pleiotropy did not introduce bias to the findings. Conclusion: A causal relationship exists between diabetes mellitus and the occurrence of pancreatic cancer. People with diabetes mellitus are at high risk of pancreatic cancer and should receive early screening. The IGF signaling pathway may be a key mediator of the effects of diabetes on pancreatic cancer pathogenesis.

DSec16 Acting in Insulin-like Peptide Producing Cells Controls Energy Homeostasis in Drosophila.

Many studies show that genetics play a major contribution to the onset of obesity. Human genome-wide association studies (GWASs) have identified hundreds of genes that are associated with obesity. However, the majority of them have not been functionally validated. SEC16B has been identified in multiple obesity GWASs but its physiological role in energy homeostasis remains unknown. Here, we use Drosophila to determine the physiological functions of dSec16 in energy metabolism. Our results showed that global RNAi of dSec16 increased food intake and triglyceride (TAG) levels. Furthermore, this TAG increase was observed in flies with a specific RNAi of dSec16 in insulin-like peptide producing cells (IPCs) with an alteration of endocrine peptides. Together, our study demonstrates that dSec16 acting in IPCs controls energy balance and advances the molecular understanding of obesity.

Bone-cartilage crosstalk informed by aging mouse bone transcriptomics and human osteoarthritis genome-wide association studies

Subchondral bone participates in crosstalk with articular cartilage to maintain joint homeostasis, and disruption of either tissue results in overall joint degeneration. Among the subchondral bone changes observed in osteoarthritis (OA), subchondral bone plate (SBP) thickening has a time-dependent relationship with cartilage degeneration and has recently been shown to be regulated by osteocytes. Here, we evaluate the effect of age on SBP thickness and cartilage degeneration in aging mice. We find that SBP thickness significantly increases by 18-months of age, corresponding temporally with increased cartilage degeneration. To identify factors in subchondral bone that may participate in bone cartilage crosstalk or OA, we leveraged mouse transcriptomic data from one joint tissue compartment – osteocyte-enriched bone – to search for enrichment with human OA in UK Biobank and Arthritis Research UK Osteoarthritis Genetics (arcOGEN) GWAS using the mouse2human (M2H, www.mouse2human.org) strategy. Genes differentially expressed in aging mouse bone are significantly enriched for human OA, showing joint site-specific (knee vs. hip) relationships, exhibit temporal associations with age, and unique gene clusters are implicated in each type of OA. Application of M2H identifies genes with known and unknown functions in osteocytes and OA development that are clinically associated with human OA. Altogether, this work prioritizes genes with a potential role in bone/cartilage crosstalk for further mechanistic study based on their association with human OA in GWAS.