Background: Our multidisciplinary team (MDT) is a large specialized team based in Semarang, Indonesia, that cares for a wide variety of pediatric and adult individuals with differences of sex development (DSD) from across Indonesia. Here, we describe our work over the last 17 years. Methods: We analyzed phenotypic, hormonal, and genetic findings from clinical records for all patients referred to our MDT during the period 2004–2020. Results: Among 1,184 DSD patients, 10% had sex chromosome DSD, 67% had 46,XY DSD, and 23% had 46,XX DSD. The most common sex chromosome anomaly was Turner syndrome (45,X) (55 cases). For patients with 46,XY DSD under-masculinization was the most common diagnosis (311 cases), and for 46,XX DSD, a defect of Müllerian development was most common (131 cases) followed by congenital adrenal hyperplasia (CAH) (116 cases). Sanger sequencing, MLPA, and targeted gene sequencing of 257 patients with 46,XY DSD found likely causative variants in 21% (55 cases), with 13 diagnostic genes implicated. The most affected gene codes for the androgen receptor. Molecular analysis identified a diagnosis for 69 of 116 patients with CAH, with 62 carrying variants in CYP21A2 including four novel variants, and 7 patients carrying variants in CYP11B1. In many cases, these genetic diagnoses influenced the clinical management of patients and their families. Conclusions: Our work has highlighted the occurrence of different DSDs in Indonesia. By applying sequencing technologies as part of our clinical care, we have delivered a number of genetic diagnoses and identified novel pathogenic variants in some genes, which may be clinically specific to Indonesia. Genetics can inform many aspects of DSD clinical management, and while many of our patients remain undiagnosed, we hope that future testing may provide answers for even more.

The development of sex and gender requires an intricate interplay of sex chromosomes, genes and genetic pathways, hormones, anatomy, psychology, and social behavior [Cools et al., 2018; Zhao and Yao, 2019; Garcia-Acero et al., 2020]. While our fundamental understanding of sex development and differentiation is still evolving, we have gained significant insights into these processes through the analysis of individuals with differences of sex development (DSD).

DSDs are conditions with atypical chromosomal, gonadal, or anatomical sex and encompasses a broad range of clinical phenotypes. Individuals born with a DSD were classified into three main groups: 46,XY DSD, 46,XX DSD, and sex chromosome DSD following a consensus meeting in Chicago [Hughes et al., 2006]. The classification of sex chromosomal DSD is subdivided into 45,X (Turner syndrome), 47,XXY (Klinefelter), 45,X/46,XY (mixed gonadal dysgenesis), and 46,XY/46,XX chimeric (ovotesticular DSD) subtypes. Although this classification has its limitations and has been subject of intense discussion, its use prevents misunderstanding and confusion in particular by delineating these conditions from other conditions such as gender dysphoria.

Patients with DSD may present with a wide range of phenotypes, from ambiguous or atypical genitalia, absence of secondary sex characteristics, primary amenorrhea, hypospadias, clitoromegaly to complex congenital malformations, such as cloacal exstrophy [Arboleda et al., 2014]. DSD cases can be rare and the rate of occurrence differs for each condition and from country to country, although their incidence is not well established [Lee et al., 2016], especially in developing countries such as Indonesia.

Indonesia is a developing country where tropical infectious diseases are closely associated with poverty and limited resources, and where noncommunicable diseases including genetic conditions are not a health care priority. In developed countries, DSD is often diagnosed in new-born babies or even before birth, and clinical management involves large interdisciplinary expert teams in a dedicated center. However, the management of DSD in Indonesia provides some unique challenges including limitations in diagnostic and treatment facilities in addition to poverty, sociocultural factors, religious and traditional beliefs [Ediati et al., 2015b].

Many patients presenting with DSD in Indonesia do not receive diagnostic evaluation, and options for medical and surgical treatment are limited. The lack of a national standardized management protocol means that DSD care is highly variable and often depends on the experience of the treating clinicians and facilities available in a particular health center/tertiary hospital.

Additionally, health professionals with less awareness of DSD find it challenging to explain the implications and management options to parents. Diagnostic facilities are usually limited to research centers at the universities or provincial referral hospitals, where some DSD genes can be detected using cytogenetics and simple molecular techniques. However, many patients will not know the cause of their DSD and current clinical management strategies depend on phenotypic appearance, imaging, chromosome analysis, and hormonal assays. In response to these challenges, a multidisciplinary team (MDT) for DSD at our center was established in 1989, and it remains the only prominent one of its kind in Indonesia, catering to a total population of >270 million.

Since clinical management of DSD depends on a correct diagnosis, genetic testing has become an essential part of the diagnostic pathway, especially where there is limited access to other tools such as detailed endocrine testing. A molecular diagnosis also provides opportunities for understanding disease mechanisms, patterns of inheritance, and planning clinical management. Advances in genetics have led to a better understanding of the molecular causes of DSD and have provided new insight into the pathophysiology with a number of genes now known to be essential for sexual differentiation in humans [Délot and Vilain, 2021]. In recent years, next-generation sequencing technologies have increased the diagnostic yield to between 43% and 69.2% of DSD cases [Abualsaud et al., 2021; Globa et al., 2022; Zidoune et al., 2022] compared to <20% reported in previous years [Arboleda et al., 2013]. Our previous study using HaloPlex Agilent Technology identified variants in 28 genes from a panel of 64 known diagnostic genes for DSD in just over 300 DSD patients with different conditions [Eggers et al., 2016]. Our earlier research on the etiological spectrum of DSD in a large cohort of underprivileged and undiagnosed patients from Indonesia found that a stepwise diagnostic approach led to a molecular or histologically established diagnosis in 29.4% of patients [Juniarto et al., 2016], indicating the importance of molecular analysis but also highlighting that other genetic causes remain to be discovered.

In this study, we reflect on the large cohort of patients with DSD seen by a MDT at a single Indonesian health center over the past 17 years. We focus on the genetic analysis and describe the distribution of gene variants among Indonesian DSD patients and the lessons learned along the way.

Patient Clinical Notes and Classification

We report here on all patients entered into the study by the MDT and the laboratory of Center for Biomedical Research (CEBIOR) of National Diponegoro University Hospital from 2004 to 2020. This includes patients referred from provincial referral hospitals/tertiary hospitals, peripheral hospitals, private hospitals, or health center by physicians (mostly pediatricians, endocrinologists, urologists, and gynecologists). Clinical notes for patients include physical descriptions such as clinical features of the external genitalia, the presence of gonads, the fusion of the labioscrotal folds, the size of the phallus, and the site of the urinary meatus on the phallus scored by the External Masculinization Score (EMS) [Ahmed et al., 2000]. Patients with EMS total scores <9 were considered as under-virilized males.

Patients were classified as suspected of androgen action disorder based on an EMS score <9, while patients with an EMS score ≥9 were classified as under-masculinization of unknown cause (UMU) [Juniarto et al., 2016]. Patients with various DSD phenotypes were referred to CEBIOR for chromosomal analysis. In 46,XY DSD patients, diagnoses were grouped into androgen action disorder, UMU, and gonadal dysgenesis (GD) based on EMS and serum concentration of LH, FSH, and testosterone levels [Juniarto et al., 2016]. The 46,XX DSD patients were grouped based on hormonal diagnosis and physical examination, mainly 17 hydroxyprogesterone for congenital adrenal hyperplasia (CAH). Some patients also had ultrasonographic or radiologic imaging.

MDT and Patient Clinical Care

The MDT consists of specialists representing various departments from the Dr. Kariadi provincial referral Hospital and the National Diponegoro University Hospital including urology, pediatric endocrinology, plastic surgery, gynecology, andrology, anesthetics, medical genetics (including a genetic counselor), psychiatry, pathology, psychology, legal medicine, and religious leaders. Discussion for 5–10 cases were carried out at bimonthly MDT meetings; their family (mainly parent) were also invited if necessary, and the genetic testing was done for the extended family. In some cases, the team referred patients or families to psychological evaluation by the psychologist member of the MDT team who had specialized training on DSD. This was most common where assessment was required to evaluate the gender identity due to patient’s doubts over their assigned gender or wishes for gender reassignment. An individual interview was conducted for patients and their parents or caregiver regarding gender role behavior, gender stereotypes, and their acceptance of DSD in the family. In pediatric cases where surgery was considered, this was discussed broadly and with the families, only proceeding with informed consent from parents. Regardless of the surgery, sex was assigned based on cytogenetics, internal and external reproductive organs, hormone levels, and psychosocial analysis, including cultural and religious issues, and in many cases, the wishes of the parents had a strong bearing on the sex choice of children. The flowchart of patient care is shown in Figure 1.

Fig. 1.

Algorithm of DSD management.

Fig. 1.

Algorithm of DSD management.

Close modal

Molecular Testing

Collaborations with Erasmus Medical Centre (Erasmus MC) in Rotterdam and Radboud University Medical Center (Radboudumc) in Nijmegen, the Netherlands, and the Murdoch Children’s Research Institute (MCRI) in Melbourne, Australia, provided genetic testing for various conditions using differing technologies.

Specifically, Radboudumc: CAH patients – Sanger sequencing and MLPA (P050-C1 kit, MRC Holland, Amsterdam, the Netherlands) of CYP21A2 and CYP11B1. Erasmus MC: CAH patients (46,XX DSD with increased 17-OHP level) – Sanger sequencing of the CYP21A2 and CYP11B1 gene. 46,XY DSD patients (with an EMS score <9) for AR, SRY, WNT4, NR3C1, and LHCGR gene Sanger sequencing analysis. MCRI: Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome – variant analysis for using microarrays and whole-exome sequencing. 46,XY DSD – targeted gene sequencing using HaloPlex Agilent Technology on patients [Eggers et al., 2016].

Our cohort consisted of 1,184 patients covering the period from 2004 to 2020. Patient numbers fluctuated each year, with lower numbers in 2014 because our laboratory moved to the main university campus and in 2020 due to limitations imposed by COVID-19 (Fig. 2). Patients were grouped into three categories (sex chromosome DSD, 46,XX DSD, and 46,XY DSD), which were subsequently divided into the appropriate subgroups (Table 1).

Fig. 2.

Number of DSD patients per year.

Fig. 2.

Number of DSD patients per year.

Close modal
Table 1.

Distribution of DSD cases during 2004–2020 according to DSD classification

No. of casesSex chromosome DSD46,XY DSD46,XX DSD
 110 (9%) 800 (68%) 274 (23%) 
No. of casesSex chromosome DSD46,XY DSD46,XX DSD
 110 (9%) 800 (68%) 274 (23%) 
Turner and its variant47,XXY Klinefelter45,X/46,XY (MGD, ovotesticular DSD)46,XX/46,XY ovotesticular DSDOthersHypospadiasGDAADUMUCAHAndrogen excessDefect of Mullerian developmentGDOthers
1,184 55 (50) 12 (10) 21 (20) 9 (8) 13 (12) 292 (36) 32 (4) 157 (20) 319 (40) 116 (42) 11 (4) 131 (48) 3 (1) 13 (5) 
Turner and its variant47,XXY Klinefelter45,X/46,XY (MGD, ovotesticular DSD)46,XX/46,XY ovotesticular DSDOthersHypospadiasGDAADUMUCAHAndrogen excessDefect of Mullerian developmentGDOthers
1,184 55 (50) 12 (10) 21 (20) 9 (8) 13 (12) 292 (36) 32 (4) 157 (20) 319 (40) 116 (42) 11 (4) 131 (48) 3 (1) 13 (5) 

Values are n (%). AAD, androgen action disorders; UMU, under-masculinization of unknown cause male under-masculinization; MGD, mixed gonadal dysgenesis; GD, gonadal dysgenesis; CAH, congenital adrenal hyperplasia.

A total of 800 patients had 46,XY DSD; in total, 152 of these patients were evaluated in Erasmus MC, while 105 patients were assessed at MCRI. Of those analyzed at MCRI, 97 were previously reported in the study by Eggers et al. [2016]. The remaining patients are still without diagnosis except for 5 patients who were found to carry causative AR variants as reported by Listyasari et al. [2019] Of those analyzed at Erasmus MC, a total of 152 cases evaluated were previously reported by Juniarto et al. [2016]; Listyasari et al. [2021] The molecular analysis of these 257 patients identified genetic variants considered diagnostic in 55 cases (21%) from 13 diagnostic DSD genes (Fig. 3). The undiagnosed cases included patients with UMU, GD, or hypospadias.

Fig. 3.

Graph showing the genes responsible for a genetic diagnosis for 46,XY DSD patients. Of a total of 257 46,XY DSD patients sequenced, 55 had variants in a DSD gene that were curated to be likely pathogenic or pathogenic. This illustrates the breakdown of these variants in terms of genes responsible.

Fig. 3.

Graph showing the genes responsible for a genetic diagnosis for 46,XY DSD patients. Of a total of 257 46,XY DSD patients sequenced, 55 had variants in a DSD gene that were curated to be likely pathogenic or pathogenic. This illustrates the breakdown of these variants in terms of genes responsible.

Close modal

Patient presented with 46,XX DSD were 274 cases. In 116 of these cases, a diagnosis of CAH was made, which was confirmed by genetic analysis for 69 cases: in 62 patients, 28 pathogenic variants were identified in the CYP21A2 gene (Table 2), while in the CYP11B1 gene (NM_000497.4), two different pathogenic variants were detected in seven cases, i.e., c.1121G>A (2 cases) and c.799G>A (5 cases). The diagnostic yield for DSD cases is available in Figure 4. The most common pathogenic variant in CAH cases is c.1069C>T p(Arg357Trp) in CYP21A2 (NM_000500.9) which was found in 49 alleles (out of 124 alleles); 15 patients were homozygous for this variant and in 2 patients this variant was part of a complex allele with several variants (Table 2) [Krone and Arlt, 2009; New et al., 2013; Baumgartner-Parzer et al., 2020]. The remaining 47 cases were diagnosed based on clinical and biochemical evaluation only. In addition to the characterized and already registered pathogenic variants in CYP21A2, we also discovered four novel variants which were not described in any variant library. In case 1 (c.1155_1157del p.[Ile386del]) has salt wasting type; case 2 (c.739G>T p.[Glu247*]) has simple virilizing type; and case 3 (c.586C>T p.[Gln196]) has salt wasting type all of which were identified at late age, have virilization of external genitalia, short stature, and hyperpigmentation. While case 4 (c.65G>A p.[Trp22*]) has salt wasting type which is recognized at early age with normal length/height-for-age and normal puberty due to good medication compliance.

Table 2.

Identified pathogenic variants in the CYP21A2 gene (NM_000500.7) in patients with CAH and their predicted phenotype based on residual enzyme activity and literature [Krone and Arlt, 2009; New et al., 2013; Baumgartner-Parzer et al., 2020]

No. of patientsAllele 1Allele 2(Assumed) inheritancePredicted phenotype based on residual enzyme activity
nucleotide changeprotein changelegacy name(predicted) enzyme activity, %nucleotide changeprotein changelegacy name(predicted) enzyme activity, %
c.293-13C>G p.? I2G 0–1 c.293-13C>G p.? I2G 0–1 Hom. SW 
c.518T>A p.(Ile173Asn) I172N 1–5 c.518T>A p.(Ile173Asn) I172N 1–5 Hom. SV 
c.518T>A p.(Ile173Asn) I172N 1–5 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW or SV 
c.955C>T p.(Gln319*) Q318X Not identified yet Het. Unknown 
15 c.1069C>T p.(Arg357Trp) R356W c.1069C>T p.(Arg357Trp) R356W Hom. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.955C>T p.(Gln319*) Q318X Comp. het. NC 
c.844G>T p.(Val282Leu) V281L 30–50 c.955C>T p.(Gln319*) Q318X Comp. het. NC 
c.92C>T p.(Pro31Leu) P30L 20–30 c.844G>T p.(Val282Leu) V281L-Q318X-R356W Comp. het. NC 
c.955C>T p.(Gln319*) 
c.1069C>T p.(Arg357Trp) 
c.92C>T p.(Pro31Leu) P30L 20–30 c.92C>T p.(Pro31Leu) P30L 20–30 Hom. NC 
c.293-13A/C>G p.? I2G 0–1 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.1069C>T p.(Arg357Trp) R356W Comp. het. NC 
c.1155_1157del p.(Ile386del) Unknown Unknown c.1069C>T p.(Arg357Trp) R356W Comp. het. Unknown 
c.1279C>T p.(Arg427Cys) R426C c.1279C>T p.(Arg427Cys) R426C Hom. SW 
c.1069C>T p.(Arg357Trp) R356W c.65G>A p.(Trp22*) Unknown 0 Comp. het. SW 
c.710T>A p.(Ile237Asn) E6 cluster or [I236N, V237E, M239K] c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.713T>A p.(Val238Glu) 
c.719T>A p.(Met240Lys) 
c.332_339del p.(Gly111fs) E3Δ8bp c.739G>T p.(Glu247*) Unknown 0 Comp. het. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.1069C>T p.(Arg357Trp) R356W Comp. het. NC 
c.60G>A p.(Trp20*) W19X c.518T>A p.(Ile173Asn) I172N 1–5 Comp. het. SW or SV 
c.1069C>T p.(Arg357Trp) R356W c.1218G>A p.(Trp406*) W405X Comp. het. SW 
c.586C>T p.(Gln196*) Unknown 0 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.1069C>T p.(Arg357Trp) R356W c.1451G>C p.(Arg484Pro) R483P Unknown Comp. het. SV 
c.844G>T p.(Val282Leu) V281L-F306 + 1nt-Q318X-R356W c.1218G>A p.(Trp406*) W405X Comp. het. SW 
c.923dup p.(Leu308fs) 
c.955C>T p.(Gln319*) 
c.1069C>T p.(Arg357Trp) 
c.293-13C>G p.? I2G 0–1 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.( ?_-107)_(939 + ?)del p.0 Deletion exon 1–7 c.(?_-107)_(939+?)del p.0 deletion exon 1–7 Hom. SW 
c.(?_-107)_(738 + ?)del p.0 Deletion exon 1–6 c.518T>A p.(Ile173Asn) I172N 1–5 Comp. het. SW or SV 
c.(?_-107)_(738 + ?)del p.0 Deletion exon 1–6 Not identified yet Het. Unknown 
c.(?_-107)_(447 + ?)del p.0 Deletion exon 1–3 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
No. of patientsAllele 1Allele 2(Assumed) inheritancePredicted phenotype based on residual enzyme activity
nucleotide changeprotein changelegacy name(predicted) enzyme activity, %nucleotide changeprotein changelegacy name(predicted) enzyme activity, %
c.293-13C>G p.? I2G 0–1 c.293-13C>G p.? I2G 0–1 Hom. SW 
c.518T>A p.(Ile173Asn) I172N 1–5 c.518T>A p.(Ile173Asn) I172N 1–5 Hom. SV 
c.518T>A p.(Ile173Asn) I172N 1–5 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW or SV 
c.955C>T p.(Gln319*) Q318X Not identified yet Het. Unknown 
15 c.1069C>T p.(Arg357Trp) R356W c.1069C>T p.(Arg357Trp) R356W Hom. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.955C>T p.(Gln319*) Q318X Comp. het. NC 
c.844G>T p.(Val282Leu) V281L 30–50 c.955C>T p.(Gln319*) Q318X Comp. het. NC 
c.92C>T p.(Pro31Leu) P30L 20–30 c.844G>T p.(Val282Leu) V281L-Q318X-R356W Comp. het. NC 
c.955C>T p.(Gln319*) 
c.1069C>T p.(Arg357Trp) 
c.92C>T p.(Pro31Leu) P30L 20–30 c.92C>T p.(Pro31Leu) P30L 20–30 Hom. NC 
c.293-13A/C>G p.? I2G 0–1 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.1069C>T p.(Arg357Trp) R356W Comp. het. NC 
c.1155_1157del p.(Ile386del) Unknown Unknown c.1069C>T p.(Arg357Trp) R356W Comp. het. Unknown 
c.1279C>T p.(Arg427Cys) R426C c.1279C>T p.(Arg427Cys) R426C Hom. SW 
c.1069C>T p.(Arg357Trp) R356W c.65G>A p.(Trp22*) Unknown 0 Comp. het. SW 
c.710T>A p.(Ile237Asn) E6 cluster or [I236N, V237E, M239K] c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.713T>A p.(Val238Glu) 
c.719T>A p.(Met240Lys) 
c.332_339del p.(Gly111fs) E3Δ8bp c.739G>T p.(Glu247*) Unknown 0 Comp. het. SW 
c.92C>T p.(Pro31Leu) P30L 20–30 c.1069C>T p.(Arg357Trp) R356W Comp. het. NC 
c.60G>A p.(Trp20*) W19X c.518T>A p.(Ile173Asn) I172N 1–5 Comp. het. SW or SV 
c.1069C>T p.(Arg357Trp) R356W c.1218G>A p.(Trp406*) W405X Comp. het. SW 
c.586C>T p.(Gln196*) Unknown 0 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.1069C>T p.(Arg357Trp) R356W c.1451G>C p.(Arg484Pro) R483P Unknown Comp. het. SV 
c.844G>T p.(Val282Leu) V281L-F306 + 1nt-Q318X-R356W c.1218G>A p.(Trp406*) W405X Comp. het. SW 
c.923dup p.(Leu308fs) 
c.955C>T p.(Gln319*) 
c.1069C>T p.(Arg357Trp) 
c.293-13C>G p.? I2G 0–1 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 
c.( ?_-107)_(939 + ?)del p.0 Deletion exon 1–7 c.(?_-107)_(939+?)del p.0 deletion exon 1–7 Hom. SW 
c.(?_-107)_(738 + ?)del p.0 Deletion exon 1–6 c.518T>A p.(Ile173Asn) I172N 1–5 Comp. het. SW or SV 
c.(?_-107)_(738 + ?)del p.0 Deletion exon 1–6 Not identified yet Het. Unknown 
c.(?_-107)_(447 + ?)del p.0 Deletion exon 1–3 c.1069C>T p.(Arg357Trp) R356W Comp. het. SW 

Bold indicates variants that have not been described before. Adherence to HGVS nomenclature rules (varnomen.hgvs.org) makes the nucleotide numbers 3 higher and the amino acid numbers 1 higher than the legacy names used in literature.

Comp. het., compound heterozygous; het., heterozygous; hom., homozygous; NC, nonclassical; SV, simple virilizing; SW, salt wasting.

Fig. 4.

Diagnostic yield of the DSD cohort in this study. a The molecular analysis of 46,XY DSD patients. From 257 patients who were tested, a genetic diagnosis was identified in 21% of cases. b In the 46,XX DSD patient cohort, from 184 patients who had genetic test, 33% were found to have a pathogenic variant in a DSD gene.

Fig. 4.

Diagnostic yield of the DSD cohort in this study. a The molecular analysis of 46,XY DSD patients. From 257 patients who were tested, a genetic diagnosis was identified in 21% of cases. b In the 46,XX DSD patient cohort, from 184 patients who had genetic test, 33% were found to have a pathogenic variant in a DSD gene.

Close modal

The remaining 158 patients with 46,XX DSD were diagnosed with unknown androgen excess disorders, defect in Mullerian development such as MRKH syndrome, GD, or Mullerian duct aplasia-renal agenesis-cervicothoracic somite dysplasia or cloacal malformation. A number of 10 patients with MRKH had genetic testing in MCRI and molecular studies for this group of disorders were not definitive, except for four cases of MRKH type 1 with clinical finding of primary amenorrhea and utero vaginal agenesis [Backhouse et al., 2019]. The other 105 patients in this category were tested for specific genes for WNT4, LHCGR, SRY, and NR3C1 using Sanger sequencing in Erasmus MC. The molecular study for this group have identified for 12 cases of WNT4 defect, 2 cases of glucocorticoid receptor defect, 2 cases of LH receptor defects, and 6 cases of SRY gene translocation as previously reported, while the remaining cases were not conclusive [Juniarto et al., 2012, 2016].

In total, 246 patients (aged 2–45 years) were referred for psychological assessment and care by the MDT. Many parents who had girls with DSD voiced concern about future infertility and asked about keeping the DSD diagnosis hidden from others, including extended family. Adolescent girls and adult women who had a variation in physical appearance (such as a muscular body, masculine behavior, deep voice) experienced more negative reactions and rejection from peers or the community compared to girls with no visible external differences in these characteristics. In most cases, patients and their parents received psychological counseling to facilitate better acceptance of their condition in the family and in society. Difficulties in dealing with societal acceptance and understanding were the most common issue reported by parents and patients, particularly when DSD characteristics were visible to others. This included discussion around educating others of their DSD condition, coping with awkward situations in the neighborhood, coping with or anticipating peer or social rejection, and marital or infertility issues. We have previously reported similar findings from 118 patients aged 6–42 years [Ediati et al., 2015a, 2015c].

Clinical care for DSD patients is challenging. Patients with DSD and their families require complex medical and psychosocial care, as well as ongoing follow-up throughout their life. Multidisciplinary care, of which psychosocial and peer support are key components, is recommended as best practice [Moran and Karkazis, 2012; Ahmed et al., 2021]. Most importantly, patient management needs to be individualized, especially for decisions related to hormonal treatment, sex assignment, surgical interventions, genetic counseling, long-term treatment monitoring, and potential fertility preservation [Ediati et al., 2017, 2018; Speiser et al., 2018].

In Indonesia, the lack of national guidelines for DSD providers means there is no standardized approach [Dessens et al., 2017]. This is a significant dilemma with serious implications for patients. A standardized network of multidisciplinary DSD teams would provide clinicians with the necessary program planning and team implementation tools.

Late presentations and delayed diagnoses of DSD are also common in Indonesia. This leads to incongruence between expressed gender and the sex assigned at birth which can cause clinically significant distress and impairment in social and emotional function [Ediati et al., 2015b]. The availability of psychological care in addition to medical care is pivotal for patients with DSD and their family. Accepting a diagnosis of DSD can be challenging for patients and families when there is little understanding of these conditions in the community. Affected individuals must make sense of their body and gender when such a discrepancy is present. Having a proper understanding of DSD and its causes can help patients and families come to terms with a diagnosis and educate the wider community about their condition to promote social acceptance. The availability of a MDT in our center has been central to providing clinical care for patients with DSD and their families. It has also raised awareness in both the medical and wider communities.

In the Indonesian clinical setting, the etiology of DSD remains largely unknown. The integration of phenotypic data, such as detailed anatomical, hormonal, histological data, and genotyping will continue to be crucial [Ahmed et al., 2022; Rey, 2022]. Currently, only some of our many patients have access to genetic sequencing technologies, but in the case of 46,XY DSD, it provided a diagnostic finding in 21% of patients. A higher genetic diagnostic rate can be achieved by increasing the availability of next-generation sequencing technology in diagnostic settings [Gomes et al., 2022]. In Indonesia, the molecular findings showed androgen insensitivity (AIS) and 5α-reductase type 2 (5ARD) deficiency as the most common causes of 46,XY DSD, both of which are often inherited. These findings provide genetic counseling opportunities, which can influence family planning and cascade testing [Juniarto et al., 2016; Listyasari et al., 2019, 2021; Marzuki et al., 2019].

In some cases, genetic testing can also provide a definitive diagnosis where traditional methods may have failed. One example is AIS and 5ARD which can be phenotypically very similar. Although these can sometimes be distinguished by hormonal testing [Veiga-Junior et al., 2012], this is not always available or accurate. In Indonesia, the DHT test required to make these diagnoses is not frequently provided in clinical settings due to cost considerations. When molecular analysis for the SRD5A2 gene is unavailable, the urinary etiocholanolone/androsterone ratio may be a relevant test for identifying 5ARD2 patients and carriers [Marzuki et al., 2021]. 5ARD is autosomal recessive, while many cases of AIS we identified had X-linked inheritance, again meaning that a genetic diagnosis may have different implications for families.

In addition to diagnostic changes, molecular sequencing of our patient cohorts also revealed variants of uncertain significance (VUS). Fortunately, functional testing has been able to resolve the causality or lack of for several of these VUS, for example, in heterozygous missense variants in GATA4 [Bergen et al., 2020], novel variants in the prodomain of BMP7 in two pairs of monozygotic concordant twins with proximal hypospadias [Bouty et al., 2019], and DHH variants in 46,XY DSD [Ayers et al., 2019]. While no diagnostic variants were found in MRKH patients (indeed, few genes are yet confirmed to cause MRKH), several candidate genes were identified that may contribute to this condition but further analysis will be required to establish whether they are causative [Backhouse et al., 2019]. Further functional analysis of VUS will be essential, especially of variants that may be unique or overrepresented in populations like Indonesia.

In our work over the past 17 years, we have seen that the yearly rate of DSD diagnoses has remained stable. However, this decreased steeply during the COVID-19 pandemic as many DSD patients/families postponed their appointments as DSD was not considered a medical emergency. However, a large group of 46,XX female patients presenting with ambiguous genitalia at birth have 21-hydroxylase deficiency (CAH). Identification of infants with CAH is urgent because if undiagnosed they may develop salt wasting, which could lead to death. Several other DSD genes cause co-morbidities. A genetic diagnosis can help monitor these or allow early intervention. For example, GATA4 variants in 46,XY DSD can contribute to congenital heart defects. Similarly, WT1 variants can underlie Wilms tumor. Understanding the pattern of inheritance is essential for genetic counseling and family planning.

For patients with DSD, clinical management by a dedicated MDT and accurate genetic diagnoses are required for optimal care. This is a challenge in developing countries, such as Indonesia, where patients face limited access to diagnostic tests such as hormonal analysis or advanced genetic analysis. We hope that by sharing the last 17 years of experience of our DSD MDT, we can spread awareness in the medical community, and work toward national guidelines for establishing this model of care across Indonesia and other developing countries.

We would like to thank the patients and their families for their participation in this study. We also thank the members of the MDT team at Dr. Kariadi Hospital/Faculty of Medicine, Diponegoro University who helped us with patient assessment.

This study protocol was reviewed and approved by the Health Research Ethics Committee, Faculty of Medicine, Diponegoro University, Semarang, Indonesia, approval number 31/H.7/KEPK/FK-RSDK/2010. Written informed consent was obtained from participants (or their parent/legal guardian/next of kin) to participate in the study.

The authors have no conflicts of interest to declare.

This study was supported by the Universitas Diponegoro WCRU Grant No. 118-03/UN7.6.1/PP/2021.

Sultana M.H. Faradz, Frank H. de Jong, Stenvert L.S. Drop, Katie Ayers, and Andrew Sinclair: substantial contributions to the conception or design of the work; drafting the work or revising it critically for important intellectual content; final approval of the version to be published; and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Nurin Aisyiyah Listyasari, Agustini Utari, Mahayu Dewi Ariani, Achmad Zulfa Juniarto, Ardy Santosa, Annastasia Ediati, Tuula K. Rinne, Dineke Westra, and Hedi C. van der Grinten: substantial contributions to the analysis, and interpretation of data for the work; drafting the work or revising it critically for important intellectual content; final approval of the version to be published; and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All authors read and approved the final manuscript.

All data generated or analyzed during this study are included in this article. Further inquiries can be directed to the corresponding author.

1.
Abualsaud
D
,
Hashem
M
,
AlHashem
A
,
Alkuraya
FS
.
Survey of disorders of sex development in a large cohort of patients with diverse Mendelian phenotypes
.
Am J Med Genet
.
2021
;
185
(
9
):
2789
800
.
2.
Ahmed
SF
,
Achermann
J
,
Alderson
J
,
Crouch
NS
,
Elford
S
,
Hughes
IA
et al
.
Society for Endocrinology UK Guidance on the initial evaluation of a suspected difference or disorder of sex development (Revised 2021)
.
Clin Endocrinol
.
2021
;
95
(
6
):
818
40
.
3.
Ahmed
SF
,
Alimusina
M
,
Batista
RL
,
Domenice
S
,
Lisboa Gomes
N
,
McGowan
R
et al
.
The use of genetics for reaching a diagnosis in XY DSD
.
Sex Dev
.
2022
16
2-3
207
224
.
4.
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
.
5.
Arboleda
VA
,
Lee
H
,
Sanchez
FJ
,
Delot
EC
,
Sandberg
DE
,
Grody
WW
et al
.
Targeted massively parallel sequencing provides comprehensive genetic diagnosis for patients with disorders of sex development
.
Clin Genet
.
2013
;
83
(
1
):
35
43
.
6.
Arboleda
VA
,
Sandberg
DE
,
Vilain
E
.
DSDs: genetics, underlying pathologies and psychosexual differentiation
.
Nat Rev Endocrinol
.
2014
;
10
:
603
15
.
7.
Ayers
K
,
van den Bergen
J
,
Robevska
G
,
Listyasari
N
,
Raza
J
,
Atta
I
et al
.
Functional analysis of novel desert hedgehog gene variants improves the clinical interpretation of genomic data and provides a more accurate diagnosis for patients with 46,XY differences of sex development
.
J Med Genet
.
2019
;
56
(
7
):
434
43
.
8.
Backhouse
B
,
Hanna
C
,
Robevska
G
,
van den Bergen
J
,
Pelosi
E
,
Simons
C
et al
.
Identification of candidate genes for mayer-rokitansky-küster-hauser syndrome using genomic approaches
.
Sex Dev
.
2019
;
13
(
1
):
26
34
.
9.
Baumgartner-Parzer
S
,
Witsch-Baumgartner
M
,
Hoeppner
W
.
EMQN best practice guidelines for molecular genetic testing and reporting of 21-hydroxylase deficiency
.
Eur J Hum Genet
.
2020
;
28
(
10
):
1341
67
.
10.
van den Bergen
JA
,
Robevska
G
,
Eggers
S
,
Riedl
S
,
Grover
SR
,
Bergman
PB
et al
.
Analysis of variants in GATA4 and FOG2/ZFPM2 demonstrates benign contribution to 46,XY disorders of sex development
.
Mol Genet Genomic Med
.
2020
;
8
(
3
):
e1095
.
11.
Bouty
A
,
Walton
K
,
Listyasari
NA
,
Robevska
G
,
Van den Bergen
J
,
Santosa
A
et al
.
Functional characterization of two new variants in the bone morphogenetic protein 7 prodomain in two pairs of monozygotic twins with hypospadias
.
J Endocr Soc
.
2019
;
3
(
4
):
814
24
.
12.
Cools
M
,
Nordenstrom
A
,
Robeva
R
,
Hall
J
,
Westerveld
P
,
Fluck
C
et al
.
Caring for individuals with a difference of sex development (DSD): a Consensus Statement
.
Nat Rev Endocrinol
.
2018
;
14
(
7
):
415
29
.
13.
Délot
EC
,
Vilain
E
.
Towards improved genetic diagnosis of human differences of sex development
.
Nat Rev Genet
.
2021
;
22
(
9
):
588
602
.
14.
Dessens
A
,
Guaragna-Filho
G
,
Kyriakou
A
,
Bryce
J
,
Sanders
C
,
Nordenskjold
A
et al
.
Understanding the needs of professionals who provide psychosocial care for children and adults with disorders of sex development
.
BMJ Paediatr Open
.
2017
;
1
:
e000132
.
15.
Ediati
A
,
Faradz
SM
,
Juniarto
AZ
,
van der Ende
J
,
Drop
SL
,
Dessens
AB
.
Emotional and behavioral problems in late-identified Indonesian patients with disorders of sex development
.
J Psychosom Res
.
2015a
79
1
76
84
.
16.
Ediati
A
,
Juniarto
AZ
,
Birnie
E
,
Drop
SL
,
Faradz
SM
,
Dessens
AB
.
Gender development in Indonesian children, adolescents, and adults with disorders of sex development
.
Arch Sex Behav
.
2015b
44
5
1339
61
.
17.
Ediati
A
,
Juniarto
AZ
,
Birnie
E
,
Drop
SLS
,
Faradz
SMH
,
Dessens
AB
.
Body image and sexuality in Indonesian adults with a disorder of sex development (DSD)
.
J Sex Res
.
2015c
52
1
15
29
.
18.
Ediati
A
,
Juniarto
AZ
,
Birnie
E
,
Okkerse
J
,
Wisniewski
A
,
Drop
S
et al
.
Social stigmatisation in late identified patients with disorders of sex development in Indonesia
.
BMJ Paediatr Open
.
2017
;
1
:
e000130
.
19.
Ediati
A
,
Verrips
GHW
,
Juniarto
AZ
,
Faradz
SMH
,
Drop
SLS
,
Dessens
AB
.
Quality of life in late-treated patients with disorders of sex development: insights for patient-centered care
.
Front Pediatr
.
2018
;
6
:
434
.
20.
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
(
1
):
243
.
21.
Garcia-Acero
M
,
Moreno-Nino
O
,
Suarez-Obando
F
,
Molina
M
,
Manotas
MC
,
Prieto
JC
et al
.
Disorders of sex development: genetic characterization of a patient cohort
.
Mol Med Rep
.
2020
;
21
(
1
):
97
106
.
22.
Globa
E
,
Zelinska
N
,
Shcherbak
Y
,
Bignon-Topalovic
J
,
Bashamboo
A
,
MсElreavey
K
.
Disorders of sex development in a large Ukrainian cohort: clinical diversity and genetic findings
.
Front Endocrinol
.
2022
;
13
:
810782
.
23.
Gomes
NL
,
Batista
RL
,
Nishi
MY
,
Lerário
AM
,
Silva
TE
,
de Moraes Narcizo
A
et al
.
Contribution of clinical and genetic approaches for diagnosing 209 index cases with 46,XY differences of sex development
.
J Clin Endocrinol Metab
.
2022
;
107
(
5
):
e1797
806
.
24.
Hughes
IA
,
Houk
CP
,
Ahmed
SF
,
Lee
PA
LWPES Consensus Group
ESPE Consensus Group
.
Consensus statement on management of intersex disorders
.
Arch Dis Child
.
2006
;
91
(
7
):
554
63
.
25.
Juniarto
AZ
,
van der Zwan
YG
,
Santosa
A
,
Ariani
MD
,
Eggers
S
,
Hersmus
R
et al
.
Drop SL: hormonal evaluation in relation to phenotype and genotype in 286 patients with a disorder of sex development from Indonesia
.
Clin Endocrinol
.
2016
;
85
(
2
):
247
57
.
26.
Juniarto
AZ
,
van der Zwan
YG
,
Santosa
A
,
Hersmus
R
,
de Jong
FH
,
Olmer
R
et al
.
Application of the new classification on patients with a disorder of sex development in Indonesia
.
Int J Endocrinol
.
2012
;
2012
:
237084
.
27.
Krone
N
,
Arlt
W
.
Genetics of congenital adrenal hyperplasia
.
Best Pract Res Clin Endocrinol Metab
.
2009
;
23
(
2
):
181
92
.
28.
Lee
PA
,
Nordenstrom
A
,
Houk
CP
,
Ahmed
SF
,
Auchus
R
,
Baratz
A
et al
.
Global disorders of sex development update since 2006: perceptions, approach and care
.
Horm Res Paediatr
.
2016
;
85
(
3
):
158
80
.
29.
Listyasari
NA
,
Juniarto
AZ
,
Robevska
G
,
Ayers
KL
,
Sinclair
AH
,
Faradz
SMH
.
Analysis of the androgen receptor (AR) gene in a cohort of Indonesian undermasculinized 46, XY DSD patients
.
Egypt J Med Hum Genet
.
2021
;
22
(
1
):
14
.
30.
Listyasari
NA
,
Robevska
G
,
Santosa
A
,
Bouty
A
,
Juniarto
AZ
,
van den Bergen
J
et al
.
Genetic analysis reveals complete androgen insensitivity syndrome in female children surgically treated for inguinal hernia
.
J Invest Surg
.
2019
1
7
.
31.
Marzuki
NS
,
Idris
FP
,
Kartapradja
H
,
Renata
S
,
Harahap
A
,
Batubara
JRL
.
Accuracy of urinary etiocholanolone/androsterone ratio as alternative to serum testosterone/dihydrotestosterone ratio for diagnosis of 5 alpha-reductase type 2 deficiency patients and carriers in Indonesia
.
Int J Endocrinol Metab
.
2021
;
19
:
e109510
.
32.
Marzuki
NS
,
Idris
FP
,
Kartapradja
HD
,
Harahap
AR
,
Batubara
JRL
.
Characterising SRD5A2 gene variants in 37 Indonesian patients with 5-alpha-reductase type 2 deficiency
.
Int J Endocrinol
.
2019
;
2019
:
7676341
.
33.
Moran
ME
,
Karkazis
K
.
Developing a multidisciplinary team for disorders of sex development: planning, implementation, and operation tools for care providers
.
Adv Urol
.
2012
;
2012
:
604135
.
34.
New
MI
,
Abraham
M
,
Gonzalez
B
,
Dumic
M
,
Razzaghy-Azar
M
,
Chitayat
D
et al
.
Genotype-phenotype correlation in 1,507 families with congenital adrenal hyperplasia owing to 21-hydroxylase deficiency
.
Proc Natl Acad Sci USA
.
2013
;
110
(
7
):
2611
6
.
35.
Rey
RA
.
Next-generation sequencing as first-line diagnostic test in patients with disorders of sex development
.
J Clin Endocrinol Metab
.
2022
;
107
(
6
):
e2628
9
.
36.
Speiser
PW
,
Arlt
W
,
Auchus
RJ
,
Baskin
LS
,
Conway
GS
,
Merke
DP
et al
.
Congenital adrenal hyperplasia due to steroid 21-hydroxylase deficiency: an endocrine society clinical practice guideline
.
J Clin Endocrinol Metab
.
2018
;
103
(
11
):
4043
88
.
37.
Veiga-Junior
NN
,
Medaets
PA
,
Petroli
RJ
,
Calais
FL
,
de Mello
MP
,
Castro
CC
et al
.
Clinical and laboratorial features that may differentiate 46,XY DSD due to partial androgen insensitivity and 5α-reductase type 2 deficiency
.
Int J Endocrinol
.
2012
;
2012
:
964876
.
38.
Zhao
F
,
Yao
HH
.
A tale of two tracts: history, current advances, and future directions of research on sexual differentiation of reproductive tracts†
.
Biol Reprod
.
2019
;
101
(
3
):
602
16
.
39.
Zidoune
H
,
Ladjouze
A
,
Chellat-Rezgoune
D
,
Boukri
A
,
Dib
SA
,
Nouri
N
et al
.
Novel genomic variants, atypical phenotypes and evidence of a digenic/oligogenic contribution to disorders/differences of sex development in a large north african cohort
.
Front Genet
.
2022
;
13
:
900574
.