Development of a comprehensive noninvasive prenatal test.

Our aim was to develop and apply a comprehensive noninvasive prenatal test (NIPT) by using high-coverage targeted next-generation sequencing to estimate fetal fraction, determine fetal sex, and detect trisomy and monogenic disease without parental genotype information. We analyzed 45 pregnancies, 40 mock samples, and eight mother-child pairs to generate 35 simulated datasets. Fetal fraction (FF) was estimated based on analysis of the single nucleotide polymorphism (SNP) allele fraction distribution. A Z-score was calculated for trisomy of chromosome 21 (T21), and fetal sex detection. Monogenic disease detection was performed through variant analysis. Model validation was performed using the simulated datasets. The novel model to estimate FF was robust and accurate (r2= 0.994, p-value < 2.2e-16). For samples with FF > 0.04, T21 detection had 100% sensitivity (95% CI: 63.06 to 100%) and 98.53% specificity (95% CI: 92.08 to 99.96%). Fetal sex was determined with 100% accuracy. We later performed a proof of concept for monogenic disease diagnosis of 5/7 skeletal dysplasia cases. In conclusion, it is feasible to perform a comprehensive NIPT by using only data from high coverage targeted sequencing, which, in addition to detecting trisomies, also make it possible to identify pathogenic variants of the candidate genes for monogenic diseases.

Critical points for an accurate human genome analysis.

Next-generation sequencing is radically changing how DNA diagnostic laboratories operate. What started as a single-gene profession is now developing into gene panel sequencing and whole-exome and whole-genome sequencing (WES/WGS) analyses. With further advances in sequencing technology and concomitant price reductions, WGS will soon become the standard and be routinely offered. Here, we focus on the critical steps involved in performing WGS, with a particular emphasis on points where WGS differs from WES, the important variables that should be taken into account, and the quality control measures that can be taken to monitor the process. The points discussed here, combined with recent publications on guidelines for reporting variants, will facilitate the routine implementation of WGS into a diagnostic setting.

The FSHD2 gene SMCHD1 is a modifier of disease severity in families affected by FSHD1.

Facioscapulohumeral muscular dystrophy type 1 (FSHD1) is caused by contraction of the D4Z4 repeat array on chromosome 4 to a size of 1-10 units. The residual number of D4Z4 units inversely correlates with clinical severity, but significant clinical variability exists. Each unit contains a copy of the DUX4 retrogene. Repeat contractions are associated with changes in D4Z4 chromatin structure that increase the likelihood of DUX4 expression in skeletal muscle, but only when the repeat resides in a genetic background that contains a DUX4 polyadenylation signal. Mutations in the structural maintenance of chromosomes flexible hinge domain containing 1 (SMCHD1) gene, encoding a chromatin modifier of D4Z4, also result in the increased likelihood of DUX4 expression in individuals with a rare form of FSHD (FSHD2). Because SMCHD1 directly binds to D4Z4 and suppresses somatic expression of DUX4, we hypothesized that SMCHD1 may act as a genetic modifier in FSHD1. We describe three unrelated individuals with FSHD1 presenting an unusual high clinical severity based on their upper-sized FSHD1 repeat array of nine units. Each of these individuals also carries a mutation in the SMCHD1 gene. Familial carriers of the FSHD1 allele without the SMCHD1 mutation were only mildly affected, suggesting a modifier effect of the SMCHD1 mutation. Knocking down SMCHD1 in FSHD1 myotubes increased DUX4 expression, lending molecular support to a modifier role for SMCHD1 in FSHD1. We conclude that FSHD1 and FSHD2 share a common pathophysiological pathway in which the FSHD2 gene can act as modifier for disease severity in families affected by FSHD1.

Noninvasive prenatal diagnosis of Huntington disease: detection of the paternally inherited expanded CAG repeat in maternal plasma.

OBJECTIVE: With a shift towards noninvasive testing, we have explored and validated the use of noninvasive prenatal diagnosis (NIPD) for Huntington disease (HD). METHODS: Fifteen couples have been included, assessing a total of n = 20 pregnancies. Fetal paternally inherited CAG repeat length was determined in total cell-free DNA from maternal plasma using a direct approach by PCR and subsequent fragment analysis. RESULTS: All fetal HD (n = 7) and intermediate (n = 3) CAG repeats could be detected in maternal plasma. Detection of repeats in the normal range (n = 10) was successful in n = 5 cases where the paternal repeat size could be distinguished from maternal repeat patterns after fragment analysis. In all other cases (n = 5), the paternal peaks coincided with the maternal peak pattern. All NIPD results were concordant with results from routine diagnostics on fetal genomic DNA from chorionic villi. CONCLUSION: In this validation study, we demonstrated that all fetuses at risk for HD could be identified noninvasively in maternal plasma. Additionally, we have confirmed results from previously described case reports that NIPD for HD can be performed using a direct approach by PCR. For future diagnostics, parental CAG profiles can be used to predict the success rate for NIPD prior to testing.

Broader spectrum of ß-thalassemia mutations in Oman: regional distribution and comparison with neighboring countries.

The objective of this study was to expand and study the molecular spectrum of ß-thalassemia (ß-thal) mutations in Oman by examining cases from seven different regions and comparing the prevalence with neighboring countries. A total of 446 cases of ß hemoglobinopathies was obtained and analyzed to determine the frequency and distribution of the different ß alleles. The molecular spectrum of ß-thal in Oman revealed the presence of 32 mutations from different origins and 11 alleles are reported for the first time in the Omani population. The wide heterogeneous spectrum of ß-thal mutations found can be associated with the history of trade and migration as well as the past domination from other countries. The presented data will facilitate the development of a comprehensive prevention strategy in Oman.

Molecular spectrum of α-globin gene defects in the Omani population.

We describe the molecular characterization of α-globin gene defects in a cohort of 634 Omani patients. A total of 21 different α gene mutations were found in 484 subjects. Overall, we identified three different large deletions, three small deletions, 11 point mutations [two on the α2 polyadenylation signal (polyA) (HBA2: c.*94A>G), and nine α chain variants], three ααα(anti 3.7) triplication, a 21 nucleotide (nt) duplication on the α1 gene and two novel (presumed) polymorphisms on the α 3.7 kb hybrid gene, namely -5 (C>T) and +46 (C>A). Of these defects, 15 have not been previously reported in the Omani population. This large heterogeneity of α-thalassemia (α-thal) observed in the Omani population could be expected in neighboring Arab countries. The high frequency of α-thal, solely or in association with ß-globin gene defects, emphasize the necessity of adding α-thal testing to pre marital programs for accurate genetic counseling.

Known and new δ-globin gene mutations and other factors influencing Hb A2 measurement in the Omani population.

Although δ-thalassemia (δ-thal) is not categorized as a severe disease, it is essential to know the molecular spectrum of the δ gene mutations frequently occurring in specific areas, particularly if these areas are characterized by a high rate of ß-thalassemia (ß-thal) such as Oman. This is because coinherited δ-globin gene defects can interfere with the basic diagnosis of a ß-thal carrier when this is based upon the measurement of the Hb A2 only. Because of that, we have investigated 33 patients with low Hb A2 levels, collected from different hospitals in Oman. Some cases had a second Hb A2 fraction, while others had only significantly lower Hb A2 levels. Among these patients, 20 did carry a δ-globin gene mutation, the rest were carriers of α thalassemia (α-thal) defects or could be iron depleted or both. In total, eight different known mutations and two novel δ variants were found. The characterization of the δ-globin gene mutation spectrum will improve carrier diagnostics and genetic counseling in the Omani population screened for ß-thal.

Successful noninvasive trisomy 18 detection using single molecule sequencing.

BACKGROUND: Noninvasive trisomy 21 detection performed by use of massively parallel sequencing is achievable with high diagnostic sensitivity and low false-positive rates. Detection of fetal trisomy 18 and 13 has been reported as well but seems to be less accurate with the use of this approach. The reduced accuracy can be explained by PCR-introduced guanine-cytosine (GC) bias influencing sequencing data. Previously, we demonstrated that sequence data generated by single molecule sequencing show virtually no GC bias and result in a more pronounced noninvasive detection of fetal trisomy 21. In this study, single molecule sequencing was used for noninvasive detection of trisomy 18 and 13. METHODS: Single molecule sequencing was performed on the Helicos platform with free DNA isolated from maternal plasma from 11 weeks of gestation onward (n = 17). Relative sequence tag density ratios were calculated against male control plasma samples and results were compared to those of previous karyotyping. RESULTS: All trisomy 18 fetuses were identified correctly with a diagnostic sensitivity and specificity of 100%. However, low diagnostic sensitivity and specificity were observed for fetal trisomy 13 detection. CONCLUSIONS: We successfully applied single molecule sequencing in combination with relative sequence tag density calculations for noninvasive trisomy 18 detection using free DNA from maternal plasma. However, noninvasive trisomy 13 detection was not accurate and seemed to be influenced by more than just GC content.

Fine-tiling array CGH to improve diagnostics for α- and ß-thalassemia rearrangements.

Implementation of multiplex ligation-dependent probe amplification (MLPA) for thalassemia causing deletions has lead to the detection of new rearrangements. Knowledge of the exact breakpoint sequences should give more insight into the molecular mechanisms underlying these rearrangements, and would facilitate the design of gap-PCRs. We have designed a custom fine-tiling array with oligonucleotides covering the complete globin gene clusters. We hybridized 27 DNA samples containing newly identified deletions and nine positive controls. We designed specific primers to amplify relatively short fragments containing the breakpoint sequence and analyzed these by direct sequencing. Results from nine positive controls showed that array comparative genomic hybridization (aCGH) is suitable to detect small and large rearrangements. We were able to locate all breakpoints to a region of approximately 2 kb. We designed breakpoint primers for 22 cases and amplification was successful in 19 cases. For 12 of these, the exact locations of the breakpoints were determined. Seven of these deletions have not been reported before. aCGH is a valuable tool for high-resolution breakpoint characterization. The combination of MLPA and aCGH has lead to relatively cheap and easy to perform PCR assays, which might be of use for laboratories as an alternative for MLPA in populations where only a limited number of specific deletions occur with high frequency.

Single molecule sequencing of free DNA from maternal plasma for noninvasive trisomy 21 detection.

BACKGROUND: Noninvasive fetal aneuploidy detection by use of free DNA from maternal plasma has recently been shown to be achievable by whole genome shotgun sequencing. The high-throughput next-generation sequencing platforms previously tested use a PCR step during sample preparation, which results in amplification bias in GC-rich areas of the human genome. To eliminate this bias, and thereby experimental noise, we have used single molecule sequencing as an alternative method. METHODS: For noninvasive trisomy 21 detection, we performed single molecule sequencing on the Helicos platform using free DNA isolated from maternal plasma from 9 weeks of gestation onwards. Relative sequence tag density ratios were calculated and results were directly compared to the previously described Illumina GAII platform. RESULTS: Sequence data generated without an amplification step show no GC bias. Therefore, with the use of single molecule sequencing all trisomy 21 fetuses could be distinguished more clearly from euploid fetuses. CONCLUSIONS: This study shows for the first time that single molecule sequencing is an attractive and easy to use alternative for reliable noninvasive fetal aneuploidy detection in diagnostics. With this approach, previously described experimental noise associated with PCR amplification, such as GC bias, can be overcome.

Non-invasive prenatal diagnosis of beta-thalassemia and sickle-cell disease using pyrophosphorolysis-activated polymerization and melting curve analysis.

OBJECTIVE: The aim of this study was to develop a pyrophosphorolysis-activated polymerization (PAP) assay for non-invasive prenatal diagnosis (NIPD) of ß-thalassemia major and sickle-cell disease (SCD). PAP is able to detect mutations in free fetal DNA in a highly contaminating environment of maternal plasma DNA. METHODS: Pyrophosphorolysis-activated polymerization primers were designed for 12 informative SNPs, genotyped by melting curve analysis (MCA) in both parents. The PAP assay was tested in a series of 13 plasma DNA samples collected from pregnant women. A retrospective NIPD was performed in a couple at risk for SCD. RESULTS: All PAP reactions were optimized and able to detect <3% target gDNA in a background of >97% wildtype gDNA. In all 13 cases, the paternal allele was detected by PAP in maternal plasma at 10 to 18 weeks of gestation. For the couple at risk, PAP showed presence of the normal paternal SNP allele in maternal plasma, which was confirmed by results of the chorionic villus sampling analysis. CONCLUSIONS: In contrast to other methods used for NIPD, the combined PAP and MCA analysis detecting the normal paternal allele is also applicable for couples at risk carrying the same mutation, provided that a previously born child is available for testing to determine the linkage to the paternal SNPs.

No haploinsufficiency but loss of heterozygosity for EXT in multiple osteochondromas.

Multiple osteochondromas (MO) is an autosomal dominant disorder caused by germline mutations in EXT1 and/or EXT2. In contrast, solitary osteochondroma (SO) is nonhereditary. Products of the EXT gene are involved in heparan sulfate (HS) biosynthesis. In this study, we investigated whether osteochondromas arise via either loss of heterozygosity (2 hits) or haploinsufficiency. An in vitro three-dimensional chondrogenic pellet model was used to compare heterozygous bone marrow-derived mesenchymal stem cells (MSCs EXT(wt/-)) of MO patients with normal MSCs and the corresponding tumor specimens (presumed EXT(-/-)). We demonstrated a second hit in EXT in five of eight osteochondromas. HS chain length and structure, in vitro chondrogenesis, and EXT expression levels were identical in both EXT(wt/-) and normal MSCs. Immunohistochemistry for HS, HS proteoglycans, and HS-dependent signaling pathways (eg, TGF-ß/BMP, Wnt, and PTHLH) also showed no differences. The cartilaginous cap of osteochondroma contained a mixture of HS-positive and HS-negative cells. Because a heterozygous EXT mutation does not affect chondrogenesis, EXT, HS, or downstream signaling pathways in MSCs, our results refute the haploinsufficiency theory. We found a second hit in 63% of analyzed osteochondromas, supporting the hypothesis that osteochondromas arise via loss of heterozygosity. The detection of the second hit may depend on the ratio of HS-positive (normal) versus HS-negative (mutated) cells in the cartilaginous cap of the osteochondroma.

Keratosis Follicularis Spinulosa Decalvans is caused by mutations in MBTPS2.

Keratosis Follicularis Spinulosa Decalvans (KFSD) is a rare genetic disorder characterized by development of hyperkeratotic follicular papules on the scalp followed by progressive alopecia of the scalp, eyelashes, and eyebrows. Associated eye findings include photophobia in childhood and corneal dystrophy. Due to the genetic and clinical heterogeneity of similar disorders, a definitive diagnosis of KFSD is often challenging. Toward identification of the causative gene we reanalyzed a large Dutch KFSD family. SNP arrays (1 M) redefined the locus to a 2.9-Mb region at Xp22.12-Xp22.11. Screening of all 14 genes in the candidate region identified MBTPS2 as the candidate gene carrying a c.1523A>G (p.Asn508Ser) missense mutation. The variant was also identified in two unrelated X-linked KFSD families and cosegregated with KFSD in all families. In symptomatic female carriers, skewed X-inactivation of the normal allele matched with increased severity of symptoms. MBTPS2 is required for cleavage of sterol regulatory element-binding proteins (SREBPs). In vitro functional expression studies of the c.1523A>G mutation showed that sterol responsiveness was reduced by half. Other missense mutations in MBTPS2 have recently been identified in patients with IFAP syndrome. We postulate that both phenotypes are in the spectrum of one genetic disorder with a partially overlapping phenotype.

Pre- and postsynaptic neuromuscular junction abnormalities in musk myasthenia.

Autoantibodies to muscle-specific kinase (MuSK) can cause myasthenia gravis (MG). The pathophysiological mechanism remains unknown. We report in vitro electrophysiological and histological studies of the neuromuscular junction in a MuSK MG patient. Low levels of presynaptic acetylcholine release and small miniature endplate potentials were found. This combination of pre- and postsynaptic abnormalities was supported by histology, revealing partially denervated postsynaptic areas, and some degeneration of postsynaptic folds. Results suggest that anti-MuSK antibodies reduce the stability of muscle-nerve contact.

Diagnostic guidelines for high-resolution melting curve (HRM) analysis: an interlaboratory validation of BRCA1 mutation scanning using the 96-well LightScanner.

Genetic analysis of BRCA1 by sequencing is often preceded by a scanning method like denaturing gradient gel electrophoresis (DGGE), protein truncation test (PTT) or DHPLC. High-resolution melting curve (HRM) analysis is a promising and economical method for high-throughput mutation scanning. The EuroGentest network (www.eurogentest.org) aims to assist with the introduction of novel technologies in the diagnostic setting. Therefore, we have performed a thorough and high-standard interlaboratory evaluation and validation of HRM, in collaboration with Idaho Technology, the manufacturer of the LightScanner (LS). Through this detailed study of 170 variants, we have generated guidelines for easy setup and implementation of HRM as a scanning technique for new genes, which are adaptable to the quality system of an individual diagnostic laboratory. This validation study includes the description of a BRCA1-specific mutation screening test using the 96-well LS. This assay comprises 40 amplicons and was evaluated using a statistically significant elaborate panel of variants and control DNA samples. All heterozygous variants were detected. Moreover, genotype analysis for nine common polymorphisms created a fast screening and detection method for these frequently occurring nonpathogenic variants. A blind study using a total of 28 patient-derived DNA samples resulted also in 100% detection and showed an average specificity of 98%, indicating a low incidence of false positives (FPs).

Intronic variants in BRCA1 and BRCA2 that affect RNA splicing can be reliably selected by splice-site prediction programs.

A large number of sequence variants identified in BRCA1 and BRCA2 cannot be distinguished as either disease-causing mutations or neutral variants. These so-called unclassified variants (UVs) include variants that are located in the intronic sequences of BRCA1 and BRCA2. The purpose of this study was to assess the use of splice-site prediction programs (SSPPs) to select intronic variants in BRCA1 and BRCA2 that are likely to affect RNA splicing. We performed in vitro molecular characterization of RNA of six intronic variants in BRCA1 and BRCA2. In four cases (BRCA1, c.81-6T>A and c.4986+5G>T; BRCA2, c.7617+2T>G and c.8754+5G>A) a deleterious effect on RNA splicing was seen, whereas the c.135-15_-12del variant in BRCA1 showed no effect on RNA splicing. In the case of the BRCA2 c.68-7T>A variant, RNA analysis was not sufficient to establish the clinical significance. Six SSPPs were used to predict whether an effect on RNA splicing was expected for these six variants as well as for 23 intronic variants in BRCA1 for which the effect on RNA splicing has been published. Out of a total of 174 predictions, 161 (93%) were informative (i.e., the wild-type splice-site was recognized). No false-negative predictions were observed; an effect on RNA splicing was always predicted by these programs. In four cases (2.5%) a false-positive prediction was observed. For DNA diagnostic laboratories, these programs are therefore very useful to select intronic variants that are likely to affect RNA splicing for further analysis.

Multiple genomic aberrations in a patient with mental retardation and hypogonadism: 45,X/46,X,psu dic(Y) karyotype, thyroid hormone receptor beta (THRB) mutation and heterozygosity for Wilson disease.

We report on multiple genomic aberrations in a patient with mental retardation. In addition, he had hypogonadism, elevated thyroid hormone levels, hearing loss, delayed speech development and mild dysmorphic features. First, we identified a mosaic karyotype, 45,X/46,X,psu dic(Y). The pseudo-dicentric Y chromosome has three short arm segments. Second, we found a germline mutation (Pro453Thr) of the thyroid hormone receptor beta (THRB) which is associated with resistance to thyroid hormone. Third, he was found to be a carrier of a heterozygous ATP7B mutation (c.2575 + 5G > C), the Wilson disease gene. Even though an array-CGH (with a density of approximately 1 Mb) did not reveal any further genomic gains or losses, we cannot exclude that all contributing factors have been identified. However, this case report shows that with increasing technological possibilities we can find more than one cause for developmental problems in a single patient. The identification of multiple causes in a single patient may complicate explaining the disorder and genetic counseling.

Search for copy number alterations in the MEFV gene using multiplex ligation probe amplification, experience from three diagnostic centres.

Familial mediterranean fever (FMF) is a hereditary autoinflammatory autosomal recessive disease caused by mutations in the MEFV gene. Despite the identification of many disease associated MEFV mutations, often the clinical diagnosis cannot be genetically confirmed. The currently used diagnostic sequencing techniques only allow the detection of point mutations, small deletions or duplications. The question as to whether larger genetic alterations are also involved in the pathophysiology of FMF remains to be answered. To address this question, we used multiplex ligation-dependent probe amplification (MLPA) on a total of 216 patients with FMF symptoms. This careful analysis revealed that not a single deletion/duplication could be detected in this large cohort of patients. This result suggests that single or multiexon MEFV gene copy number changes do not contribute substantially, if at all, to the MEFV mutation spectrum.

Identification of copy number variants associated with BPES-like phenotypes.

Blepharophimosis-Ptosis-Epicanthus inversus syndrome (BPES) is a well-characterized rare syndrome that includes an eyelid malformation associated with (type I) or without premature ovarian failure (type II). Patients with typical BPES have four major characteristics: blepharophimosis, ptosis, epicanthus inversus and telecanthus. Mutations in the FOXL2 gene, encoding a forkhead transcription factor, are responsible for the majority of both types of BPES. However, many patients with BPES-like features, i.e., having at least two major characteristics of BPES, have an unidentified cause. Here, we report on a group of 27 patients with BPES-like features, but without an identified genetic defect in the FOXL2 gene or flanking region. These patients were analyzed with whole-genome high-density arrays in order to identify copy number variants (CNVs) that might explain the BPES-like phenotype. In nine out of 27 patients (33%) CNVs not previously described as polymorphisms were detected. Four of these patients displayed psychomotor retardation as an additional clinical characteristic. In conclusion, we demonstrate that BPES-like phenotypes are frequently caused by CNVs, and we emphasize the importance of whole-genome copy number screening to identify the underlying genetic causes of these phenotypes.

A novel (Leu183Pro-)mutation in the HFE-gene co-inherited with the Cys282Tyr mutation in two unrelated Dutch hemochromatosis patients.

We describe a novel heterozygous mutation in exon 3 of the HFE-gene that was co-inherited with Cys282Tyr in two unrelated Dutch men both presenting a classical form of hereditary hemochromatosis. Heterozygosity for this mutation was also found in one out of 100 healthy controls of Dutch descent. This c.548T>C mutation converts a leucine to a proline residue at position 183 in the alpha2-helix of the HFE-protein (Leu183Pro). Standard bioinformatics analysis shows that the mutation is likely to disturb the HFE interaction with TfR1. This disrupting role of the mutation in the iron regulatory pathway is further corroborated by the familial co-occurrence of the observed compound heterozygosity with increased serum iron parameters. Haplotype analysis strongly suggests that this novel mutation arose from a common ancestor in the distant past. These findings may have implications for HFE-testing of iron overloaded heterozygous Cys282Tyr-patients of Northern European origin and their relatives.