Next Article in Journal
The Alteration of Intestinal Microbiota Profile and Immune Response in Epinephelus coioides during Pathogen Infection
Next Article in Special Issue
Resolving Differential Diagnostic Problems in von Willebrand Disease, in Fibrinogen Disorders, in Prekallikrein Deficiency and in Hereditary Hemorrhagic Telangiectasia by Next-Generation Sequencing
Previous Article in Journal
The Aging Bowel Dysfunction and Elderly Vulnerability towards COVID-19 Infection
Previous Article in Special Issue
The Effect of Synonymous Single-Nucleotide Polymorphisms on an Atypical Cystic Fibrosis Clinical Presentation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Rapid PCR-Free Next-Generation Sequencing Method for the Detection of Copy Number Variations in Prenatal Samples

1
Department of Obstetrics and Gynecology, Chinese Academy of Medical Sciences, Peking Union Medical College, Peking Union Medical College Hospital, Beijing 100730, China
2
Berry Genomics Corporation, Beijing 102200, China
*
Authors to whom correspondence should be addressed.
Equal contribution to the work.
Life 2021, 11(2), 98; https://doi.org/10.3390/life11020098
Submission received: 14 December 2020 / Revised: 15 January 2021 / Accepted: 25 January 2021 / Published: 28 January 2021

Abstract

:
Next-generation sequencing (NGS) is emerging as a new method for the detection of clinically significant copy number variants (CNVs). In this study, we developed and validated rapid CNV-sequencing (rCNV-seq) for clinical application in prenatal diagnosis. Low-pass whole-genome sequencing was performed on PCR libraries prepared from amniocyte genomic DNA. From 10–40 ng of input DNA, PCR-free libraries consistently produced sequencing data with high unique read mapping ratios, low read redundancy, low coefficient of variation for all chromosomes and high genomic coverage. In validation studies, reliable and accurate CNV detection using PCR-free-based rCNV-seq was demonstrated for a range of common trisomies and sex chromosome aneuploidies as well as microdeletion and duplication syndromes. In reproducibility studies, CNV copy number and genomic intervals closely matched those defined by chromosome microarray analysis. Clinical testing of genomic DNA samples from 217 women referred for prenatal diagnosis identified eight samples (3.7%) with known chromosome disorders. We conclude that PCR-free-based rCNV-seq is a sensitive, specific, reproducible and efficient method that can be used in any NGS-based diagnostic laboratory for detection of clinically significant CNVs.

1. Introduction

There are over 100 recurrent chromosome disorders that affect the human population [1]. These syndromes are caused by copy number variations (CNVs) that include both whole chromosome and segmental aneuploidies that arise during gametogenesis or in the preimplantation period of embryo development [2]. While the large majority of fetal CNVs lead to early miscarriage [3,4], a small proportion are developmentally competent and are compatible with a livebirth outcome. Approximately 15% of all congenital abnormalities are associated with a pathogenic CNV [5].
The most common whole chromosome aneuploidies of the newborn are trisomies T21 (Down’s syndrome), T18 (Edward’s syndrome) and T13 (Patau’s syndrome) and the sex chromosome aneuploidies 45,XO (Turner’s syndrome), 47,XXX (Triple X syndrome), 47,XXY (Klinefelter’s syndromes) and 47,XYY (Jacob’s syndrome) [6]. On the other hand, segmental aneuploidies are associated with a variety of microdeletion and microduplication syndromes (MMS). While the disease prevalence of different MMS varies from around 1 in 3000 to 1 in 25,000 [7], as a collective group, they are relatively common with 22q11.2 microdeletions (DiGeorge Syndrome) being the most frequent MMS seen in the newborn [8].
In prenatal diagnosis practice, women with a high fetal risk for a chromosome abnormality are usually referred for invasive testing by amniocentesis at 15–16 weeks gestation or by chorionic villous sampling at 10–11 weeks gestation. High risk is generally inferred by either an abnormal maternal serum screening score, advanced maternal age or presence of a soft ultrasound marker or an ultrasound structural abnormality [5,9]. Invasive chromosome testing for aneuploidy is usually routinely performed by karyotyping, which has a CNV resolution of around 5–10 Mb in size [5]. For identification of smaller CNVs, chromosome microarray analysis (CMA) using high density SNP arrays is the most widely used method [10,11]. For both karyotyping and CMA, amniocytes are normally cultured to generate sufficient cells for analysis, and results are generally available after 2 weeks. More recently, next-generation sequencing (NGS) has been developed as an alternative method for detection of CNVs, with a detection resolution down to 0.1 Mb [12,13,14], which is sufficient to detect the vast majority of MMS. In a recent large cohort study of patients referred for invasive testing, we demonstrated that CNV-seq applied to amniocyte DNA samples can provide an increased yield of pathogenic CNVs compared to karyotyping [15]. More recently, in a study of over 1000 women referred for invasive testing, low-pass whole-genome sequencing performed similarly to CMA for detection of clinically significant CNVs [16].
In the field of molecular diagnostics, CMA is still the preferred methodology for identification of clinically significant CNVs in prenatal samples. To advance the clinical application of NGS, improvements in the sensitivity, specificity, reproducibility and versatility of NGS-based CNV detection methods are still needed before this technology can be generally considered for routine chromosome testing. One essential part of the NGS workflow is the library preparation method, with most CNV-seq tests currently using a low DNA input (50–200 ng) and a PCR step to amplify genomic fragments. As an alternative approach, we developed and validated a PCR-free-based rapid CNV-seq (rCNV-seq) method suitable for analysis of low nanograms amounts of genomic DNA from uncultured amniocytes. In a prospective study of 217 patients referred for chromosome testing, we show that rCNV-seq of the amniotic cell genomic DNA can reliably and accurately detect clinically significant CNVs.

2. Materials and Methods

2.1. Patient Samples

All blood and amniocentesis samples were collected at the Prenatal Care Unit of Peking Union Medical College Hospital (PUMCH). The research study was approved by the Ethics Committee for Drug Clinical Trials in PUMCH (approval number KS2019136), and written informed consent was obtained from each patient. For validation studies, peripheral blood samples (2 mL) were taken from patients with known chromosome disorders and genomic DNA (gDNA) extracted using AxyPrep Mag Tissue-Blood gDNA kit (Axygen, Corning NY, USA). For patients with a suspected chromosome abnormality due to indications such as advanced maternal age, an abnormal maternal serum screening result or a structural abnormality revealed by ultrasound, invasive testing was performed on amniocentesis samples (10 mL of amniotic fluid) obtained at 15–16 weeks gestation. The fluid was centrifuged, and gDNA from the amniocyte cell pellet was purified using the AxyPrep Mag Tissue-Blood gDNA kit. Prior to CNV testing, amniocyte DNA was checked for maternal cell contamination (MCC) using an STR-based semi quantitative PCR assay [17]; MCC levels of <5% were considered acceptable for clinical testing.

2.2. Construction of NGS Libraries

An overview of the three library construction methods used in this study is shown in Figure 1. All genomic samples for library construction were quantified using Qubit 3.0 (Invitrogen, Waltham, MA, USA). For PCR-free-frag library construction used by rCNV-seq (Method 1), gDNA (10–40 ng) was initially treated by dsDNA fragmentase at 37 °C (NEBNext dsDNA Fragmentase, New England Biolabs, Ipswich, MA, USA) to produce smaller derivatives with an average size of ~200 bps. Prepared gDNA fragments were then end-repaired, A-tailed and then ligated with barcoded sequencing adaptors using a proprietary DNA repair kit (KR2000, Berry Genomics, Beijing, China) to generate libraries for sequencing. For PCR-free-soni (Method 2) and PCR-soni (Method 3) used for commercial CNV-seq, gDNA (1 µg, PCR-free-soni; 100–200 ng PCR-soni) was sheared by sonification and fragments of 350 bps size selected on agarose gels using TruSeq DNA PCR-Free Low Throughput Library Prep Kit (Illumina, San Diego, CA, USA) and TruSeq Nano DNA Low Throughput Library Prep Kit (Illumina, San Diego, CA, USA), respectively (Figure 1).

2.3. Copy Number Variation Sequencing and Data Analysis

Single end sequencing was performed on the NextSeq CN500 platform (Illumina, San Diego, CA, USA) with a run time of 6.5 h to generate approximately 5 million raw 45 bp reads per sample. Raw reads were then edited to remove artificial adaptor sequences, and the true 36 bp genome sequences were then mapped to the hg19 reference genome using the Burrows and Wheeler algorithm [18]. On average, approximately 2.8–3.2 million reads were uniquely mapped for data analysis. Reads were allocated to 20 kb bins along the length of each chromosome, and CNVs were identified from 24 chromosome copy number (CN) plots, as previously described [13,19]. Duplications were defined as CN >2.8, deletions CN <1.2, disomy (1.8 < CN < 2.2), mosaic trisomy (2.2 < CN < 2.8) and mosaic monosomy (1.2 < CN < 1.8).
Genomic variant databases including DGV (Database of Genomic Variants), http://projects.tcag.ca/variation), OMIM (Online Mendelian Inheritance in Man) (http://www.omim.org), PubMed (http://www.ncbi.nlm.nih.gov/pubmed) and UCSC (University of California, Santa Cruz) (http://genome.ucsc.edu/, hg19) were used as a reference source of CNVs. The pathogenicity of detected CNVs was assessed following ACMG guidelines [20].

2.4. Karyotyping

Cultured amniocytes were karyotyped by standard procedures [21]. Cytogenetic analysis of Giemsa-stained metaphase spreads was performed at a resolution of 320 bands.

2.5. Chromosome Microarray Analysis

CMA was performed using the CytoScan™ HD Array Kit (ThermoFisher Scientific, Waltham, MA, USA) according to the recommended protocol. The array contains more than 2.6 million SNPs and can detect copy number changes across the genome at a resolution of 25–50 kb. Genomic DNA samples were labeled and hybridized to the array according to the manufacturer’s recommended protocol. Fluorescence signals were scanned using the GeneChip scanner (ThermoFisher Scientific, Waltham, MA, USA) and chromosome copy number changes called by Applied Biosystems™ GeneChip Command Console Software (version 3.2.2).

3. Results

3.1. High Performance of PCR-Free-Based rCNV-Seq

The performance of PCR-free-frag-based libraries (rCNV-seq) and two control PCR-free-soni and PCR-soni libraries (CNV-seq) (Figure 1) was bioinformatically assessed for 17 replicate normal female gDNA samples extracted from postnatal blood samples (Figure 2). The percentage of uniquely mapped reads was significantly higher for PCR-free-frag (rCNV-seq) compared to PCR-free-soni (p < 0.0001) and PCR-soni (p < 0.0001) libraries. Both PCR-free-frag and PCR-free-soni libraries had a significantly lower read redundancy rate compared to the PCR-soni libraries (p < 0.0001). The coefficient of variation (CV) achieved with PCR-free-frag libraries was also lower across all chromosomes and for each individual chromosome, compared to PCR-free-soni and PCR-soni libraries. The median read sequencing depth, indicative of increased bin read numbers, was also higher for PCR-free-frag compared to PCR-free-soni and PCR-soni libraries. Further, the Guanidine/Cytosine (GC) content of the sequencing reads was higher for PCR-free-frag than PCR-free-soni. Lastly, the Q30 values were significantly higher for PCR-free-frag compared to PCR-free-soni (p < 0.0001) and PCR-soni (p < 0.0001). Based on evaluation of these key Quality Control (QC) sequencing indicators, rCNV-seq using PCR-free-frag libraries showed the highest performance values.
We further assessed the impact of input gDNA amount on library yield, unique mapping ratio and redundancy for PCR-free rCNV-seq (Figure 3). Mean DNA library concentrations were low (<200 p mole) at input DNA amounts <10 ng. Higher mean library yields suitable for sequencing were obtained with DNA input levels ranging from 50–800 ng. Despite different DNA inputs (1–800 ng), the unique mapping ratio was relatively stable at 60–65%. In terms of read redundancy, lower and higher DNA inputs were associated with a slightly higher ratio. Based on these assessment criteria, rCNV-seq provided high-quality sequencing data at input DNA amounts of ≥10 ng.

3.2. Validation of rCNV-Seq Using PCR-Free-Frag Library Preparations

Our newly designed rCNV-seq method based on PCR-free-frag libraries was further evaluated for the ability to detect known chromosome disorders, previously detected by CMA. These included single samples identified with trisomies T21, T18 and T13 and sex chromosome aneuploidies (SCAs) 45,XO, 47,XXX, 47,XXY and 47,XYY. In addition, we also selected MMS samples including three cases each of Cri du Chat and William–Beuren syndrome, two cases each of Wolf–Hirschhorn and Di George syndrome and one case each of Prader-Willi/Angelman, Smith–Magenis and Miller–Dieker syndrome. Further, we also included four additional samples with variants of uncertain significance (VOUS) involving CNVS < 1 Mb in size. In these experiments, we used a low DNA input of 40 ng as the starting template for PCR-free-frag library construction. Following rCNV-seq analysis of all samples, the specific CNV that was originally identified by CMA was also detected by rCNV-Seq (Table 1, Figure 4). In addition, the CNV intervals defined closely matched those defined by CMA and the expected CN prediction of one for deletions and three for duplications was observed. Taken together, our rCNV-seq protocol was highly sensitive for detecting the underlying causative CNVs associated with these MMS.
In detailed reproducibility experiments using 10 replicate samples of 10ng for library construction, we further evaluated the reliability of rCNV-seq to correctly detect the precise CNV interval, benchmarking it against the CNV interval previously defined by CMA (Figure 5). For 45,XO, one copy of chromosome X was clearly deleted (CN = 1) in all replicates. Likewise, for T21, the whole q arm of chromosome 21 was clearly duplicated (CN = 3) in all replicates. For the 5p15.33p15.1 (17.3 Mb), 17p13.3p13.2 (6.0 Mb), 17p12p11.2 (4.9 Mb), 15q11.2q13 (4.8 Mb), 22q11.21 (2.6 Mb), 7q11.23 (1.36 Mb) and 8p23.2 (0.83 Mb) deletions and the 11q14.3 duplication (0.72 Mb) CNVs, all were detected in the 10 replicates and the mapped interval closely matched the coordinates of the CMA-defined interval (within 93–100%).

3.3. Diagnostic Performance of rCNV-Seq Using Prenatal Samples

To assess the diagnostic performance of rCNV-seq, 40 ng of amniocyte gDNA was analyzed from 217 pregnant women referred for chromosome testing for a variety of different clinical indications (Table 2). MCC checks confirmed that all amniocyte genomic DNA samples were >99% fetal DNA. In eight samples, rCNV-seq detected whole chromosome fetal aneuploidies, including T21 (n = 6), T18 (n = 1) and 47,XXX (n = 1). These aneuploidies were confirmed by follow up karyotyping (Figure 6). There were no samples carrying pathogenic fetal CNVs associated with an MMS. However, eight samples were identified with small non-pathogenic duplications (0.96–1.81 Mb) that were classified as VOUS (Figure 7).

4. Discussion

In this study, we developed rCNV-seq based on a PCR library free sequencing methodology for detection of clinically significant CNVs in prenatal samples. The optimized rCNV-seq method accommodates input DNA sheared enzymatically without additional size selection steps and then combines the end repair and dA (adenine)-tailing into one step, removing the need for a PCR process and making library preparation much faster and easier (Figure 1). Equally important, library construction can be completed in a single tube, making the workflow more conducive to automation. In validation studies, rCNV-seq was highly sensitive and specific for the detection of common whole chromosome aneuploidies such as T21, T18, T13, 45,XO, 47,XXY, 47,XXX and 47,XYY as well as for segmental aneuploidies associated with different types of MMS. Detection of CNVs was highly reproducible, even when gDNA amounts as low as 10 ng were used to construct PCR libraries. In clinical studies of amniocentesis prenatal samples at risk for a fetal CNV, we also showed the ability of rCNV-seq to correctly identify both pathogenic and VOUS CNVs.
One of the main drivers for reliable and accurate calling of CNVs is the quality and depth of the sequencing data that are binned across each chromosome. In this regard, rCNV-seq based on PCR-free libraries demonstrated improved performance over current CNV-seq methods, showing an average Q30 for the sequencing data as high as 94%. Compared to PCR, PCR-free-frag libraries achieved a higher proportion of uniquely mapped reads (fewer duplicate reads), a lower read redundancy and a smaller CV for all chromosomes, which are important parameters to improve the signal-to-noise ratio and allow more accurate detection of CNVs and their correct copy number. In addition, rCNV-seq exhibited some additional advantages for implementation of the method into molecular diagnostic laboratories with NGS capacity. Firstly, PCR-free libraries can be generated in around 2 h compared to 4–5 h for current PCR-free and PCR-dependent commercial methods. Thus, from sample receipt in the morning and with a sequencing run time of 6.5 h on the Next-Seq platform, it is possible to generate same day results and improves the overall efficiency of laboratory workflow and staff time. Second, the method can generate representative libraries from a wide range of input DNA amounts, improving versatility for the analysis of different sample types and compromised samples where the gDNA is limiting.
Benchmarking against CMA, rCNV-seq demonstrated a high degree of accuracy and reproducibility. For the majority of CNVs tested, there was a high concordance between the two methods for correctly calling the CNV interval and copy number. There were some exceptions where the CNV interval varied up to 33% for the genome location called. These differences are probably related to the inherent limitations of both techniques, whereby CMA has reduced probe coverage for some genome positions, whereas rCNV-seq has reduced coverage in highly repetitive regions. While high sensitivity and specificity is an important parameter for calling CNVs, reproducibility is also a key factor for defining the resolution of CNV detection, since each rCNV-seq analysis is based on a set of randomly mapped sequencing reads mapped across the 24 chromosomes. Previous studies have shown that while CNVs as small as 0.1 Mb can be detected by CNV-seq [13], there were no data provided on the reproducibility of detection. In this study, using gDNA inputs of 10 ng, high reproducibility of CNV detection was demonstrated for CNVs as large as 31.9 Mb to as small as 0.7 Mb. Nonetheless, further studies on a range of smaller CNVs are still needed to determine the minimum resolution of rCNV-seq for a single-pass diagnostic test.
Based on our findings, rCNV-seq is also suitable as an alternative to CMA for routine prenatal testing where a gDNA sample is generally available from cultured amniocytes. However, there are also clinical situations whereby rCNV-seq has some clear advantages over CMA. Firstly, some women are referred for an invasive chromosome test close to the legal limit for termination of pregnancy, which is usually 20 weeks in most countries [22]. Secondly, the amniocentesis samples collected can be substandard, lacking sufficient fetal cells for analysis and, occasionally, there can be cell culture failure [23]. In these two scenarios, it may not be possible to use CMA methods, because they rely on cultured amniocytes (one to two weeks) and a minimum of 0.5–1 µg of input DNA for analysis. In these urgent clinical situations, it will be possible to use rCNV-seq to achieve a rapid and accurate diagnosis using either low amounts of gDNA from uncultured amniocytes or even the available cell free fetal DNA, which has been previously shown to be an ideal NGS template for CNV detection in pregnancies with abnormal fetal ultrasound structural abnormalities [24].

5. Conclusions

PCR-free-based rCNV-seq is a rapid, sensitive, specific and efficient method that can be used in any NGS-based diagnostic laboratory to achieve a reliable and accurate CNV diagnosis. Further studies on a large number of amniocentesis samples are needed to define the true clinical utility and versatility of rCNV-seq for prenatal chromosome testing.

Author Contributions

Conceptualization, X.Z., J.L.; methodology X.C., X.Z.; prospective clinical study validation, X.Z., Y.J., Q.Q., N.H.; formal analysis, X.C., C.L., M.X.; data curation, X.C., X.Z., Y.J., Q.Q., N.H.; writing—review and editing, D.S.C., J.L.; project administration and funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the CAMS Innovation Fund for Medical Sciences (CIFMS), grant number 2016-I2M-1-008; 2020-I2M-C&T-B-046.

Institutional Review Board Statement

The research study was approved by the Ethics Committee for Drug Clinical Trials in PUMCH (approval number KS2019136).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

X.C., C.L., M.X., and D.S.C. are employees of Berry Genomics. None have stock or bond holdings.

References

  1. Mikhail, F.M. Copy number variations and human genetic disease. Curr. Opin. Pediatr. 2014, 26, 646–652. [Google Scholar] [CrossRef]
  2. Daughtry, B.L.; Chavez, S.L. Chromosomal instability in mammalian pre-implantation embryos: Potential causes, detection methods, and clinical consequences. Cell Tissue Res. 2016, 363, 201–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bagheri, H.; Mercier, E.; Qiao, Y.; Stephenson, M.D.; Rajcan-Separovic, E. Genomic characteristics of miscarriage copy number variants. Mol. Hum. Reprod. 2015, 21, 655–661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chen, Y.; Bartanus, J.; Liang, D.; Zhu, H.; Breman, A.M.; Smith, J.L.; Wang, H.; Ren, Z.; Patel, A.; Stankiewicz, P.; et al. Characterization of chromosomal abnormalities in pregnancy losses reveals critical genes and loci for human early development. Hum. Mutat. 2017, 38, 669–677. [Google Scholar] [CrossRef] [PubMed]
  5. Wapner, R.J.; Martin, C.L.; Levy, B.; Ballif, B.C.; Eng, C.M.; Zachary, J.M.; Savage, M.; Platt, L.D.; Saltzman, D.; Grobman, W.A.; et al. Chromosomal microarray versus karyotyping for prenatal diagnosis. N. Engl. J. Med. 2012, 367, 2175–2184. [Google Scholar] [CrossRef] [Green Version]
  6. Moorthie, S.; Blencowe, H.; Darlison, M.W.; Gibbons, S.; Lawn, J.E.; Mastroiacovo, P.; Morris, J.K.; Modell, B.; Congenital Disorders Expert Group. Chromosomal disorders: Estimating baseline birth prevalence and pregnancy outcomes worldwide. J. Community Genet. 2018, 9, 377–386. [Google Scholar] [CrossRef]
  7. Goldenberg, P. An Update on Common Chromosome Microdeletion and Microduplication Syndromes. Pediatr. Ann. 2018, 47, e198–e203. [Google Scholar] [CrossRef]
  8. Lopez-Rivera, E.; Liu, Y.P.; Verbitsky, M.; Anderson, B.R.; Capone, V.P.; Otto, E.A.; Yan, Z.; Mitrotti, A.; Martino, J.; Steers, N.J.; et al. Genetic Drivers of Kidney Defects in the DiGeorge Syndrome. N. Engl. J. Med. 2017, 376, 742–754. [Google Scholar] [CrossRef] [Green Version]
  9. Levy, B.; Wapner, R. Prenatal diagnosis by chromosomal microarray analysis. Fertil. Steril. 2018, 109, 201–212. [Google Scholar] [CrossRef] [Green Version]
  10. Schaaf, C.P.; Wiszniewska, J.; Beaudet, A.L. Copy number and SNP arrays in clinical diagnostics. Annu. Rev. Genom. Hum. Genet. 2011, 12, 25–51. [Google Scholar] [CrossRef] [Green Version]
  11. Srebniak, M.I.; Diderich, K.E.; Joosten, M.; Govaerts, L.C.; Knijnenburg, J.; de Vries, F.A.; Boter, M.; Lont, D.; Knapen, M.F.; de Wit, M.C.; et al. Prenatal SNP array testing in 1000 fetuses with ultrasound anomalies: Causative, unexpected and susceptibility CNVs. Eur. J. Hum. Genet. 2016, 24, 645–651. [Google Scholar] [CrossRef] [Green Version]
  12. Xie, C.; Tammi, M.T. CNV-seq, a new method to detect copy number variation using high-throughput sequencing. BMC Bioinform. 2009, 10, 80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liang, D.; Peng, Y.; Lv, W.; Deng, L.; Zhang, Y.; Li, H.; Yang, P.; Zhang, J.; Song, Z.; Xu, G.; et al. Copy number variation sequencing for comprehensive diagnosis of chromosome disease syndromes. J. Mol. Diagn. 2014, 16, 519–526. [Google Scholar] [CrossRef] [PubMed]
  14. Dong, Z.; Xie, W.; Chen, H.; Xu, J.; Wang, H.; Li, Y.; Wang, J.; Chen, F.; Choy, K.W.; Jiang, H. Copy-Number Variants Detection by Low-Pass Whole-Genome Sequencing. Curr. Protoc. Hum. Genet. 2017, 94, 8.17.1–8.17.16. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, J.; Chen, L.; Zhou, C.; Wang, L.; Xie, H.; Xiao, Y.; Zhu, H.; Hu, T.; Zhang, Z.; Zhu, Q.; et al. Prospective chromosome analysis of 3429 amniocentesis samples in China using copy number variation sequencing. Am. J. Obstet. Gynecol. 2018, 219, 287.e1–287.e18. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, H.; Dong, Z.; Zhang, R.; Chau, M.H.K.; Yang, Z.; Tsang, K.Y.C.; Wong, H.K.; Gui, B.; Meng, Z.; Xiao, K.; et al. Low-pass genome sequencing versus chromosomal microarray analysis: Implementation in prenatal diagnosis. Genet. Med. 2020, 22, 500–510. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, S.; Song, L.; Cram, D.S.; Xiong, L.; Wang, K.; Wu, R.; Liu, J.; Deng, K.; Jia, B.; Zhong, M.; et al. Traditional karyotyping vs copy number variation sequencing for detection of chromosomal abnormalities associated with spontaneous miscarriage. Ultrasound Obstet. Gynecol. 2015, 46, 472–477. [Google Scholar] [CrossRef] [Green Version]
  18. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Y.; Chen, Y.; Tian, F.; Zhang, J.; Song, Z.; Wu, Y.; Han, X.; Hu, W.; Ma, D.; Cram, D.; et al. Maternal mosaicism is a significant contributor to discordant sex chromosomal aneuploidies associated with noninvasive prenatal testing. Clin. Chem. 2014, 60, 251–259. [Google Scholar] [CrossRef]
  20. Riggs, E.R.; Andersen, E.F.; Cherry, A.M.; Kantarci, S.; Kearney, H.; Patel, A.; Raca, G.; Ritter, D.I.; South, S.T.; Thorland, E.C.; et al. Technical standards for the interpretation and reporting of constitutional copy-number variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet. Med. 2020, 22, 245–257. [Google Scholar] [CrossRef] [Green Version]
  21. Howe, B.; Umrigar, A.; Tsien, F. Chromosome preparation from cultured cells. J. Vis. Exp. 2014, e50203. [Google Scholar] [CrossRef] [Green Version]
  22. Rayburn, W.F.; LaFerla, J.J. Second-trimester pregnancy termination for genetic abnormalities. J. Reprod. Med. 1982, 27, 584–588. [Google Scholar] [PubMed]
  23. Alfirevic, Z.; Navaratnam, K.; Mujezinovic, F. Amniocentesis and chorionic villus sampling for prenatal diagnosis. Cochrane Database Syst. Rev. 2017, 9, CD003252. [Google Scholar] [CrossRef] [PubMed]
  24. Qi, Q.; Lu, S.; Zhou, X.; Yao, F.; Hao, N.; Yin, G.; Li, W.; Bai, J.; Li, N.; Cram, D.S. Copy number variation sequencing-based prenatal diagnosis using cell-free fetal DNA in amniotic fluid. Prenat. Diagn. 2016, 36, 576–583. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overview of next-generation sequencing (NGS) library construction workflows. The overall time required to generate PCR-free-frag, PCR-free-soni and PCR-soni libraries was 2.1, 3.75 and 4.75 h, respectively.
Figure 1. Overview of next-generation sequencing (NGS) library construction workflows. The overall time required to generate PCR-free-frag, PCR-free-soni and PCR-soni libraries was 2.1, 3.75 and 4.75 h, respectively.
Life 11 00098 g001
Figure 2. Comparison of key sequencing QC data to evaluate the performance of NGS library construction methods. (a) Ratio of unique read mapping. (b) Ratio of read redundancy. (c) Coefficient of variation (CV) (SD/mean) across the whole genome. (d) CV for each chromosome. (e) The uniformity of genome sequence coverage measured by median read depth. (f) Unique GC ratio of sequencing reads. (g) Q30 (chance that the base call is incorrect is 1 in 1000).
Figure 2. Comparison of key sequencing QC data to evaluate the performance of NGS library construction methods. (a) Ratio of unique read mapping. (b) Ratio of read redundancy. (c) Coefficient of variation (CV) (SD/mean) across the whole genome. (d) CV for each chromosome. (e) The uniformity of genome sequence coverage measured by median read depth. (f) Unique GC ratio of sequencing reads. (g) Q30 (chance that the base call is incorrect is 1 in 1000).
Life 11 00098 g002
Figure 3. Performance of rCNV sequencing using different amounts of DNA inputs. (a) Library yield. (b) Ratio of unique read mapping and read redundancy. (c) CV. (d) Unique GC ratio. (e) Median read depth.
Figure 3. Performance of rCNV sequencing using different amounts of DNA inputs. (a) Library yield. (b) Ratio of unique read mapping and read redundancy. (c) CV. (d) Unique GC ratio. (e) Median read depth.
Life 11 00098 g003
Figure 4. Chromosome plots produced from rCNV-seq analysis of samples with known CNVs. Y-axis, copy number and X-axis, copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Figure 4. Chromosome plots produced from rCNV-seq analysis of samples with known CNVs. Y-axis, copy number and X-axis, copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Life 11 00098 g004
Figure 5. Reproducibility of rCNV-seq for detection of CNVs, copy number (CN) and genome interval coordinates. The top panel describes how the CNV intervals detected by rCNV-seq were compared to CMA results. The bottom panel provides the level of CNV concordance (%) between rCNV-seq and CMA.
Figure 5. Reproducibility of rCNV-seq for detection of CNVs, copy number (CN) and genome interval coordinates. The top panel describes how the CNV intervals detected by rCNV-seq were compared to CMA results. The bottom panel provides the level of CNV concordance (%) between rCNV-seq and CMA.
Life 11 00098 g005
Figure 6. Confirmation of whole chromosome aneuploidies by karyotyping. A. Karyotyping. Red arrows indicate chromosome gains. B. rCNV-seq results. The chromosome copy number gains determined by rCNV-seq were concordant with karyotyping. rCNV-seq plots. Y-axis copy number and X-axis mean copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Figure 6. Confirmation of whole chromosome aneuploidies by karyotyping. A. Karyotyping. Red arrows indicate chromosome gains. B. rCNV-seq results. The chromosome copy number gains determined by rCNV-seq were concordant with karyotyping. rCNV-seq plots. Y-axis copy number and X-axis mean copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Life 11 00098 g006
Figure 7. rCNV-seq detection of CNVs in 217 clinical samples. All CNVs identified were classified as VOUS. Y-axis copy number and X-axis mean copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Figure 7. rCNV-seq detection of CNVs in 217 clinical samples. All CNVs identified were classified as VOUS. Y-axis copy number and X-axis mean copy number of the sequencing bins (gray dots). Blue line signifies the mean copy along the length of the chromosome. Red line indicates regions of repetitive DNA, and the black box denotes the centromere.
Life 11 00098 g007
Table 1. Validation of rCNV-seq for detection of known CNVs.
Table 1. Validation of rCNV-seq for detection of known CNVs.
Chromosome DisordersCMArCNV-Seq (PCR-Free-Frag)Interpretation
Common trisomies
T21 ǂ47XN,+2147XN,+21Concordant
T1847XN,+1847XN,+18Concordant
T13 ǂ47XN,+1347XN,+13Concordant
Sex chromosome aneuploidies (SCAs)
45,XO ǂ45,XO45,XOConcordant
47,XXY ǂ47,XXY47,XXYConcordant
47,XXX47,XXX47,XXXConcordant
47,XYY47,XYY47,XYYConcordant
Microdeletion/microduplication syndromes (MMS)
Cri du Chat syndrome (5p15.3-p15.2)
Sample 1 ǂ5p15.33p13.3(113576-31928290) ×1, 31.81 Mb5p15.33-p13.3(del, 31.90 Mb)Concordant
Sample 25p15.2p15.33(151,737-14,756,030) ×1, 14.6 Mb5p15.33-p15.2(del, 17.30 Mb)Concordant
Sample 35p13.3p15.33(151,737-33,120,547) ×1, 32.97 Mb5p15.33-p13.3(del, 33.14 Mb)Concordant
Williams–Beuren syndrome (7q11.23)
Sample 1 ǂ7q11.23(72653992-74146927) ×1, 1.49 Mb7q11.23(del, 1.32 Mb)Concordant
Sample 2 ǂ7q11.23(72549979-74374748) ×1, 1.82 Mb7q11.23(del, 1.36 Mb)Concordant
Sample 37q11.23(72,726,578-74,139,390) ×1, 1.41 Mb7q11.23(del, 1.39 Mb)Concordant
Wolf–Hirschhorn syndrome (4p16.3)
Sample 1 ǂ4p16.1p16.3(71,552-8,510,870) ×1, 8.44 Mb4p16.3-p16.1(del, 8.48 Mb)Concordant
Sample 24p16.1p16.3(143,413-8,368,180) ×1, 8.22 Mb4p16.3-p16.1(del, 8.35 Mb)Concordant
DiGeorge syndrome (22q11.2)
Sample 1 ǂ22q11.21(18648855-21800471) ×1, 3.15 Mb22q11.21(del, 2.49 Mb)Concordant
Sample 222q11.21(18,919,942-21,440,514) ×1, 2.52 Mb22q11.21(del, 2.55 Mb)Concordant
Prader-Willi/Angelman syndrome (15q11-q13)
Sample 1 ǂ15q11.2q13.1(23290787-28560664) ×1, 5.27 Mb15q11.2-q13.1(del, 4.84 Mb)Concordant
Smith–Magenis syndrome (17p11.2)
Sample 117p11.2p12(15,816,892-20,193,169) ×1, 4.38 Mb17p12-p11.2(del, 4.92 Mb)Concordant
Miller–Dieker syndrome (17p13.3-p13.2)
Sample 1 ǂ17p13.3p13.2(1,972,209-5,945,876) ×1, 4.0 Mb17p13.3p13.2(del, 5.95 Mb)Concordant
VOUS (8p23.2)
Sample 1 ǂ8p23.2(3,427,306-4,275,392) ×1, 0.8 Mb8p23.2(del, 0.83 Mb)Concordant
VOUS (11q14.3)
Sample 111q14.3(91,347,506-92,085,142) ×3, 0.7 Mb11q14.3(dup, 0.73 Mb)Concordant
VOUS (17p13.3)
Sample 117p13.3(1,004,599-1,518,383) ×3, 0.51 Mb17p13.3(dup,0.68 Mb)Concordant
VOUS (5p12-p11)
Sample 1 ǂ5p12p11(45,455,695-46,242,541) ×3, 0.78 Mb5p12-p11(dup, 0.85 Mb)Concordant
ǂ Amniotic fluid.
Table 2. Clinical performance of rCNV-seq for detection of CNVs in 217 amniotic fluid samples.
Table 2. Clinical performance of rCNV-seq for detection of CNVs in 217 amniotic fluid samples.
Primary Indication for Chromosome TestingSample (n)rCNV-Seq Result
AbnormalNormal
Advanced Maternal Age (≥35 Years)1524 (T21)148
Abnormal serum screening result (<35 years)412 (T21 and 47,XXX)39
Previously given birth to child with a chromosomal disorder81 (T18)7
Structural abnormality by ultrasound71 (T21)6
One parent is a carrier of a chromosomal disorder202
Advanced maternal age (≥35 years) and secondary indication #707
Total2178209
# Secondary indications included abnormal serum screening result, previously given birth to child with a chromosomal disorder, structural abnormality by ultrasound or one parent is a carrier of a chromosomal disorder.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, X.; Chen, X.; Jiang, Y.; Qi, Q.; Hao, N.; Liu, C.; Xu, M.; Cram, D.S.; Liu, J. A Rapid PCR-Free Next-Generation Sequencing Method for the Detection of Copy Number Variations in Prenatal Samples. Life 2021, 11, 98. https://doi.org/10.3390/life11020098

AMA Style

Zhou X, Chen X, Jiang Y, Qi Q, Hao N, Liu C, Xu M, Cram DS, Liu J. A Rapid PCR-Free Next-Generation Sequencing Method for the Detection of Copy Number Variations in Prenatal Samples. Life. 2021; 11(2):98. https://doi.org/10.3390/life11020098

Chicago/Turabian Style

Zhou, Xiya, Xiangbin Chen, Yulin Jiang, Qingwei Qi, Na Hao, Chengkun Liu, Mengnan Xu, David S. Cram, and Juntao Liu. 2021. "A Rapid PCR-Free Next-Generation Sequencing Method for the Detection of Copy Number Variations in Prenatal Samples" Life 11, no. 2: 98. https://doi.org/10.3390/life11020098

APA Style

Zhou, X., Chen, X., Jiang, Y., Qi, Q., Hao, N., Liu, C., Xu, M., Cram, D. S., & Liu, J. (2021). A Rapid PCR-Free Next-Generation Sequencing Method for the Detection of Copy Number Variations in Prenatal Samples. Life, 11(2), 98. https://doi.org/10.3390/life11020098

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop