- 1Department of Endocrinology, Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 2Research Centre for Clinical Medicine, Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
Background: Androgen insensitivity syndrome (AIS) is a rare X-linked genetic disease and one of the causes of 46,XY disorder of sexual development. The unstraightforward diagnosis of AIS and the gender assignment dilemma still make a plague for this disorder due to the overlapping clinical phenotypes.
Methods: Peripheral blood mononuclear cells (PBMCs) of partial AIS (PAIS) patients and healthy controls were separated, and RNA-seq was performed to investigate transcriptome variance. Then, tissue-specific gene expression, functional enrichment, and protein–protein interaction (PPI) network analyses were performed; and the key modules were identified. Finally, the RNA expression of differentially expressed genes (DEGs) of interest was validated by quantitative real-time PCR (qRT-PCR).
Results: In our dataset, a total of 725 DEGs were captured, with functionally enriched reproduction and immune-related pathways and Gene Ontology (GO) functions. The most highly specific systems centered on hematologic/immune and reproductive/endocrine systems. We finally filtered out CCR1, PPBP, PF4, CLU, KMT2D, GP6, and SPARC by the key gene clusters of the PPI network and manual screening of tissue-specific gene expression. These genes provide novel insight into the pathogenesis of AIS in the immune system or metabolism and bring forward possible molecular markers for clinical screening. The qRT-PCR results showed a consistent trend in the expression levels of related genes between PAIS patients and healthy controls.
Conclusion: The present study sheds light on the molecular mechanisms underlying the pathogenesis and progression of AIS, providing potential targets for diagnosis and future investigation.
Introduction
Androgen insensitivity syndrome (AIS; OMIM#300068), one of the common causes of 46,XY disorder/difference of sex development (DSD), is estimated to affect one in 20,000–100,000 live births and accounts for 40%–80% of 46,XY DSD patients (1–4). This disorder is divided into three categories, ranging from complete feminization with external genitalia [complete AIS (CAIS)] to various degrees of undervirilization [partial AIS (PAIS)] and to male infertility and/or gynecomastia [mild AIS (MAIS)]. AIS is an X-linked recessive disease caused by inactivating mutations in the androgen receptor (AR) gene (Xq11-q12), resulting in complete or partial resistance to the physiological effects of androgen in 46,XY individuals (5).
AR is a nuclear receptor that belongs to the steroid hormone receptor family. By binding to AR, androgens, such as testosterone (T) and dihydrotestosterone (DHT), could translocate into the nucleus of androgen target cells and bind to androgen response elements (AREs), which would regulate the expression of target genes (6). The AR-mediated effects of androgens are well known to be normally key elements for achieving appropriate internal and external male sex differentiation in both embryonic and secondary sexual development. In addition, normal AR function enables many target tissues to respond to androgen signaling, including the development and maintenance of bone and muscle mass, production and maturation of immune cells and thymocytes, and regulation of metabolism (7). The expression of various target genes is also controlled by AR, which could regulate cell growth, biosynthetic processes, and metabolism. Moreover, the repair of DNA damage could be affected by androgen signaling by regulating gene expression and transcription-related processes (8). Therefore, defects in AR could compromise the virilization process, exhibiting variably impaired masculinization of external genitalia, spermatogenesis arrest, gynecomastia, sparseness or absence of pubic and axillary hair, increased risk of gonadal tumors, and abnormal sex hormone levels (9). However, this disorder also leads to abnormal symptoms of other systems in males: slightly increased height as compared with normal females, disorder of skeletal development, inherited X-linked neurodegenerative disease, increased incidence of insulin resistance and cardiovascular risk, and risk of infection and autoimmune diseases were reported in CAIS and aromatase knockout (ARKO) mice (7, 10–13). However, except for AR mutation, the mechanisms underlying the genetic etiology of AIS have rarely been investigated.
The clinical diagnosis of AIS, especially PAIS, still represents a demanding challenge due to the phenotypes overlapping with 46,XY DSDs caused by other etiologies. Current evaluations rely mainly on clinical and biochemical features, karyotype analysis, AR genotype analysis, and functional studies of external genital skin-derived fibroblasts to detect the expression and binding capacity of AR receptors and androgen-regulated apolipoprotein D (APOD) (1). Pathogenic mutations of AR themselves result in lasting changes in gene expression in genital fibroblasts (14). However, the invasive procedure restricts the widespread clinical application of direct evidence derived from external genital skin fibroblasts. Moreover, in some patients with typical AIS phenotypes and impaired function of AR in the study of fibroblasts, no AR mutation was identified (14).
In this study, we performed comprehensive transcriptomic profiling of peripheral blood mononuclear cells (PBMCs) from genotype-proven PAIS patients and normal controls via RNA-seq analyses and focused on the differentially expressed genes (DEGs) of pubertal sex differentiation and fertility in males. The mechanism of AR defects was investigated from the potential genes that had never been studied in connection to AIS previously. Our results, linking the clinical features and molecular basis, contribute to the understanding of the mechanisms underlying the pathogenesis of AIS and provide new insights for convenient and economical measures that are helpful for clinical screening of AIS.
Materials and Methods
Subjects
The peripheral blood of three PAIS patients and three healthy male volunteers (age-matched) was obtained from Shanghai Ninth People’s Hospital (Shanghai Jiao Tong University School of Medicine, Shanghai, China). All the patients presented with micropenis, hypospadias, cryptorchid, and elevated serum testosterone levels. The external masculinization score (EMS) index was applied to describe the degree of undervirilization of the patients. Genetic analysis was performed by targeted next-generation sequencing (NGS), and two missense mutations (p.V686A and p.R761W) and a very rare synonymous mutation p.S889S [an exonic splicing mutation (15)] of the AR gene were identified in three AIS patients (Table 1). The variations were judged as pathogenic according to American College of Medical Genetics and Genomics (ACMG) guidelines, consistent with the classical phenotype, including micropenis, hypospadias, bifid scrotum, and cryptorchid. All subjects provided written informed consent. The Ethics Committee of the Shanghai Ninth People’s Hospital reviewed and approved this study.
Isolation of Peripheral Blood Mononuclear Cells and RNA-seq Library Preparation
Five hundred microliters of saturated ethylenediaminetetraacetic acid (EDTA) was diluted with 50 ml of phosphate-buffered saline (PBS). The separated anticoagulated plasma samples were transferred into 50-ml centrifuge tubes and mixed with EDTA–PBS up to 50 ml. Along the walls of new 50-ml centrifuge tubes, the plasma–EDTA–PBS mixtures were carefully added to 15 ml of Lymphoprep (STEMCELL Technologies, Vancouver, Canada), a density gradient medium. PBMCs were centrifuged through the Lymphoprep layer during centrifugation at 1,200 g for 20 min at 4°C and then washed with EDTA–PBS. Total RNA was extracted from PBMCs by using TRIzol reagent (Invitrogen, USA) for RNA-seq. Single-end reads of 50 bp were used to sequence the libraries on the Illumina HiSeq platform (Igenecode Technology). For data analysis, Cufflinks (version 2.2.1) was used to annotate the fragment per kilobase of exons per million mapped reads (FPKM) against GENCODE (Release M22, GRCm38. p6).
Differential Expression Analysis
To identify the DEGs between PAIS patients and normal males, gene expression profiles were compared by using the R-Bioconductor package limma. Student’s t-tests were used to calculate p-values. The DEGis were screened out by the following selection criteria: |log2 fold change (FC) > 1|, and the threshold value of the p-value was 0.05.
Functional Enrichment Analysis of Differentially Expressed Genes
The functional categories Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses of DEGs were performed by using Metascape (http://metascape.org/gp/index.html) (16). Functional enrichment research on large scales is often conducted through GO analysis, in which gene functions are categorized as biological process (BP), molecular function (MF), and cellular component (CC). DEGs were assigned to specific pathways using KEGG pathway analysis, which enabled the construction of reaction, molecular interaction, and relationship networks (17). The following criteria in Metascape were used to obtain a significant difference: minimum count of 3, enrichment factor of >1.5, and pick selective GO clusters. p < 0.05 was considered statistically significant.
Protein–Protein Interaction Network Analysis and Gene Cluster Identification
The Search Tool for the Retrieval of Interacting Genes (STRING; http://string-db.org) (version 10.0) online database was used to predict the protein–protein interaction (PPI) network between proteins encoded by the top 236 DEGs (|log2 FC > 1.5|) and other proteins. PPI pairs with a combined score >0.4 were extracted. Subsequently, a PPI network was generated using Cytoscape Software (http://cytoscape.org/). We used the Molecular Complex Detection (MCODE) plugin to perform gene network clustering analysis, which identified the most significant modules within PPI networks with the standard degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and maximum depth = 100.
Hub Gene Selection and Analysis
The degree of each protein node was calculated, and the top 10 genes identified as hub genes were picked according to the CytoHubba plugin of Cytoscape. The Cytoscape plugin Biological Networks Gene Oncology tool (BiNGO) (version 3.0.3) was used to visualize the BP analysis of the top five hub genes (degree >10) (18). The color depth of nodes represents the corrected p-value of ontologies, and the size of nodes indicates the numbers of genes that are contained in the ontologies. p < 0.01 indicated a significant difference.
Tissue-Specific Gene Expression Analysis
The tissue-specific expression of the DEGs was analyzed by using the online resource BioGPS (http://biogps.org). The inclusion criteria were as follows: 1) expression level of >10 times the median for transcripts that mapped to a single tissue, and 2) the second highest expression level was no more than 1/3 as high as the highest level (19).
Transcription Factor Analysis
Epigenetic landscape in silico detection (LISA) was used for transcription factor (TF) analysis (20). The upregulated and downregulated DEGs were entered into the gene lists to analyze the key TFs contributing to DEGs. TFs ranking lists and volcano plots were created automatically.
Quantitative Real-Time PCR
Reverse transcription was conducted using the PrimeScript RT Master Mix Kit (Takara, Dalian, China). The mRNA level was assessed using TB Green Premix Ex Taq Kit (Takara, Dalian, China) following the instructions. Primers are available in Supplementary Table 4. The relative expression of mRNA was calculated by the 2−ΔΔCt method with normalization to GAPDH.
Results
Differential Expression Analysis
The clinical features, serum hormones, and AR mutations of the three cases are summarized in Table 1. To determine the transcriptomic variance among the PAIS patients and normal controls, RNA-seq was carried out on each sample to obtain the transcriptome profiles (Figure 1A). Following the elimination of adapter sequences and low-quality reads, signals of 14,318 transcripts were obtained from the six samples. Groups were well discriminated from patients and normal in a principal component analysis (PCA) model, as shown in Figure 1B. Based on the criteria of |log2 FC > 1| and p < 0.05, we observed that PBMCs of patients exhibited significantly different gene expression patterns (Figure 1C). In Figure 2, according to |log2 FC|, we chose the top 100 DEGs to visualize their expression patterns, chromosomal locations, and logarithmic adjusted p-value shown in the inner layer. Compared with those in normal individuals, a total of 725 DEGs, including 495 upregulated and 230 downregulated genes, were identified in PAIS patients (Supplementary Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g001.jpg)
Figure 1 Experiment process and principal component analysis (PCA) and volcano plot of differentially expressed genes (DEGs). (A) Schematic diagram of experiment procedure. (B) PCA of two groups. (C) DEGs in partial androgen insensitivity syndrome (PAIS) patients and healthy controls samples are shown in the volcano plot, with blue dots representing significantly downregulated genes and red dots representing significantly upregulated genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g002.jpg)
Figure 2 Circos plot of expression patterns and chromosomal positions of top 100 differentially expressed genes (DEGs). The outer circle represents chromosomes, and lines coming from each gene point to their specific chromosomal locations. The group of patients is represented in the inner circular heatmaps. According to |log2 fold change|, the top four upregulated genes (red) and the top four downregulated genes (blue) are connected with red and blue lines in the center of the Circos plot.
Functional Enrichment Analysis of Differentially Expressed Genes
According to the p-value <0.05, 15 and five KEGG pathways of the DEGs with upregulation and downregulation respectively were screened out, and the pie chart was drawn according to the enrichment score (Figures 3A, B). The enriched signaling pathways with upregulated DEGs included platelet activation, proteoglycans in cancer, and chemokine signaling pathways. The highly expressed genes were related to sex differentiation and cell development signaling pathways, such as phagosome, cAMP, notch, and adherens junction. The downregulated DEGs were significantly enriched in ribosome, TGF-β. For GO, the top 13 of these terms based on their p-value are shown in chord plots (Figure 4). The DEGs were enriched mainly in BPs, with upregulated DEGs specifically centered on platelet activation, actin polymerization or depolymerization, regulation of cell-substrate junction organization, platelet degranulation, regulation of neuron projection arborization, cell junction assembly, etc. The majority of the downregulated DEGs were related to DNA repair, gene expression, meiotic nuclear division, and apoptosis. Regarding MF, actin binding, calcium ion transmembrane transporter activity, and kinase activity were the most remarkably enriched GO terms in the upregulated DEGs, while structural constituent of ribosome, cyclin-dependent protein serine/threonine kinase regulator activity, and rRNA binding were mainly enriched in the downregulated DEGs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g003.jpg)
Figure 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of all differentially expressed genes (DEGs) and biological process analysis of the hub genes (degree >10). (A) KEGG pathways of upregulated DEGs. (B) KEGG pathways of downregulated DEGs. (C) The biological process analysis of hub genes (degree >10). The color depth of nodes refers to the corrected p-value of ontologies. The size of nodes refers to the numbers of genes that are involved in the ontologies.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g004.jpg)
Figure 4 Gene Ontology (GO) analyses of all differentially expressed genes (DEGs). (A) GO enrichment of the upregulated DEGs. (B) GO enrichment of the downregulated DEGs.
Protein–Protein Interaction Network Analysis and Gene Cluster Identification
According to the criteria of |log2 FC > 1.5|, 182 nodes and 123 edges comprised the PPI network; the network had an interaction score >0.4 and PPI enrichment p-value < 1.91e−07 based on the STRING online database (Figure 5A). We identified four key modules with upregulated genes by MCODE to further analyze the PPI network of DEGs (Figures 5B–E, Table 2). Genes in the first gene module with the highest score were primarily related to gene silencing and platelet- and immune system-related functions. Furthermore, the results showed that genes in the other three modules were enriched mainly in methylation, cell adhesion synaptic transmission, and inflammatory response.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g005.jpg)
Figure 5 Protein–protein interaction (PPI) network and the key modules of top 236 differentially expressed genes (DEGs) (|log2 (fold change) >1.5|). (A) Cytoscape network visualization of the 182 nodes and 123 edges that were obtained with interaction scores >0.4 according to the STRING online database. The nodes denote proteins, and the edges correspond to the interactions between two proteins. Red, blue, and yellow colors represent upregulated, downregulated, and hub genes, respectively (B–E) Four key modules were identified by MCODE.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-t002.jpg)
Table 2 MCODE was used to process the data downloaded from the STRING to further mine gene clusters.
Hub Gene Selection and Analysis
In the PPI network, the top 10 genes that were evaluated by connectivity degree were considered hub genes for subsequent analysis, including PF4, GP1BA, ITGA2B, GP9, PPBP, SPARC, GP6, ITGB3, CTTN, and KMT2D. The BP analysis of the top five hub genes (degree >10) is shown in Figure 3C, and these genes were found to be principally involved in coagulation, cell adhesion, cytolysis, cell division, and immune response. Many hub genes, such as GP1BA, GP9, ITGA2B, ITGB3, PPBP, and PF4, were found to be key genes enriched in the chemokine signaling pathway, focal adhesion, extracellular matrix (ECM)–receptor interaction, and PI3K–Akt signaling pathway within each breed.
Tissue-Specific Gene Expression Analysis
BioGPS identified 75 genes among the top 234 DEGs (cutoff of |log2 FC > 1.5|) that were expressed in a particular tissue or organ system (Supplementary Table 2). The most highly represented systems were hematologic/immune and reproductive/endocrine systems, which accounted for 57.3% of the total tissue-specific expression genes. The neurologic and muscle systems had a relatively low percentage, which was 18.6% (14/75) and 8% (6/75) of tissue-specific genes, respectively, while the urinary and digestive systems had similar levels of enrichment (2.6%, 2/75). In addition, the circulatory and fat systems had the lowest enrichment levels (1.3%, 1/75) (Table 3).
General Analysis of Transcription Factors
TFs bind to specific DNA sequences to regulate gene expression activation or repression. Coactivators and TFs interact to control the spatial and temporal expression of genes (21). LISA TF analysis was performed to identify the TFs that caused differential gene expression in PAIS. Following inputting selected genes for analysis, key TFs were subsequently generated (Figure 6). The top 10 upregulated and four downregulated candidates may participate in the genesis and progression of PAIS, including RNF2, FOXA1, TAL1, BRD4, GLI2, 5HMC, PBX3, RING1 and TALE, SUZ12, MED1, and MEF2C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g006.jpg)
Figure 6 Key transcription factors (TFs) generated by Lisa analysis. The top 10 TFs are displayed in orange and marked with gene symbols beside them. Key TFs potentially contributing to the differentially expressed genes (DEGs) in upregulated and downregulated genes.
Potential Genes Related to Disorder of Sex Development
The results revealed that the gene expression profiles were significantly disturbed in patients with PAIS compared with normal controls. This disturbance is especially visible in genes closely related to sex development (Figure 7). In our analysis, we distinguished three downregulated DEGs (MEI1, CHD5, and MORN2) in connection to spermatogenesis, spermatid development, and spermatid differentiation (Supplementary Table 3). ACTR3C, ITGA7, and ATF4 were involved in Sertoli cell function, whereas HSD17B3 and CYP1B1 participated in testosterone metabolic processes. Androgen signaling pathway including important transcriptional coactivators and corepressors, such as NCOA2, ARID1A, DAB2, NCOR2, ROCK2, and NR2C2, also showed significant differences. Several genes were associated mostly with sex differentiation, aging, and cell development. In addition, we also observed significant transcriptome differences in autophagy between patients and controls.
Identification of Genes of Interest
MCODE analysis detected key clusters consisting of 22 genes (Table 2). Combined with the tissue-specific gene expression analysis, we identified genes of interest that were particularly upregulated in the hematologic/immune system (CCR1, PPBP, and PF4) and the metabolism/endocrine system (GP9 and SPARC). Additionally, we manually identified two additional genes of interest (CLU and KMT2D) that are potentially involved in the pathogenesis of PAIS by the GeneCards database. The abbreviations, names, and functions for these 22 genes of interest are shown in Table 4.
Validation of the Differentially Expressed Genes of Interest in Partial Androgen Insensitivity Syndrome Patients
The selected biomarkers were validated in PBMCs of PAIS using qRT-PCR analysis. Consistent with the results of the mRNA microarray, the expression level of CCR1 was obviously higher in PAIS patients than in controls. However, the expression levels of PPBP, PF4, CLU, KMT2D, GP6, and SPARC showed the same trends as the bioinformatics analysis results, although there was no significant statistical difference between the two groups due to the limited number of samples (Figure 8). Notably, the expression levels of the above target genes were the lowest in the PBMCs of patient 1 with synonymous mutation p.S889S, suggesting relatively mild impairment of this variation with a small fraction of normal AR transcripts. In addition, the expression levels of CYP1B1, HSD17B3, and MEI1, which are involved in the androgen signaling pathway, were also changed (Figure 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/731107/fendo-12-731107-HTML-r1/image_m/fendo-12-731107-g008.jpg)
Figure 8 RNA expression of genes of interest were measured in partial androgen insensitivity syndrome (PAIS) and healthy samples. RNA expression of CCR1, PPBP, PF4, CLU, KMT2D, GP6, SPARC, CYP1B1, HSD17B3, and MEI1 was measured. p-Values were calculated using a two-sided unpaired Student’s t-test. *p < 0.05; ns, not significant.
Discussion
AIS is a rare endocrine disorder falling within the category of 46,XY DSD caused by mutations within the AR gene located on the X chromosome long arm and is classified as CAIS, PAIS, or MAIS (22). PAIS, a relatively frequent type, has a variable clinical presentation depending on the degree to which the external genitalia responds to androgens (23, 24). Psychological distress occurred due to the vague gender identity of the patients (25). The risk of germ cell tumors seems to be higher in patients with PAIS than in patients with CAIS (26). Therefore, it is urgent to reveal potential markers that are helpful for accurate diagnosis and prompt treatment. However, effective biomarkers for diagnosis remain scarce, as well as the reasonable predictors for androgen therapy, which actually are important for sex assignment accordingly (27). To our knowledge, for the first time, we investigated the variance of gene expression profiles of PBMCs in genotyping proven PAIS by RNA-seq and identified potential hub genes as candidate predicted factors for the disorder.
AR and the androgen–AR signaling pathway play significant roles in sexually dimorphic development and the function of male reproductive and non-reproductive organs, such as sex development, spermatogenesis, regulation of growth and metabolism, muscle mass/function, bone metabolism, and certain aspects of the immune and cardiovascular system (28). More attention has been drawn to the impaired function of AR in sexual development, but defects in other regulatory processes or systems have long been neglected. Our results are consistent with previous limited studies, providing reasonable support from different aspects. In this study, we identified 725 DEGs, including 495 upregulated and 230 downregulated genes, in the PBMCs of PAIS patients compared with age-matched 46,XY controls. The upregulated genes screened out by combined analysis were mostly enriched in platelet activation, proteoglycans in cancer, chemokine signaling pathway, phagosome, cAMP, notch, and adherent junction, while the downregulated genes were enriched mainly in ribosome and TGF-β signaling pathway. Notably, proteoglycans and their modifying enzymes have been implicated in tumorigenesis in numerous cancers, including germ cell tumors (29–31). Moreover, previous studies have reported that adhesion-related genes, which are involved in cytoskeleton remodeling, actin dynamics, integrin signaling, and focal adhesion formation, were upregulated after suppression of the action of both testosterone and follicle-stimulating hormone (FSH) by facilitating aggregation/clumping of round spermatids together and formation of giant multinucleated germ cells (32, 33). Cell adhesion was also related to sexual immaturity (34). Additionally, phagosome, cAMP, TGF-β, and notch signaling pathways play crucial roles in sexual differentiation and the maintenance of gonad function (35–38). Furthermore, genes from the key modules were related mostly to immune function, including the chemokine signaling pathway, focal adhesion, ECM–receptor interaction, and PI3K–Akt signaling pathway. GO enrichment analysis showed that reduced expression was primarily enriched in DNA repair, translation, and meiotic division processes, which are closely correlated with development and spermatogenesis.
According to our data, the top TFs of upregulated genes (e.g., RNF2, FOXA1, BRD4, Gli2, 5HMC, and PBX3) and downregulated genes (e.g., TALE and MEF2C) were identified, most of which have been indicated to participate in regulating sex development (39–42). FOXA1 has been reported to act as a pioneer factor that opens up compact chromatin to facilitate the recruitment of AR. The expression of FOXA1 was observed to cooperate and imbalance the AR cistrome in regulating AR signaling in prostate cancer cells, which may form a mechanism compensating for the functional defects of AR, even in PBMCs with a low abundance of AR expression (43). Recent studies have demonstrated that Gli2 is a main regulatory gene mediating the hedgehog signaling pathway, which not only is a crucial factor for initial development but also plays a vital role in coordination with androgen signaling and in sexually dimorphic development of the external genitalia (44). Therefore, the alteration of Gli2 may represent a regulatory or compensatory mechanism in this disease, underlying the mechanism of hypospadias. Unlike other TFs, MEF2C plays a crucial role in myogenesis of cardiac, skeletal, and smooth muscle and is essential for the differentiation of muscle (45, 46). Current literature on PBX3 is known to be a functionally significant TF in a range of cancers, including prostate cancer, and linked to shorter overall survival; and the miR inhibitor of PBX3, miR-let-7d, is regulated by androgen (47). Using our preliminary analysis, we prioritized TFs that were likely to play a role in the pathophysiology of PAIS.
Then, BioGPS showed that hematologic/immune and reproductive/endocrine factors were significantly dysregulated, elucidating the common occurrence of immune dysfunction and sexual abnormalities in PAIS patients. The overexpression of immune system genes in the patient transcriptome may be a possible early indication of autoimmune disorders that could be involved in the onset of autoimmune diseases seen in AIS patients, including celiac disease (48) and primary and secondary Sjögren’s syndrome (49). Similarly, a higher prevalence of autoimmune diseases has been suggested in patients with male hypogonadism (50). Previous studies have shown that androgens have a suppressive effect on the immune system and play a key role in the development of immune cells (7). The TGF-β pathway was proven to regulate immune tolerance and homeostasis against self- and microbe antigens (51), also appearing in our functional enrichment results.
Combined BioGPS and PPI key module analyses discovered a number of novel genes that have been previously not been known to be related to this disorder. Three selected genes of interest (CCR1, PPBP, and PF4) were associated with the hematologic/immune system, among which PPBP and PF4 served as the hub genes. In an original study, CCR1 was shown to be a G protein-coupled receptor expressed by various leukocyte types that leads to inflammation and activation of robust cell infiltration, so it is considered to be an appropriate target for autoimmune and inflammatory therapeutics (52). Intriguingly, PPBP (CXCL7) and PF4 (CXCL4) belong to the CXC chemokine family, which plays a pivotal role in the development or progression of rheumatoid arthritis (RA), antiphospholipid syndrome, and systemic lupus erythematosus (53, 54). The level of PPBP was found to be significantly higher in vasculitis patients and female RA patients than in ordinary people (55, 56). PF4 protein functions as a chemoattractant for multiple types of cells, as well as an inhibitor of T-cell proliferation, angiogenesis, and hematopoiesis. Thus, we speculated that certain interrelations exist between AIS and the immune system. However, clinical data and basic research on these issues are scarce, and the immune response varies in knockout mice (T-ARKO or Tfm), which may not be comparable with AIS patients because of orchidectomy and supplementation with E2. In summary, whether AIS patients are susceptible to autoimmune diseases, inflammation, and infectious diseases should be studied further.
Impaired AR signaling can be responsible for an increased risk of metabolic diseases (10, 57, 58). The genes in reproduction/endocrine metabolism we found were CLU, KMT2D, GP6, and SPARC. CLU, thought to be involved in spermatogenesis, encoding clusterin, may increase in non-capacitated sperm to play a protective role in mice and bovines (59, 60). The protein encoded by KMT2D is histone H3 lysine 4-methyltransferase, which regulates muscle and adipose tissue development. Heterozygous Kmt2d+/− mice demonstrated improved obesity, lipid accumulation, glucose tolerance, and insulin sensitivity (61). GP6, as a new acute coronary syndrome marker encoding platelet glycoprotein VI, is essential for the formation of arterial thrombosis (62). SPARC, a profibrotic protein secreted by adipocytes, is implicated in the pathogenesis of inflammation, dyslipidemia, increased cardiovascular risk, obesity, and type 2 diabetes (63). These data indicated that dysregulation of KMT2D, GP6, and SPARC might be responsible for the increased risk of metabolic syndrome, diabetes, and cardiovascular diseases in AIS.
By manual screening of these DEGs, we showed that some robust genes were disordered in male reproductive development, such as spermatogenesis, sex differentiation, autophagy, aging, and androgen signaling pathways. Among these genes, CYP1B1, an enzyme belonging to the cytochrome P450 superfamily, has roles in the masculinization process that are regulated by androgens. During genital tubercle masculinization, CYP1B1 may function to metabolize androgens, such as testosterone. Consistent with Tanase-Nakao et al. (64), who reported in the transcriptomics of external genital fibroblasts of PAIS patients, CYP1B1 was also found to be elevated in the PBMCs of patients. As a potential androgen target gene, CYP1B1 may gain intriguing insights from the current study for new biomarkers in the diagnosis of AIS.
HSD17B3 encodes hydroxysteroid 17-beta dehydrogenase 3, which catalyzes the conversion of the more biologically active testosterone in the testis and is required for normal fetal development of the male genitalia. HSD17B3 expression has been reported in the Sertoli cells of the gonads in CAIS patients, which preserves the characteristics of the fetal gonad in the postpubertal testis in androgen synthesis (65). For the upregulation under elevated androgen levels in PAIS patients, it could be assumed that the expression of HSD17B3 might be modulated by the AR-targeted gene but not the product testosterone. The downregulated genes CHD5 and MEI1, encoding chromodomain helicase DNA binding protein 5, which regulates the histone-to-protamine chromatin remodeling process, and meiosis inhibitor protein 1, which affects male meiosis, respectively, play a spermatogenetic role (66, 67).
Previous studies aimed to identify potential biomarkers for the process of germ cell differentiation by using testicular tissue from CAIS patients (68) and mined predictors of genital development and AR function in genital skin fibroblasts from PAIS patients (64). A few genes have been identified as candidates to decipher androgen-dependent normal and abnormal genital development by genital fibroblasts (14, 69). ApoD has not been identified as a DEG due to its low abundance, which may result from the specific expression profile in different tissues. The limitation of the study is the small sample size, only included in genotype-proven PAIS patients. Mutation-negative AIS should be investigated as a separate group in future work (70).
In our study, we first performed a system-wide study of PAIS using PBMCs, with different compositions of DEGs in the results due to the specific expression profile of different tissues. However, the common DEGs (CYP1B1 and HSD17B3) could also be identified in PBMCs and in genital skin fibroblasts of PAIS patients, which could be easy biomarkers for diagnosis. Our approach provides valuable information for reproductive-related molecular biomarkers for screening and casts light on the understanding of AR and the possible pathogenesis of AIS.
Data Availability Statement
The original contributions presented in the study are publicly available. These data can be found here: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA767615.
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethical Committee of the Shanghai Ninth People’s Hospital, China. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YP performed the data analysis. YP and HZ contributed to the writing and revision of this manuscript. BH, YX, and LX enrolled the patients and sorted out the data. HS and JQ conceived and designed the experiments and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the National Natural Science Foundation of China (grant no. 81873652), the Natural Science Foundation of Shanghai (grant no. 21ZR1438300), and Clinical Research Project of Multi-Disciplinary Team, Shanghai Ninth People’s Hospital, Shanghai Jiaotong University School of Medicine (201903).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.731107/full#supplementary-material
References
1. Hughes IA, Davies JD, Bunch TI, Pasterski V, Mastroyannopoulou K, MacDougall J. Androgen Insensitivity Syndrome. Lancet (2012) 380(9851):1419–28. doi: 10.1016/S0140-6736(12)60071-3
2. Costagliola G, Cosci O, Masini B, Baldinotti F, Caligo MA, Tyutyusheva N, et al. Disorders of Sexual Development With XY Karyotype and Female Phenotype: Clinical Findings and Genetic Background in a Cohort From a Single Centre. J Endocrinol Invest (2021) 44(1):145–51. doi: 10.1007/s40618-020-01284-8
3. Berglund A, Johannsen TH, Stochholm K, Viuff MH, Fedder J, Main KM, et al. Incidence, Prevalence, Diagnostic Delay, and Clinical Presentation of Female 46, XY Disorders of Sex Development. J Clin Endocrinol Metab (2016) 101(12):4532–40. doi: 10.1210/jc.2016-2248
4. Audi L, Fernandez-Cancio M, Carrascosa A, Andaluz P, Torán N, Piró C, et al. Novel (60%) and Recurrent (40%) Androgen Receptor Gene Mutations in a Series of 59 Patients With a 46,XY Disorder of Sex Development. J Clin Endocrinol Metab (2010) 95(4):1876–88. doi: 10.1210/jc.2009-2146
5. Mendoza N, Rodriguez-Alcalá C, Motos MA, Salamanca A. Androgen Insensitivity Syndrome: An Update on the Management of Adolescents and Young People. J Pediatr Adolesc Gynecol (2017) 30(1):2–8. doi: 10.1016/j.jpag.2016.08.013
6. Wilson S, Qi J, Filipp FV. Refinement of the Androgen Response Element Based on ChIP-Seq in Androgen-Insensitive and Androgen-Responsive Prostate Cancer Cell Lines. Sci Rep (2016) 6:32611. doi: 10.1038/srep32611
7. Chang C, Yeh S, Lee SO, Chang TM. Androgen Receptor (AR) Pathophysiological Roles in Androgen-Related Diseases in Skin, Bone/Muscle, Metabolic Syndrome and Neuron/Immune Systems: Lessons Learned From Mice Lacking AR in Specific Cells. Nucl Recept Signal (2013) 11:e001. doi: 10.1621/nrs.11001
8. Jividen K, Kedzierska KZ, Yang CS, Szlachta K, Ratan A, Paschal BM. Genomic Analysis of DNA Repair Genes and Androgen Signaling in Prostate Cancer. BMC Cancer (2018) 18(1):960. doi: 10.1186/s12885-018-4848-x
9. Tyutyusheva N, Mancini I, Baroncelli GI, D'Elios S, Peroni D, Meriggiola MC, et al. Complete Androgen Insensitivity Syndrome: From Bench to Bed. Int J Mol Sci (2021) 22(3):1264. doi: 10.3390/ijms22031264
10. Dati E, Baroncelli GI, Mora S, Russo G, Baldinotti F, Parrini D, et al. Body Composition and Metabolic Profile in Women With Complete Androgen Insensitivity Syndrome. Sex Dev (2009) 3(4):188–93. doi: 10.1159/000228719
11. Marcus R, Leary D, Schneider DL, Shane E, Favus M, Quigley CA. The Contribution of Testosterone to Skeletal Development and Maintenance: Lessons From the Androgen Insensitivity Syndrome. J Clin Endocrinol Metab (2000) 85(3):1032–7. doi: 10.1210/jcem.85.3.6428
12. Danilovic DL, Correa PH, Costa EM, Melo KF, Mendonca BB, Arnhold IJ. Height and Bone Mineral Density in Androgen Insensitivity Syndrome With Mutations in the Androgen Receptor Gene. Osteoporos Int (2007) 18(3):369–74. doi: 10.1007/s00198-006-0243-6
13. La Spada AR, Wilson EM, Lubahn DB, Harding AE, Fischbeck KH. Androgen Receptor Gene Mutations in X-Linked Spinal and Bulbar Muscular Atrophy. Nature (1991) 352(6330):77–9. doi: 10.1038/352077a0
14. Holterhus PM, Deppe U, Werner R, Richter-Unruh A, Bebermeier JH, Wünsch L, et al. Intrinsic Androgen-Dependent Gene Expression Patterns Revealed by Comparison of Genital Fibroblasts From Normal Males and Individuals With Complete and Partial Androgen Insensitivity Syndrome. BMC Genomics (2007) 8:376. doi: 10.1186/1471-2164-8-376
15. Hellwinkel OJ, Holterhus PM, Struve D, Marschke C, Homburg N, Hiort O. A Unique Exonic Splicing Mutation in the Human Androgen Receptor Gene Indicates a Physiologic Relevance of Regular Androgen Receptor Transcript Variants. J Clin Endocrinol Metab (2001) 86(6):2569–75. doi: 10.1210/jcem.86.6.7543
16. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
17. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res (2000) 28(1):27–30. doi: 10.1093/nar/28.1.27
18. Maere S, Heymans K, Kuiper M. BiNGO: A Cytoscape Plugin to Assess Overrepresentation of Gene Ontology Categories in Biological Networks. Bioinformatics (2005) 21(16):3448–9. doi: 10.1093/bioinformatics/bti551
19. Wang H, Zhu H, Zhu W, Xu Y, Wang N, Han B, et al. Bioinformatic Analysis Identifies Potential Key Genes in the Pathogenesis of Turner Syndrome. Front Endocrinol (Lausanne) (2020) 11:104. doi: 10.3389/fendo.2020.00104
20. Qin Q, Fan J, Zheng R, Wan C, Mei S, Wu Q, et al. Lisa: Inferring Transcriptional Regulators Through Integrative Modeling of Public Chromatin Accessibility and ChIP-Seq Data. Genome Biol (2020) 21(1):32. doi: 10.1186/s13059-020-1934-6
21. Ulasov AV, Rosenkranz AA, Sobolev AS. Transcription Factors: Time to Deliver. J Control Release (2018) 269:24–35. doi: 10.1016/j.jconrel.2017.11.004
22. Hughes IA, Werner R, Bunch TI, Pasterski V, Mastroyannopoulou K, MacDougall J, et al. Androgen Insensitivity Syndrome. Semin Reprod Med (2012) 30(5):432–42. doi: 10.1016/S0140-6736(12)60071-3
23. Ahmed SF, Khwaja O, Hughes IA. The Role of a Clinical Score in the Assessment of Ambiguous Genitalia. BJU Int (2000) 85(1):120–4. doi: 10.1046/j.1464-410x.2000.00354.x
24. Lee SW, Kwak DS, Jung IS, Kwak JH, Park JH, Hong SM, et al. Partial Androgen Insensitivity Syndrome Presenting With Gynecomastia. Endocrinol Metab (Seoul) (2015) 30(2):226–30. doi: 10.3803/EnM.2015.30.2.226
25. Bouvattier C, Mignot B, Lefèvre H, Morel Y, Bougnères P. Impaired Sexual Activity in Male Adults With Partial Androgen Insensitivity. J Clin Endocrinol Metab (2006) 91(9):3310–5. doi: 10.1210/jc.2006-0218
26. Cools M, Drop SL, Wolffenbuttel KP, Oosterhuis JW, Looijenga LH. Germ Cell Tumors in the Intersex Gonad: Old Paths, New Directions, Moving Frontiers. Endocr Rev (2006) 27(5):468–84. doi: 10.1210/er.2006-0005
27. Mongan NP, Tadokoro-Cuccaro R, Bunch T, Hughes IA. Androgen Insensitivity Syndrome. Best Pract Res Clin Endocrinol Metab (2015) 29(4):569–80. doi: 10.1016/j.beem.2015.04.005
28. Shukla GC, Plaga AR, Shankar E, Gupta S. Androgen Receptor-Related Diseases: What Do We Know? Andrology (2016) 4(3):366–81. doi: 10.1111/andr.12167
29. Theocharis AD, Skandalis SS, Tzanakakis GN, Karamanos NK. Proteoglycans in Health and Disease: Novel Roles for Proteoglycans in Malignancy and Their Pharmacological Targeting. FEBS J (2010) 277(19):3904–23. doi: 10.1111/j.1742-4658.2010.07800.x
30. Rosen SD, Lemjabbar-Alaoui H. Sulf-2: An Extracellular Modulator of Cell Signaling and a Cancer Target Candidate. Expert Opin Ther Targets (2010) 14(9):935–49. doi: 10.1517/14728222.2010.504718
31. Fuster MM, Esko JD. The Sweet and Sour of Cancer: Glycans as Novel Therapeutic Targets. Nat Rev Cancer (2005) 5(7):526–42. doi: 10.1038/nrc1649
32. O’Donnell L, Pratis K, Wagenfeld A, Gottwald U, Müller J, Leder G, et al. Transcriptional Profiling of the Hormone-Responsive Stages of Spermatogenesis Reveals Cell-, Stage-, and Hormone-Specific Events. Endocrinology (2009) 150(11):5074–84. doi: 10.1210/en.2009-0755
33. Kavarthapu R, Anbazhagan R, Raju M, Morris CT, Pickel J, Dufau ML. Targeted Knock-in Mice With a Human Mutation in GRTH/DDX25 Reveals the Essential Role of Phosphorylated GRTH in Spermatid Development During Spermatogenesis. Hum Mol Genet (2019) 28(15):2561–72. doi: 10.1093/hmg/ddz079
34. Li Y, Li J, Fang C, Shi L, Tan J, Xiong Y, et al. Genome-Wide Differential Expression of Genes and Small RNAs in Testis of Two Different Porcine Breeds and at Two Different Ages. Sci Rep (2016) 6:26852. doi: 10.1038/srep26852
35. Itman C, Mendis S, Barakat B, Loveland KL. All in the Family: TGF-Beta Family Action in Testis Development. Reproduction (2006) 132(2):233–46. doi: 10.1530/rep.1.01075
36. Chen H, Wang Y, Ge R, Zirkin BR. Leydig Cell Stem Cells: Identification, Proliferation and Differentiation. Mol Cell Endocrinol (2017) 445:65–73. doi: 10.1016/j.mce.2016.10.010
37. Lasala C, Schteingart HF, Arouche N, Bedecarrás P, Grinspon RP, Picard JY, et al. SOX9 and SF1 Are Involved in Cyclic AMP-Mediated Upregulation of Anti-Mullerian Gene Expression in the Testicular Prepubertal Sertoli Cell Line SMAT1. Am J Physiol Endocrinol Metab (2011) 301(3):E539–547. doi: 10.1152/ajpendo.00187.2011
38. Tang EI, Lee WM, Cheng CY. Coordination of Actin- and Microtubule-Based Cytoskeletons Supports Transport of Spermatids and Residual Bodies/Phagosomes During Spermatogenesis in the Rat Testis. Endocrinology (2016) 157(4):1644–59. doi: 10.1210/en.2015-1962
39. Yokobayashi S, Liang CY, Kohler H, Nestorov P, Liu Z, Vidal M, et al. PRC1 Coordinates Timing of Sexual Differentiation of Female Primordial Germ Cells. Nature (2013) 495(7440):236–40. doi: 10.1038/nature11918
40. Hong SH, You JS. SOX9 Is Controlled by the BRD4 Inhibitor JQ1 via Multiple Regulation Mechanisms. Biochem Biophys Res Commun (2019) 511(4):746–52. doi: 10.1016/j.bbrc.2019.02.135
41. de Souza AF, Pieri NCG, Roballo KCS, Bressan FF, Casals JB, Ambrósio CE, et al. Dynamics of Male Canine Germ Cell Development. PloS One (2018) 13(2):e0193026. doi: 10.1371/journal.pone.0193026
42. Sonnet W, Rezsohazy R, Donnay I. Characterization of TALE Genes Expression During the First Lineage Segregation in Mammalian Embryos. Dev Dyn (2012) 241(11):1827–39. doi: 10.1002/dvdy.23873
43. Jin HJ, Zhao JC, Wu L, Kim J, Yu J. Cooperativity and Equilibrium With FOXA1 Define the Androgen Receptor Transcriptional Program. Nat Commun (2014) 5:3972. doi: 10.1038/ncomms4972
44. Miyagawa S, Matsumaru D, Murashima A, Omori A, Satoh Y, Haraguchi R, et al. The Role of Sonic Hedgehog-Gli2 Pathway in the Masculinization of External Genitalia. Endocrinology (2011) 152(7):2894–903. doi: 10.1210/en.2011-0263
45. Juszczuk-Kubiak E, Flisikowski K, Wicinska K. Nucleotide Sequence and Variations of the Bovine Myocyte Enhancer Factor 2C (MEF2C) Gene Promoter in Bos Taurus Cattle. Mol Biol Rep (2011) 38(2):1269–76. doi: 10.1007/s11033-010-0226-8
46. Dodou E, Xu S-M, Black BL. Mef2c Is Activated Directly by Myogenic Basic Helix-Loop-Helix Proteins During Skeletal Muscle Development In Vivo. Mech Dev (2003) 120(9):1021–32. doi: 10.1016/S0925-4773(03)00178-3
47. Morgan R, Pandha HS. PBX3 in Cancer. Cancers (Basel) (2020) 12(2):431. doi: 10.3390/cancers12020431
48. Walker MD, Zylberberg HM, Green PHR, Katz MS. Endocrine Complications of Celiac Disease: A Case Report and Review of the Literature. Endocr Res (2019) 44(1-2):27–45. doi: 10.1080/07435800.2018.1509868
49. Sullivan DA, Sullivan BD, Evans JE, Schirra F, Yamagami H, Liu M, et al. Androgen Deficiency, Meibomian Gland Dysfunction, and Evaporative Dry Eye. Ann N Y Acad Sci (2002) 966:211–22. doi: 10.1111/j.1749-6632.2002.tb04217.x
50. Jiménez-Balderas FJ, Tápia-Serrano R, Fonseca ME, Arellano J, Beltrán A, Yáñez P, et al. High Frequency of Association of Rheumatic/Autoimmune Diseases and Untreated Male Hypogonadism With Severe Testicular Dysfunction. Arthritis Res (2001) 3(6):362–7. doi: 10.1186/ar328
51. Sanjabi S, Oh SA, Li MO. Regulation of the Immune Response by TGF-Beta: From Conception to Autoimmunity and Infection. Cold Spring Harb Perspect Biol (2017) 9(6):a022236. doi: 10.1101/cshperspect.a022236
52. Cheng JF, Jack R. CCR1 Antagonists. Mol Divers (2008) 12(1):17–23. doi: 10.1007/s11030-008-9076-x
53. Patsouras MD, Sikara MP, Grika EP, Moutsopoulos HM, Tzioufas AG, Vlachoyiannopoulos PG. Elevated Expression of Platelet-Derived Chemokines in Patients With Antiphospholipid Syndrome. J Autoimmun (2015) 65:30–7. doi: 10.1016/j.jaut.2015.08.001
54. Xia W, Wu J, Deng FY, Wu LF, Zhang YH, Guo YF, et al. Integrative Analysis for Identification of Shared Markers From Various Functional Cells/Tissues for Rheumatoid Arthritis. Immunogenetics (2017) 69(2):77–86. doi: 10.1007/s00251-016-0956-4
55. Karatoprak C, Uyar S, Abanonu GB, Pehlevan SM, Okuroglu N, Demirtunc R. The Levels of Beta-Thromboglobulin in Female Rheumatoid Arthritis Patients as Activation Criteria. Rheumatol Int (2013) 33(5):1229–32. doi: 10.1007/s00296-012-2511-4
56. Wen M, Qu X, Zhu Y, Tian L, Shen Z, Yang X, et al. Exploration of Novel Biomarkers in Vasculitis by Integrated Bioinfomatic Approaches. J Comput Biol (2019) 26(12):1448–57. doi: 10.1089/cmb.2019.0176
57. Lin HY, Xu Q, Yeh S, Wang RS, Sparks JD, Chang C. Insulin and Leptin Resistance With Hyperleptinemia in Mice Lacking Androgen Receptor.pdf. Diabetes (2005) 54(6):1717–25. doi: 10.2337/diabetes.54.6.1717
58. Gava G, Mancini I, Orsili I, Bertelloni S, Alvisi S, Seracchioli R, et al. Bone Mineral Density, Body Composition and Metabolic Profiles in Adult Women With Complete Androgen Insensitivity Syndrome and Removed Gonads Using Oral or Transdermal Estrogens. Eur J Endocrinol (2019) 181(6):711–8. doi: 10.1530/EJE-19-0383
59. Saewu A, Kadunganattil S, Raghupathy R, Kongmanas K, Diaz-Astudillo , Hermo L, et al. Clusterin in the Mouse Epididymis: Possible Roles in Sperm Maturation and Capacitation. Reproduction (2017) 154(6):867–80. doi: 10.1530/REP-17-0518
60. Kumar S, Deb R, Singh U, Ganguly I, Mandal DK, Tyagi S, et al. Bovine Circadian Locomotor Output Cycles Kaput (CLOCK) and Clusterin (CLU) mRNA Quantitation in Ejaculated Crossbred Bull Spermatozoa. Reprod Domest Anim (2015) 50(3):505–9. doi: 10.1111/rda.12522
61. Froimchuk E, Jang Y, Ge K. Histone H3 Lysine 4 Methyltransferase KMT2D. Gene (2017) 627:337–42. doi: 10.1016/j.gene.2017.06.056
62. Bigalke B, Langer H, Geisler T, Lindemann S, Gawaz M. Platelet Glycoprotein VI: A Novel Marker for Acute Coronary Syndrome. Semin Thromb Hemost (2007) 33(2):179–84. doi: 10.1055/s-2007-969032
63. Kos K, Wilding JP. SPARC: A Key Player in the Pathologies Associated With Obesity and Diabetes. Nat Rev Endocrinol (2010) 6(4):225–35. doi: 10.1038/nrendo.2010.18
64. Tanase-Nakao K, Mizuno K, Hayashi Y, Kojima Y, Hara M, Matsumoto KI, et al. Dihydrotestosterone Induces Minor Transcriptional Alterations in Genital Skin Fibroblasts of Children With and Without Androgen Insensitivity. Endocr J (2019) 66(4):387–93. doi: 10.1507/endocrj.EJ18-0494
65. Jarzabek K, Koda M, Chrusciel M, Kanczuga-Koda L, Sobczynska-Tomaszewska A, Rahman NA, et al. Features of the Fetal Gonad in Androgen Synthesis in the Postpubertal Testis Are Preserved in Complete Androgen Insensitivity Syndrome Due to a Novel Genetic Splice Site Donor Variant in Androgen Receptor Gene Intron 1. J Steroid Biochem Mol Biol (2019) 193:105420. doi: 10.1016/j.jsbmb.2019.105420
66. Li W, Wu J, Kim SY, Zhao M, Hearn SA, Zhang MQ, et al. Chd5 Orchestrates Chromatin Remodelling During Sperm Development. Nat Commun (2014) 5:3812. doi: 10.1038/ncomms4812
67. Grelon M, Gendrot G, Vezon D, Pelletier G. The Arabidopsis MEI1 Gene Encodes a Protein With Five BRCT Domains That Is Involved in Meiosis-Specific DNA Repair Events Independent of SPO11-Induced DSBs. Plant J (2003) 35(4):465–75. doi: 10.1046/j.1365-313X.2003.01820.x
68. Malcher A, Jedrzejczak P, Stokowy T, Monem S, Nowicka-Bauer K, Zimna A, et al. Novel Mutations Segregating With Complete Androgen Insensitivity Syndrome and Their Molecular Characteristics. Int J Mol Sci (2019) 20(21):5418. doi: 10.3390/ijms20215418
69. Holterhus PM, Hiort O, Demeter J, Brown PO, Brooks JD. Differential Gene-Expression Patterns in Genital Fibroblasts of Normal Males and 46,XY Females With Androgen Insensitivity Syndrome: Evidence for Early Programming Involving the Androgen Receptor. Genome Biol (2003) 4(6):R37. doi: 10.1186/gb-2003-4-6-r37
Keywords: androgen insensitivity syndrome, RNA transcriptome, differentially expressed genes, reproduction, immune, metabolism
Citation: Peng Y, Zhu H, Han B, Xu Y, Liu X, Song H and Qiao J (2021) Identification of Potential Genes in Pathogenesis and Diagnostic Value Analysis of Partial Androgen Insensitivity Syndrome Using Bioinformatics Analysis. Front. Endocrinol. 12:731107. doi: 10.3389/fendo.2021.731107
Received: 26 June 2021; Accepted: 25 October 2021;
Published: 18 November 2021.
Edited by:
Ani Amelia Zainuddin, National University of Malaysia, MalaysiaReviewed by:
Hao Chen, Nantong University, ChinaRafael Loch Batista, Johns Hopkins University, United States
Copyright © 2021 Peng, Zhu, Han, Xu, Liu, Song and Qiao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Qiao, qiaoj2001@126.com; Huaidong Song, huaidong_s1966@163.com
†These authors have contributed equally to this work