Non-targeted metabolomics revealed novel links between serum metabolites and primary ovarian insufficiency: a Mendelian randomization study

Background

Primary ovarian insufficiency (POI) is a common clinical endocrine disorder with a high heterogeneity in both endocrine hormones and etiological phenotypes. However, the etiology of POI remains unclear. Herein, we unraveled the causality of genetically determined metabolites (GDMs) on POI through Mendelian randomization (MR) study with the overarching goal of disclosing underlying mechanisms.

Methods

Genetic links with 486 metabolites were retrieved from GWAS data of 7824 European participants as exposures, while GWAS data concerning POI were utilized as the outcome. Via MR analysis, we selected inverse-variance weighted (IVW) method for primary analysis and several additional MR methods (MR-Egger, weighted median, and MR-PRESSO) for sensitivity analyses. MR-Egger intercept and Cochran's Q statistical analysis were conducted to assess potential heterogeneity and pleiotropy. In addition, genetic variations in the key target metabolite were scrutinized further. We conducted replication, meta-analysis, and linkage disequilibrium score regression (LDSC) to reinforce our findings. The MR Steiger test and reverse MR analysis were utilized to assess the robustness of genetic directionality. Furthermore, to deeply explore causality, we performed colocalization analysis and metabolic pathway analysis.

Results

Via IVW methods, our study identified 33 metabolites that might exert a causal effect on POI development. X-11437 showed a robustly significant relationship with POI in four MR analysis methods (P IVW=0.0119; P weighted-median =0.0145; PMR-Egger =0.0499; PMR-PRESSO =0.0248). Among the identified metabolites, N-acetylalanine emerged as the most significant in the primary MR analysis using IVW method, reinforcing its pivotal status as a serum biomarker indicative of an elevated POI risk with the most notable P-value (P IVW=0.0007; PMR-PRESSO =0.0022). Multiple analyses were implemented to further demonstrate the reliability and stability of our deduction of causality. Reverse MR analysis did not provide evidence for the causal effects of POI on 33 metabolites. Colocalization analysis revealed that some causal associations between metabolites and POI might be driven by shared genetic variants.

Conclusion

By incorporating genomics with metabolomics, this study sought to offer a comprehensive analysis in causal impact of serum metabolome phenotypes on risks of POI with implications for underlying mechanisms, disease screening and prevention.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

ROS-induced oxidative stress is a major contributor to sperm cryoinjury

STUDY QUESTION

What is the mechanism behind cryoinjury in human sperm, particularly concerning the interplay between reactive oxygen species (ROS) and autophagy, and how does it subsequently affect sperm fate?

SUMMARY ANSWER

The freeze-thaw operation induces oxidative stress by generating abundant ROS, which impairs sperm motility and activates autophagy, ultimately guiding the sperm toward programmed cell death such as apoptosis and necrosis, as well as triggering premature capacitation.

WHAT IS KNOWN ALREADY

Both ROS-induced oxidative stress and autophagy are thought to exert an influence on the quality of frozen-thawed sperm.

STUDY DESIGN, SIZE, DURATION

Overall, 84 semen specimens were collected from young healthy fertile males, with careful quality evaluation. The specimens were split into three groups to investigate the ROS-induced cryoinjury: normal control without any treatment, sperm treated with 0.5 mM hydrogen peroxide (H2O2) for 1 h, and sperm thawed following cryopreservation. Samples from 48 individuals underwent computer-assisted human sperm analysis (CASA) to evaluate sperm quality in response to the treatments. Semen samples from three donors were analyzed for changes in the sperm proteome after H2O2 treatment, and another set of samples from three donors were analyzed for changes following the freeze-thaw process. The other 30 samples were used for fluorescence-staining and western blotting.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Sperm motility parameters, including progressive motility (PR %) and total motility (PR + NP %), were evaluated using the CASA system on a minimum of 200 spermatozoa. The proteomic profiles were determined with label-free mass spectrometry (MS/MS) and protein identification was performed via ion search against the NCBI human database. Subsequently, comprehensive bioinformatics was applied to detect significant proteomic changes and functional enrichment. Fluorescence-staining and western blot analyses were also conducted to confirm the proteomic changes on selected key proteins. The ROS level was measured using 2',7'-dichlorodihydrofluorescein diacetate labeling and the abundance of bioactive mitochondria was determined by evaluating the inner mitochondrial membrane potential (MMP) level. Molecular behaviors of sequestosome-1 (p62 or SQSTM1) and microtubule-associated proteins 1A/1B light chain 3 (LC3) were monitored to evaluate the state of apoptosis in human sperm. Fluorescent probes oxazole yellow (YO-PRO-1) and propidium iodide (PI) were utilized to monitor programmed cell death, namely apoptosis and necrosis. Additionally, gradient concentrations of antioxidant coenzyme Q10 (CoQ10) were introduced to suppress ROS impacts on sperm.

MAIN RESULTS AND THE ROLE OF CHANCE

The CASA analysis revealed a significant decrease in sperm motility for both the H2O2-treatment and freeze-thaw groups. Fluorescence staining showed that high ROS levels were produced in the treated sperm and the MMPs were largely reduced. The introduction of CoQ10 at concentrations of 20 and 30 μM resulted in a significant rescue of progressive motility (P < 0.05). The result suggested that excessive ROS could be the major cause of sperm motility impairment, likely by damaging mitochondrial energy generation. Autophagy was significantly activated in sperm when they were under oxidative stress, as evidenced by the upregulation of p62 and the increased conversion of LC3 as well as the upregulation of several autophagy-related proteins, such as charged multivesicular body protein 2a, mitochondrial import receptor subunit TOM22 homolog, and WD repeat domain phosphoinositide-interacting protein 2. Additionally, fluorescent staining indicated the occurrence of apoptosis and necrosis in both H2O2-treated sperm and post-thaw sperm. The cell death process can be suppressed when CoQ10 is introduced, which consolidates the view that ROS could be the major contributor to sperm cryoinjury. The freeze-thaw process could also initiate sperm premature capacitation, demonstrated by the prominent increase in tyrosine phosphorylated proteins, verified with anti-phosphotyrosine antibody and immunofluorescence assays. The upregulation of capacitation-related proteins, such as hyaluronidase 3 and Folate receptor alpha, supported this finding.

LARGE SCALE DATA

The data underlying this article are available in the article and its online supplementary material.

LIMITATIONS, REASONS FOR CAUTION

The semen samples were obtained exclusively from young, healthy, and fertile males with progressive motility exceeding 60%, which might overemphasize the positive effects while possibly neglecting the negative impacts of cryoinjury. Additionally, the H2O2 treatment conditions in this study may not precisely mimic the oxidative stress experienced by sperm after thawing from cryopreservation, potentially resulting in the omission of certain molecular alterations.

WIDER IMPLICATIONS OF THE FINDINGS

This study provides substantial proteomic data for a comprehensive and deeper understanding of the impact of cryopreservation on sperm quality. It will facilitate the design of optimal protocols for utilizing cryopreserved sperm to improve applications, such as ART, and help resolve various adverse situations caused by chemotherapy, radiotherapy, and surgery.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by grants from the Major Innovation Project of Research Institute of National Health Commission (#2022GJZD01-3) and the National Key R&D Program of China (#2018YFC1003600). All authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Gene Expression Profiling Reveals Potential Players of Sex Determination and Asymmetrical Development in Chicken Embryo Gonads

Despite the notable progress made in recent years, the understanding of the genetic control of gonadal sex differentiation and asymmetrical ovariogenesis in chicken during embryonic development remains incomplete. This study aimed to identify potential key genes and speculate about the mechanisms associated with ovary and testis development via an analysis of the results of PacBio and Illumina transcriptome sequencing of embryonic chicken gonads at the initiation of sexual differentiation (E4.5, E5.5, and E6.5). PacBio sequencing detected 328 and 233 significantly up-regulated transcript isoforms in females and males at E4.5, respectively. Illumina sequencing detected 95, 296 and 445 DEGs at E4.5, E5.5, and E6.5, respectively. Moreover, both sexes showed asymmetrical expression in gonads, and more DEGs were detected on the left side. There were 12 DEGs involved in cell proliferation shared between males and females in the left gonads. GO analysis suggested that coagulation pathways may be involved in the degradation of the right gonad in females and that blood oxygen transport pathways may be involved in preventing the degradation of the right gonad in males. These results provide a comprehensive expression profile of chicken embryo gonads at the initiation of sexual differentiation, which can serve as a theoretical basis for further understanding the mechanism of bird sex determination and its evolutionary process.

PIP kinases: A versatile family that demands further therapeutic attention

Phosphoinositides are membrane-localized phospholipids that regulate a plethora of essential cellular processes. These lipid signaling molecules are critical for cell homeostasis and therefore their levels are strictly regulated by the coordinated action of several families of lipid kinases and phosphatases. In this review, we provide a focused perspective on the phosphatidylinositol phosphate kinase (PIPK) family and the three subfamilies that compose it: Type I PIPKs or phosphatidylinositol-4-phosphate 5-kinases (PI4P5Ks), Type II PIPKs or phosphatidylinositol-5-phosphate 4-kinases (PI5P4Ks), and Type III PIPKs or phosphatidylinositol-3-phosphate 5-kinases (PIKfyve). Each subfamily is responsible for catalyzing a hydroxyl phosphorylation on specific phosphoinositide species to generate a double phosphorylated lipid, therefore regulating the levels of both substrate and product. Here, we summarize our current knowledge about the functions and regulation of each PIPK subfamily. Further, we highlight the roles of these kinases in various in vivo genetic models and give an overview of their involvement in multiple pathological conditions. The phosphoinositide field has been long focused on targeting PI3K signaling, but growing evidence suggests that it is time to draw attention to the other phosphoinositide kinases. The discovery of the involvement of PIPKs in the pathogenesis of multiple diseases has prompted substantial efforts to turn these enzymes into pharmacological targets. An increasingly refined knowledge of the biology of PIPKs in a variety of in vitro and in vivo models will facilitate the development of effective approaches for therapeutic intervention with the potential to translate into meaningful clinical benefits for patients suffering from cancer, immunological and infectious diseases, and neurodegenerative disorders.

The role of PIP5K1A in cancer development and progression

Malignant tumors are formed via a pathological process of uncontrolled cell division that seriously endangers human physical and mental health. The PI3K/AKT signaling pathway plays an important role in the occurrence and development of various cancers. As a lipid kinase, PIP5K1A acts on the upstream of the PI3K/AKT signaling pathway and has a variety of biological functions, including cell differentiation, cell migration, and sperm development. An increasing number of studies have shown that the overexpression of PIP5K1A promotes the growth, invasion, and migration of cancer cells, and the inhibition of PIP5K1A can effectively hinder tumor progression. These findings imply that PIP5K1A are potential markers and targets for cancers. The aim of this study was to systemically review the structure and function of PIP5K1A, the relationship between PIP5K1A and tumors and the potential therapeutic value of PIP5K1A inhibitors in cancer. PIP5K1A affects the occurrence and progression of many tumors and will provide new strategies for cancer diagnosis and treatment.

Expression pattern and the roles of phosphatidylinositol phosphatases in testis

Phosphoinositides (PIs) are relatively rare lipid components of the cellular membranes. Their homeostasis is tightly controlled by specific PI kinases and phosphatases. PIs play essential roles in cellular signaling, cytoskeletal organization, and secretory processes in various diseases and normal physiology. Gene targeting experiments strongly suggest that in mice with deficiency of several PI phosphatases such as Pten, Mtmrs, Inpp4b, and Inpp5b, spermatogenesis is affected, resulting in partial or complete infertility. Similarly, in men, loss of several of the PIP phosphatases is observed in infertility characterized by the lack of mature sperm. Using available gene expression databases, we compare expression of known PI phosphatases in various testicular cell types, infertility patients, and mouse age-dependent testicular gene expression, and discuss their potential roles in testis physiology and spermatogenesis.

Whole-genome homozygosity mapping reveals candidate regions affecting bull fertility in US Holstein cattle

Background

Achieving rapid genetic progress while maintaining adequate genetic diversity is one of the main challenges facing the dairy industry. The increase in inbreeding can be used to monitor the loss of genetic diversity. Inbreeding tends to increase the proportion of homozygous loci, some of which cause homozygosity of recessive alleles that results in reduced performance. This phenomenon is known as inbreeding depression and tends to be most prominent on fitness-related traits, such as male fertility. Traditionally, inbreeding has been monitored using pedigree information, or more recently, genomic data. Alternatively, it can be quantified using runs of homozygosity (ROH), i.e., contiguous lengths of homozygous genotypes observed in an individual’s chromosome.

Results

The objective of this study was to evaluate the association between ROH and sire conception rate. ROH were evaluated using 268 k genetic markers in 11,790 US Holstein bulls. Interestingly, either the sum, mean, or maximum length of ROH were negatively associated with bull fertility. The association analysis between ROH and sire fertility was performed comparing 300 high-fertility vs. 300 low-fertility bulls. Both the average and sum of ROH length were higher in the low-fertility group. The enrichment of ROH regions in bulls with low fertility was assessed using a Fisher’s exact test. Nine regions were significantly enriched in low-fertility compared to high-fertility bulls. Notably, these regions harbor genes that are closely related to sperm biology and male fertility, including genes exclusively or highly expressed in testis.

Conclusions

The results of this study can help not only to manage inbreeding in genomic selection programs by designing custom mating schemes, but also to better understand the mechanisms underlying male fertility in dairy cattle.

Interrogating the protein interactomes of RAS isoforms identifies PIP5K1A as a KRAS-specific vulnerability

In human cancers, oncogenic mutations commonly occur in the RAS genes KRAS, NRAS, or HRAS, but there are no clinical RAS inhibitors. Mutations are more prevalent in KRAS, possibly suggesting a unique oncogenic activity mediated by KRAS-specific interaction partners, which might be targeted. Here, we determine the specific protein interactomes of each RAS isoform by BirA proximity-dependent biotin identification. The combined interactomes are screened by CRISPR-Cas9 loss-of-function assays for proteins required for oncogenic KRAS-dependent, NRAS-dependent, or HRAS-dependent proliferation and censored for druggable proteins. Using this strategy, we identify phosphatidylinositol phosphate kinase PIP5K1A as a KRAS-specific interactor and show that PIP5K1A binds to a unique region in KRAS. Furthermore, PIP5K1A depletion specifically reduces oncogenic KRAS signaling and proliferation, and sensitizes pancreatic cancer cell lines to a MAPK inhibitor. These results suggest PIP5K1A as a potential target in KRAS signaling for the treatment of KRAS-mutant cancers.

RAS isoforms are frequently mutated in cancer but their inhibition remains challenging. By comparing the protein interactomes of the highly similar isoforms HRAS, NRAS and KRAS, the authors here identify PIP5K1A as a KRAS-specific interactor and a target to inhibit KRAS-driven cell growth.

Regulation of HGF-induced hepatocyte proliferation by the small GTPase Arf6 through the PIP2-producing enzyme PIP5K1A

HGF and its receptor c-Met are critical molecules in various biological processes. Others and we have previously shown that the small GTPase Arf6 plays a pivotal role in HGF signaling in hepatocytes. However, the molecular mechanism of how Arf6 regulates HGF signaling is unclear. Here, we show that Arf6 plays an important role in HGF-stimulated hepatocyte proliferation and liver regeneration through the phosphatidylinositol 4,5-bisphosphate (PIP₂)-producing enzyme PIP5K1A. We find that knockdown of Arf6 and PIP5K1A in HepG2 cells inhibits HGF-stimulated proliferation, Akt activation, and generation of phosphatidylinositol 3,4,5-trisphosphate (PIP₃) and its precursor PIP₂. Interestingly, PIP5K1A is recruited to c-Met upon HGF stimulation in an Arf6 activity-dependent manner. Finally, we show that hepatocyte proliferation and liver regeneration after partial hepatectomy are suppressed in Pip5k1a knockout mice. These results provide insight into the molecular mechanism for HGF-stimulated hepatocyte proliferation and liver regeneration: Arf6 recruits PIP5K1A to c-Met and activates it upon HGF stimulation to produce PIP₂ and subsequently PIP₃, which in turn activates Akt to promote hepatocyte proliferation, thereby accelerating liver regeneration after liver injury.

Activation of the Small G Protein Arf6 by Dynamin2 through Guanine Nucleotide Exchange Factors in Endocytosis

The small G protein Arf6 and the GTPase dynamin2 (Dyn2) play key roles in clathrin-mediated endocytosis (CME). However, their functional relationship remains obscure. Here, we show that Arf6 functions as a downstream molecule of Dyn2 in CME. Wild type of Dyn2 overexpressed in HeLa cells markedly activates Arf6, while a GTPase-lacking Dyn2 mutant does not. Of the Arf6-specific guanine nucleotide exchange factors, EFA6A, EFA6B, and EFA6D specifically interact with Dyn2. Furthermore, overexpression of dominant negative mutants or knockdown of EFA6B and EFA6D significantly inhibit Dyn2-induced Arf6 activation. Finally, overexpression of the binding region peptide of EFA6B for Dyn2 or knockdown of EFA6B and EFA6D significantly suppresses clathrin-mediated transferrin uptake. These results provide evidence for a novel Arf6 activation mechanism by Dyn2 through EFA6B and EFA6D in CME in a manner dependent upon the GTPase activity of Dyn2.

The role of PI3K/AKT-related PIP5K1α and the discovery of its selective inhibitor for treatment of advanced prostate cancer

Nitrogen-containing heterocyclic compounds are an important class of molecules that are commonly used for the synthesis of candidate drugs. Phosphatidylinositol-4-phosphate 5-kinase-α (PIP5Kα) is a lipid kinase, similar to PI3K. However, the role of PIP5K1α in oncogenic processes and the development of inhibitors that selectively target PIP5K1α have not been reported. In the present study we report that overexpression of PIP5K1α is associated with poor prognosis in prostate cancer and correlates with an elevated level of the androgen receptor. Overexpression of PIP5K1α in PNT1A nonmalignant cells results in an increased AKT activity and an increased survival, as well as invasive malignant phenotype, whereas siRNA-mediated knockdown of PIP5K1α in aggressive PC-3 cells leads to a reduced AKT activity and an inhibition in tumor growth in xenograft mice. We further report a previously unidentified role for PIP5K1α as a druggable target for our newly developed compound ISA-2011B using a high-throughput KINOMEscan platform. ISA-2011B was discovered during our synthetic studies of C-1 indol-3-yl substituted 1,2,3,4-tetrahydroisoquinolines via a Pictet-Spengler approach. ISA-2011B significantly inhibits growth of tumor cells in xenograft mice, and we show that this is mediated by targeting PIP5K1α-associated PI3K/AKT and the downstream survival, proliferation, and invasion pathways. Further, siRNA-mediated knockdown of PIP5K1α exerts similar effects on PC3 cells as ISA-2011B treatment, significantly inhibiting AKT activity, increasing apoptosis and reducing invasion. Thus, PIP5K1α has high potential as a drug target, and compound ISA-2011B is interesting for further development of targeted cancer therapy.