|Home | About | Journals | Submit | Contact Us | Français|
Accumulating evidence suggests that genomic structural variations, particularly copy number variations (CNV), are a common occurrence in humans that may bear phenotypic consequences for living individuals possessing the variant. While precise estimates vary, large-scale karyotypic abnormalities are present in 6–12% of stillbirths (SB). However, due to inherent limitations of conventional cytogenetics, the contribution of genomic aberrations to stillbirth is likely underrepresented. High-resolution copy number variant (CNV) analysis by genomic array-based profiling may overcome such limitations.
Prospectively acquired SB cases >22 weeks underwent classification of "unexplained" stillbirth by Wigglesworth and Aberdeen criteria after extensive testing and rigorous multidisciplinary audit. Genome-wide analysis was conducted using high resolution Illumina SNP arrays (HumanCNV370 Duo) on placental and fetal samples. Potential alternate detection methods were completed by one or more of three independent means (quantitative PCR, Illumina1M or Agilent105K CGH arrays).
In our cohort of 54 stillbirths, 29 met strict unexplained criteria. Among these, we identified 24 putative novel CNVs. Subsequent interrogation detected18 of 24 CNVs (75%) in placental samples, 8 of which were also confirmed in available fetal samples; none were present in maternal blood.
We describe the potential of whole-genome placental profiling to identify small genomic imbalances which might contribute to a small proportion of well-characterized, unexplained stillbirths.
It is increasingly evident that genomic structural variations commonly occur in humans (1–2). These structural variations are comprised of inversions, translocations, deletions, duplications, or more complex combinations of such events. Deletions or duplications result in a change in the number of copies of a genomic segment and are referred to as copy number variations (CNVs). It has been suggested that structural variations such as CNVs contribute more to the genomic difference between individuals than single nucleotide polymorphisms (SNPs) (3). At least one potentially beneficial reproductive phenotype, increased fecundity, has been shown to be associated with a structural variation (4). However, a much larger number of deleterious effects (genomic disorders) have been shown to be associated with imbalances in array-based copy number analyses (5–8).
Recent estimates from a survey of 25,655 stillbirths after the 20th week of pregnancy supports an incidence of 6.2 per 1000 livebirths, or approximately 1/160 (9, 10). Identifiable causes of stillbirth are generally ascribed to one of three categories: fetal anomalies, placentation or cord morbidities, or maternal comorbidities. Accurate determination of the relative contribution of these and other factors, such as environmental exposures, to stillbirth is hampered by inaccurate fetal death records (11, 12) and incomplete evaluation of such potential comorbidities. As such, the cause of as great as 60% of stillbirths remains classified as “unknown” (13–16). The poorly understood epidemiology of stillbirths recently prompted the National Institute of Health to identify it as a high priority research focus, resulting in the Stillbirth Collaborative Research Network (13, 14). Because the empiric recurrence risk of stillbirth approximates 3–5% (10), identification of the etiology in index cases is essential to better stratify recurrence estimates, potentially reduce parental anxiety, and diminish future unnecessary and expensive antenatal testing in subsequent pregnancies (10, 4).
Previous studies to identify genomic abnormalities in stillbirth samples have primarily employed G-banded karyotyping and fluorescence in situ hybridization (FISH). G-banded chromosomal karyotyping at standard 450–600 band resolutions detects aneuploidy and chromosomal segmental abnormalities larger than 3 to 5 megabases (Mb). FISH can detect smaller-size abnormalities in the range of 50–150 Kb, but it is limited by the small number of probes (maximum ~5) that can be tested simultaneously requiring an informed selection of the specific probes to be used. Using these methods, cytogenetic abnormalities have been described in a small percentage (5–13%) of all stillbirths (13, 17, 18). Thus, array based approaches to identifying genomic abnormalities in stillbirths might be expected to have two main advantages over cytogenetic methods. First, depending on genomic coverage and density of the probes on the array, they can offer higher resolution analysis resulting in detection of smaller abnormalities to increase the test resolution by orders of magnitude. Secondly, array based approaches require preparation of DNA from tissue samples, avoiding the need for successful cell culturing. Other investigators have employed array comparative genomic hybridization (aCGH) in the genetic evaluation of anomalous stillbirths (19) as well as instances of fetal loss from 10 to 20 weeks of gestation (20). To date, employment of aCGH in the majority (60%) of instances of stillbirth (i.e., unexplained, non-anomalous stillbirths) has not been undertaken.
Given that CNVs are a common occurrence in humans and that stillbirths are generally only examined using low resolution cytogenetic methods if at all, it seems likely that some proportion of stillbirths ascribed to unknown causes may be due to CNVs. We hypothesized that in a well-described population-based cohort, a small but definable number of cases of unexplained stillbirth would be associated with small, cytogenetically undetected copy number imbalances identifiable by array-based platforms. Our overall objective was to estimate the genome-wide prevalence of copy number changes in truly unexplained stillbirths.
Following IRB approval, placental samples were prospectively acquired from stillbirths (n=33) and control livebirths (n=10) between January 2005 and November 2007 collected at the Policlinico Hospital of Modena, Italy. The diagnosis of stillbirth was based on the World Health Organization criteria defined as fetal death >22 weeks of gestation, or birth weight ≥500g if the gestational age was unknown (7, 29). The IRB protocol allowing the index study with diagnostic evaluation and sample repository with capacity for future use on stillborn infants was approved and granted by the Ministry of Education, University and Research. Subsequent IRB approval was obtained thereafter for genetics and DNA-based research.
All stillbirths underwent evaluation according to a standardized protocol. Following initial evaluation, a panel of 2 obstetricians, 1 neonatologist and 1 pathologist reviewed the cases and classified them based on the Wigglesworth and Aberdeen stillbirth classification systems (21). The cause of stillbirth was assigned only when the majority (>2 of 4) of the clinician panel agreed. Complete parameters for evaluation entailed maternal serologies, expansive microbiology, and fetal karyotype. For maternal serologic testing, TORCH (HIV, toxoplasmosis, rubella, and CMV serologies), maternal serum parvovirus titers (IgG and IgM), maternal measures of thyroid function (TSH, free T3 and T4 in reflex), indirect Coombs, and testing for inherited and acquired thrombophilias (factor V Leiden, PT20210, APLS, anti-cardiolipin antibody, and protein C and S) were systematically obtained following enrolment. Gram staining, acid-fast bacilli (AFB), and aerobic and anaerobic cultures on vaginal, placental, fetal swabs (pre and pos-auricular, auricular, and nasal-pharyngeal cavity), and on fetal blood drawn by intra-cardiac fetal puncture prior to full autopsy. Fetal karyotype to a 600 G-banding resolution was performed on amniocentesis samples or fetal blood, cell cultures from placental biopsies, skin samples, or fascia lata. Finally, complete fetal autopsy was performed by a perinatal pathologist and a detailed placental gross and histological exam was undertaken. All data was recorded, and reviewed in a systematic and unbiased manner by the panel described as above.
DNA isolation was attempted from frozen tissues and/or blood samples employing Gentra Puregene DNA purification kit (Qiagen, Valencia, CA) in accordance with the manufacturer’s protocol and was completed as needed until high purity was reached (Quant-iT PicoGreen dsDNA Genomic, Perkin Elmer, Waltham, MA).
The Illumina CNV370-Duo BeadChip (Illumina Corporation, San Diego, CA) was initially used to perform genome wide copy number profiling on DNA from placental samples. A subset of placental samples in which CNVs were identified was also interrogated using the Illumina 1M-Duo Beadchip (Illumina Corporation, San Diego, CA) following the same methodology. For these assays, a single-base pair extension reaction was performed, followed by staining with fluorescent dyes. The resulting products were imaged with the BeadArray Reader. Allele detection and genotype calling were performed using the BeadStudio software (Illumina Corporation, San Diego, CA) based on human genome build NCBI 36/UCSC hg18. CNVs were identified using the cnvPartition BeadStudio plugin. CNVs identified by cnvPartition that were smaller the 1000bp were removed from consideration. A total of 10 live born controls and 33 stillbirths were interrogated using the array. Of the 33 stillbirth arrays, 3 had a call rate less than 95% and were removed from subsequent analysis.
Three micrograms of genomic DNA from placentas and three micrograms of reference human female or male DNA prepared from single human placentas were labeled and hybridized to the Human Genome CGH Microarray Kit 105A (Agilent Technologies, Inc, Santa Clara, CA) in a rotating oven at 65ºC for 40 hours, according to the manufacturer’s recommended protocol with the following modifications. Briefly, DNA was digested with RsaI and AluI for 2 hours at 37ºC, followed by heat inactivation at 65ºC for 20 minutes. DNA was then labeled with the Agilent Genomic DNA Labeling Kit PLUS using random primers. Experimental samples were labeled with Cyanine-5 dUTP; reference samples were labeled with Cyanine-3 dUTP for two hours at 37ºC. After heat inactivation and sample clean-up, equal amounts of experimental and reference DNA were mixed, heat-denatured and blocked for 30 minutes at 37ºC with Cot-1 DNA in blocking buffer. After blocking, the mixture was applied to the arrays and hybridized in a rotating oven (20 rpm) at 65ºC for 40 hours. Slides were then washed and dried using Agilent’s oligo-array wash buffer 1 and 2, followed by one acetonitrile wash, drying and stabilization wash. Slides were scanned on an Agilent DNA microarray scanner model G2565BA and data extraction was performed with Agilent G2567AA Feature Extraction 9.1 software using the CGH014693 design file as template for automated gridding and the CGH-v4_91 protocol to assign spot-intensity values and ratios to each extraction set. Data were visualized and analyzed based on human genome build NCBI 36/UCSC hg18 with Agilent’s CGHanalytics and displayed as log2 ratios using following settings: 1-fold cut-off, ADM-2 aberration algorithm with threshold 10.0 and a 2 Mb moving-average window.
CNVs identified in stillbirths were compared to CNVs identified in controls and/or present in The Centre for Applied Genomics (TCAG) Database of Genomic Variants (http://projects.tcag.ca/variation/) based on human genome build NCBI 36/UCSC hg18. The TCAG contains CNVs identified in normal (live) individuals from a large number of studies. Stillbirth CNVs completely overlapped by CNVs in live born controls or in the TCAG database were removed from further consideration.
Placental DNA sample qPCR analysis was performed for all 24 identified CNVs in Wigglesworth classification 2 cases. Fetal DNA sample qPCR detection was performed for 19 of 24 CNVs, because fetal samples were not available for cases 46 and 48. For cases not identified as Wiggleworths 2, qPCR analysis was performed on 5 of 6 identified CNVs in both placental and fetal samples. Maternal qPCR was performed on 13 of 24 CNVs for which samples were available. For validation, twenty ng of genomic DNA was amplified using PCR primers designed with Primer Express software (Applied Biosystems, Foster City, CA; Supplementary table). All reactions were done in triplicate. Each target region was amplified using the Biorad iQ5 optical system (Biorad Laboratories, Inc., Hercules, CA) for quantitative real-time PCR (SYBR-Green). Data were normalized against an endogenous reference gene (GAPDH) and a melting curve analysis was performed to verify PCR product specificity. Relative gene copy number was determined by the comparative threshold cycle method (Δ ΔCt) after standard curves with serial dilutions were obtained for each amplification reaction. To obtain ratios and generate graphs of copy numbers, we calculated the ΔCT of threshold cycles of normal female DNA [ΔCT(HF)] minus the ΔCT of threshold cycles of G17 DNA [ΔCT(G17)] and the ΔCT of threshold cycles of normal male DNA [ΔCT(HM)] minus the ΔCT of threshold cycle of G17 DNA [ΔCT(G17)]. As a control, we also calculated the ΔCT of threshold cycles of normal female DNA [ΔCT(HF)] minus the ΔCT of threshold cycles of normal male DNA [ΔCT(HM)]. Ratios for each amplification product and comparison were then calculated as 2ΔΔCT. qPCR 2-ΔΔCT ratios in the following ranges were considered to confirm the array CNV call: homozygous deletion 0–0.25; heterozygous deletion 0.25–0.75; duplication >1.25.
In order to identify stillbirth cases with unknown etiology, cases were examined and classified as described in Methods. The examination included fetal or placental G-banded karyotyping in order to identify cases with microscopically visible chromosomal abnormalities. 33 cases classified as Wigglesworth (21) 2 (unexplained antepartum fetal death) were chosen for this study, and 29 cases were ultimately successfully interrogated by array (Table I). These cases were further subdivided into stillborn fetuses with normal growth (n=17, ≥10th% adjusted for gestational age by birth weight and gender; 33, 34) and intrauterine growth restriction (n=12; <10th% adjusted for gestational age; 33, 34) cases. In addition to examining stillbirths classified as unexplained antepartum fetal death, we also interrogated 4 placental samples from stillbirths with a Wigglesworth classification of 3, 5, or 6 (Table II). As a control, 10 placental samples from liveborn infants were also interrogated by array and will be the subject of future publications. Suffice it to say, none of the large (>1 Mb) noted deletions nor duplications (Table III and andIV)IV) were present in the normal liveborn population (data not shown).
CNVs were initially identified in the 29 placental samples using the Illumina CNV370-Duo Beadchip arrays. CNVs with genomic coordinates that completely overlapped with those CNVs identified in liveborn controls (n=10) or contained in The Centre for Applied Genomics Database of Genomic Variants (http://projects.tcag.ca/variation/) were removed from further consideration. This resulted in the identification of 24 novel CNVs in 17 samples (Tables III and andIV)IV) representing 59% of samples interrogated. The CNVs were comprised of 11 deletions (Table III) and 13 duplications (Table IV) and ranged in size from 8,984 bp to 2,886,321 bp.
We further measured our ability to detect CNVs in non-W2 expanded Wigglesworth classification stillbirths (W5 and W6), by examining stillborns affected by either non-syndromic anomalies (case 2) or common perinatal infections (cases 7, 22, and 47). For example and as detailed in Table II, we interrogated case 2 (W6 classification) because although cranial asymmetry, hypotelorism, hydrops, and a hydrothorax was noted, these findings did not associate with a microscopically detected cytogenetic abnormality (karyotypically normal). Cases 22 and 47 were interrogated because although CMV, listeria, and parvovirus B19 serologies were positive, there were no associated abnormalities such as hydrops, microcephaly, nor fetal calcifications (Table II). Thus in each of these instances we were curious whether we could detect copy number imbalances with our ascribed methodologies. Our detection of novel deletions and duplications was akin to the W2 schema (Table V).
At present, it is unclear what platforms represent “optimal” methods to detect CNVs and, what methods, if any, actually serve to “validate” detected CNVs. In order to examine reproducibility of our findings quantitative PCR and two other platforms were employed. Specifically, the subset of placental samples with identified CNVs on the CNV370-Duo were also interrogated using the more expansive Illumina 1M-Duo Beadchip and the lower density Agilent 105K Comparative Genomic Hybridization Array.
Identified CNVs were examined by qPCR in both the interrogated placental samples and the associated fetal samples with the exception of cases 46 and 48 for which no fetal samples were available (Tables III and andIV).IV). Results of qPCR assays on the placental samples detected 18 of the 24 total identified CNVs (75%). When this was limited to only larger CNVs (>200 Kb) this was noted to be 96% (deletions and duplications; Tables III and andIV)IV) with all placental deletions (8/8 or 100%; Table III) detected on the CNV370-Duo also detected by qPCR.
Of the CNV deletions detected by qPCR on placental DNA, all (6/6, or 100%; Table III) were also confirmed by qPCR assays on fetal DNA samples in the 7 of 8 cases where fetal samples were available. In contrast, among the larger duplications detected by qPCR on placental DNA, only 2/7 (29%) were also detected in fetal tissue. This is likely due to the presence of probable maternal contamination, as validated by our agreement between available maternal samples and placental samples with discordance of fetal tissue (Table IV). These are therefore most likely to represent so-called benign duplications.
With respect to determination of CNV zygosity, three of the deletions were identified by cnvPartition as being homozygous. QPCR identified two of these putative homozygous deletions (sample 3, chr6:33134224–33143208 and sample 48, chr12:18693321–18703198) but the Ct values were consistent with heterozygous deletions rather than homozygous (ΔΔCt>0.25) and the remaining putative homozygous deletion (sample 45, chr11:41755318–41788391) was not detected by qPCR.
Placental samples with identified CNVs were also interrogated employing the Illumina 1M-Duo Beadchip in order to assess reproducibility of results using similar platforms, with the latter being of markedly higher density. The Illumina 1M-Duo Beadchip (1.2M) differs from the CNV370-Duo Beadchip (370K) only in the number of SNPs interrogated by the array. Of the 17 samples in which CNVs were identified using the CNV370-Duo, 6 samples (representing 9 identified CNVs) were also interrogated using the 1M-Duo Beadchip. The CNV370-Duo and 1M-Duo Beadchips agreed on 6 of the 9 (67%) CNV calls. In a comparison of only those CNVs also confirmed by qPCR, the CNV370-Duo and 1M-Duo Beadchips agreed on 5 the 6 (83%) confirmed CNVs. When we limited this analysis only to larger CNVs >50 Kb, 6 of 7 (87%) were detected.
Finally, placental samples were also interrogated using the less dense Agilent 105K Comparative Genomic Hybridization (CGH) Array. The Agilent 105K CGH array differs from the CNV370-Duo Beadchip in both the number and genomic location of probes, as well as basic analysis principles. The Agilent aCGH platforms rely on a comparative hybridization of a differentially labeled interrogated DNA sample and a reference DNA sample, with copy number aberrations derived from the comparative strength of the hybridization signal from each call sample. The Illumina platforms, however, rely on a single sample hybridization, whereby copy number information is derived from combined analysis of absence of SNP heterozygosity and signal strength data at various oligonucleotides in aggregation. Of the 17 samples in which CNVs were identified using the CNV370-Duo Beadchip, 6 samples, representing 12 identified CNVs, were also interrogated using the Agilent 105K CGH Array. A total of 3 samples, representing 4 identified CNVs, were interrogated using all three array platforms. All 4 CNVs interrogated by all three platforms were confirmed by qPCR. As reported in Tables III and andIV,IV, the vast majority of CNVs detected on either Illumina platform but not the Agilent platform were due to lack of probe coverage on the latter, less dense array. When we again limit our analysis to CNV deletions 100% of the cases would demonstrate a deletion regardless of the platform employed (Table III). However, this was not true for duplications (Table IV).
Interestingly and as shown in Table V, in cases 2 and 22 all duplications and deletions detected by the Illumina CNV370Duo were detected by qPCR of both placental and fetal tissue (6/6; 100%) but only the deletions were detected by the Agilient 105K array (Table V).
Deletions affecting genes are most likely to play a causative role in stillbirth. A total of 4 deletions confirmed by qPCR on the placental and fetal samples, with the exception of case 48 for which there was no fetal sample, affected genes (Table VI).
By far the most compelling evidence for an identifiable deletion with a potential causative role for the stillbirth was observed for case 4 in which 25 contiguous genes were deleted (Table VI).
Other less notable and less expansive deletion occurred. For example in case 3 this involved deletion of the HLA-DPA1 3′ UTR. The potential mechanism from this deletion resulting in dysregulation of gene expression might arise from the concomitant observations that conserved miRNAs have been identified that target HLA-DPA1 (MicroCosm Targets http://www.ebi.ac.uk/enright-srv/microcosm/htdocs/targets/v5/), and miRNA regulation occurs primarily through binding to the 3′ UTR. The genes deleted in case 48 are all olfactory receptor encoding genes and as such seem unlikely to play a causative role even though deletions of olfactory receptors have also been identified in spontaneous abortions between 10–20 weeks (15). As a third example, a deletion occurring in an intron of CADM2 in case 33 could play a causative role as such a deletion can result in the gene being incorrectly spliced. An additionally likely possibility for cases 33 and 48 are that these variants are in fact not de novo but may arise from the maternal germline. However, the most compelling case for a likely de novo and contributing variant to stillbirth was observed in case 4.
The available clinical data on case 4 indicated that there were no fetal abnormalities or maternal comorbidities, and prior karyotype was reported as normal. Array based copy number analysis resulted in the identification of a 2.89 Mbp heterozygous deletion on chromosome 9 in this case (Figure 1), which was detected by all array platforms (Tables III, ,IV,IV, ,VI)VI) and confirmed by qPCR assays on the placenta and fetal tissue, with the high qPCR ΔΔCt in the fetal blood sample most likely resulting from sample contamination by maternal blood (Figure 1B). The deletion resulted in the loss of 25 RefSeq Genes (Table 3 and Figure 1C). The large size and high gene content of this deletion strongly suggests that it likely played a causative role in the stillbirth.
We have demonstrated for the first time that genome-wide array-based profiling of placentas from well-characterized, systematically classified unexplained stillbirths enables the detection of CNVs that were not observed by cytogenetic methods. These findings suggest that small genomic imbalances might contribute to a small proportion of non-anamalous, unexplained stillbirths. While in this study the likelihood of truly significant imbalances was low, it was not negligible (Table IV) and occurred in a carefully controlled research setting where karyotypes had been already performed.
Previous studies to identify genomic abnormalities in stillbirth samples have primarily employed G-banded karyotyping and fluorescence in situ hybridization (FISH). Using these methods, cytogenetic abnormalities have been described in a small percentage (5–13%) of all stillbirths (13, 17, 18). While in two studies this number does not demonstrably increase when fetal anatomic anomalies are present (17, 18), in a 3rd study chromosomal abnormalities were significantly higher (p<0.001) and observed in up to 38% of dysmorphic stillbirths (13). Karyotyping and FISH are generally accompanied by a high failure rate due to the inability to culture cells from fetal demise tissue samples (16%; 17). Presently, two types of arrays are available for genome-wide analysis of copy number imbalances: single nucleotide polymorphism (SNP) based arrays and array comparative genomic hybridization without SNP analysis.
CNVs detectable by arrays can be quite large as demonstrated by the 2.89 Mbp deletion in case 4, and yet would still be under the threshold of detection by broadly utilized traditional cytogenetic methods. Thus, the higher resolution provided by arrays suggests that they can be used as an alternative or at least a complementary method to standard karyotype analysis and FISH for examining genomic abnormalities in unexplained stillbirths. While the cost of clinical array-based copy number analysis is still higher than that of standard karyotyping, it is expected to continue to decrease as development will improve automation to facilitate high-throughput sample processing, data analysis, and reliable result reporting (22). We propose that the nominal difference in cost between the two approaches is well balanced by the higher resolution of genomic abnormalities in stillbirths provided by an array. In addition, besides its capacity to detect clinically significant copy number changes, array based analysis also avoids the low rate of detection of cytogenetic abnormalities resulting from failure of fetal or placental cell growth in vitro (17). A third advantage is that results can be more rapidly available because there is no need for precedent cell culture (23).
There are some inherent limitations to our study. First, we did not have uniformly available maternal DNA for the establishment of the de novo status of described CNVs. Moreover, in no instance was paternal DNA available. Optimal characterization of de novo status and likely contributory role of detected CNVs would require comparison to both parental genomes. Second, the samples were not initially collected with the purpose of array-based copy number analysis in mind, and thus in not all cases were we able to obtain stringent quality fetal DNA needed for confirmation. We believe that protocols optimized for maintenance of DNA integrity (including immediate DNA extraction) could have potentially allowed for a higher detection rate as well as greater ability for confirmation. Third, there are a number of inherent limitations to the array based technology itself. For example, CGH arrays cannot detect balanced rearrangements and are unreliable in detection of triploidy (24, 25). That said, de novo truly balanced rearrangements are unlikely to be a cause of stillbirth.
Despite these limitations, the strengths of our study are several-fold. First, we employed a number of practical applications in our protocol design. For example, we chose to use arrays to interrogate placental samples in order to examine the applicability of using placental DNA as a substitute for assaying DNA from fetal samples for a number of clinically practical reasons. In some cases direct fetal DNA assay is possible because fascia lata (or other fetal tissue) can be acceptably obtained. However, in other instances the parents may not be comfortable with consenting for direct fetal tissue sampling because they have concerns about tissue extraction on the fetus (14, 26). In other instances, extreme maceration in cases of longer interval fetal demise samples may necessitate utilization of placental tissue as the only viable alternative (27).
A second strength to our study was the strong concordance between results from placental and fetal DNA with respect to detection of deletions. While this was not true for detection of CNV duplications, there are likely several reasons. First, relatively small (<1 Mb) duplications are of unlikely significance in our sample set since the overwhelming majority detected in placenta were also detected in maternal blood but not in the fetus (Table IV). Second, of the four cases where CNV deletions (cases 4 and 48) and duplications (cases 37 and 48) were detected by multiple modalities (Illumina and/or Agilent platforms and/or qPCR) we observed concordance between placenta and fetus. The exception to this was case 37, which could not be tested secondary to the unavailability of fetal DNA.
Alternatively, CNVs detected solely in placental DNA and not in fetal tissue DNA may also be of interest if they represent genomic abnormities that are deleterious to the placenta, thereby contributing to the stillbirth. Chromosomal aneuploidy differences between placental and fetal genomes, known as confined placental mosaicism, have been observed by others to be associated with stillbirths (28). Thus, it seems possible that smaller scale CNV differences between the placenta and fetus may play a similar role and it would be of great value to interrogate both placenta and fetal tissue in future analyses of normal pregnancies and those complicated by stillbirth or other prenatal complications. A potential pertinent example of this is with case 37 which demonstrated a duplication of chr9 at the PAPP-A region. Lower circulating levels in PAPP-A on maternal serum screening have been associated with miscarriage (30); CNV duplications may result in either increased or decreased expression depending upon the genomic disruption (1).
A third strength of our study is our utilization of multiple measures for detection of structural variants. qPCR detection of putative CNVs can be time consuming and costly and does not necessarily imply “validation”. Alternately, confirmation using a second array platform of similar or higher resolution (i.e., the Illumina 1M but not Agilent 105K) might offer an alternative to qPCR especially in cases where multiple CNVs must be validated in a single sample because those CNVs can be validated by a single array platform rather than multiple qPCR assays. Warren et al. (20) initially identified CNVs in spontaneous abortions between 10–20 weeks using a low resolution BAC based array and then validated identified CNVs using a higher resolution CGH array. We initially identified CNVs using an array with 370K SNP probes and attempted to additionally detect those found using a lower resolution 105K CGH array. The low validation rate by the second array of 30% for CNVs also confirmed by qPCR is largely due to the lower resolution of the second array compared to the first. Validation by a second array platform would be an alternative to qPCR if the resolution is similar or greater on the second array. This will become more important as higher resolution arrays are used, with an inherently greater chance of detecting more than one CNV. Furthermore, it would allow for detection of additional CNVs that may be found by one type of array technology, but not by another (i.e., SNP based versus aCGH). This would be rare for large deletions and duplications (as illustrated by Case 4) but could be important for small yet clinically significant copy number deletions.
A fourth strength of this study is its robust characterization of stillbirths as “unexplained”, with further sub-classification by virtue of gestational age and gender-adjusted fetal growth. While others (19) interrogated 15 anomalous and non-anomalous stillbirths with aCGH, this is the first such report of determination of structural genetic variation in systematically acquired and well-characterized unexplained fetal deaths. By benefitting from the robust experience of a panel of experts in the characterization of these stillbirths, we avoided potential bias by not including anomalous stillbirths in the W2 criterion.
A fifth strength of our study is its potential for broader clinical implications. Unexplained stillbirths are clinically problematic, as they comprise the vast majority of stillbirths and occur most notably in the third trimester. The ability to utilize newer genome-wide technologies in a relatively inexpensive manner in this population would constitute an efficacious use of available resources to obtain a definitive diagnosis and thereby minimize parental stress and grief. For example, consider a hypothetical couple who present with a 39 week non-anomalous, unexplained stillbirth. In the course of evaluation, an array was completed which demonstrated a novel or known clinically significant 2.9 Mbp deletion which encompassed 34 genes. Both parents are then interrogated, found not to contain the structural variant, thus confirming an association between the fetal demise and a spontaneous de novo event. In any future pregnancies, the couple could choose to employ CVS or amniocentesis for targeted detection of the same structural variant. If the structural variant is not present,, then expensive other antenatal testing (i.e., once or twice weekly non-stress tests through the third trimester) might be further avoided. Larger multicenter studies will likely be needed to confirm that such an approach will be cost-effective and clinically beneficial.
Finally, such larger studies on stillbirth samples may reveal new copy number imbalances, associated with an important proportion of cases that have not previously been ascertained because we lacked the technologies to do so. Such discoveries have clearly resulted from postnatal application of array-based copy number analysis in children with congenital anomalies, dysmorphic features, and developmental delays. This will not only result in better understanding of causes of stillbirth, but also in improved diagnosis (ACOG Committee Opinion 440, 2009). Array based copy number analysis has now replaced karyotyping as the first-line diagnostic test in pediatric genetic analyses. Similar developments are predicted for prenatal diagnosis and the evaluation of fetal demise.
In sum, we have employed multiple modalities in detecting and analyzing copy number variants in unexplained stillbirths. We have observed that in at least one instance, this yielded a likely deletion potentially causative of stillbirth which was not detected with cytogenetic analysis. While the overall yield may be conceived of as low and is most likely to be of significance when CNV deletions are detected, it bears mention that under conventional protocols such detection would avoid the emotionally and clinically troublesome classification of this fetal death as having been “unknown”. Further studies with protocols optimized for such advancing genetic technologies are essential in this common perinatal disorder.
This work was funded by the NIH New Innovator Award to K.A.T. (DP2OD001500-01).