Positive predictive values and outcomes for uninformative cell-free DNA tests: An Italian multicentric Cytogenetic and cytogenomic Audit of diagnOstic testing (ICARO study)

OBJECTIVES

To establish the positive predictive values (PPV) of cfDNA testing based on data from a nationwide survey of independent clinical cytogenetics laboratories.

METHODS

Prenatal diagnostic test results obtained by Italian laboratories between 2013 and March 2020 were compiled for women with positive non-invasive prenatal tests (NIPT), without an NIPT result, and cases where there was sex discordancy between the NIPT and ultrasound. PPV and other summary data were reviewed.

RESULTS

Diagnostic test results were collected for 1327 women with a positive NIPT. The highest PPVs were for Trisomy (T) 21 (624/671, 93%) and XYY (26/27, 96.3%), while rare autosomal trisomies (9/47, 19.1%) and recurrent microdeletions (8/55, 14.5%) had the lowest PPVs. PPVs for T21, T18, and T13 were significantly higher when diagnostic confirmation was carried out on chorionic villi (97.5%) compared to amniotic fluid (89.5%) (p < 0.001). In 19/139 (13.9%), of no result cases, a cytogenetic abnormality was detected. Follow-up genetic testing provided explanations for 3/6 cases with a fetal sex discordancy between NIPT and ultrasound.

CONCLUSIONS

NIPT PPVs differ across the conditions screened and the tissues studied in diagnostic testing. This variability, issues associated with fetal sex discordancy, and no results, illustrate the importance of pre- and post-test counselling.

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.

Cerebral folate transporter deficiency syndrome in three siblings: Why genetic testing for developmental and epileptic encephalopathies should be performed early and include the FOLR1 gene

Cerebral folate transporter deficiency syndrome, caused by FOLR-1 mutations is characterized by late infantile onset, severe developmental regression, epilepsy, and leukodystrophy. An extremely low concentration of 5-methyltetrahydrofolate in the cerebrospinal fluid provides a crucial clue to its diagnosis and is a treatment target. Oral or intravenous folinic acid (5-formyltetrahydrofolate) administration improves clinical symptoms and brain magnetic resonance imaging (MRI) findings. We describe three siblings carrying a novel homozygous FOLR1 nonsense mutation, that were referred due to intractable epilepsy and progressive neurological decline. Brain MRI showed hypomyelination and cerebellar atrophy. Folinic acid (oral and intravenous) supplementation, initiated after over 15 years illness, has failed to result in any sizeable clinical or neurophysiological improvement. Cerebral folate transport deficiency bears overlapping clinical features with many severe developmental encephalopathies. It is crucial to recognize FOLR1 signs and establish an early clinical and molecular diagnosis in order to provide timely folinic acid treatment and improve outcome.

Lynch syndrome and Muir-Torre phenotype associated with a recurrent variant in the 3'UTR of the MSH6 gene

A MSH6 3'UTR variant (c.*23_26dup) was found in 13 unrelated families consulted for Lynch/Muir-Torre Syndrome. This variant, which is very rare in the genomic databases, was absent in healthy controls and strongly segregated with the disease in the studied pedigrees. All tumors were defective for MSH2/MSH6/MSH3 proteins expression, but only MSH2 somatic pathogenic mutations were found in 5 of the 12 sequenced tumors. Moreover, we had no evidence of MSH6 transcript decrease in carriers, whereas MSH2 transcript was downregulated. Additional evaluations performed in representative carriers, including karyotype, arrayCGH and Linked-Reads whole genome sequencing, failed to evidence any MSH2 germline pathogenic variant. Posterior probability of pathogenicity for MSH6 c.*23_26dup was obtained from a multifactorial analysis incorporating segregation and phenotypic data and resulted >0.999, allowing to classify the variant as pathogenic (InSiGHT Class 5). Carriers shared a common haplotype involving MSH2/MSH6 loci, then a cryptic disease-associated variant, linked with MSH6 c.*23_26dup, cannot be completely excluded. Even if it is not clear whether the MSH6 variant is pathogenic per se or simply a marker of a disease-associated MSH2/MSH6 haplotype, all data collected on patients and pedigrees prompted us to manage the variant as pathogenic and to offer predictive testing within these families.

SETD5 Gene Haploinsufficiency in Three Patients With Suspected KBG Syndrome

Mendelian disorders of the epigenetic machinery (MDEMs), also named chromatin modifying disorders, are a broad group of neurodevelopmental disorders, caused by mutations in functionally related chromatin genes. Mental retardation autosomal dominant 23 (MRD23) syndrome, due to SETD5 gene mutations, falls into this group of disorders. KBG syndrome, caused by ANKRD11 gene haploinsufficiency, is a chromatin related syndrome not formally belonging to this category. We performed high resolution array CGH and trio-based WES on three molecularly unsolved patients with an initial KBGS clinical diagnosis. A de novo deletion of 116 kb partially involving SETD5 and two de novo frameshift variants in SETD5 were identified in the patients. The clinical re-evaluation of the patients was consistent with the molecular findings, though still compatible with KBGS due to overlapping phenotypic features of KBGS and MRD23. Careful detailed expert phenotyping ascertained some facial and physical features that were consistent with MRD23 rather than KBGS. Our results provide further examples that loss-of-function pathogenic variants in genes encoding factors shaping the epigenetic landscape, lead to a wide phenotypic range with significant clinical overlap. We recommend that clinicians consider SETD5 gene haploinsufficiency in the differential diagnosis of KBGS. Due to overlap of clinical features, careful and detailed phenotyping is important and a large gene panel approach is recommended in the diagnostic workup of patients with a clinical suspicion of KBGS.

High-resolution array-CGH analysis on 46,XX patients affected by early onset primary ovarian insufficiency discloses new genes involved in ovarian function

STUDY QUESTION

Can high resolution array-CGH analysis on a cohort of women showing a primary ovarian insufficiency (POI) phenotype in young age identify copy number variants (CNVs) with a deleterious effect on ovarian function?

SUMMARY ANSWER

This approach has proved effective to clarify the role of CNVs in POI pathogenesis and to better unveil both novel candidate genes and pathogenic mechanisms.

WHAT IS KNOWN ALREADY

POI describes the progression toward the cessation of ovarian function before the age of 40 years. Genetic causes are highly heterogeneous and despite several genes being associated with ovarian failure, most of genetic basis of POI still needs to be elucidated.

STUDY DESIGN, SIZE, DURATION

The current study included 67 46,XX patients with early onset POI (<19 years) and 134 control females recruited between 2012 and 2016 at the Medical Cytogenetics and Molecular Genetics Lab, IRCCS Istituto Auxologico Italiano.

PARTICIPANTS/MATERIALS, SETTING, METHODS

High resolution array-CGH analysis was carried out on POI patients’ DNA. Results of patients and female controls were analyzed to search for rare CNVs. All variants were validated and subjected to a gene content analysis and disease gene prioritization based on the present literature to find out new ovary candidate genes. Case-control study with statistical analysis was carried out to validate our approach and evaluate any ovary CNVs/gene enrichment. Characterization of particular CNVs with molecular and functional studies was performed to assess their pathogenic involvement in POI.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified 37 ovary-related CNVs involving 44 genes with a role in ovary in 32 patients. All except one of the selected CNVs were not observed in the control group. Possible involvement of the CNVs in POI pathogenesis was further corroborated by a case-control analysis that showed a significant enrichment of ovary-related CNVs/genes in patients (P = 0.0132; P = 0.0126). Disease gene prioritization identified both previously reported POI genes (e.g. BMP15, DIAPH2, CPEB1, BNC1) and new candidates supported by transcript and functional studies, such as TP63 with a role in oocyte genomic integrity and VLDLR which is involved in steroidogenesis.

LARGE SCALE DATA

ClinVar database (http://www.ncbi.nlm.nih.gov/clinvar/); accession numbers SCV000787656 to SCV000787743.

LIMITATIONS, REASONS FOR CAUTION

This is a descriptive analysis for almost all of the CNVs identified. Inheritance studies of CNVs in some non-familial sporadic cases was not performed as the parents’ DNA samples were not available. Addionally, RT-qPCR analyses were carried out in few cases as RNA samples were not always available and the genes were not expressed in blood.

WIDER IMPLICATIONS OF THE FINDINGS

Our array-CGH screening turned out to be efficient in identifying different CNVs possibly implicated in disease onset, thus supporting the extremely wide genetic heterogeneity of POI. Since almost 50% of cases are negative rare ovary-related CNVs, array-CGH together with next generation sequencing might represent the most suitable approach to obtain a comprehensive genetic characterization of POI patients.

STUDY FUNDING/COMPETING INTEREST(S)

Supported by Italian Ministry of Health grants ‘Ricerca Corrente’ (08C203_2012) and ‘Ricerca Finalizzata’ (GR-2011-02351636, BIOEFFECT) to IRCCS Istituto Auxologico Italiano.

Testing single/combined clinical categories on 5110 Italian patients with developmental phenotypes to improve array-based detection rate

Background

Chromosomal microarray analysis (CMA) is nowadays widely used in the diagnostic path of patients with clinical phenotypes. However, there is no ascertained evidence to date on how to assemble single/combined clinical categories of developmental phenotypic findings to improve the array‐based detection rate.

Methods

The Italian Society of Human Genetics coordinated a retrospective study which included CMA results of 5,110 Italian patients referred to 17 genetics laboratories for variable combined clinical phenotypes.

Results

Non‐polymorphic copy number variants (CNVs) were identified in 1512 patients (30%) and 615 (32%) present in 552 patients (11%) were classified as pathogenic. CNVs were analysed according to type, size, inheritance pattern, distribution among chromosomes, and association to known syndromes. In addition, the evaluation of the detection rate of clinical subgroups of patients allowed to associate dysmorphisms and/or congenital malformations combined with any other single clinical sign to an increased detection rate, whereas non‐syndromic neurodevelopmental signs and non‐syndromic congenital malformations to a decreased detection rate.

Conclusions

Our retrospective study resulted in confirming the high detection rate of CMA and indicated new clinical markers useful to optimize their inclusion in the diagnostic and rehabilitative path of patients with developmental phenotypes.

Molecular Etiology Disclosed by Array CGH in Patients With Silver-Russell Syndrome or Similar Phenotypes

Introduction: Silver–Russell syndrome (SRS) is an imprinting disorder primarily caused by genetic and epigenetic aberrations on chromosomes 11 and 7. SRS is a rare growth retardation disorder often misdiagnosed due to its heterogeneous and non-specific clinical features. The Netchine–Harbison clinical scoring system (NH-CSS) is the recommended tool for differentiating patients into clinical SRS or unlikely SRS. However, the clinical diagnosis is molecularly confirmed only in about 60% of patients, leaving the remaining substantial proportion of SRS patients with unknown genetic etiology.

Materials and Methods: A cohort of 34 Italian patients with SRS or SRS-like features scored according to the NH-CSS and without any SRS-associated (epi)genetic alterations was analyzed by high-resolution array-based comparative genomic hybridization (CGH) in order to identify potentially pathogenic copy number variants (CNVs).

Results and Discussion: In seven patients, making up 21% of the initial cohort, five pathogenic and two potentially pathogenic CNVs were found involving distinct genomic regions either previously associated with growth delay conditions (1q24.3-q25.3, 17p13.3, 17q22, and 22q11.2-q11.22) and with SRS spectrum (7p12.1 and 7p15.3-p14.3) or outlined for the first time (19q13.42), providing a better definition of reported and as yet unreported SRS overlapping syndromes. All the variants involve genes with a defined role in growth pathways, and for two genes mapping at 7p, IGF2BP3 and GRB10, the association with SRS turns out to be reinforced. The deleterious effect of the two potentially pathogenic variants, comprising GRB10 and ZNF331 genes, was explored by targeted approaches, though further studies are needed to validate their pathogenic role in the SRS etiology. In conclusion, we reconfirm the utility of performing a genome-wide scan to achieve a differential diagnosis in patients with SRS or similar features and to highlight novel chromosome alterations associated with SRS and growth retardation disorders.

9q34.3 microduplications lead to neurodevelopmental disorders through EHMT1 overexpression

Both copy number losses and gains occur within subtelomeric 9q34 region without common breakpoints. The microdeletions cause Kleefstra syndrome (KS), whose responsible gene is EHMT1. A 9q34 duplication syndrome (9q34 DS) had been reported in literature, but it has never been characterized by a detailed molecular analysis of the gene content and endpoints. To the best of our knowledge, we report on the first patient carrying the smallest 9q34.3 duplication containing EHMT1 as the only relevant gene. We compared him with 21 reported patients described here as carrying 9q34.3 duplications encompassing the entire gene and extending within ~ 3 Mb. By surveying the available clinical and molecular cytogenetic data, we were able to discover that similar neurodevelopmental disorders (NDDs) were shared by patient carriers of even very differently sized duplications. Moreover, some facial features of the 9q34 DS were more represented than those of KS. However, an accurate in silico analysis of the genes mapped in all the duplications allowed us to support EHMT1 as being sufficient to cause a NDD phenotype. Wider patient cohorts are needed to ascertain whether the rearrangements have full causative role or simply confer the susceptibility to NDDs and possibly to identify the cognitive and behavioral profile associated with the increased dosage of EHMT1.

13q mosaic deletion including RB1 associated to mild phenotype and no cancer outcome - case report and review of the literature

Background

The 13q deletion syndrome is a rare chromosome disorder associated with wide phenotypic spectrum, which is related to size and location of the deleted region and includes intellectual disability, growth retardation, craniofacial dysmorphisms, congenital malformations, and increased risk of retinoblastoma.

Case presentation

Here, we report on a teenage boy with a mild phenotype characterized by obesity, hyperactivity, dysphagia, dysgraphia, sleep disturbance, and minor dysmorphic features (round face, bushy eyebrows, and stubby hands). Array Comparative Genomic Hybridization on blood identified a mosaic 13q14.13-13q31.1 deletion, with a mosaicism rate around 40%, which was confirmed by quantitative PCR and interphase Fluorescent In Situ Hybridization (iFISH) on both blood genomic DNA and cultured/uncultured blood lymphocytes, respectively. Conversely, karyotype analysis on blood estimated a mosaicism rate of 24% and iFISH on buccal smears revealed a borderline value of 0.4%, suggesting the absence of 13q deletion in this cell line.

Conclusions

The comparison with previous patients carrying similar deletions informed that the proband clinical presentation is the mildest reported to date, thus supporting the burden of mosaicism in modulating the phenotype also in case of large chromosomal rearrangements. Characterization of further cases by in-depth mosaicism rate in tissues with different embryonic origins might contribute in the future to a better definition of genotype-phenotype correlation, including tumor risk.

A balanced reciprocal translocation t(10;15)(q22.3;q26.1) interrupting ACAN gene in a family with proportionate short stature

PURPOSE

Few examples of the involvement of a single gene in idiopathic short stature have been described until now. Our aim was to identify the causative gene of proportionate short stature in a large family showing co-segregation of the phenotype with the reciprocal translocation t(10;15)(q22;q24).

METHODS

FISH mapping was carried out with BACs and long-range PCR probes to identify the smallest genomic regions harboring the translocation breakpoints. Real-Time RT-PCR was performed in blood after pre-amplification of target genes cDNA.

RESULT

The affected family members presented with a final height of between - 2.41 and - 4.18 SDS and very mild skeletal dysmorphisms. Growth rates of the proband and of her cousin, whose childhood and pre-pubertal bone age corresponded to the chronological age, showed a poor growth spurt during treatment with rhGH. However, their adult height was greater than that of their untreated mothers, suggesting efficacy of GH therapy. Breakpoint mapping revealed that the translocation t(10;15)(q22.3;q26.1) disrupts, on 15q, the ACAN gene at intron 1, decreasing its transcriptional expression.

CONCLUSIONS

This is the first description of a chromosome rearrangement disrupting ACAN and leading to its haploinsufficiency. ACAN loss of function should be considered a potential underpinning of short patients who display a poor growth spurt and belong to families with autosomal dominant segregation of proportionate short stature. Besides this core phenotype, literature review suggests that advanced bone age, early onset osteochondritis dissecans, osteoarthritis, intervertebral disc disease as well as craniofacial dysmorphisms can be important suggestive phenotypes in affected families.

Constitutional de novo deletion of the FBXW7 gene in a patient with focal segmental glomerulosclerosis and multiple primitive tumors

Multiple primary malignant neoplasms are rare entities in the clinical setting, but represent an important issue in the clinical management of patients since they could be expression of a genetic predisposition to malignancy. A high resolution genome wide array CGH led us to identify the first case of a de novo constitutional deletion confined to the FBXW7 gene, a well known tumor suppressor, in a patient with a syndromic phenotype characterized by focal segmental glomerulosclerosis and multiple primary early/atypical onset tumors, including Hodgkin's lymphoma, Wilms tumor and breast cancer. Other genetic defects may be associated with patient's phenotype. In this light, constitutional mutations at BRCA1, BRCA2, TP53, PALB2 and WT1 genes were excluded by performing sequencing and MLPA analysis; similarly, we ruled out constitutional abnormalities at the imprinted 11p15 region by methylation specific -MLPA assay. Our observations sustain the role of FBXW7 as cancer predisposition gene and expand the spectrum of its possible associated diseases.

Familial intragenic duplication of ANKRD11 underlying three patients of KBG syndrome

BACKGROUND

KBG syndrome, a rare autosomal disorder characterised by distinctive craniofacial and skeletal features and developmental delay, is caused by haploinsufficiency of the ANKRD11 gene.

RESULTS

Here we describe two siblings with multiple symptoms characteristic of KBG and their mother with a milder phenotype. In the siblings, array-based comparative genomic hybridization (array CGH) identified an intragenic microduplication affecting ANKRD11 that was absent from the parents' array CGH profiles. Microsatellite analysis revealed the maternal origin of the rearrangement and interphase fluorescent in situ hybridization (i-FISH) experiments identified the rearrangement in low-level mosaicism in the mother. Molecular characterisation of the duplication allele demonstrated the presence of two mutant ANKRD11 transcripts containing a premature stop codon and predicting a truncated non-functional protein.

CONCLUSIONS

Similarly to deletions and point mutations, this novel pathogenetic rearrangement causes haploinsufficiency of ANKRD11, resulting in KBG syndrome.

New case of trichorinophalangeal syndrome-like phenotype with a de novo t(2;8)(p16.1;q23.3) translocation which does not disrupt the TRPS1 gene

BACKGROUND

Trichorhinophalangeal syndrome (TRPS) is a rare autosomal dominant genetic disorder characterised by distinctive craniofacial and skeletal abnormalities. TRPS is generally associated with mutations in the TRPS1 gene at 8q23.3 or microdeletions of the 8q23.3-q24.11 region. However, three deletions affecting the same chromosome region and a familial translocation t(8;13) co-segregating with TRPS, which do not encompass or disrupt the TRPS1 gene, have been reported. A deregulated expression of TRPS1 has been hypothesised as cause of the TRPS phenotype of these patients.

CASE PRESENTATION

We report the clinical and molecular characterisation of a 57-year-old Caucasian woman carrying the t(2;8)(p16.1;q23.3) de novo balanced translocation. The proband presented with peculiar clinical features (severe craniofacial dysmorphism, alopecia universalis, severe scoliosis, mitral valve prolapse, mild mental impairment and normal growth parameters) that partially overlap with TRPS I. Mutational and array CGH analyses ruled out any genetic defect affecting TRPS1 or genomic alteration at the translocation breakpoint or elsewhere in the genome. Breakpoint mapping excluded disruption of TRPS1, and revealed that the chromosome 8q23.3 breakpoint was located within the IVS10 of the long intergenic non-coding RNA LINC00536, at approximately 300 kb from the TRPS1 5' end. Conversely, the 2p16.1 breakpoint mapped within a LINE sequence, in a region that lacks transcriptional regulatory elements. As a result of the translocation, nucleotide base pair additions and deletions were detected at both breakpoint junction fragments, and an evolutionarily conserved VISTA enhancer element from 2p16.1 was relocated at approximately 325 kb from the TRPS1 promoter.

CONCLUSIONS

We suggest that the disruption of the genomic architecture of cis regulatory elements downstream the TRPS1 5' region, combined with the translocation of a novel enhancer element nearby TRPS1, might be the pathogenetic mechanism underpinning the proband's phenotype. The clinical and genetic characterisation of the present subject allowed us to make a genetic diagnosis in the context of a known syndrome, contributing to a better comprehension of the complex transcriptional regulation of TRPS1 and TRPS ethiopathogenesis.