Targeted whole exome sequencing and Drosophila modelling to unveil the molecular basis of primary ovarian insufficiency.

STUDY QUESTION: Can a targeted whole exome sequencing (WES) on a cohort of women showing a primary ovarian insufficiency (POI) phenotype at a young age, combined with a study of copy number variations, identify variants in candidate genes confirming their deleterious effect on ovarian function? SUMMARY ANSWER: This integrated approach has proved effective in identifying novel candidate genes unveiling mechanisms involved in POI pathogenesis. WHAT IS KNOWN ALREADY: POI, a condition occurring in 1% of women under 40 years of age, affects women's fertility leading to a premature loss of ovarian reserve. The genetic causes of POI are highly heterogeneous and several determinants contributing to its prominent oligogenic inheritance pattern still need to be elucidated. STUDY DESIGN, SIZE, DURATION: WES screening for pathogenic variants of 41 Italian women with non-syndromic primary and early secondary amenorrhoea occurring before age 25 was replicated on another 60 POI patients, including 35 French and 25 American women, to reveal statistically significant shared variants. PARTICIPANTS/MATERIALS, SETTING, METHODS: The Italian POI patients' DNA were processed by targeted WES including 542 RefSeq genes expressed or functioning during distinct reproductive or ovarian processes (e.g. DNA repair, meiosis, oocyte maturation, folliculogenesis and menopause). Extremely rare variants were filtered and selected by means of a Fisher Exact test using several publicly available datasets. A case-control Burden test was applied to highlight the most significant genes using two ad-hoc control female cohorts. To support the obtained data, the identified genes were screened on a novel cohort of 60 Caucasian POI patients and the same case-control analysis was carried out. Comparative analysis of the human identified genes was performed on mouse and Drosophila melanogaster by analysing the orthologous genes in their ovarian phenotype, and two of the selected genes were fruit fly modelled to explore their role in fertility. MAIN RESULTS AND THE ROLE OF CHANCE: The filtering steps applied to search for extremely rare pathogenic variants in the Italian cohort revealed 64 validated single-nucleotide variants/Indels in 59 genes in 30 out of 41 screened women. Burden test analysis highlighted 13 ovarian genes as being the most enriched and significant. To validate these findings, filtering steps and Burden analysis on the second cohort of Caucasian patients yielded 11 significantly enriched genes. Among them, AFP, DMRT3, MOV10, FYN and MYC were significant in both patient cohorts and hence were considered strong candidates for POI. Mouse and Drosophila comparative analysis evaluated a conserved role through the evolution of several candidates, and functional studies using a Drosophila model, when applicable, supported the conserved role of the MOV10 armitage and DMRT3 dmrt93B orthologues in female fertility. LARGE SCALE DATA: The datasets for the Italian cohort generated during the current study are publicly available at ClinVar database (http://www.ncbi.nlm.nih.gov/clinvar/): accession numbers SCV001364312 to SCV001364375. LIMITATIONS, REASONS FOR CAUTION: This is a targeted WES analysis hunting variants in candidate genes previously identified by different genomic approaches. For most of the investigated sporadic cases, we could not track the parental inheritance, due to unavailability of the parents' DNA samples; in addition, we might have overlooked additional rare variants in novel candidate POI genes extracted from the exome data. On the contrary, we might have considered some inherited variants whose clinical significance is uncertain and might not be causative for the patients' phenotype. Additionally, as regards the Drosophila model, it will be extremely important in the future to have more mutants or RNAi strains available for each candidate gene in order to validate their role in POI pathogenesis. WIDER IMPLICATIONS OF THE FINDINGS: The genomic, statistical, comparative and functional approaches integrated in our study convincingly support the extremely heterogeneous oligogenic nature of POI, and confirm the maintenance across the evolution of some key genes safeguarding fertility and successful reproduction. Two principal classes of genes were identified: (i) genes primarily involved in meiosis, namely in synaptonemal complex formation, asymmetric division and oocyte maturation and (ii) genes safeguarding cell maintenance (piRNA and DNA repair pathways). STUDY FUNDING/COMPETING INTEREST(S): This work was supported by Italian Ministry of Health grants 'Ricerca Corrente' (08C621_2016 and 08C924_2019) provided to IRCCS Istituto Auxologico Italiano, and by 'Piano Sostegno alla Ricerca' (PSR2020_FINELLI_LINEA_B) provided by the University of Milan; M.P.B. was supported by Telethon-Italy (grant number GG14181). There are no conflicts of interest.

High-resolution microarray analysis unravels complex Xq28 aberrations in patients and carriers affected by X-linked blue cone monochromacy.

The human X chromosome contains ∼ 1600 genes, about 15% of which have been associated with a specific genetic condition, mainly affecting males. Blue cone monochromacy (BCM) is an X-linked condition caused by a loss-of-function of both the OPN1LW and OPN1MW opsin genes. The cone opsin gene cluster is composed of 2-9 paralogs with 99.8% sequence homology and is susceptible to deletions, duplications, and mutations. Current diagnostic tests employ polymerase chain reaction (PCR)-based technologies; however, alterations remain undetermined in 10% of patients. Furthermore, carrier testing in females is limited or unavailable. High-resolution X chromosome-targeted CGH microarray was applied to test for rearrangements in males with BCM and female carriers from three unrelated families. Pathogenic alterations were revealed in all probands, characterized by sequencing of the breakpoint junctions and quantitative real-time PCR. In two families, we identified a novel founder mutation that consisted of a complex 3-kb deletion that embraced the cis-regulatory locus control region and insertion of an additional aberrant OPN1MW gene. The application of high-resolution X-chromosome microarray in clinical diagnosis brings significant advantages in detection of small aberrations that are beyond the resolution of clinically available aCGH analysis and which can improve molecular diagnosis of the known conditions and unravel previously unrecognized X-linked diseases.

Novel Inactivating Mutation of the FSH Receptor in Two Siblings of Indian Origin With Premature Ovarian Failure.

CONTEXT: Inactivating FSH receptor (FSHR) mutations can affect ovarian function, resulting in variable clinical presentations ranging from primary amenorrhea to premature menopause. FSHR mutations have been largely reported in the Finnish population, but in patients of Asian Indian descent, the incidence of FSHR mutations is extremely rare. CASE DESCRIPTION: Two female siblings of Indian descent were diagnosed with primary ovarian failure and hypergonadotropic hypogonadism. The daughters were the result of a consanguineous marriage between second cousins. A combination of comparative genomic hybridization plus single nucleotide polymorphism array and whole exome sequencing was conducted on the family to identify potential causative genetic variants. CONCLUSION: Both daughters were found to have a novel pathogenic variant in FSHR (c.1253T>G, p.Ile418Ser), inherited as an autosomal recessive trait from heterozygous parents. This loss of function mutation is located in exon 10 of FSHR affecting the second transmembrane helix of the FSHR protein. The transmembrane domain of FSHR is highly conserved across species and is involved in signal transduction. The FSHR c.1253T>G variant is next to a known pathogenic variant, rs12190966 (c.1255G>A, p.Ala419Thr), previously reported in a Finnish woman with primary amenorrhea.

Application of chromosomal microarray in the evaluation of abnormal prenatal findings.

We performed karyotype and array comparative genomic hybridization (aCGH) analyses on 177 prenatal samples, including 162 (92%) samples from fetuses with sonographic anomalies. Overall 12 fetuses (6.8%) had abnormal karyotype and 42 (23.7%) fetuses had abnormal microarray results: 20 (11.3%) with pathogenic copy number variations (CNVs), 16 with CNVs of uncertain clinical significance, 4 with CNVs establishing carrier status for recessive, X-linked, or susceptibility to late onset dominant disease, and two CNVs with pseudomosaicism due to in vitro cultural artifacts. For 23 pregnancies (13%), aCGH contributed important new information. Our results highlight the interpretation challenges associated with CNVs of unclear significance, incidental findings, as well as technical aspects. Array CGH analysis significantly improved the detection of genomic imbalances in prenatal diagnosis of pregnancies with structural birth defects.

A Case of Agonadism, Skeletal Malformations, Bicuspid Aortic Valve, and Delayed Development with a 16p13.3 Duplication Including GNG13 and SOX8 Upstream Enhancers: Are Either, Both or Neither Involved in the Phenotype?

We report a female patient with delayed growth and development, skeletal and cardiac defects, and a male XY sex chromosome complement with early failure of gonad development. SRY sequencing was normal. Array comparative genome hybridization (CGH) analysis revealed a gain in copy number in the subtelomeric region of the short arm of chromosome 16, encompassing a region of approximately 560 kb in size including GNG13 which may be involved in ovarian development. The proximal breakpoint of the duplication maps about 18 kb upstream of SOX8 and involves evolutionary conserved regulatory elements. SOX8, like SOX9, is a transcription factor expressed in many tissues, including neural crest, nervous system, muscle, cartilage, adrenal gland, kidney, and testis. There was no increase in GNG13 or SOX8 expression in the patient's lymphoblastoid line. It is possible that an alteration of SOX8 or/and GNG13 expression is responsible for the multiple congenital anomalies and sex reversal in our patient.

Identification of critical regions for clinical features of distal 10q deletion syndrome.

Array comparative genomic hybridization studies were performed to further characterize cytogenetic abnormalities found originally by karyotype and fluorescence in situ hybridization in five clinical cases of distal 10q deletions, including several with complex cytogenetic rearrangements and one with a partial male-to-female sex-reversal phenotype. These results have enabled us to narrow the previously proposed critical regions for the craniofacial, urogenital, and neuropsychiatric disease-related manifestations associated with distal 10q deletion syndrome. Furthermore, we propose that haploinsufficiency of the DOCK1 gene may play a crucial role in the pathogenesis of the 10q deletion syndrome. We hypothesize that alteration of DOCK1 and/or other genes involved in regulation and signaling of multiple pathways can explain the wide range of phenotypic variability between patients with similar or identical cytogenetic abnormalities.

Complex chromosome 17p rearrangements associated with low-copy repeats in two patients with congenital anomalies.

Recent molecular cytogenetic data have shown that the constitution of complex chromosome rearrangements (CCRs) may be more complicated than previously thought. The complicated nature of these rearrangements challenges the accurate delineation of the chromosomal breakpoints and mechanisms involved. Here, we report a molecular cytogenetic analysis of two patients with congenital anomalies and unbalanced de novo CCRs involving chromosome 17p using high-resolution array-based comparative genomic hybridization (array CGH) and fluorescent in situ hybridization (FISH). In the first patient, a 4-month-old boy with developmental delay, hypotonia, growth retardation, coronal synostosis, mild hypertelorism, and bilateral club feet, we found a duplication of the Charcot-Marie-Tooth disease type 1A and Smith-Magenis syndrome (SMS) chromosome regions, inverted insertion of the Miller-Dieker lissencephaly syndrome region into the SMS region, and two microdeletions including a terminal deletion of 17p. The latter, together with a duplication of 21q22.3-qter detected by array CGH, are likely the unbalanced product of a translocation t(17;21)(p13.3;q22.3). In the second patient, an 8-year-old girl with mental retardation, short stature, microcephaly and mild dysmorphic features, we identified four submicroscopic interspersed 17p duplications. All 17 breakpoints were examined in detail by FISH analysis. We found that four of the breakpoints mapped within known low-copy repeats (LCRs), including LCR17pA, middle SMS-REP/LCR17pB block, and LCR17pC. Our findings suggest that the LCR burden in proximal 17p may have stimulated the formation of these CCRs and, thus, that genome architectural features such as LCRs may have been instrumental in the generation of these CCRs.

Two novel translocation breakpoints upstream of SOX9 define borders of the proximal and distal breakpoint cluster region in campomelic dysplasia.

The semilethal skeletal malformation syndrome campomelic dysplasia (CD) with or without XY sex reversal is caused by mutations within the SOX9 gene on 17q24.3 or by chromosomal aberrations (translocations, inversions or deletions) with breakpoints outside the SOX9 coding region. The previously published CD translocation breakpoints upstream of SOX9 fall into two clusters: a proximal cluster with breakpoints between 50-300 kb and a distal cluster with breakpoints between 899-932 kb. Here, we present clinical, cytogenetic and molecular data from two novel CD translocation cases. Case 1 with karyotype 46,XY,t(1;17)(q42.1;q24.3) has characteristic symptoms of CD, including mild tibial bowing, cryptorchidism and hypospadias. By standard fluorescence in situ hybridization (FISH) and by high-resolution fiber FISH, the 17q breakpoint was mapped 375 kb from SOX9, defining the centromeric border of the proximal breakpoint cluster region. Case 2 with karyotype 46,X,t(Y;17)(q11.2;q24.3) has the acampomelic form of CD and complete XY sex reversal. By FISH and somatic cell hybrid analysis, the 17q breakpoint was mapped 789 kb from SOX9, defining the telomeric border of the distal breakpoint cluster region. We discuss the structure of the 1 Mb cis-control region upstream of SOX9 and the correlation between the position of the 14 mapped translocation breakpoints with respect to disease severity and XY sex reversal.

A natural history of a child with monosomy 5p syndrome (Cat-cry/Cri-du-chat syndrome) during the 18 years of follow-up.

A record of a natural history of a long-term case study devoted to monosomy 5p (Cat-cry/Cri-du-chat) syndrome has been described rarely. Knowledge on the range of the changes in phenotype attributable to advancing age can be useful in clinical diagnosis of monosomy 5p at the different developmental stages, including adolescence, as well in prognosis for genetic counseling. In this case a detailed analysis of the morphologic phenotype in a girl with del(5)(p13.3) observed from 4 months to 18 years of age is reported. The comparative analysis of the girl's phenotype in different developmental stages has revealed that microcephaly, flat occipital region, face asymmetry, wide spaced palpebral fissures, epicanthic folds, small mouth fissure, thin mucous lip, small and low set ears and short IV metacarpals has not changed with advancing age. However, facial asymmetry was more evident, frontal tubers were less prominent, nasal root and back became prominent nasal back became elongated, the subnasal region was shorter and marked malocclusion appeared.

Interstitial deletion of 10p and atrial septal defect in DiGeorge 2 syndrome.

We present molecular genetic investigations of a 4-year-old boy with craniofacial dysmorphism and developmental delay. Trivial mitral and tricuspid regurgitation without gross structural abnormality was diagnosed by echocardiography. High-resolution chromosome analysis revealed an interstitial deletion, del(10)(p12.1p12.32). To characterize the deletion size and breakpoints, we performed fluorescence in situ hybridization analysis using 27 BAC clones. Our data demonstrate an approximately 5.5 Mb deletion del(10)(p12.1p12.31). Surprisingly, the BAC clone RP11-56H7 that contains NEBL, an apparent downstream gene of the cardiogenic transcription factor HAND2 previously shown to be deleted in the patients with DiGeorge 2 syndrome and 10p13 deletion, was deleted in our patient with 10p12.1-p12.31 deletion. In addition, we provide clinical data and results of molecular analysis for a patient with multiple congenital anomalies including Ebstein's anomaly, kidney malformations, and 10p13-p14 deletion. We also reviewed 19 patients with congenital heart defects and deletions involving 10p and propose that atrial septal defect (ASD) is a common cardiac anomaly associated with DiGeorge 2 syndrome. Based on genotype-phenotype analysis of published patients and those reported herein, we propose an approximately 1.0 Mb critical region between loci D10S547 and D10S2176 in 10p14 to be associated with ASD. Considering that septal defects are the most frequent congenital heart anomaly, we suggest that further investigations in the 10p critical region are important to identify gene(s) responsible for this common birth defect.

Omphalocele in trisomy 3q: further delineation of phenotype.

We report a case of a patient with omphalocele, dysmorphic features, and mild developmental delay associated with a chromosomal aberration. Chromosome studies showed that the propositus carries a maternally derived unbalanced translocation der(4)t(3;4)(q27.3;q32.3), resulting in trisomy for region 3q27.3-->qter and monosomy for 4q32.3-->qter. Because the association between dup3q and omphalocele has been reported in several cases, we analyzed the data on 93 previously reported patients with partial trisomy of the long arm of chromosome 3 and compared the clinical features between the cases. The imbalance of chromosome 3 in the patient was further defined by fluorescence in situ hybridization (FISH) studies using bacterial artificial chromosome (BAC) clones. BAC clone RP11-171N2 was identified as a breakpoint-spanning clone in the patient and his mother. Based on our comparative analysis, we have delineated that the smallest region of overlap (SRO) associated with omphalocele is from BAC 171N2 to 3qter. We hypothesize that the SRO contains a gene(s) important in normal abdominal wall development and is of potential interest for further investigation.