- Research
- Open access
- Published:
Application and insights of targeted next-generation sequencing in a large cohort of 46,XY disorders of sex development in Chinese
Biology of Sex Differences volume 15, Article number: 73 (2024)
Abstract
Purpose
46,XY disorders of sex development (46,XY DSD) are characterized by incomplete masculinization of genitalia with reduced androgenization. Accurate clinical management remains challenging, especially based solely on physical examination. Targeted next-generation sequencing (NGS) with known pathogenic genes provides a powerful tool for diagnosis efficiency. This study aims to identify the prevalent genetic variants by targeted NGS technology and investigate the diagnostic rate in a large cohort of 46,XY DSD patients, with most of them presenting atypical phenotypes.
Methods
Two different DSD panels were developed for sequencing purposes, targeting a cohort of 402 patients diagnosed with 46,XY DSD, who were recruited from the Department of Urology at Children’s Hospital, Zhejiang University School of Medicine (Hangzhou, China). The detailed clinical characteristics were evaluated, and peripheral blood was collected for targeted panels to find the patients’ variants. The clinical significance of these variants was annotated according to American College of Medical Genetics and Genomics (ACMG) guidelines.
Results
A total of 108 variants across 42 genes were found in 107 patients, including 46 pathogenic or likely pathogenic variants, with 45.7%(21/46) being novel. Among these genes, SRD5A2, AR, FGFR1, LHCGR, NR5A1, CHD7 were the most frequently observed. Besides, we also detected some uncommon causative genes like SOS1, and GNAS. Oligogenic variants were also identified in 9 patients, including several combinations PROKR2/FGFR1/CYP11B1, PROKR2/ATRX, PROKR2/AR, FGFR1/LHCGR/POR, FGFR1/NR5A1, GATA4/NR5A1, WNT4/AR, MAP3K1/FOXL2, WNT4/AR, and SOS1/FOXL2.
Conclusion
The overall genetic diagnostic rate was 11.2%(45/402), with an additional 15.4% (62/402) having variants of uncertain significance. Additionally, trio/duo patients had a higher genetic diagnostic rate (13.4%) compared to singletons (8.6%), with a higher proportion of singletons (15.1%) presenting variants of uncertain significance. In conclusion, targeted gene panels identified pathogenic variants in a Chinese 46,XY DSD cohort, expanding the genetic understanding and providing evidence for known pathogenic genes’ involvement.
Plain English summary
46,XY disorders of sex development (46,XY DSD) are conditions where individuals don’t fully develop male genitalia due to reduced androgen hormones. Diagnosing these conditions based only on physical exams is difficult. This study used advanced genetic testing called targeted next-generation sequencing (NGS) to identify common genetic variations in a large group of 46,XY DSD patients, many of whom had unusual symptoms. We examined 402 patients with DSD and a 46,XY karyotype, focusing on 142 candidate genes related to sex development. We found genetic variations in 107 patients, including 45 that were likely responsible for their condition. Some of these variations were new discoveries. The most commonly affected genes were SRD5A2, AR, FGFR1, LHCGR, NR5A1, CHD7. We also found that some patients had variations in multiple genes, suggesting complex genetic causes. Overall, we were able to diagnose 11.2% of patients based on our genetic testing, with another15.4% having uncertain results. Patients tested as a trio or duo (with their parents) had a higher diagnosis rate than those tested alone. This study helps expand our understanding of the genetic factors behind 46,XY DSD in the Chinese population.
Highlights
Two gene panels were designed, comprising 142 and 271 candidate genes associated with sex development, to sequence 402 patients with 46,XY DSD.
A total of 108 variants across 42 genes were identified in 107 patients, with 46 classified as pathogenic or likely pathogenic, including several novel variants. SRD5A2, AR, FGFR1, LHCGR, NR5A1, CHD7 were among the most frequently observed genes with variants.
The overall genetic diagnostic rate was 11.2%, with an additional 15.4% of patients having variants of uncertain significance, and oligogenic variants were detected in 9 patients, suggesting complex genetic interactions in 46,XY DSD.
Patients with trio/duo genetic information had a higher genetic diagnostic rate (13.4%) compared to singletons (8.6%), underscoring the importance of parental genetic data.
Phenotypic variability and potential genetic heterogeneity underscore the complexity of 46,XY DSD, highlighting the necessity for further research and multi-center collaboration.
Introduction
Disorders of sexual development (DSDs) are a diverse group of congenital complex conditions in which the development of chromosomal, gonadal, or anatomical sex is atypical [1]. The exact incidence of DSD is not clear; but when considering all congenital genital anomalies, it can be as high as 1 in 200 to 300 births [2]. Signs of DSD can manifest at birth with atypical external genital anatomy, in childhood with bilateral inguinal hernias, in puberty with atypical secondary sex characteristics, or in adulthood with infertility [3]. Individuals with DSD often experience challenges related to self-esteem and gender identity [4]. Additionally, there is a significantly higher risk of gonadal tumors in patients with DSD compared to the general population [5]. The distress associated with overall health concerns, as well as the shame and fear of social stigma for individuals and their families, often lead to social withdrawal and isolation, consequently bring a profound psychopathological burden on them. Therefore, early diagnosis and intervention are necessary for the medical management of DSD. However, due to the complexity and heterogeneity of these conditions, clinical management of DSD can be challenging, interdisciplinary teams are needed.
According to the DSD international consensus, DSDs are classified into 46,XY, 46,XX, and sex chromosome DSD [1]. Among them, 46,XY DSD is the most common and complex type. It is characterized by female or ambiguous external genitalia such as hypospadias, cryptorchidism, micropenis, and testicular/gonadal dysgenesis, resulting from incomplete virilization with or without the presence of Mullerian structures [6]. Considering the aetio-pathogenesis of process of fetal sex differentiation, 46,XY DSD may be endocrine-related or not [7]. Non-endocrine disorders, such as isolated hypospadias, can arise from abnormal morphogenesis of testicular primordia [7]. On the other hand, endocrine-related 46,XY DSD can be attributed to early-onset gonadal failure, specifically affecting Leydig or Sertoli cell functions, or defects in male hormone at target tissues during embryogenesis, puberty, and adulthood [3]. Impaired secretion of anti-mullerian hormone (AMH) from Stertoli cells will lead to the failure of Mullerian ducts regression. Similarly, decreased androgen production such as testosterone (T) or dihydrotestosterone (DHT) from Leydig cells can impact the development of Wolffian ducts, influencing the formation of epidymides, vasa deferentia, and seminal vesicles [8].
Mutational and functional analyses of patients with DSDs and mouse models have revealed a large number of genes involved in male differentiation [9].These genes can be categorized according to their functions as follows: (1) genes related to the development and diferentiation of gonadal, such as SRY, NR5A1, WT1, MAP3K1, GATA4, and MAMLD1. (2) genes associated with hormone synthesis and action, including steroid hormone and androgen, such as AR, SRD5A2, AKR1C2, CHD7, LHCGR, ANOS1, and CYP11A1. In addition to single gene mutations, copy number variations (CNVs) have also been shown to be associated with DSD. For instance, deletions of 9p24.3 and 10q26.1, as well as duplications of Xp21.2 and 1p35, have been identified [10].
The clinical phenotype of 46,XY DSD patients often presents similar manifestations, making accurate diagnose challenging based solely on physical examination. However, having knowledge of the pathogenic genes involved in these patients can enhance diagnosis efficiency, potentially reducing the need for additional expensive biochemical and radiological assessments. The decreasing costs of Targeted Next-Generation Sequencing (NGS) have rendered this approach increasingly viable in clinical practice. The primary objective of this study was to identify gene variants in a large cohort of 46,XY DSD patients and specifically investigate any unique gene variants or hotspots within the Chinese population, and reveal the prevalent variants of 46,XY DSD by comparing with previous studies. By doing so, we aimed to contribute substantial evidence to facilitate early diagnosis and intervention for affected individuals.
Methods
Patient recruitment and clinical assessment
We conducted a retrospective analysis involving 402 patients suspected of DSD at the Department of Urology, Children’s Hospital, Zhejiang University School of Medicine, between 1/2017 and 12/2020, for further investigation. All patients were subjected to detailed clinical and genital examinations, family history pedigree analysis, and the assessment of associated abnormalities. The patient inclusion criteria were as follows: (1) patients with 46,XY karyotype and (2) patients with external genital malformation or gonadal dysplasia, including female external genitalia, clitoromegaly, ambiguous external genitalia, perineal hypospadias, cryptorchidism, and micropenis.
The study was ethically approved by the Human Subjects Committees of Children’s Hospital, Zhejiang University School of Medicine(2018-IRB-076). All procedures performed in studies involving human participants were in accordance with the ethical standards of the Ethical Committee. An informed consent was obtained from the patients or their guardians. Written informed consent was obtained from the patients, the patients’ parents, or their legal guardians.
Targeted gene panel
Two gene panels were designed for screening 46,XY DSD, comprising 142 candidate genes (M014) and 271 genes (KY043), respectively. Genes included in these panels were selected based on the knowledge of DSD sourced from PubMed, OMIM, and other genetic testing registry databases. These genes were classified into several categories, including synthesis or action of androgen, development and differentiation of gonadal, synthesis and activation of Steroid hormone, syndromic disorders, development and differentiation of germ cell, and others) (Supplementary Tables 3 & 4).
DNA library preparation
Genomic DNA was extracted from peripheral blood leukocytes using the Qiagen DNA Blood kit (Qiagen, Dusseldorf, Germany) from each sample. A minimum of 3 µg DNA was fragmented to an average size of 180 bp using a Bioruptor sonicator (Diagenode). Paired-end sequencing libraries were then prepared using a DNA sample prep reagent set 1 (NEBNext). The library preparation was followed by the recommended protocols from Illumina, the process involved end repair, adapter ligation, and polymerase chain reaction (PCR) amplication.
Targeted genes enrichment and sequencing
The exon regions and exon-intron boundaries of the target genes were captured using GenCap customized DSD Kit (MyGenostics Inc. Beijing, China) according to previously described methods [11]. In brief, 1 µg DNA library was mixed with Buffer BL and GenCap DSD probe (MyGenostics, Beijing, China). The mixture was then heated at 95 °C for 7 min, followed by 65 °C for 2 min in a PCR machine. Subsequently, 23 µl of the 65 °C pre-warmed Buffer HY (MyGenostics Inc, Beijing, China) was added, and the mixture was maintained at 65 °C with the PCR lid heat on for 22 h for hybridization. After that, 50 µl MyOne beads (Life Technology) underwent washing in 500µL 1X binding buffer three times and were then resuspended in 80 µl 1X binding buffer. Following this, 64 µl of 2X binding buffer was added to the hybrid mixture, and the entire mixture was transferred into a tube containing 80 µl MyOne beads. The mixure was rotated for 1 h on a rotator. The beads were then washed with WB1 buffer at room temperature for 15 min once and with WB3 buffer at 65 °C for 15 min three times. Afterward, the bound DNA was eluted with Buffer Elute. The eluted DNA was finally amplified for 15 cycles using the following program: 98 °C for 30 s (1 cycle); 98 °C for 25 s, 65 °C for 30 s, 72 °C for 30 s (15 cycles); 72 °C for 5 min (1 cycle). The PCR product was purified using SPRI beads (Beckman Coulter) according to manufacturer’s protocol. The enriched libraries were sequenced on Illumina HiSeq X ten sequencer with paired read of 150 bp.
Variant calling and annotation
Paired-end sequence reads in fastq files were filtered to remove low-quality reads by using Fastp [12] with default parameters. After ensuring quality control, the clean reads were aligned to the UCSC hg19 human reference genome with Burrows-Wheeler Aligner Maximal Exact Match (BWA-MEM) [13]. Duplicated reads were marked with Picard tools [14]. Variants of Single Nucleotide Variations (SNVs) and Insertions/deletions (Indels) were called using Genome Analysis Toolkit (GATK 4.1.7.0) [15] HaplotypeCaller to generate gVCF files for joint genotyping. All samples were jointly genotyped and variants were selected with GATK Selectvariants for hard filtering with the following criteria; (a) Quality by depth < 2; (b) Mapping Quality < 40; (c) approximate read depth < 6; (d) phred-scaled p-value using Fisher’s exact test to detect strand bias (FS) > 60 for SNVs and > 200 for Indels. ANNOVAR [16] was used to annotate variant function, in silico predictions of deleteriousness, and variant populations frequencies.
Pathogenicity analysis of the variants
Variants were classified as pathogenic (P), likely pathogenic (LP), variant of uncertain significance (VUS), likely benign, or benign according to the American College of Medical Genetics and genomics guidelines [17]. Pathogenic and likely pathogenic variants were prioritized if indicated in Clinvar, Human Gene Mutation Database (HGMD, 2022.12), or had been reported to be associated with DSD. The remaining variants were filtered and prioritized based on MAF less than 0.1% in all population datasets, and function annotation of missense, inframe, frameshift indels, canonical splice site, or nonsense variants. The pathogenicity of variants was further calibrated using Varsome [18]. Variants identified as pathogenic, likely pathogenic, or VUS were verified by Sanger sequencing and further confirmation of origin with available parental samples, benign and likely benign variants were discarded.
Results
Sequencing quality and cohort characteristics
The Average sequencing depth on target for panel M014 and KY043 is 265× (range from 127× to 660×) and 461× (range from 165× to 1316×), respectively (Supplementary Fig. 1A). More than 76% of targeted regions were covered with greater than 20× for each sample in M014 panel, while more than 90% of targeted regions were covered with greater than 20× for samples in KY043 panel (Supplementary Fig. 1B).
Total of 402 patients with a 46,XY karyotype and genetic testing result have been recruited for analysis. Among them, 185 are singletons, 30 are duo samples, and 187 are trio-based samples. The mean age at diagnosis was 3.4 ± 3 years old, specifically, 54 patients were first assessed at minipubertal age (< 6 months), 294 patients at prepubertal age (6 months-9 years), 54 patients at pubertal age (9–18 years). Majority of patients are Han Chinese, and 7 of them are minority Chinese.
As listed in Table 1, the patients presented with a variety of DSD phenotypes. Seventeen patients have a determined etiologic diagnosis. Among the remaining patients with unknown etiology, the most common clinical manifestation was the micropenis and hypospadias, accounted for 50.6% (205/402) and 22.5% (91/402), respectively. Three of those 402 patients were raised as females because of ambiguous genitalia and need for male sex assignment after genetic diagnosis and karyotype (supplementary Table 1).
Variants analysis
After filtering by the ACMG guideline, a total of 108 variants involved 42 genes were positively detected in 107 patients (Supplementary Table 1). These genes were distributed all around the chromosomes (Fig. 1A). Furthermore, variants were mainly enriched in genes associated with the synthesis or action of androgen (28.6%), development and differentiation of gonadal (28.6%), and synthesis and activation of steroid hormone (19.1%) (Fig. 1B). For the 108 variants, 14.8% (16/108) were pathogenic variants, 27.8% (30/108) were likely pathogenic variants, and the remaining 57.4% (62/108) were VUS (Fig. 1C). Among these variants, missense variants were the most common type accounting for 79.6% (86/108), followed by frameshift variants accounting for 7.4% (8/108), and nonsense variants for 6.5% (7/108) (Fig. 1D). Additionally, four splicing variants and two nonframeshift insertion/deletion variants were identified. One patient had a VUS synonymous variant compound with a pathogenic missense variant in SRD5A2 gene (Fig. 1D, supplementary Table 1).
The clinical and molecular characteristics including main identified variants, clinical features, karyotypes, and transmission of positive cases are summarized in Table 1. And we visualized patients- genes information in patients in Fig. 2A. Information of the negative cases and the details of the evidence level applied appropriately are provided in Supplementary Table S2. SRD5A2 was the most frequent gene, including 2 novel variants and 7 reported variants (Fig. 2B) occurring in 15 patients (Table 2).
Oligogenic variants
In addition to monogenic variants, oligogenic variants were also identified in several patients. After excluding the variants without interaction or combinations with low pathogenicity by ORVAL(Oligogenic Resource for Variant AnaLysis (ibsquare.be), we finally identified 7 patients with digenic variants and 2 patients with trigenic variants, involving various gene combinations and including PROKR2/FGFR1/CYP11B1, PROKR2/ATRX, PROKR2/AR, FGFR1/LHCGR/POR, FGFR1/NR5A1, GATA4/NR5A1, WNT4/AR, MAP3K1/FOXL2, and SOS1/FOXL2 (Table 3). All 9 patients presented with phenotypes of micropenis, cryptorchidism, or hypospadias.
Molecular genetic diagnosis
An overall diagnostic rate was 11.2% (45/402) with pathogenic and likely pathogenic variants. There were another 15.4% of patients with VUS variants (Fig. 3A). Among those 45 positive cases, the genetic diagnostic rate was over 50% in patients who have determined etiologic diagnoses with 5α-reductase 2 deficiency and GD. For those with clinically unknown etiology, the number of patients with isolated phenotype of micropenis and hypospadias was 205 and 91, respectively; and the diagnostic rate was 6.8% (14/205) and 14.3% (13/91) (Fig. 3B). Among those trio/duo patients, 13.4% (29 of 217) of patients have been determined as genetic diagnostic patients with P/LP variants; while in singletons,8.6% (16 of 185) of patients were with genetic diagnostic findings. The proportion of VUS in singleton patients (28 of 185, 15.1%) was similar to that in trios (33 of 217, 15.2%) (Fig. 3C).
Discussion
This study conducted a retrospective analysis in patients devolved to 46,XY DSD from 2017 to 2020, and two panels were designed for investigating genetic variants. In contrast to other research, our targeted panel encompasses a broader spectrum of candidate genes, including those crucial for sex development regulation. However, we have still missed some latest genes associated with DSD, which is a flaw in our design. As a result, 108 variants involved 42 genes were identified in 107 patients, with 46 of these variants being pathogenic or likely pathogenic. The overall diagnostic rate was 11.2% (45/402) and excluded the carriers pathogenic or likely pathogenic variants. Among these variants, SRD5A2 and AR were the most frequently affected genes, followed by FGFR1, LHCGR, NR5A1 and CHD7. These findings are basically consistent with previous studies conducted in China [19,20,21].
Notably, in our study, a majority of patients exhibited various DSD-related phenotype without a clear etiology. This observation is likely attributed to the fact that many patients opt for genetic diagnosis only when the cause of the disease is not readily apparent. However, obtaining a genetic diagnosis is of importance, not only for elucidating the etiology but also for guiding patient management, including considerations related to potential gender development, assessment of adrenal and gonadal function, and evaluating the risk of sexual adenocarcinoma, associated morbidity, and long-term outcomes [22, 23]. Previous research has indicated that the risk of germ cell tumors in patients with gonadal dysgenesis can reach up to 30%, while in AIS patients, it is approximately 15% [24]. Therefore, it is important for 46,XY DSD patients to identify the etiology as early as possible.
5α-reductase 2 deficiency is an autosomal recessive inherited disease caused by homozygous or compound heterozygous variants of SRD5A2 gene. In our study, we identified 12 different SRD5A2 variants in 17 unrelated 46,XY DSD patients. Of these, only three variants were novel, one of which was a synonymous variant found in DSD-0177 patient, who also harbored another SRD5A2 variant p.R227Q. Synonymous variants are typically regarded as benign, yet they may exert pathogenic effects by affecting mRNA stability [25]. Further functional assays are needed to ascertain the pathogenicity of such variants. Notably, positions 6, 196, 203, 227, 235 and 246 are recognized hotspots within the SRD5A2 gene [26, 27]. Among our patient cohort, SRD5A2 R227Q was the most prevalent variant, followed by R246Q. In a previous study involving 190 individuals from diverse regions in China, the prevalence of R227Q in south China was found to be 62.62% [28]. In our cohort, approximately 64.7% (11/17) of patients carried R227Q, with 5 patients being homozygosity and 6 being compound heterozygosity. Nonetheless, no distinct genotype–phenotype correlations were observed among these patients. Interestingly, two patients (DSD-0249 and DSD-0278) diagnosed with micropenis carry only one heterozygous variant of SRD5A2 gene, with no additional pathogenic variants identified, which cannot explain the phenotype. This is likely due to the panel’s incapacity to detect copy number variations (CNVs) or intronic region variants. Whole Exome Sequencing (WES) or Whole Genome Sequencing (WGS) would be valuable to uncover additional SRD5A2 variants or the other causative genetic factors in these patients. DSD-0251 was provisionally diagnosed with 5α-reductase 2 deficiency based on the assessment of T/DHT level, but no SRD5A2 variant was detected. And endocrinological tests are not always reliable for diagnosing 5α-reductase 2 deficiency because T/DHT ratios following hCG stimulation vary according to the age and the severity of the enzyme defects [29, 30]. The diagnostic sensitivity of T/ DHT ratio of 10 is around 78%, but the specificity is only 72% [31]. Moreover, not all patients with 5-α reductase 2 deficiency diagnosed by clinical and hormonal findings were found to carry pathogenic variants in SRD5A2 [30, 32]. Compared to those with genetically diagnosed 5-α reductase 2 deficiency, the average age is much younger in undiagnosed patients, at 1.3 ± 2.8 years old [30], consistent with our patient who was at 4 month old at diagnosis. This highlights the essential role of genetic testing for definitive diagnosis, especially in infants with clinical-hormonal diagnosis.
AR variants have primarily been associated with AIS. The AR protein consists of four functional domains: N-terminal domain (NTD), DNA binding domain (DBD), the C-terminal ligand binding domain (LBD), and a hinge region (HR) linking LBD and DBD [33]. In our cohort, 7 variants including two novel variants and five reported variants were identified. Notably, p.V731M has been recurrently observed as a somatic variant in prostate cancer patients [34,35,36], suggesting a potential association with an increased risk of testicular germ cell tumor (TGCT) in AIS patients [37]. While G444A was previously identified in colorectal cancer patients as a germline variant [38], it was first detected in a DSD patient in our study. Additionally, M887V and H886Y, previously reported in AIS patients [39, 40], were detected in DSD-0077, presenting ambiguous genitalia, and DSD-0292, displaying hypospadias. However, there were no AR or other gene variants detected in our three patients diagnosed with AIS. Interestingly, most patients in our cohort carrying AR variants exhibited no additional signs of androgen insensitivity beyond hypospadias and micropenis. This suggests that AR variants may be related to milder phenotypes, consistent with findings from previous studies [41, 42]. Although R841H has been identified as a hotspot variant in a Chinese cohort with a frequency of 5.6% (3/54) [20], it was not detect in our cohort.
CHD7 is associated with CHARGE syndrome and hypogonadotropic hypogonadism-5 with or without anosmia [43]. Loss-of-function pathogenic variants in the CHD7 gene account for approximately 65–70% of CHARGE syndrome cases [44]. However, not all patients carrying with CHD7 variants develop classical CHARGE syndrom phenotype. In a study of 40 patients with idiopathic hypogonadotropic hypogonadism (IHH), two adult patients with cryptorchidism, abnormal testicular, and/or abnormal penis were identified to carry CHD7 variants [45]. Another patient with middle hypospadias and right cryptorchidism and without other malformations was also found to carry a missense variant of CHD7, although the correlation between the phenotype and the variant could not be established [46]. Patients carry CHD7 variants may also present with microphallus, cryptorchidism, and hypospadias without other malformations [21]. There is also a missense variant in CHD7 reported in a patient with pure 46,XY gonadal dysgenesis in clinvar (VCV001344532.1 ). In our study, we identified CHD7 variants in 3 patients with the phenotype of micropenis, cryptorchidism, or gonadal dysgenesis, no further clinical phenoypes were recorded. All these reports suggested that patients with CHD7 variants may also present atypical CHARGE syndrome. However, the variant p.G1684S has been widely associated with CHARGE syndrome [47,48,49,50] .The other two variants are novel truncating CHD7 variants reporting as de novo inheritance, a form of which was frequently identified in CHARGE syndrome than the other atypical CHARGE features [51]. Further clinical evaluations are needed for these patients. Unfortunately, due to lost follow-up for these three patients, we could not further assess the clinical manifestations. Patients with IHH typically exhibit no symptoms before puberty and it is particularly difficult to evaluate the hypothalamic-pituitary‐gonadal (H‐P‐G) axis inchildren [21], making it easy for their condition to be overlooked during childhood [52, 53]. In our cohort, all patients with variants associated with IHH exhibited no other malformations. These patients may require regular follow-up in the future. Furthermore, genetic testing is helpful for monitoring the prospective puberty.
In our cohort, only three MAP3K1 variants were identified and all were VUS, which is consistent with findings in other regions of China [20] and Korea [54]. However, this incidence is much lower than what has been reported in Caucasians, where the prevalence of MAP3K1 gene variants ranges from 15 to 20% [55]. These disparities suggest the genetic heterogeneity in European and Asian populations. In recent years, GATA4 variants have been identified in patients with congenital heart disease(CHD) with various complications. GATA4 variants have also been found in patients with 46,XY DSD due to impaired testis formation with or without CHD [56]. In this cohort, we observed three patients carrying a heterozygous GATA4 variant P163S. Among them, one exhibited symptoms of atrial septal defect in addition to DSD, while the others did not show any CHD- associated symptoms. The GATA4 P163S variant has been previously reported and linked to tetralogy of Fallot [57]. The presence of multiple phenotypes among GATA4 variants carriers may be attributed to incomplete penetrance, variable interactions with partner genes, and oligogenic mechanisms. Nonetheless, it is important to note that patients without CHD phenotypes also require diligent monitoring of their cardiac health in the future.
In this study, we also identified several relatively uncommon causative genes associated with DSD, including TP63 and GNAS. Notably, TP63 has not been previously reported in DSD patients, with only a limited number of studies indicating its association with ovarian development in TP63 knockout mice [58, 59]. GNAS is a complex imprinted gene locus and encodes multiple transcripts sharing exons 2–13 but with alternative first exons [60]. Genetic defects affecting GNAS can cause several human diseases such as Pseudopseudo hypoparathyroidism (PHP, MIM:612462, 603233, 103580 ), end-organ resistance to parathyroid hormone (PTH, MIM: 612463), progressive osseous heteroplasia (POH, MIM: 166350). Somatic variants can also found in multiple benign and malignant tumors and McCune-Albright syndrome (MIM: 174800). The GNAS variant p.G142X we identified is located in the first exon of XLas (NM_080425.3) transcript, which is primarily expressed in neuroendocrine tissues and essential for normal fetal growth and development [61, 62]. This variant was initially reported in a patient with high bone mass, unclosed cranial suture, persistent hypophosphatemia, and elevated parathyroid hormone (PTH) levels [61]. However, genital development was not assessed in this family. The variant was inherited from the patient’s father and transmitted to his daughter. Both the proband and his father experienced bone fractures, while his 12-year-old daughter showed no symptoms. Both of our patients, DSD-0266 and DSD-0322, carry the same nonsense variant p.G142X. Recently, patient DSD-0266 presented with bone and joint instability and growth retardation (April,2024). He had elevated alkaline phosphatase (ALP) levels and an estradiol level of 124 pmol/L, but normal with thyroid function tests. Unfortunately, patient DSD-0322 was lost to follow-up. The cases of our patients, along with the previous reported cases, expand the spectrum of clinical manifestations associated with GNAS pathogenic variants. Interestingly, the absence of symptoms in the girl carrying this variant suggests the p.G142X variant in XLas may affect males exclusively. Further investigations will be conducted to explore the role of this variant.
In addition to monogenic variants, it is crucial to recognize the potential involvement of oligogenic variants in cases where patients exhibit negative results. Previous studies underscore that fewer than 50% of individuals with DSD receive a definitive genetic diagnosis, prompting the exploration of oligogenic diseases [56]. We have identified 9 patients with multiple variants with potential interaction and disease causing pathogenicity. These patients present with DSD phenotypes such as micropenis, cryptorchidism, hypospadias, or combination of these conditions. However, even when carrying the same genetic variant, the phenotype of each patient varies, suggesting other genetic modifiers may influence the unique phenotype of each patient. Furthermore, it is suggested that the greater the number of affected genes and variants an individual harbors, the more severe or variable the phenotype is likely to be [63, 64]. We found four patients had known pathogenic variants in PROKR2 combined with other DSD gene variants and one patient with monoallelic PROKR2 variant. Among them, three patients habored the most common PROKR2 W178S variant in the Chinese CHH population [65]. Although all three present with micropenis, patient DSD-0115 with ATRX P1548L had right cryptorchidism, laterla hernia, and hyposmia, while DSD-0260 with both ATRX and AR variants didn’t show any other phenotypes in recorded. Interestingly, patient DSD-0110 had both heterozygous PROKR2 W178S variant and known causative homozygous SRD5A2 p.R227Q variant, but no interaction was found in these two variants and the patient presented micropenis phenotype only. In contrast, the patient DSD-0016 with a monoallelic PROKR2 variant presented with more severe phenotypes, including bilateral cryptorchidism, growth retardation, and facial dysmorphia. This is likely due to undiscovered genetic variants outside the panels we tested, which have pleiotropic roles that interact with the monoallelic PROKR2 variant. The correlation between digenicity and severe phenotype was not always consistent [66]. This aligns with findings indicating that members of digenic families often had milder or similar phenotypes compared to probands carrying variants in the FGFR1 gene and other genetic variants [66]. Our patients did not exhibit significant phenotypic differences whether they carried monoallelic FGFR1 variants alone or in combination with other heterozygous variants in both autosomal dominant and recessive genes. This observation was also noted in patients carrying NR5A1 variants alone or combined with other genetic variants. Nevertheless, the genotype-phenotype correlation in 46,XY DSD may largely depend on the nature of the secondary hit, explaining the broad spectrum of phenotypes [67]. However, we couldn’t establish a phenotype-genotype correlation, several factors may contribute to this, including the limited sample size with diverse genetic combinations hindering correlation efforts, and the absence of comprehensive clinical descriptions for all patients.
Limitation
In our study, we observed a diagnostic rate of approximately 11.2%, and the rate is higher in trio/due patients (13.4%) than in singletons (8.6%). This rate is much lower than that observed in other studies conducted on the Chinese populations or other populations, where diagnostic rate exceeded 40%, especially in cohort using WES [19, 20, 68, 69]. Despite utilizing two panels with slightly different numbers of genes, our findings were limited by the absence of the latest DSD-related genes in the panels’ lists. WES which currently is the most effective and cost-efficient methods for identifying novel causes in undiagnosed DSD patients is the preferable tool. Less stringent patients’ selection criteria and the lack of clear differentiation of phenotype in patients may also limit our findings. Most of other groups are focusing on specific and severe phenotypes of 46,XY DSD such as 5α-reductase 2 deficiency, GD, AIS, HH[28, 70–72]. A multidisciplinary team and longitudinal follow-up should be conducted in recording physical examinations and biochemical (hormonal) evaluations. Combing these with genetic testing will improve the accuracy of genetic diagnoses. Additionally, it’s worth noting that we did not detect DNA copy number variations which may also contribute the phenotypes of DSD [73] due to the limitations of our gene panel. Furthermore, due to the lack of further validation, patients identified with VUS could not be conclusively diagnosed especially in singletons.
Conclusion
In summary, we utilized the targeted gene panels to identify potential pathogenic gene variants in a large cohort of Chinese 46,XY DSD patients. This study represents the general investigation of such a large cohort of 46,XY DSD in eastern China, contributing to our understanding of the gene spectrum, mutation spectrum, and phenotypic spectrum of 46,XY DSD in this region. However, it’s important to acknowledge potential limitations in our study, particularly since all patients were from a single center and the gene panels used were relatively outdated. Therefore, future research should involve multi-center collaboration to further advance our knowledge in this field.
Data availability
No datasets were generated or analysed during the current study.
References
Hughes IA, Houk C, Ahmed SF, Lee PA, Group LC, Group EC. Consensus statement on management of intersex disorders. Arch Dis Child. 2006;91:7554–63. https://doi.org/10.1136/adc.2006.098319.
Nordenvall AS, Frisen L, Nordenstrom A, Lichtenstein P, Nordenskjold A. Population based nationwide study of hypospadias in Sweden, 1973 to 2009: incidence and risk factors. J Urol. 2014;191:3:783–9. https://doi.org/10.1016/j.juro.2013.09.058.
Wisniewski AB, Batista RL, Costa EMF, Finlayson C, Sircili MHP, Denes FT, et al. Management of 46,XY Differences/Disorders of Sex Development (DSD) throughout life. Endocr Rev. 2019;40:6:1547–72. https://doi.org/10.1210/er.2019-00049.
Lampalzer U, Briken P, Schweizer K. Psychosocial care and support in the field of intersex/diverse sex development (dsd): counselling experiences, localisation and needed improvements. Int J Impot Res. 2021;33:2228–42. https://doi.org/10.1038/s41443-021-00422-x.
Lu L, Luo F, Wang X. Gonadal tumor risk in pediatric and adolescent phenotypic females with disorders of sex development and Y chromosomal constitution with different genetic etiologies. Front Pediatr. 2022;10:856128. https://doi.org/10.3389/fped.2022.856128.
Mendonca BB, Domenice S, Arnhold IJ, Costa EM. 46,XY disorders of sex development (DSD). Clin Endocrinol (Oxf). 2009;70. https://doi.org/10.1111/j.1365-2265.2008.03392.x. :2:173 – 87.
Rey RA, Grinspon RP. Normal male sexual differentiation and aetiology of disorders of sex development. Best Pract Res Clin Endocrinol Metab. 2011;25:2221–38. https://doi.org/10.1016/j.beem.2010.08.013.
Ono M, Harley VR. Disorders of sex development: new genes, new concepts. Nat Rev Endocrinol. 2013;9:2:79–91. https://doi.org/10.1038/nrendo.2012.235.
Delot EC, Vilain E. Towards improved genetic diagnosis of human differences of sex development. Nat Rev Genet. 2021;22:9:588–602. https://doi.org/10.1038/s41576-021-00365-5.
Ahmed SF, Bashamboo A, Lucas-Herald A, McElreavey K. Understanding the genetic aetiology in patients with XY DSD. Br Med Bull. 2013;106:67–89. https://doi.org/10.1093/bmb/ldt008.
Luo X, Wang R, Sun Y, Qiu W, Lu D, Wang Y, et al. Deep intronic PAH variants explain missing heritability in Hyperphenylalaninemia. J Mol Diagn. 2023;25:5. https://doi.org/10.1016/j.jmoldx.2023.02.001.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i17884–90. https://doi.org/10.1093/bioinformatics/bty560.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:14:1754–60. https://doi.org/10.1093/bioinformatics/btp324.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:16. https://doi.org/10.1093/bioinformatics/btp352.
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinf. 2013;43:11. https://doi.org/10.1002/0471250953.bi1110s43.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:16e164. https://doi.org/10.1093/nar/gkq603.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:5405–24. https://doi.org/10.1038/gim.2015.30.
Kopanos C, Tsiolkas V, Kouris A, Chapple CE, Albarca Aguilera M, Meyer R, et al. VarSome: the human genomic variant search engine. Bioinformatics. 2019;35:11. https://doi.org/10.1093/bioinformatics/bty897.
Wang H, Zhang L, Wang N, Zhu H, Han B, Sun F, et al. Next-generation sequencing reveals genetic landscape in 46, XY disorders of sexual development patients with variable phenotypes. Hum Genet. 2018;137:3265–77. https://doi.org/10.1007/s00439-018-1879-y.
Yu BQ, Liu ZX, Gao YJ, Wang X, Mao JF, Nie M, et al. Prevalence of gene mutations in a Chinese 46,XY disorders of sex development cohort detected by targeted next-generation sequencing. Asian J Androl. 2021;23(1):69–73. https://doi.org/10.4103/aja.aja_36_20.
Zhang B, Song Y, Li W, Gong C. Variant analysis of the chromodomain helicase DNA-binding protein 7 in pediatric disorders of sex development. Pediatr Investig. 2019;3:131–8. https://doi.org/10.1002/ped4.12111.
Alimussina M, Diver LA, McGowan R, Ahmed SF. Genetic testing of XY newborns with a suspected disorder of sex development. Curr Opin Pediatr. 2018;30:4. https://doi.org/10.1097/MOP.0000000000000644.
Yatsenko SA, Witchel SF. Genetic approach to ambiguous genitalia and disorders of sex development: what clinicians need to know. Semin Perinatol. 2017;41:4. https://doi.org/10.1053/j.semperi.2017.03.016.
Pleskacova J, Hersmus R, Oosterhuis JW, Setyawati BA, Faradz SM, Cools M, et al. Tumor risk in disorders of sex development. Sex Dev. 2010;4:4–5. https://doi.org/10.1159/000314536.
Sauna ZE, Kimchi-Sarfaty C. Understanding the contribution of synonymous mutations to human disease. Nat Rev Genet. 2011;12:10. https://doi.org/10.1038/nrg3051.
Batista RL, Mendonca BB. Integrative and Analytical Review of the 5-Alpha-reductase type 2 Deficiency Worldwide. Appl Clin Genet. 2020;13:83–96. https://doi.org/10.2147/TACG.S198178.
Han B, Cheng T, Zhu H, Yu J, Zhu WJ, Song HD, et al. Genetic analysis of 25 patients with 5alpha-Reductase Deficiency in Chinese Population. Biomed Res Int. 2020;2020:1789514. https://doi.org/10.1155/2020/1789514.
Gui B, Song Y, Su Z, Luo FH, Chen L, Wang X, et al. New insights into 5alpha-reductase type 2 deficiency based on a multi-centre study: regional distribution and genotype-phenotype profiling of SRD5A2 in 190 Chinese patients. J Med Genet. 2019;56:10685–92. https://doi.org/10.1136/jmedgenet-2018-105915.
Palui R, Ravichandran L, Kamalanathan S, Chapla A, Sahoo J, Narayanan N, et al. Clinical, hormonal, and genetic spectrum of 46 XY disorders of sexual development (DSD) patients. Indian J Pediatr. 2024. https://doi.org/10.1007/s12098-024-05144-8.
Akcan N, Uyguner O, Bas F, Altunoglu U, Toksoy G, Karaman B, et al. Mutations in AR or SRD5A2 genes: clinical findings, endocrine pitfalls, and genetic features of children with 46,XY DSD. J Clin Res Pediatr Endocrinol. 2022;14:2153–71. https://doi.org/10.4274/jcrpe.galenos.2022.2021-9-19.
Hughes IA. Consequences of the Chicago DSD Consensus: a personal perspective. Horm Metab Res. 2015;47:5:394–400. https://doi.org/10.1055/s-0035-1545274.
Mendonca BB, Inacio M, Costa EM, Arnhold IJ, Silva FA, Nicolau W, et al. Male pseudohermaphroditism due to steroid 5alpha-reductase 2 deficiency. Diagnosis, psychological evaluation, and management. Med (Baltim). 1996;75:2:64–76. https://doi.org/10.1097/00005792-199603000-00003.
Davey RA, Grossmann M. Androgen receptor structure, function and Biology: from bench to Bedside. Clin Biochem Rev. 2016;37:1:3–15. https://www.ncbi.nlm.nih.gov/pubmed/27057074.
Newmark JR, Hardy DO, Tonb DC, Carter BS, Epstein JI, Isaacs WB, et al. Androgen receptor gene mutations in human prostate cancer. Proc Natl Acad Sci U S A. 1992;89:14. https://doi.org/10.1073/pnas.89.14.6319.
Lv SD, Wang HY, Yu XP, Zhai QL, Wu YB, Wei Q, et al. Integrative molecular characterization of Chinese prostate cancer specimens. Asian J Androl. 2020;22(2):162–8. https://doi.org/10.4103/aja.aja_36_19.
Torquato S, Pallavajjala A, Goldstein A, Toro PV, Silberstein JL, Lee J, et al. Genetic alterations detected in cell-free DNA are Associated with Enzalutamide and Abiraterone Resistance in Castration-resistant prostate Cancer. JCO Precis Oncol. 2019;3. https://doi.org/10.1200/PO.18.00227.
Lanciotti L, Cofini M, Leonardi A, Bertozzi M, Penta L, Esposito S. Different clinical presentations and management in complete androgen insensitivity syndrome (CAIS). Int J Environ Res Public Health. 2019;16:7. https://doi.org/10.3390/ijerph16071268.
Yang J, Zhao S, Su J, Liu S, Wu Z, Ma W, et al. Comprehensive genomic profiling reveals prognostic signatures and insights into the molecular landscape of colorectal cancer. Front Oncol. 2023;13:1285508. https://doi.org/10.3389/fonc.2023.1285508.
Bouvattier C, Mignot B, Lefevre H, Morel Y, Bougneres P. Impaired sexual activity in male adults with partial androgen insensitivity. J Clin Endocrinol Metab. 2006;91:9:3310–5. https://doi.org/10.1210/jc.2006-0218.
Su L, Cheng J, Yin X, Liu G, Lu Z, Sheng H, et al. Clinical and molecular characteristics in 15 patients with androgen receptor gene mutations from South China. Andrologia. 2017;49:10. https://doi.org/10.1111/and.12763.
Kalfa N, Philibert P, Werner R, Audran F, Bashamboo A, Lehors H, et al. Minor hypospadias: the tip of the iceberg of the partial androgen insensitivity syndrome. PLoS ONE. 2013;8:4e61824. https://doi.org/10.1371/journal.pone.0061824.
Ea V, Bergougnoux A, Philibert P, Servant-Fauconnet N, Faure A, Breaud J, et al. How far should we explore Hypospadias? Next-generation sequencing Applied to a large cohort of Hypospadiac patients. Eur Urol. 2021;79:4507–15. https://doi.org/10.1016/j.eururo.2020.12.036.
Reddy NC, Majidi SP, Kong L, Nemera M, Ferguson CJ, Moore M, et al. CHARGE syndrome protein CHD7 regulates epigenomic activation of enhancers in granule cell precursors and gyrification of the cerebellum. Nat Commun. 2021;12:15702. https://doi.org/10.1038/s41467-021-25846-3.
van Ravenswaaij-Arts C, Martin DM. New insights and advances in CHARGE syndrome: diagnosis, etiologies, treatments, and research discoveries. Am J Med Genet C Semin Med Genet. 2017;175:4:397–406. https://doi.org/10.1002/ajmg.c.31592.
Wang Y, Jiang W, Xia X. Clinical and molecular features of 40 Chinese patients with idiopathic hypogonadotropic hypogonadism. Transl Androl Urol. 2023;12:9. https://doi.org/10.21037/tau-23-225.
Zhang W, Shi J, Zhang C, Jiang X, Wang J, Wang W, et al. Identification of gene variants in 130 Han Chinese patients with hypospadias by targeted next-generation sequencing. Mol Genet Genomic Med. 2019;7:8e827. https://doi.org/10.1002/mgg3.827.
Janssen N, Bergman JE, Swertz MA, Tranebjaerg L, Lodahl M, Schoots J, et al. Mutation update on the CHD7 gene involved in CHARGE syndrome. Hum Mutat. 2012;33:8:1149–60. https://doi.org/10.1002/humu.22086.
Husu E, Hove HD, Farholt S, Bille M, Tranebjaerg L, Vogel I, et al. Phenotype in 18 Danish subjects with genetically verified CHARGE syndrome. Clin Genet. 2013;83:2. https://doi.org/10.1111/j.1399-0004.2012.01884.x.
Bademci G, Cengiz FB, Foster Ii J, Duman D, Sennaroglu L, Diaz-Horta O, et al. Variations in multiple syndromic deafness genes mimic non-syndromic hearing loss. Sci Rep. 2016;6:31622. https://doi.org/10.1038/srep31622.
Marcos S, Sarfati J, Leroy C, Fouveaut C, Parent P, Metz C, et al. The prevalence of CHD7 missense versus truncating mutations is higher in patients with Kallmann syndrome than in typical CHARGE patients. J Clin Endocrinol Metab. 2014;99:E102138–43. https://doi.org/10.1210/jc.2014-2110.
Balasubramanian R, Choi JH, Francescatto L, Willer J, Horton ER, Asimacopoulos EP, et al. Functionally compromised CHD7 alleles in patients with isolated GnRH deficiency. Proc Natl Acad Sci U S A. 2014;111:50:17953–8. https://doi.org/10.1073/pnas.1417438111.
Stamou MI, Georgopoulos NA. Kallmann syndrome: phenotype and genotype of hypogonadotropic hypogonadism. Metabolism. 2018;86:124–34. https://doi.org/10.1016/j.metabol.2017.10.012.
Swee DS, Quinton R. Congenital Hypogonadotrophic Hypogonadism: Minipuberty and the case for neonatal diagnosis. Front Endocrinol (Lausanne). 2019;10:97. https://doi.org/10.3389/fendo.2019.00097.
Kim JH, Kang E, Heo SH, Kim GH, Jang JH, Cho EH, et al. Diagnostic yield of targeted gene panel sequencing to identify the genetic etiology of disorders of sex development. Mol Cell Endocrinol. 2017;444:19–25. https://doi.org/10.1016/j.mce.2017.01.037.
Ostrer H. Disorders of sex development (DSDs): an update. J Clin Endocrinol Metab. 2014;99:5:1503–9. https://doi.org/10.1210/jc.2013-3690.
Eggers S, Sadedin S, van den Bergen JA, Robevska G, Ohnesorg T, Hewitt J, et al. Disorders of sex development: insights from targeted gene sequencing of a large international patient cohort. Genome Biol. 2016;17:1243. https://doi.org/10.1186/s13059-016-1105-y.
Zhao Z, Zhan Y, Chen W, Ma X, Sheng W, Huang G. Functional analysis of rare variants of GATA4 identified in Chinese patients with congenital heart defect. Genesis. 2019;57:11–2. https://doi.org/10.1002/dvg.23333.
Deutsch GB, Zielonka EM, Coutandin D, Weber TA, Schafer B, Hannewald J, et al. DNA damage in oocytes induces a switch of the quality control factor TAp63alpha from dimer to tetramer. Cell. 2011;144:4566–76. https://doi.org/10.1016/j.cell.2011.01.013.
Kim DA, Suh EK. Defying DNA double-strand break-induced death during prophase I meiosis by temporal TAp63alpha phosphorylation regulation in developing mouse oocytes. Mol Cell Biol. 2014;34:8. https://doi.org/10.1128/MCB.01223-13.
Turan S, Bastepe M. GNAS Spectrum of disorders. Curr Osteoporos Rep. 2015;13(3):146–58. https://doi.org/10.1007/s11914-015-0268-x.
Chen X, Meng Y, Tang M, Wang Y, Xie Y, Wan S, et al. A paternally inherited non-sense variant c.424G > T (p.G142*) in the first exon of XLalphas in an adult patient with hypophosphatemia and osteopetrosis. Clin Genet. 2020;97:5. https://doi.org/10.1111/cge.13734.
Richard N, Molin A, Coudray N, Rault-Guillaume P, Juppner H, Kottler ML. Paternal GNAS mutations lead to severe intrauterine growth retardation (IUGR) and provide evidence for a role of XLalphas in fetal development. J Clin Endocrinol Metab. 2013;98:E91549–56. https://doi.org/10.1210/jc.2013-1667.
Mazen I, Abdel-Hamid M, Mekkawy M, Bignon-Topalovic J, Boudjenah R, El Gammal M, et al. Identification of NR5A1 mutations and possible digenic inheritance in 46,XY gonadal dysgenesis. Sex Dev. 2016;10:3147–51. https://doi.org/10.1159/000445983.
Sarfati J, Guiochon-Mantel A, Rondard P, Arnulf I, Garcia-Pinero A, Wolczynski S, et al. A comparative phenotypic study of Kallmann syndrome patients carrying monoallelic and biallelic mutations in the prokineticin 2 or prokineticin receptor 2 genes. J Clin Endocrinol Metab. 2010;95:2. https://doi.org/10.1210/jc.2009-0843.
Wang Y, Qin M, Fan L, Gong C. Correlation Analysis of Genotypes and phenotypes in Chinese Male Pediatric patients with congenital hypogonadotropic hypogonadism. Front Endocrinol (Lausanne). 2022;13:846801. https://doi.org/10.3389/fendo.2022.846801.
Sykiotis GP, Plummer L, Hughes VA, Au M, Durrani S, Nayak-Young S, et al. Oligogenic basis of isolated gonadotropin-releasing hormone deficiency. Proc Natl Acad Sci U S A. 2010;107:34:15140–4. https://doi.org/10.1073/pnas.1009622107.
Fluck CE, Audi L, Fernandez-Cancio M, Sauter KS, Martinez de LaPiscina I, Castano L, et al. Broad Phenotypes of Disorders/Differences of Sex Development in MAMLD1 patients through oligogenic disease. Front Genet. 2019;10:746. https://doi.org/10.3389/fgene.2019.00746.
YF Xu ea. New insights from unbiased panel and whole-exome sequencing in a large Chinese Cohort with Disorders of Sex Developmen. Eur Soc Endocrinol. 2019; https://doi.org/10.1530/EJE-19-0111..
Xu Y, Wang Y, Li N, Yao R, Li G, Li J, et al. New insights from unbiased panel and whole-exome sequencing in a large Chinese cohort with disorders of sex development. Eur J Endocrinol. 2019;181:3:311 − 23; https://doi.org/10.1530/EJE-19-0111.
Fan L, Song Y, Polak M, Li L, Ren X, Zhang B, et al. Clinical characteristics and genotype-phenotype correlations of 130 Chinese children in a high-homogeneity single-center cohort with 5alpha-reductase 2 deficiency. Mol Genet Genomic Med. 2020;8:10:e1431; https://doi.org/10.1002/mgg3.1431.
Xie QG, Luo P, Xia K, Li ZQ, Xu Z, Su C, et al. 46,XY disorders of sex development: the use of NGS for prevalent variants. Hum Genet. 2022;141:12:1863-73; https://doi.org/10.1007/s00439-022-02465-6.
Tang Y, Chen Y, Wang J, Zhang Q, Wang Y, Xu Y, et al. Clinical characteristics and genetic expansion of 46,XY disorders of sex development children in a Chinese prospective study. Endocr Connect. 2023;12:10; https://doi.org/10.1530/EC-23-0029.
White S, Ohnesorg T, Notini A, Roeszler K, Hewitt J, Daggag H, et al. Copy number variation in patients with disorders of sex development due to 46,XY gonadal dysgenesis. PLoS One. 2011;6:3:e17793; https://doi.org/10.1371/journal.pone.0017793.
Acknowledgements
AcknowledgmentsWe would like to thank all the participants in this study and all the reviewers who participated in the review.
Funding
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Author information
Authors and Affiliations
Contributions
Study conception was designed by HY-C, GJ-C, X-Y, and L-Y. Material preparation, data collection, and analysis were performed by HY-C, GJ-C, FX-L, Y-H, LF-Z and YJ-Z . Variants prediction and verification were completed by HY-C, ZY-J. the first draft of the manuscript was written by HY-C and GJ-C. X-Y and L-Y critically revised the manuscript for important intellectual content. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the local IRB (2018-IRB-076) and informed consent was obtained prior to recruitment.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, H., Chen, G., Li, F. et al. Application and insights of targeted next-generation sequencing in a large cohort of 46,XY disorders of sex development in Chinese. Biol Sex Differ 15, 73 (2024). https://doi.org/10.1186/s13293-024-00648-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13293-024-00648-6