MATERIAL AND METHODS: All the information for CYP1B1 missense variants was retrieved from the dbSNP database. Seven different tools, namely: SIFT, PolyPhen-2, PROVEAN, SNAP2, PANTHER, PhD-SNP, and Predict-SNP, were used for functional annotation, and two packages, which were I-Mutant 2.0 and MUpro, were used to predict the effect of the variants on protein stability. A phylogenetic conservation analysis using deleterious variants was performed by the ConSurf server. The 3D structures of the wild-type and mutants were generated using the I-TASSER tool, and a 50 ns molecular dynamic simulation (MDS) was executed using the GROMACS webserver to determine the stability of mutants compared to the native protein. Co-expression, protein-protein interaction (PPI), gene ontology (GO), and pathway analyses were additionally performed for the CYP1B1 in-depth study.

RESULTS: All the retrieved data from the dbSNP database was subjected to functional, structural, and phylogenetic analysis. From the conducted analyses, a total of 19 high-risk variants (P52L, G61E, G90R, P118L, E173K, D291G, Y349D, G365W, G365R, R368H, R368C, D374N, N423Y, D430E, P442A, R444Q, F445L, R469W, and C470Y) were screened out that were considered to be deleterious to the CYP1B1 gene. The phylogenetic analysis revealed that the majority of the variants occurred in highly conserved regions. The MD simulation analysis exhibited that all mutants' average root mean square deviation (RMSD) values were higher compared to the wild-type protein, which could potentially cause CYP1B1 protein dysfunction, leading to the severity of the disease. Moreover, it has been discovered that CYP1A1, VCAN, HSD17B1, HSD17B2, and AKR1C3 are highly co-expressed and interact with CYP1B1. Besides, the CYP1B1 protein is primarily involved in the metabolism of xenobiotics, chemical carcinogenesis, the retinal metabolic process, and steroid hormone biosynthesis pathways, demonstrating its multifaceted and important roles.

METHODS: We utilized a combination of chromosome microarray, exome sequencing, and genome sequencing with structural variant and trio analysis to investigate a cohort of 41 predominantly sporadic cases.

RESULTS: We identified likely causative variants in 54% (22/41) of cases, including 51% (19/37) of sporadic cases and 75% (3/4) of cases initially referred as familial ASD. Two-thirds of sporadic cases were found to have heterozygous variants, which in most cases were de novo. Approximately one-third (7/22) of genetic diagnoses were found in rarely reported or recently identified ASD genes including PXDN, GJA8, COL4A1, ITPR1, CPAMD8, as well as the new phenotypic association of Axenfeld-Rieger anomaly with a homozygous ADAMTS17 variant. The remainder of the variants were in key ASD genes including FOXC1, PITX2, CYP1B1, FOXE3, and PAX6.