Identifying potential biomarkers for non-obstructive azoospermia using WGCNA and machine learning algorithms

Objective

The cause and mechanism of non-obstructive azoospermia (NOA) is complicated; therefore, an effective therapy strategy is yet to be developed. This study aimed to analyse the pathogenesis of NOA at the molecular biological level and to identify the core regulatory genes, which could be utilised as potential biomarkers.

Methods

Three NOA microarray datasets (GSE45885, GSE108886, and GSE145467) were collected from the GEO database and merged into training sets; a further dataset (GSE45887) was then defined as the validation set. Differential gene analysis, consensus cluster analysis, and WGCNA were used to identify preliminary signature genes; then, enrichment analysis was applied to these previously screened signature genes. Next, 4 machine learning algorithms (RF, SVM, GLM, and XGB) were used to detect potential biomarkers that are most closely associated with NOA. Finally, a diagnostic model was constructed from these potential biomarkers and visualised as a nomogram. The differential expression and predictive reliability of the biomarkers were confirmed using the validation set. Furthermore, the competing endogenous RNA network was constructed to identify the regulatory mechanisms of potential biomarkers; further, the CIBERSORT algorithm was used to calculate immune infiltration status among the samples.

Results

A total of 215 differentially expressed genes (DEGs) were identified between NOA and control groups (27 upregulated and 188 downregulated genes). The WGCNA results identified 1123 genes in the MEblue module as target genes that are highly correlated with NOA positivity. The NOA samples were divided into 2 clusters using consensus clustering; further, 1027 genes in the MEblue module, which were screened by WGCNA, were considered to be target genes that are highly correlated with NOA classification. The 129 overlapping genes were then established as signature genes. The XGB algorithm that had the maximum AUC value (AUC=0.946) and the minimum residual value was used to further screen the signature genes. IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers. This 5 biomarker model had the highest AUC value, of up to 0.982, compared to other single biomarker models; additionally, the results of this biomarker model were verified in the validation set.

Conclusions

As IL20RB, C9orf117, HILS1, PAOX, and DZIP1 have been determined to possess the strongest association with NOA, these five genes could be used as potential therapeutic targets for NOA patients. Furthermore, the model constructed using these five genes, which possessed the highest diagnostic accuracy, may be an effective biomarker model that warrants further experimental validation.

The highly conserved FOXJ1 target CFAP161 is dispensable for motile ciliary function in mouse and Xenopus

Cilia are protrusions of the cell surface and composed of hundreds of proteins many of which are evolutionary and functionally well conserved. In cells assembling motile cilia the expression of numerous ciliary components is under the control of the transcription factor FOXJ1. Here, we analyse the evolutionary conserved FOXJ1 target CFAP161 in Xenopus and mouse. In both species Cfap161 expression correlates with the presence of motile cilia and depends on FOXJ1. Tagged CFAP161 localises to the basal bodies of multiciliated cells of the Xenopus larval epidermis, and in mice CFAP161 protein localises to the axoneme. Surprisingly, disruption of the Cfap161 gene in both species did not lead to motile cilia-related phenotypes, which contrasts with the conserved expression in cells carrying motile cilia and high sequence conservation. In mice mutation of Cfap161 stabilised the mutant mRNA making genetic compensation triggered by mRNA decay unlikely. However, genes related to microtubules and cilia, microtubule motor activity and inner dyneins were dysregulated, which might buffer the Cfap161 mutation.

Genome-Wide Association Study Identifies Genomic Loci of Sex Determination and Gonadosomatic Index Traits in Large Yellow Croaker (Larimichthys crocea)

Larimichthys crocea is one of the traditional marine culture fishes in China, widely distributed in South China Sea, East Sea, and southern Yellow Sea. Sex dimorphism is evident in this species that females present a substantial growth strength than males, suggesting breeding females could obtain more economic benefits in L. crocea aquaculture industry. With the continuous expansion of aquaculture industry, both identifying sex-associated genome region and understanding the genetic basis underlying gonad differentiation and development matter to not only sex control aquaculture but also breeding industry. Thus, genome-wide association analysis (GWAS) of sex determination was conducted with a random breeding population of 905 individuals (including 463 females and 442 males) by ddRAD sequencing. For sex determination, 21 significant single nucleotide polymorphisms (SNPs) in chromosome (Chr) 22 were identified. Surrounding these SNPs, we founded 14 candidate genes, including dmrt1, dmrt3, and piwil2, fam102a, and odf2. The sex-associated region was narrowed down further to 2.4 Mb on Chr22 through F_st scanning and insertion-deletion (InDel) analysis. Besides, 3 SNPs in the supposed sex-determining region on Chr22 were identified as highly associated with gonad differentiation through GWAS on gonadosomatic index (GSI) in 350 males and 231 females. Because of the significant difference of GSI between females and males of L. crocea, GWAS on GSI of different genders was also conducted independently. Finally, we identified a SNP in Chr18 showing genome-wide significant association with male GSI (MGSI) and three genes axl, cyp2a10, and cyp2g1 involved in the gonadal development regulation process of aromatase. Overall, this study explored the genetic basis of sex determination mechanism and provided novel insights into gonad differentiation and development, offering solid genetic support for sex control breeding, marker-assisted selection, and marine resources conservation.

Androglobin gene expression patterns and FOXJ1-dependent regulation indicate its functional association with ciliogenesis

Androglobin (ADGB) represents the latest addition to the globin superfamily in metazoans. The chimeric protein comprises a calpain domain and a unique circularly permutated globin domain. ADGB expression levels are most abundant in mammalian testis, but its cell-type-specific expression, regulation, and function have remained unexplored. Analyzing bulk and single-cell mRNA-Seq data from mammalian tissues, we found that—in addition to the testes—ADGB is prominently expressed in the female reproductive tract, lungs, and brain, specifically being associated with cell types forming motile cilia. Correlation analysis suggested coregulation of ADGB with FOXJ1, a crucial transcription factor of ciliogenesis. Investigating the transcriptional regulation of the ADGB gene, we characterized its promoter using epigenomic datasets, exogenous promoter-dependent luciferase assays, and CRISPR/dCas9-VPR-mediated activation approaches. Reporter gene assays revealed that FOXJ1 indeed substantially enhanced luciferase activity driven by the ADGB promoter. ChIP assays confirmed binding of FOXJ1 to the endogenous ADGB promoter region. We dissected the minimal sequence required for FOXJ1-dependent regulation and fine mapped the FOXJ1 binding site to two evolutionarily conserved regions within the ADGB promoter. FOXJ1 overexpression significantly increased endogenous ADGB mRNA levels in HEK293 and MCF-7 cells. Similar results were observed upon RFX2 overexpression, another key transcription factor in ciliogenesis. The complex transcriptional regulation of the ADGB locus was illustrated by identifying a distal enhancer, responsible for synergistic regulation by RFX2 and FOXJ1. Finally, cell culture studies indicated an ADGB-dependent increase in the number of ciliated cells upon overexpression of the full-length protein, confirming a ciliogenesis-associated role of ADGB in mammals.

Motile cilia genetics and cell biology: big results from little mice

Our understanding of motile cilia and their role in disease has increased tremendously over the last two decades, with critical information and insight coming from the analysis of mouse models. Motile cilia form on specific epithelial cell types and typically beat in a coordinated, whip-like manner to facilitate the flow and clearance of fluids along the cell surface. Defects in formation and function of motile cilia result in primary ciliary dyskinesia (PCD), a genetically heterogeneous disorder with a well-characterized phenotype but no effective treatment. A number of model systems, ranging from unicellular eukaryotes to mammals, have provided information about the genetics, biochemistry, and structure of motile cilia. However, with remarkable resources available for genetic manipulation and developmental, pathological, and physiological analysis of phenotype, the mouse has risen to the forefront of understanding mammalian motile cilia and modeling PCD. This is evidenced by a large number of relevant mouse lines and an extensive body of genetic and phenotypic data. More recently, application of innovative cell biological techniques to these models has enabled substantial advancement in elucidating the molecular and cellular mechanisms underlying the biogenesis and function of mammalian motile cilia. In this article, we will review genetic and cell biological studies of motile cilia in mouse models and their contributions to our understanding of motile cilia and PCD pathogenesis.

Scalable phylogenetic profiling using MinHash uncovers likely eukaryotic sexual reproduction genes

Phylogenetic profiling is a computational method to predict genes involved in the same biological process by identifying protein families which tend to be jointly lost or retained across the tree of life. Phylogenetic profiling has customarily been more widely used with prokaryotes than eukaryotes, because the method is thought to require many diverse genomes. There are now many eukaryotic genomes available, but these are considerably larger, and typical phylogenetic profiling methods require at least quadratic time as a function of the number of genes. We introduce a fast, scalable phylogenetic profiling approach entitled HogProf, which leverages hierarchical orthologous groups for the construction of large profiles and locality-sensitive hashing for efficient retrieval of similar profiles. We show that the approach outperforms Enhanced Phylogenetic Tree, a phylogeny-based method, and use the tool to reconstruct networks and query for interactors of the kinetochore complex as well as conserved proteins involved in sexual reproduction: Hap2, Spo11 and Gex1. HogProf enables large-scale phylogenetic profiling across the three domains of life, and will be useful to predict biological pathways among the hundreds of thousands of eukaryotic species that will become available in the coming few years. HogProf is available at https://github.com/DessimozLab/HogProf.

Genes that are involved in the same biological process tend to co-evolve. This property is exploited by the technique of phylogenetic profiling, which identifies co-evolving (and therefore likely functionally related) genes through patterns of correlated gene retention and loss in evolution and across species. However, conventional methods to computing and clustering these correlated genes do not scale with increasing numbers of genomes. HogProf is a novel phylogenetic profiling tool built on probabilistic data structures. It allows the user to construct searchable databases containing the evolutionary history of hundreds of thousands of protein families. Such fast detection of coevolution takes advantage of the rapidly increasing amount of genomic data publicly available, and can uncover unknown biological networks and guide in-vivo research and experimentation. We have applied our tool to describe the biological networks underpinning sexual reproduction in eukaryotes.

The FOXJ1 target Cfap206 is required for sperm motility, mucociliary clearance of the airways and brain development

Cilia are complex cellular protrusions consisting of hundreds of proteins. Defects in ciliary structure and function, many of which have not been characterised molecularly, cause ciliopathies: a heterogeneous group of human syndromes. Here, we report on the FOXJ1 target gene Cfap206, orthologues of which so far have only been studied in Chlamydomonas and Tetrahymena In mouse and Xenopus, Cfap206 was co-expressed with and dependent on Foxj1 CFAP206 protein localised to the basal body and to the axoneme of motile cilia. In Xenopus crispant larvae, the ciliary beat frequency of skin multiciliated cells was enhanced and bead transport across the epidermal mucociliary epithelium was reduced. Likewise, Cfap206 knockout mice revealed ciliary phenotypes. Electron tomography of immotile knockout mouse sperm flagella indicated a role in radial spoke formation reminiscent of FAP206 function in Tetrahymena Male infertility, hydrocephalus and impaired mucociliary clearance of the airways in the absence of laterality defects in Cfap206 mutant mice suggests that Cfap206 may represent a candidate for the subgroup of human primary ciliary dyskinesias caused by radial spoke defects.

A framework for high-resolution phenotyping of candidate male infertility mutants: from human to mouse

Male infertility is a heterogeneous condition of largely unknown etiology that affects at least 7% of men worldwide. Classical genetic approaches and emerging next-generation sequencing studies support genetic variants as a frequent cause of male infertility. Meanwhile, the barriers to transmission of this disease mean that most individual genetic cases will be rare, but because of the large percentage of the genome required for spermatogenesis, the number of distinct causal mutations is potentially large. Identifying bona fide causes of male infertility thus requires advanced filtering techniques to select for high-probability candidates, including the ability to test causality in animal models. The mouse remains the gold standard for defining the genotype-phenotype connection in male fertility. Here, we present a best practice guide consisting of (a) major points to consider when interpreting next-generation sequencing data performed on infertile men, and, (b) a systematic strategy to categorize infertility types and how they relate to human male infertility. Phenotyping infertility in mice can involve investigating the function of multiple cell types across the testis and epididymis, as well as sperm function. These findings will feed into the diagnosis and treatment of male infertility as well as male health broadly.

Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries)

Maximizing the number of offspring born per female is a key functionality trait in commercial- and/or subsistence-oriented livestock enterprises. Although the number of offspring born is closely associated with female fertility and reproductive success, the genetic control of these traits remains poorly understood in sub-Saharan Africa livestock. Using selection signature analysis performed on Ovine HD BeadChip data from the prolific Bonga sheep in Ethiopia, 41 candidate regions under selection were identified. The analysis revealed one strong selection signature on a candidate region on chromosome X spanning BMP15, suggesting this to be the primary candidate prolificacy gene in the breed. The analysis also identified several candidate regions spanning genes not reported before in prolific sheep but underlying fertility and reproduction in other species. The genes associated with female reproduction traits included SPOCK1 (age at first oestrus), GPR173 (mediator of ovarian cyclicity), HB-EGF (signalling early pregnancy success) and SMARCAL1 and HMGN3a (regulate gene expression during embryogenesis). The genes involved in male reproduction were FOXJ1 (sperm function and successful fertilization) and NME5 (spermatogenesis). We also observed genes such as PKD2L2, MAGED1 and KDM3B, which have been associated with diverse fertility traits in both sexes of other species. The results confirm the complexity of the genetic mechanisms underlying reproduction while suggesting that prolificacy in the Bonga sheep, and possibly African indigenous sheep is partly under the control of BMP15 while other genes that enhance male and female fertility are essential for reproductive fitness.

Role of primary cilia in non-dividing and post-mitotic cells

The essential role of primary (non-motile) cilia during the development of multi-cellular tissues and organs is well established and is underlined by severe disease manifestations caused by mutations in cilia-associated molecules that are collectively termed ciliopathies. However, the role of primary cilia in non-dividing and terminally differentiated, post-mitotic cells is less well understood. Although the prevention of cells from re-entering the cell cycle may represent a major chore, primary cilia have recently been linked to DNA damage responses, autophagy and mitochondria. Given this connectivity, primary cilia in non-dividing cells are well positioned to form a signaling hub outside of the nucleus. Such a center could integrate information to initiate responses and to maintain cellular homeostasis if cell survival is jeopardized. These more discrete functions may remain undetected until differentiated cells are confronted with emergencies.

1700012B09Rik, a FOXJ1 effector gene active in ciliated tissues of the mouse but not essential for motile ciliogenesis

In humans and mice, motile cilia occur on the surface of the embryonic ventral node, on respiratory and ependymal epithelia and in reproductive organs where they ensure normal left-right asymmetry of the organism, mucociliary clearance of airways, homeostasis of the cerebrospinal fluid and fertility. The genetic programme for the formation of motile cilia, thus critical for normal development and health, is switched on by the key transcription factor FOXJ1. In previous microarray screens for murine FOXJ1 effectors, we identified candidates for novel factors involved in motile ciliogenesis, including both genes that are well conserved throughout metazoa and beyond, like FOXJ1 itself, and genes without overt homologues outside higher vertebrates. Here we examine one of the novel murine FOXJ1 effectors, the uncharacterised 1700012B09Rik whose homologues appear to be restricted to higher vertebrates. In mouse embryos and adults, 1700012B09Rik is predominantly expressed in motile ciliated tissues in a FOXJ1-dependent manner. 1700012B09RIK protein localises to basal bodies of cilia in cultured cells. Detailed analysis of 1700012B09Rik^lacZ knock-out mice reveals no impaired function of motile cilia or non-motile cilia. In conclusion, this novel FOXJ1 effector is associated mainly with motile cilia but - in contrast to other known FOXJ1 targets - its putative ciliary function is not essential for development or health in the mouse, consistent with a late emergence during evolution of motile ciliogenesis.