Abstract
Objective: To examine the performance of chromosome-selective sequencing of cell-free (cf) DNA in maternal blood for assessment of fetal sex chromosome aneuploidies. Methods: This was a case-control study of 177 stored maternal plasma samples, obtained before fetal karyotyping at 11-13 weeks of gestation, from 59 singleton pregnancies with fetal sex chromosome aneuploidies (45,X, n = 49; 47,XXX, n = 6; 47,XXY, n = 1; 47,XYY, n = 3) and 118 with euploid fetuses (46,XY, n = 59; 46,XX, n = 59). Digital analysis of selected regions (DANSR™) on chromosomes 21, 18, 13, X and Y was performed and the fetal-fraction optimized risk of trisomy evaluation (FORTE™) algorithm was used to estimate the risk for non-disomic genotypes. Performance was calculated at a risk cut-off of 1:100. Results: Analysis of cfDNA provided risk scores for 172 (97.2%) samples; 4 samples (45,X, n = 2; 46,XY, n = 1; 46,XX, n = 1) had an insufficient fetal cfDNA fraction for reliable testing and 1 case (47,XXX) failed laboratory quality control metrics. The classification was correct in 43 (91.5%) of 47 cases of 45,X, all 5 of 47,XXX, 1 of 47,XXY and 3 of 47,XYY. There were no false-positive results for monosomy X. Discussion: Analysis of cfDNA by chromosome-selective sequencing can correctly classify fetal sex chromosome aneuploidy with reasonably high sensitivity.
Introduction
Prenatal screening for fetal aneuploidies has traditionally focused on trisomy 21, and more recently on trisomies 18 and 13. Sex chromosome aneuploidies, including monosomy X (45,X), Klinefelter syndrome (47,XXY or 48,XXYY), Triple X syndrome (47,XXX), and 47,XYY, with a combined prevalence of 1:500 are more common than the major trisomies [1,2,3]. Although most cases of sex chromosome aneuploidies are generally mild without intellectual disability, some have a well-established phenotype that can include physical abnormalities, learning delays and infertility [1,2,3]. It may therefore be desirable to some parents that these conditions could be diagnosed prenatally with the option of pregnancy termination [4,5,6,7,8,9,10]. However, the traditional methods of screening for trisomies, including maternal age, maternal serum biochemical testing and ultrasound examination of the fetus, are not effective in detecting sex chromosome aneuploidies, except cases of Turner syndrome presenting with cystic hygromas.
The introduction of cell-free (cf) DNA analysis in maternal blood has now made it possible to screen not only for trisomies 21, 18 and 13, but also potentially for other aneuploidies, including sex chromosome aneuploidies. Chromosome-selective sequencing of cfDNA, referred to as digital analysis of selected regions (DANSR™), combined with an algorithm (fetal-fraction optimized risk of trisomy evaluation; FORTE™) that accounts for a priori risk factors and fetal fraction, has been shown in a variety of clinical settings to provide an accurate risk assessment for trisomies 21, 18 and 13 [11,12,13,14,15,16,17,18]. This same approach lends itself to possible expansion of the testing platform to include evaluation of chromosomes X and Y. The objective of this case-control study is to report the clinical performance of chromosome-selective sequencing of cfDNA in maternal blood and the FORTE algorithm for the assessment of fetal sex chromosome aneuploidies.
Methods
Study Population
This was a case-control study of 177 stored maternal plasma samples from 59 singleton pregnancies with fetal sex chromosome aneuploidies (45,X, n = 49; 47,XXX, n = 6; 47,XXY, n = 1; 47,XYY, n = 3) and 118 with euploid fetuses (46,XY, n = 59; 46,XX, n = 59). No cases of fetal mosaicism were included. All samples were collected before invasive testing for fetal karyotyping at 11-13 weeks of gestation because screening by the combined test of maternal age, fetal nuchal translucency (NT) thickness and serum free β-hCG and PAPP-A indicated an increased risk for fetal trisomies [19]. Gestational age was determined from the measurement of the fetal crown-rump length [20].
Maternal venous blood (10 ml) collected before chorionic villous sampling in ethylene-diamine-tetraacetic acid, EDTA, BD Vacutainer™ tubes (Becton Dickinson UK Ltd., Oxford, UK) was processed within 15 min of collection. The samples were centrifuged at 2,000 g for 10 min to separate the plasma from packed cells and buffy coat, and subsequently at 16,000 g for 10 min to further separate cell debris. Plasma samples (2 ml each) were divided into 0.5-ml aliquots in separate eppendorf tubes, which were labeled with a unique patient identifier and stored at -80°C until subsequent analysis. The study was approved by the UK NHS Research Ethics Committee and all patients gave written informed consent to participate.
We searched our database and selected 59 cases of sex chromosome aneuploidies with available stored plasma and then selected samples from 59 pregnancies with euploid male fetuses and 59 with euploid female fetuses. Plasma samples (4 tubes of 0.5 ml per patient) from selected cases were sent overnight on dry ice from London to the USA for cfDNA analysis (Harmony™ Prenatal Test, Ariosa Diagnostics Inc., San Jose, Calif., USA). The information given to the laboratory for each case was: the patient unique identifier, maternal age, maternal weight, method of conception and date of blood collection. The laboratory was informed that the samples were collected from pregnancies with sex chromosome aneuploidies and euploid controls, but not the karyotype of individual samples.
Laboratory Analysis
cfDNA was extracted from maternal plasma, chromosome-selective sequencing with DANSR was carried out and analyzed using a FORTE algorithm as previously described [11,12]. All laboratory personnel were blinded to the identity and origin of plasma samples.
DANSR on chromosomes 21, 18, 13, X and Y was performed and a FORTE algorithm was used to estimate risk for non-disomic genotypes [11,12]. For assessment of risk for sex chromosome aneuploidies, the original DANSR assays were expanded to include 32 regions on chromosome Y and 599 on chromosome X. Similarly, the FORTE algorithm used to assess trisomies 21, 18 and 13 [11,12] was adapted to assess sex chromosome aneuploidies. Models of observing proportions of median Y or X assay counts relative to median assay counts on chromosomes 13, 18, and 21 were based on normal or truncated normal (in the case of no Y presence) distributions using standard deviations estimated by Monte Carlo simulations. All models assuming Y chromosome copy counts of 0, 1 or 2, and X chromosome copy counts of 1, 2 or 3 were constructed and combined to consider joint probabilities representing genotypes of monosomy X, XX, XY, XXX, XXY, XYY and XXYY, making the assumption that observed X and Y proportions are independent.
For example, in the case of 45,X the model can be described by the formula:
P (pX|X = 1, f) P (pY|Y = 0, f) P (X = 1, Y = 0|MA, GA),
where pX is the proportion of X assay counts, pY is the proportion of Y assay counts, f is the observed fetal fraction, and P (X = 1, Y = 0|MA, GA) is the prevalence of monosomy X for the given maternal age (MA) and gestational age (GA). Values derived from these models were normalized using the standard Bayesian approach of summing values across all possible genotypes to obtain a final score for each genotype. Non-disomic genotypes were reported when their score was above 1 in 10,000, otherwise the disomic genotype with the highest score was reported as long as it was greater than 99%. A risk cut-off of 1 in 100 for non-disomic genotypes was used for calculation of the detection rate (DR) and false-positive rate (FPR).
Results
Maternal and pregnancy characteristics of the study population are summarized in table 1. In the euploid group, the median maternal age (36 years), NT thickness (2.8 mm) and serum-free β-hCG (1.4 multiples of the normal medium; MoM) were higher and serum PAPP-A (0.7 MoM) was lower than the respective values of 31 years, 1.8 mm, 1.0 and 1.0 MoM in more than 87,000 pregnancies which had first trimester combined screening in our center [21]. In monosomy X, the maternal age and serum-free β-hCG were similar to our screened population, but the median fetal NT (8.3 mm) was very much higher and serum PAPP-A (0.4 MoM) was lower; the NT was ≥3.5 mm in 45 (91.8%) of the 49 cases. In the group with other sex chromosome aneuploidies, compared to our screened population, median maternal age and fetal NT were higher, but serum-free β-hCG and PAPP-A were lower.
Analysis of cfDNA provided risk scores for 172 (97.2%) of the samples. Four samples (45,X, n = 2; 46,XY, n = 1; 46,XX, n = 1) had an insufficient fetal cfDNA fraction for reliable testing and 1 case (47,XXX) failed laboratory quality control metrics.
In the 172 cases with risk scores, the median fetal fraction was 12.0% (range 4.8-32.0%). The performance of screening by cfDNA analysis in the assessment of sex chromosome aneuploidies is summarized in table 2. The classification was correct in 43 (91.5%) of 47 cases of 45,X, all 5 with 47,XXX, 1 of 47,XXY and 3 of 47,XYY. In 115 of the 116 euploid pregnancies with results, the classifications were correct, but in 1 case of 46,XX the classification was 47,XXX with a risk score of 55 in 100.
The distribution of estimated risk for monosomy X by cfDNA analysis in maternal blood from pregnancies with 45,X and euploid fetuses, plotted against fetal fraction is shown in figure 1. There were 39 cases with a risk of >99 in 100, and 1 each with a risk of 91 in 100, 90 in 100, 80 in 100, 2.8 in 100, 1.6 in 100, 1 in 285, 1 in 2,500, and 1 in 7,692 and <1 in 10,000. Consequently, 4 of the 47 cases were classified as screen negative because the risk was below the pre-specified cut-off of 1 in 100; 3 were classified as normal female and 1 as normal male. In the euploid pregnancies there were no false-positive cases for monosomy X.
Discussion
The findings of this study demonstrate the feasibility of chromosome-selective sequencing of cfDNA in maternal blood for assessment of risk for fetal monosomy X and other sex chromosome aneuploidies; 91.5% of cases of 45,X were identified at an FPR of 0%. Although all cases with 47,XYY, 47,XXY and 47,XXX were correctly classified, the number of cases examined is too small for definite conclusions to be drawn for the performance of screening for these aneuploidies.
Several studies have reported on the feasibility of cfDNA analysis in the diagnosis of 45,X and other sex chromosome aneuploidies, but only four studies have examined more than 10 cases. Bianchi et al. [22] carried out massively parallel shotgun sequencing (MPSS) in 20 cases of 45,X and 462 non 45,X. In the 45,X group, 15 were classified correctly, 1 was classified as non-45,X and 4 were unclassified, giving a DR of 75.0%. In the non-45,X, 416 were classified correctly, 1 was classified as 45,X and 45 were unclassified, giving an FPR of 0.2%. It is of interest that in the case of trisomies 21, 18 and 13 the rate of unclassified result was only 0.4-1.4%, compared to 9.9% for 45,X. Similarly, Mazloom et al. [23] carried out MPSS in 21 cases of 45,X and 390 non-45,X. The result from cfDNA analysis was unclassified in 21 (5.1%) cases (45,X, n = 3; non-45,X, n = 18) and the overall DR and FPR were 81.0 and 0.3%, respectively. A proof of principle study investigating the feasibility of developing a genome-wide assay of cfDNA in maternal plasma capable of replacing invasive testing for major aneuploidies, which used MPSS at 20 million reads per sample in 176 pregnancies, correctly classified all 15 cases of 45,X at an FPR of 0% [24]. Samango-Sprouse et al. [25] used targeted sequencing and allelic ratio analysis of SNPs covering chromosomes 21, 18, 13, X and Y, to examine 13 cases of 45,X and 188 non-45,X. In 14 (7.0%) cases, including 1 case of 45,X, the samples did not pass quality control. In those with results, the DR and FPR for 45,X were 91.7 and 0%, respectively. The DR of monosomy X from the combined data of these four studies was 87% (59 of 68) at FPR of 0.2% (2 of 1,188) [22,23,24,25].
The poorer performance of cfDNA analysis in screening for X chromosome aneuploidy, as well as trisomy 13 compared to trisomies 21 and 18 using MPSS technology could, at least in part, be due to the highly variable amplification of chromosomes X and 13 because of lower guanosine-cytosine content [26,27,28,29]. This would not be a factor for the DANSR methodology as loci selection avoids regions of the chromosome where this is an issue.
An additional factor that complicates cfDNA screening for sex chromosome aneuploidies is the high incidence of maternal and fetal mosaicism for these aneuploidies. For example, fetal mosaicism can account for up to 50% of sex aneuploidy cases [2,3,30,31]. This could result in a quantitative assessment of the X and Y chromosome from cfDNA analysis which is not concordant with the fetal karyotype as assessed by either CVS or amniocentesis. In the case of maternal sex chromosome mosaicism, whether previously known or unknown at the time of testing, the assay might result in counts of X chromosome fragments predicted as a fetal sex chromosome aneuploidy when in fact the fetus is euploid. In normal females there is an age-related loss of an X chromosome in white blood cells [32]. Moreover, an unknown maternal karyotype of 47,XXX could confound assay results because up to 90% of such women are not aware that they have a third X chromosome [3]. A limitation of our study was that we did not perform maternal karyotyping to determine if there were any cases of mosaicism.
Conventional prenatal screening has never directly sought to uncover fetal sex chromosome aneuploidies, and their detection was coincidental in pregnancies undergoing invasive testing following screening for trisomy 21 [33,34]. The objective of screening for trisomy 21 over the last 4 decades has been to increase the DR and decrease the rate of unnecessary invasive tests. Such decrease in invasive testing would inevitably reduce the coincidental detection of sex chromosome aneuploidies other than monosomy X. In patients presenting with very high fetal NT during the first trimester or cystic hygromas/hydrops during the second trimester, it could be argued that the investigation of choice would be invasive testing for fetal karyotype evaluation, including a sub-chromosomal analysis with microarray [35], rather than cfDNA analysis for assessment of risk for 45,X.
This study has confirmed the feasibility of cfDNA analysis in the assessment of risk for fetal monosomy X and other sex chromosome aneuploidies. Whether or not we should be screening the pregnant population for these aneuploidies, in light of their wide yet mild phenotypic spectrum and the possibility of uncovering a previously unknown aneuploidy in the mother, warrants further consideration of both individual patient preferences and the clinical utility of such an endeavor. Another important factor to be considered in expanding the indications of cfDNA testing from screening for trisomies to include sex chromosome and other aneuploidies is the decrease in efficiency of screening because of the likely increase in cumulative FPR.
Acknowledgments
The study was supported by a grant from The Fetal Medicine Foundation (UK Charity No. 1037116). Analysis of samples was performed at their own expense by Ariosa Diagnostics (San Jose, Calif., USA).