Implication of androgen receptor gene dysfunction in human Mllerian duct anomalies

Background: Mullerian duct anomalies (MDAs) are congenital developmental disorders exhibiting as a variety of malformations of female reproductive tract. The identified etiology of MDAs is limited. The present study aimed to unravel the underlying genetic causes of MDAs. Methods: Rare variants in androgen receptor (AR) were called from the cohort consists of patients with MDAs and underwent whole exome sequencing (WES) at Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing, China. Sanger sequencing was used to confirm the causative genetic mutations. In silico analysis were used to classify the pathogenicity of each variant. Molecular modeling and simulations were conducted to investigate the conformational changes between the wild-type (WT) and mutant proteins. Results: A total of 3 rare heterozygous variants in AR from the MDAs cohort in our institution were identified, with unknown effects. All variants were missense mutations, including c.173A > T, c.558C > A and c.1208C > T, and were absent or rare in East Asian populations in Genome Aggregation Database and the Exome Aggregation Consortium Database. According to the American College of Medical Genetics and Genomics guidelines, c.1208C > T variant was classified as likely pathogenic, while the other two were variants of uncertain significance. During molecular dynamics simulations, WT and mutant proteins all reached stable status according to rootmean-square variance. Values of radius of gyration showed that Q58L and S186R protein would be more compact than WT, while the structure of A403V became looser. Despite, in comparison with WT, the number of hydrogen bonds increased in Q58L, while decreased in the other two variants. Furthermore, the solventaccessible surface area diminished in Q58L and A403V while enlarged in S186R proteins, when compared with WT. Conclusions: To our knowledge, this is the first report regarding the association of AR mutation and MDAs. The identification of these variants, predicted to damage the structure and function of AR protein, not only expanded the mutational spectrum of causative genes of MDAs, but provide novel molecular genetic reference for future studies.