Oligogenic analysis across broad phenotypes of 46,XY differences in sex development associated with NR5A1/SF-1 variants: findings from the international SF1next study.

Background: Oligogenic inheritance has been suggested as a possible mechanism to explain the broad phenotype observed in individuals with differences of sex development (DSD) harbouring NR5A1/SF-1 variants.

Methods: We investigated genetic patterns of possible oligogenicity in a cohort of 30 individuals with NR5A1/SF-1 variants and 46,XY DSD recruited from the international SF1next study, using whole exome sequencing (WES) on family trios whenever available. WES data were analysed using a tailored filtering algorithm designed to identify rare variants in DSD and SF-1-related genes. Identified variants were subsequently tested using the Oligogenic Resource for Variant Analysis (ORVAL) bioinformatics platform for a possible combined pathogenicity with the individual NR5A1/SF-1 variant.

Findings: In 73% (22/30) of the individuals with NR5A1/SF-1 related 46,XY DSD, we identified one to seven additional variants, predominantly in known DSD-related genes, that might contribute to the phenotype. We found identical variants in eight unrelated individuals with DSD in DSD-related genes (e.g., TBCE, FLNB, GLI3 and PDGFRA) and different variants in eight genes frequently associated with DSD (e.g., CDH23, FLNB, GLI2, KAT6B, MYO7A, PKD1, SPRY4 and ZFPM2) in 15 index cases. Our study also identified combinations with NR5A1/SF-1 variants and variants in novel candidate genes.

Interpretation: These findings highlight the complex genetic landscape of DSD associated with NR5A1/SF-1, where in several cases, the use of advanced genetic testing and filtering with specific algorithms and machine learning tools revealed additional genetic hits that may contribute to the phenotype.

Funding: Swiss National Science Foundation and Boveri Foundation Zurich.