Recurrent neural network for predicting absence of heterozygosity from low pass WGS with ultra-low depth.

BACKGROUND: The absence of heterozygosity (AOH) is a kind of genomic change characterized by a long contiguous region of homozygous alleles in a chromosome, which may cause human genetic disorders. However, no method of low-pass whole genome sequencing (LP-WGS) has been reported for the detection of AOH in a low-pass setting of less than onefold. We developed a method, termed CNVseq-AOH, for predicting the absence of heterozygosity using LP-WGS with ultra-low sequencing data, which overcomes the sparse nature of typical LP-WGS data by combing population-based haplotype information, adjustable sliding windows, and recurrent neural network (RNN). We tested the feasibility of CNVseq-AOH for the detection of AOH in 409 cases (11 AOH regions for model training and 863 AOH regions for validation) from the 1000 Genomes Project (1KGP). AOH detection using CNVseq-AOH was also performed on 6 clinical cases with previously ascertained AOHs by whole exome sequencing (WES). RESULTS: Using SNP-based microarray results as reference (AOHs detected by CNVseq-AOH with at least a 50% overlap with the AOHs detected by chromosomal microarray analysis), 409 samples (863 AOH regions) in the 1KGP were used for concordant analysis. For 784 AOHs on autosomes and 79 AOHs on the X chromosome, CNVseq-AOH can predict AOHs with a concordant rate of 96.23% and 59.49% respectively based on the analysis of 0.1-fold LP-WGS data, which is far lower than the current standard in the field. Using 0.1-fold LP-WGS data, CNVseq-AOH revealed 5 additional AOHs (larger than 10 Mb in size) in the 409 samples. We further analyzed AOHs larger than 10 Mb, which is recommended for reporting the possibility of UPD. For the 291 AOH regions larger than 10 Mb, CNVseq-AOH can predict AOHs with a concordant rate of 99.66% with only 0.1-fold LP-WGS data. In the 6 clinical cases, CNVseq-AOH revealed all 15 known AOH regions. CONCLUSIONS: Here we reported a method for analyzing LP-WGS data to accurately identify regions of AOH, which possesses great potential to improve genetic testing of AOH.

Endoscopic ultrasound avoids adverse events in high probability choledocholithiasis patients with a negative computed tomography.

BACKGROUND: The current guideline recommends patients who meet high probability criteria for choledocholithiasis to receive endoscopic retrograde cholangiopancreatography (ERCP). However, adverse events can occur during ERCP. Our goal is to determine whether endoscopic ultrasound (EUS) before ERCP can avoid unnecessary ERCP complications, especially in patients with a negative CT scan. METHODS: A total of 604 patients with high probability of choledocholithiasis were screened and 104 patients were prospectively enrolled. Patients with malignant biliary obstruction, altered GI anatomy, and choledocholithiasis on CT scan were excluded. Among them, 44 patients received EUS first, and ERCP if choledocholithiasis present (EUS-first group). The other 60 patients received ERCP directly (ERCP-first group). The baseline characteristics, presence of choledocholithiasis, and complications were compared between groups. All patients were followed for 3 months to determine the difference in recurrent biliary event rate. Cost-effectiveness was compared between the two strategies. RESULTS: There was no marked difference in age, sex, laboratory data, presenting with pancreatitis, and risk factors for choledocholithiasis. Overall, 51 patients (49.0%) had choledocholithiasis, which did not justify the risk of direct ERCP. In the EUS-first group, 27 (61.4%) ERCP procedures were prevented. The overall complication rate was significantly lower in the EUS-first group compared to the ERCP-fist group (6.8% vs. 21.7%, P = 0.04). The number-needed-to-treat to avoid one unnecessary adverse event was 6.71. After a 3-month follow-up, the cumulative recurrence biliary event rates were similar (13.6% vs. 15.0%, P = 0.803). EUS-first strategy was more cost-effective than the ERCP-first strategy (mean cost 2322.89$ vs. 3175.63$, P = 0.002). CONCLUSIONS: In high-probability choledocholithiasis patients with a negative CT, the EUS-first strategy is cost-effective, which can prevent unnecessary ERCP procedures and their complications.

Noninvasive prenatal diagnosis of cobalamin C (cblC) deficiency through target region sequencing of cell-free DNA in maternal plasma.

OBJECTIVE: This study aimed to validate the feasibility of haplotype-based noninvasive prenatal diagnosis (NIPD) of cobalamin C (cblC) deficiency. METHOD: This method includes three steps: First, targeted sequencing was performed on 21 families affected by cblC deficiency (including the couples and probands). Second, parental haplotypes linked with the pathogenic variant were determined using the genotypes of trios. Then, the fetal haplotypes were inferred through a parental haplotype assisted hidden Markov model (HMM). The NIPD results were confirmed by using the invasive procedures. RESULTS: Twenty-one fetal genotypes were successfully inferred by NIPD including three compound heterozygotes with cblC deficiency, nine heterozygote carriers of cblC deficiency, and nine normal fetuses. The NIPD results were confirmed using the invasive procedures with 100% concordant rate. CONCLUSION: This result has shown that haplotype-based NIPD of cblC deficiency has high concordant rate and indicated potential clinical utility as a pregnancy diagnosis method for high-risk carrier couples.

Noninvasive prenatal diagnosis for Duchenne muscular dystrophy based on the direct haplotype phasing.

OBJECTIVE: We aimed to investigate the validity of noninvasive prenatal diagnosis (NIPD) based on direct haplotype phasing without the proband or other family members and its feasibility for clinical application in the case of Duchenne muscular dystrophy (DMD). METHODS: Thirteen singleton-pregnancy families affected by DMD were recruited. The pathogenic variants in the pregnant females have been identified by multiplex ligation-dependent probe amplification (MLPA). We resolved maternal haplotypes for each family by performing targeted linked-read sequencing of their high molecular weight DNA, respectively. Then, we integrated the maternal haplotypes and the targeted sequencing results of maternal plasma DNA to infer the fetal haplotype and the DMD gene variant status. The fetal genotypes were further validated by using chorionic villus sampling. RESULTS: The method of directly resolving maternal haplotype through targeted linked-read sequencing was smoothly performed in 12 participated families, but one failed (F11). The predicted variant status of 12 fetuses was correct, which had been confirmed by invasive prenatal diagnosis. CONCLUSION: Direct haplotyping of NIPD based on linked-read sequencing for DMD is accurate.

Accuracy and depth evaluation of clinical low pass genome sequencing in the detection of mosaic aneuploidies and CNVs.

BACKGROUND: Low-pass genome sequencing (LP GS) has shown distinct advantages over traditional methods for the detection of mosaicism. However, no study has systematically evaluated the accuracy of LP GS in the detection of mosaic aneuploidies and copy number variants (CNVs) in prenatal diagnosis. Moreover, the influence of sequencing depth on mosaicism detection of LP GS has not been fully evaluated. METHODS: To evaluate the accuracy of LP GS in the detection of mosaic aneuploidies and mosaic CNVs, 27 samples with known aneuploidies and CNVs and 1 negative female sample were used to generate 6 simulated samples and 21 virtual samples, each sample contained 9 different mosaic levels. Mosaic levels were simulated by pooling reads or DNA from each positive sample and the negative sample according to a series of percentages (ranging from 3 to 40%). Then, the influence of sequencing depth on LP GS in the detection of mosaic aneuploidies and CNVs was evaluated by downsampling. RESULTS: To evaluate the accuracy of LP GS in the detection of mosaic aneuploidies and CNVs, a comparative analysis of mosaic levels was performed using 6 simulated samples and 21 virtual samples with 35 M million (M) uniquely aligned high-quality reads (UAHRs). For mosaic levels > 30%, the average difference (detected mosaic levels vs. theoretical mosaic levels) of 6 mosaic CNVs in simulated samples was 4.0%, and the average difference (detected mosaic levels vs. mosaic levels of Y chromosome) of 6 mosaic aneuploidies and 15 mosaic CNVs in virtual samples was 2.7%. Furthermore, LP GS had a higher detection rate and accuracy for the detection of mosaic aneuploidies and CNVs of larger sizes, especially mosaic aneuploidies. For depth evaluation, the results of LP GS in downsampling samples were compared with those of LP GS using 35 M UAHRs. The detection sensitivity of LP GS for 6 mosaic aneuploidies and 15 mosaic CNVs in virtual samples increased with UAHR. For mosaic levels > 30%, the total detection sensitivity reached a plateau at 30 M UAHRs. With 30 M UAHRs, the total detection sensitivity was 99.2% for virtual samples. CONCLUSIONS: We demonstrated the accuracy of LP GS in mosaicism detection using simulated data and virtual samples, respectively. Thirty M UAHRs (single-end 35 bp) were optimal for LP GS in the detection of mosaic aneuploidies and most mosaic CNVs larger than 1.48 Mb (Megabases) with mosaic levels > 30%. These results could provide a reference for laboratories that perform clinical LP GS in the detection of mosaic aneuploidies and CNVs.

Test development, optimization and validation of a WGS pipeline for genetic disorders.

BACKGROUND: With advances in massive parallel sequencing (MPS) technology, whole-genome sequencing (WGS) has gradually evolved into the first-tier diagnostic test for genetic disorders. However, deployment practice and pipeline testing for clinical WGS are lacking. METHODS: In this study, we introduced a whole WGS pipeline for genetic disorders, which included the entire process from obtaining a sample to clinical reporting. All samples that underwent WGS were constructed using polymerase chain reaction (PCR)-free library preparation protocols and sequenced on the MGISEQ-2000 platform. Bioinformatics pipelines were developed for the simultaneous detection of various types of variants, including single nucleotide variants (SNVs), insertions and deletions (indels), copy number variants (CNVs) and balanced rearrangements, mitochondrial (MT) variants, and other complex variants such as repeat expansion, pseudogenes and absence of heterozygosity (AOH). A semiautomatic pipeline was developed for the interpretation of potential SNVs and CNVs. Forty-five samples (including 14 positive commercially available samples, 23 laboratory-held positive cell lines and 8 clinical cases) with known variants were used to validate the whole pipeline. RESULTS: In this study, a whole WGS pipeline for genetic disorders was developed and optimized. Forty-five samples with known variants (6 with SNVs and Indels, 3 with MT variants, 5 with aneuploidies, 1 with triploidy, 23 with CNVs, 5 with balanced rearrangements, 2 with repeat expansions, 1 with AOHs, and 1 with exon 7-8 deletion of SMN1 gene) validated the effectiveness of our pipeline. CONCLUSIONS: This study has been piloted in test development, optimization, and validation of the WGS pipeline for genetic disorders. A set of best practices were recommended using our pipeline, along with a dataset of positive samples for benchmarking.

Next generation sequencing is a highly reliable method to analyze exon 7 deletion of survival motor neuron 1 (SMN1) gene.

Spinal muscular atrophy (SMA) is one of the most common and severe genetic diseases. SMA carrier screening is an effective way to identify couples at risk of having affected children. Next-generation sequencing (NGS)-based expanded carrier screening could detect SMN1 gene copy number without extra experiment and with high cost performance. However, its performance has not been fully evaluated. Here we conducted a systematic comparative study to evaluate the performance of three common methods. 478 samples were analyzed with multiplex ligation probe amplification (MLPA), real-time quantitative polymerase chain reaction (qPCR) and NGS, simultaneously. Taking MLPA-based results as the reference, for 0 copy, 1 copy and ≥ 2 copy SMN1 analysis with NGS, the sensitivity, specificity and precision were all 100%. Using qPCR method, the sensitivity was 100%, 97.52% and 94.30%, respectively; 98.63%, 95.48% and 100% for specificity; and 72.72%, 88.72% and 100% for precision. NGS repeatability was higher than that of qPCR. Moreover, among three methods, NGS had the lowest retest rate. Thus, NGS is a relatively more reliable method for SMN1 gene copy number detection. In expanded carrier screening, compared with the combination of multiple methods, NGS method could reduce the test cost and simplify the screening process.

Performance characterization of PCR-free whole genome sequencing for clinical diagnosis.

ABSTRACT: To evaluate the performance of polymerase chain reaction (PCR)-free whole genome sequencing (WGS) for clinical diagnosis, and thereby revealing how experimental parameters affect variant detection.Five NA12878 samples were sequenced using MGISEQ-2000. NA12878 samples underwent WGS with differing deoxyribonucleic acid (DNA) input and library preparation protocol (PCR-based vs PCR-free protocols for library preparation). The depth of coverage and genotype quality of each sample were compared. The performance of each sample was measured for sensitivity, coverage of depth and breadth of coverage of disease-related genes, and copy number variants. We also developed a systematic WGS pipeline (PCR-free) for the analysis of 11 clinical cases.In general, NA12878-2 (PCR-free WGS) showed better depth of coverage and genotype quality distribution than NA12878-1 (PCR-based WGS). With a mean depth of ∼40×, the sensitivity of homozygous and heterozygous single nucleotide polymorphisms (SNPs) of NA12878-2 showed higher sensitivity (>99.77% and >99.82%) than NA12878-1, and positive predictive value exceeded 99.98% and 99.07%. The sensitivity and positive predictive value of homozygous and heterozygous indels for NA12878-2 (PCR-free WGS) showed great improvement than NA128878-1. The breadths of coverage for disease-related genes and copy number variants are slightly better for samples with PCR-free library preparation protocol than the sample with PCR-based library preparation protocol. DNA input also influences the performance of variant detection in samples with PCR-free WGS. All the 19 previously confirmed variants in 11 clinical cases were successfully detected by our WGS pipeline (PCR free).Different experimental parameters may affect variant detection for clinical WGS. Clinical scientists should know the range of sensitivity of variants for different methods of WGS, which would be useful when interpreting and delivering clinical reports.

Noninvasive prenatal diagnosis of monogenic disorders based on direct haplotype phasing through targeted linked-read sequencing.

BACKGROUND: Though massively parallel sequencing has been widely applied to noninvasive prenatal screen for common trisomy, the clinical use of massively parallel sequencing to noninvasive prenatal diagnose monogenic disorders is limited. This study was to develop a method for directly determining paternal haplotypes for noninvasive prenatal diagnosis of monogenic disorders without requiring proband's samples. METHODS: The study recruited 40 families at high risk for autosomal recessive diseases. The targeted linked-read sequencing was performed on high molecular weight (HMW) DNA of parents using customized probes designed to capture targeted genes and single-nucleotide polymorphisms (SNPs) distributed within 1Mb flanking region of targeted genes. Plasma DNA from pregnant mothers also underwent targeted sequencing using the same probes to determine fetal haplotypes according to parental haplotypes. The results were further confirmed by invasive prenatal diagnosis. RESULTS: Seventy-eight parental haplotypes of targeted gene were successfully determined by targeted linked-read sequencing. The predicted fetal inheritance of variant was correctly deduced in 38 families in which the variants had been confirmed by invasive prenatal diagnosis. Two families were determined to be no-call. CONCLUSIONS: Targeted linked-read sequencing method demonstrated to be an effective means to phase personal haplotype for noninvasive prenatal diagnosis of monogenic disorders.

Characterizing sensitivity and coverage of clinical WGS as a diagnostic test for genetic disorders.

BACKGROUND: Due to its reduced cost and incomparable advantages, WGS is likely to lead to changes in clinical diagnosis of rare and undiagnosed diseases. However, the sensitivity and breadth of coverage of clinical WGS as a diagnostic test for genetic disorders has not been fully evaluated. METHODS: Here, the performance of WGS in NA12878, the YH cell line, and the Chinese trios were measured by assessing their sensitivity, PPV, depth and breadth of coverage using MGISEQ-2000. We also compared the performance of WES and WGS using NA12878. The sensitivity and PPV were tested using the family-based trio design for the Chinese trios. We further developed a systematic WGS pipeline for the analysis of 8 clinical cases. RESULTS: In general, the sensitivity and PPV for SNV/indel detection increased with mean depth and reached a plateau at an ~ 40X mean depth using down-sampling samples of NA12878. With a mean depth of 40X, the sensitivity of homozygous and heterozygous SNPs of NA12878 was > 99.25% and > 99.50%, respectively, and the PPV was 99.97% and 98.96%. Homozygous and heterozygous indels showed lower sensitivity and PPV. The sensitivity and PPV were still not 100% even with a mean depth of ~ 150X. We also observed a substantial variation in the sensitivity of CNV detection across different tools, especially in CNVs with a size less than 1 kb. In general, the breadth of coverage for disease-associated genes and CNVs increased with mean depth. The sensitivity and coverage of WGS (~ 40X) was better than WES (~ 120X). Among the Chinese trios with an ~ 40X mean depth, the sensitivity among offspring was > 99.48% and > 96.36% for SNP and indel detection, and the PPVs were 99.86% and 97.93%. All 12 previously validated variants in the 8 clinical cases were successfully detected using our WGS pipeline. CONCLUSIONS: The current standard of a mean depth of 40X may be sufficient for SNV/indel detection and identification of most CNVs. It would be advisable for clinical scientists to determine the range of sensitivity and PPV for different classes of variants for a particular WGS pipeline, which would be useful when interpreting and delivering clinical reports.

Noninvasive prenatal testing of α-thalassemia and ß-thalassemia through population-based parental haplotyping.

BACKGROUND: Noninvasive prenatal testing (NIPT) of recessive monogenic diseases depends heavily on knowing the correct parental haplotypes. However, the currently used family-based haplotyping method requires pedigrees, and molecular haplotyping is highly challenging due to its high cost, long turnaround time, and complexity. Here, we proposed a new two-step approach, population-based haplotyping-NIPT (PBH-NIPT), using α-thalassemia and ß-thalassemia as prototypes. METHODS: First, we deduced parental haplotypes with Beagle 4.0 with training on a large retrospective carrier screening dataset (4356 thalassemia carrier screening-positive cases). Second, we inferred fetal haplotypes using a parental haplotype-assisted hidden Markov model (HMM) and the Viterbi algorithm. RESULTS: With this approach, we enrolled 59 couples at risk of having a fetus with thalassemia and successfully inferred 94.1% (111/118) of fetal alleles. We confirmed these alleles by invasive prenatal diagnosis, with 99.1% (110/111) accuracy (95% CI, 95.1-100%). CONCLUSIONS: These results demonstrate that PBH-NIPT is a sensitive, fast, and inexpensive strategy for NIPT of thalassemia.

NGS-based spinal muscular atrophy carrier screening of 10,585 diverse couples in China: a pan-ethnic study.

In this study, we performed a spinal muscular atrophy carrier screening investigation with NGS-based method. First, the validation for NGS-based method was implemented in 2255 samples using real-time PCR. The concordance between the NGS-based method and real-time PCR for the detection of SMA carrier and patient were up to 100%. Then, we applied this NGS-based method in 10,585 self-reported normal couples (34 Chinese ethnic groups from 5 provinces in South China) for SMA carrier screening. The overall carrier frequency was 1 in 73.8 (1.4%). It varied substantially between ethnic groups, highest in Dai ethnicity (4.3%), and no significant difference was found between five provinces. One couple was detected as carriers with an elevated risk of having an SMA affected baby. The distribution of SMN1:SMN2 genotype was also revealed in this study. Among the individuals with normal phenotype, the exon 7 copy-number ratio of SMN1 to SMN2 proved the gene conversion between them. With NGS-based method, we investigated SMA carrier status in Chinese population for the first time, and our results demonstrated that it is a promising alternative for SMA carrier screening and could provide data support and reference for future clinical application.

Whole Genome Sequencing in the Evaluation of Fetal Structural Anomalies: A Parallel Test with Chromosomal Microarray Plus Whole Exome Sequencing.

Whole genome sequencing (WGS) is a powerful tool for postnatal genetic diagnosis, but relevant clinical studies in the field of prenatal diagnosis are limited. The present study aimed to prospectively evaluate the utility of WGS compared with chromosomal microarray (CMA) and whole exome sequencing (WES) in the prenatal diagnosis of fetal structural anomalies. We performed trio WGS (≈40-fold) in parallel with CMA in 111 fetuses with structural or growth anomalies, and sequentially performed WES when CMA was negative (CMA plus WES). In comparison, WGS not only detected all pathogenic genetic variants in 22 diagnosed cases identified by CMA plus WES, yielding a diagnostic rate of 19.8% (22/110), but also provided additional and clinically significant information, including a case of balanced translocations and a case of intrauterine infection, which might not be detectable by CMA or WES. WGS also required less DNA (100 ng) as input and could provide a rapid turnaround time (TAT, 18 ± 6 days) compared with that (31 ± 8 days) of the CMA plus WES. Our results showed that WGS provided more comprehensive and precise genetic information with a rapid TAT and less DNA required than CMA plus WES, which enables it as an alternative prenatal diagnosis test for fetal structural anomalies.

Haplotype-Based noninvasive prenatal diagnosis for duchenne muscular dystrophy: A pilot study in South China.

OBJECTIVE: To explore the accuracy and feasibility of noninvasive prenatal diagnosis (NIPD) for Duchenne Muscular Dystrophy (DMD) based on the haplotype approach. METHODS: We recruited singleton pregnancies at-risk of DMD at 12-25 weeks of gestation from 17 families who all had a proband child affected by DMD. We have identified the pathogenic mutations in probands and their mothers by multiplex ligation-dependent probe amplification (MLPA). To construct parental haplotypes, we performed captured sequencing on genomic DNA from parents and probands. The integration analysis of parental haplotypes and targeted sequencing results of maternal plasma DNA were used to infer the fetal haplotype and genotypes in the DMD gene. Fetal DMD genotypes were further confirmed by invasive prenatal diagnosis. RESULT: We have successfully performed the haplotype-based NIPD in all recruited families. Ten fetuses were identified as normal, including four female and six male fetuses. Four female fetuses were carriers, and the other three male fetuses were affected by DMD with exons 49-52 deletion, exons 8-37 deletion and c.628 G > T mutation, respectively. The results of NIPD were consistent with those of invasive diagnosis. CONCLUSION: Haplotype-based NIPD for DMD by targeted sequencing is promising and has the potential for clinical application.

Pilot study of expanded carrier screening for 11 recessive diseases in China: results from 10,476 ethnically diverse couples.

Expanded carrier screening (ECS) has been demonstrated to increase the detection rate of carriers compared with traditional tests. The aim of this study was to assess the potential value of ECS for clinical application in Southern China, a region with high prevalence of thalassemia and with diverse ethnic groups, and to provide a reference for future implementations in areas with similar population characteristics. A total of 10,476 prenatal/preconception couples from 34 self-reported ethnic groups were simultaneously tested and analyzed anonymously for 11 Mendelian disorders using targeted next-generation sequencing. Overall, 27.49% of individuals without self-reported family history of disorders were found to be carriers of at least 1 of the 11 conditions, and the carrier frequency varied greatly between ethnic groups, ranging from 4.15% to 81.35%. Furthermore, 255 couples (2.43%) were identified as carrier couples at an elevated risk having an affected baby, sixty-five of which would not have been identified through the existing screening strategy, which only detects thalassemia. The modeled risk of fetuses being affected by any of the selected disorders was 531 per 100,000 (95% CI, 497-567 per 100,000). Our data demonstrate the feasibility of ECS, and provide evidence that ECS is a promising alternative to traditional one-condition screening strategies. The lessons learned from this experience should be applicable for other countries or regions with diverse ethnic groups.

Haplotype-based Noninvasive Prenatal Diagnosis of Hyperphenylalaninemia through Targeted Sequencing of Maternal Plasma.

Here we developed a haplotype-based noninvasive prenatal diagnosis method for hyperphenylalaninemia (HPA) and demonstrated its accuracy and feasibility during early pregnancy. Capture sequencing was performed on genomic DNA from parents and probands using customized hybridization probes targeting highly heterozygous single-nucleotide polymorphisms located within the 1 M region flanking phenylalanine hydroxylase (PAH) and 6-pyruvoyltetrahydropterin (PTS) and its coding region to determine the parental haplotypes and linkage to pathogenic mutations. Maternal plasma DNA obtained at 12-20 weeks of gestation was also subjected to targeted sequencing to deduce the fetal haplotypes based on the parental haplotypes. The fetal genotypes were further validated by invasive prenatal diagnosis. Haplotype-based noninvasive prenatal testing was successfully performed in 13 families. Five fetuses were identified to harbor bi-allelic pathogenic variants of PAH, four fetuses were carriers of one heterozygous PAH variant, three fetuses were normal, and the fetus of the 6-pyruvoyl tetrahydrobiopterin synthase family was identified as normal. The fetal genotypes at two gestational weeks from the same PAH family were identical. All results were consistent with the prenatal diagnosis based on amniotic fluid. Haplotype-based noninvasive prenatal testing for HPA through targeted sequencing is accurate and feasible during early gestation.

A Pilot Study of Noninvasive Prenatal Diagnosis of Alpha- and Beta-Thalassemia with Target Capture Sequencing of Cell-Free Fetal DNA in Maternal Blood.

AIMS: Thalassemia is a dangerous hematolytic genetic disease. In south China, ∼24% Chinese carry alpha-thalassemia or beta-thalassemia gene mutations. Given the fact that the invasive sampling procedures can only be performed by professionals in experienced centers, it may increase the risk of miscarriage or infection. Thus, most people are worried about the invasive operation. As such, a noninvasive and accurate prenatal diagnosis is needed for appropriate genetic counseling for families with high risks. Here we sought to develop capture probes and their companion analysis methods for the noninvasive prenatal detection of deletional and nondeletional thalassemia. MATERIALS AND METHODS: Two families diagnosed as carriers of either beta-thalassemia gene or Southeast Asian deletional alpha-thalassemia gene mutation were recruited. The maternal plasma and amniotic fluid were collected for prenatal diagnosis. Probes targeting exons of the genes of interest and the highly heterozygous SNPs within the 1Mb flanking region were designed. The target capture sequencing was performed with plasma DNA from the pregnant woman and genomic DNA from the couples and their children. Then the parental haplotype was constructed by the trios-based strategy. The fetal haplotype was deduced from the parental haplotype with a hidden Markov model-based algorithm. RESULTS: The fetal genotypes were successfully deduced in both families noninvasively. The noninvasively constructed haplotypes of both fetuses were identical to the invasive prenatal diagnosis results with an accuracy rate of 100% in the target region. CONCLUSION: Our study demonstrates that the effective noninvasive prenatal diagnosis of alpha-thalassemia and beta-thalassemia can be achieved with the targeted capture sequencing and the haplotype-assisted analysis method.

Noninvasive prenatal diagnosis of 21-Hydroxylase deficiency using target capture sequencing of maternal plasma DNA.

Here, we aimed to validate a noninvasive method using capture sequencing for prenatal diagnosis of congenital adrenal hyperplasia due to 21-Hydroxylase deficiency (21-OHD). Noninvasive prenatal diagnosis (NIPD) of 21-OHD was based on 14 plasma samples collected from 12 families, including four plasma sample collected during the first trimester. Targeted capture sequencing was performed using genomic DNA from the parents and child trios to determine the pathogenic and wild-type alleles associated with the haplotypes. Maternal plasma DNA was also sequenced to determine the fetal inheritance of the allele using hidden Markov model-based haplotype linkage analysis. The effect of fetal DNA fraction and sequencing depth on the accuracy of NIPD was investigated. The lower limit of fetal DNA fraction was 2% and the threshold mean sequence depth was 38, suggesting potential advantage if used in early gestation. The CYP21A2 genotype of the fetus was accurately determined in all the 14 plasma samples as early as day 1 and 8 weeks of gestation. Results suggest the accuracy and feasibility of NIPD of 21-OHD using a small target capture region with a low threshold for fetal DNA fraction and sequence depth. Our method is cost-effective and suggests diagnostic applications in clinical practice.

Noninvasive prenatal diagnosis for X-linked disease by maternal plasma sequencing in a family of Hemophilia B.

OBJECTIVE: To apply a Hidden Markov Model to test Hemophilia B in a fetus by maternal plasma sequencing only employing proband and maternal haplotypes. CASE REPORT: A family at risk for Hemophilia B was recruited in this study. We performed genetic diagnosis on the proband using our targeted capture system (containing F9 gene coding region, highly heterozygous SNPs and a 13-kb chromosome Y specific region), and revealed a causative F9 gene mutation (c.190T>C, p.Cys64Arg). Maternal plasma cell-free DNA obtained at 8 weeks of gestation was targeted-captured and sequenced using the customized system. The fetus inherited the F9 (c.190T>C, p.Cys64Arg) mutation according to the Hidden Markov Model. The mother continued the pregnancy. CONCLUSIONS: This study is the first report of a haplotype-based approach in NIPD of Hemophilia B. With further evaluation, this method might be useful for NIPD of Hemophilia B and for other X-linked single-gene disorders.

Non-Invasive Prenatal Diagnosis of Lethal Skeletal Dysplasia by Targeted Capture Sequencing of Maternal Plasma.

BACKGROUND: Since the discovery of cell-free foetal DNA in the plasma of pregnant women, many non-invasive prenatal testing assays have been developed. In the area of skeletal dysplasia diagnosis, some PCR-based non-invasive prenatal testing assays have been developed to facilitate the ultrasound diagnosis of skeletal dysplasias that are caused by de novo mutations. However, skeletal dysplasias are a group of heterogeneous genetic diseases, the PCR-based method is hard to detect multiple gene or loci simultaneously, and the diagnosis rate is highly dependent on the accuracy of the ultrasound diagnosis. In this study, we investigated the feasibility of using targeted capture sequencing to detect foetal de novo pathogenic mutations responsible for skeletal dysplasia. METHODOLOGY/PRINCIPAL FINDINGS: Three families whose foetuses were affected by skeletal dysplasia and two control families whose foetuses were affected by other single gene diseases were included in this study. Sixteen genes related to some common lethal skeletal dysplasias were selected for analysis, and probes were designed to capture the coding regions of these genes. Targeted capture sequencing was performed on the maternal plasma DNA, the maternal genomic DNA, and the paternal genomic DNA. The de novo pathogenic variants in the plasma DNA data were identified using a bioinformatical process developed for low frequency mutation detection and a strict variant interpretation strategy. The causal variants could be specifically identified in the plasma, and the results were identical to those obtained by sequencing amniotic fluid samples. Furthermore, a mean of 97% foetal specific alleles, which are alleles that are not shared by maternal genomic DNA and amniotic fluid DNA, were identified successfully in plasma samples. CONCLUSIONS/SIGNIFICANCE: Our study shows that capture sequencing of maternal plasma DNA can be used to non-invasive detection of de novo pathogenic variants. This method has the potential to be used to facilitate the prenatal diagnosis of skeletal dysplasia.