A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders

preprint OA: closed CC-BY-NC-ND-4.0
📄 Open PDF Full text JSON View at publisher
⚙ AI-generated deep summary by qwen3.7-flash, 2026-09-09 · read from full text ⓘ

This study developed and clinically validated a novel non-invasive prenatal testing method using ultra-sensitive amplicon-based sequencing of cell-free DNA to detect monogenic autosomal recessive disorders. The researchers applied this approach to a cohort of 57 patients with known fetal genotypes, achieving high sensitivity and specificity for identifying Sickle Cell Anaemia without requiring prior knowledge of the paternal genotype. The methodology utilizes internal controls and bias factor correction to calculate allelic imbalance probabilities based on fetal fraction and read counts. Relevance to endometriosis: The paper does not explicitly discuss endometriosis or adenomyosis; it was included in the corpus via a keyword match in the upstream search index.

Read from the paper's body, not the abstract. Not a substitute for reading the paper. No clinical advice. How this works

Abstract

Non-invasive prenatal testing (NIPT) to date is used in the clinic primarily to detect foetal aneuploidy. Few studies so far have focused on the detection of monogenic autosomal recessive disorders where mother and foetus carry the same mutation. In particular, NIPT is currently not available for the detection of Sickle Cell Anaemia (SCA), the most common monogenic disorder world-wide and the most common indication for invasive prenatal testing in high-income countries. Here, we report the clinical validation of a novel diagnostic approach that combines ultra-sensitive amplicon-based sequencing of cell-free DNA (cfDNA) with internal controls and bias factor correction to calculate the probability for the presence of allelic imbalance from maternal plasma without prior knowledge of the paternal genotype. Identification of the foetal genotype was determined using a hierarchical probabilistic model based on the relative number of reads from the sequencing, along with the foetal fraction. NIPT was performed on a cohort of 57 patients, all of whom had previously undergone invasive prenatal testing so that the foetal genotype was known. Overall, NIPT demonstrated 100% sensitivity and negative predictive value for foetal fractions higher than 0.5%, and 100% specificity and positive predictive value for foetal fractions higher than or equal to 4%. Our methodology can be used as a safe, non-invasive screening tool in any clinical scenarios where early prenatal diagnosis of SCA or other recessive disorders is important.
Full text 31,901 characters · extracted from preprint-html · click to expand
A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders Anthony Cutts , Dimitrios V. Vavoulis , Mary Petrou , Frances Smith , Barnaby Clark , Shirley Henderson , Anna Schuh doi: https://doi.org/10.1101/635342 Anthony Cutts 1 Molecular Diagnostic Centre, Department of Oncology, University of Oxford , Oxford, UK 2 NIHR Oxford Biomedical Research Centre, University of Oxford , Oxford, UK 3 Department of Oncology, University of Oxford , Oxford, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Dimitrios V. Vavoulis 1 Molecular Diagnostic Centre, Department of Oncology, University of Oxford , Oxford, UK 2 NIHR Oxford Biomedical Research Centre, University of Oxford , Oxford, UK 3 Department of Oncology, University of Oxford , Oxford, UK 4 Wellcome Centre for Human Genetics, University of Oxford , Oxford, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Mary Petrou 6 Institute of Women’s Health, University College London , London, UK . 7 Haemoglobinopathy Genetics Service, University College London Hospitals NHS Trust , London, UK . Find this author on Google Scholar Find this author on PubMed Search for this author on this site Frances Smith 8 Molecular Diagnostics, Viapath LLP at Kings College Hospital , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Barnaby Clark 9 Red Cell Biology, King’s College London , London, UK 10 Precision Medicine, King’s College Hospital , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shirley Henderson 11 Genomics England , London, UK . 12 NHS England , London, UK . Find this author on Google Scholar Find this author on PubMed Search for this author on this site Anna Schuh 1 Molecular Diagnostic Centre, Department of Oncology, University of Oxford , Oxford, UK 2 NIHR Oxford Biomedical Research Centre, University of Oxford , Oxford, UK 3 Department of Oncology, University of Oxford , Oxford, UK 5 Department of Hematology, Oxford University Hospitals NHS Trust , Oxford, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: anna.schuh{at}oncology.ox.ac.uk Abstract Full Text Info/History Metrics Preview PDF Abstract Non-invasive prenatal testing (NIPT) to date is used in the clinic primarily to detect foetal aneuploidy. Few studies so far have focused on the detection of monogenic autosomal recessive disorders where mother and foetus carry the same mutation. In particular, NIPT is currently not available for the detection of Sickle Cell Anaemia (SCA), the most common monogenic disorder world-wide and the most common indication for invasive prenatal testing in high-income countries. Here, we report the clinical validation of a novel diagnostic approach that combines ultra-sensitive amplicon-based sequencing of cell-free DNA (cfDNA) with internal controls and bias factor correction to calculate the probability for the presence of allelic imbalance from maternal plasma without prior knowledge of the paternal genotype. Identification of the foetal genotype was determined using a hierarchical probabilistic model based on the relative number of reads from the sequencing, along with the foetal fraction. NIPT was performed on a cohort of 57 patients, all of whom had previously undergone invasive prenatal testing so that the foetal genotype was known. Overall, NIPT demonstrated 100% sensitivity and negative predictive value for foetal fractions higher than 0.5%, and 100% specificity and positive predictive value for foetal fractions higher than or equal to 4%. Our methodology can be used as a safe, non-invasive screening tool in any clinical scenarios where early prenatal diagnosis of SCA or other recessive disorders is important. Main Since the discovery of circulating cell-free DNA (cfDNA) in the plasma of pregnant women in the late 1990s 1 , 2 , its potential for prenatal diagnosis has been the focus of intensive technological innovation. Screening for chromosomal abnormalities is now introduced in several countries including the UK 3 . However, the so-called ‘combined test’ (which utilises ultrasound scanning to measure foetal nuchal translucency, maternal age and blood tests to measure pregnancy-associated plasma protein A and free beta-human chorionic gonadotrophin) remains the first-stage test. Women at high-risk of carrying an affected baby are offered NIPT and, if this returns a positive result, they are given the option of confirmatory invasive testing by amniocentesis. Technologies detecting dominant de-novo mutations or dominant mutations in the father from cfDNA have been previously described 4 . For recessive disorders characterised by compound heterozygosity (e.g. cystic fibrosis), the same technology can be applied if the mutations in both parents are different. Where mutations are the same, methods that utilise the linkage of disease-causing mutations to paternal SNPs have been proposed 5 , 6 . However, this approach requires testing a previously born child for linkage to the paternal SNPs and it is subject to errors due to recombination events. An alternative approach avoiding the need for complex family work-up, would be to measure small allelic imbalances caused by foetal DNA in the maternal circulation. However, this requires precise allele quantification and, so far, rapid NIPT from cfDNA for conditions such as SCA, where the mother and father carry the same mutation, remains elusive. SCA is an autosomal recessive disease characterised by a single base-pair substitution in the beta globin gene. Due to the protective effect of the mutation against malaria, carrier frequencies in sub-Saharan Africa are 20% or higher. Over 224,200 infants are born annually with SCA worldwide 7 , including at least 1,000 in the US, making SCA by far the most common monogenic disease indication for invasive prenatal testing (IPT) in high-income countries. Invasive prenatal diagnosis by amniocentesis or chorionic villus sampling is costly and carries a 1-2% risk of miscarriage. Only one study so far has reported the use of cfDNA for the detection of allelic imbalance in SCA using digital droplet-PCR 8 and, to our knowledge, no clinical services currently offer the test worldwide. Here, we describe the development and validation of a highly sensitive and specific next-generation sequencing (NGS) approach for the non-invasive prenatal diagnosis (NIPD) of SCA. Our methodology does not require the paternal genotype and combines optimised PCR of the affected locus (including sequencing and PCR error correction) with precise estimation of the foetal cell fraction and different internal controls ( Figure 1 ). The raw assay data are analysed using a bespoke statistical methodology ( Figure 2 ), which estimates the respective foetal disease status. We calibrated and, subsequently, tested our method on 57 patients, achieving very high sensitivity (94%) and specificity (88%) ( Figure 3 and Table 1 ). The method is also applicable to other autosomal recessive and autosomal dominant diseases. View this table: View inline View popup Download powerpoint Table 1. Training and test cohorts used in this study and predictions from application of the NIPT. False positives are indicated in yellow. The false negative is indicated in red. Download figure Open in new tab Figure 1: Overview of the next-generation sequencing. Foetal fraction determination was based on methylation status described previously 9 . Library preparation was performed in triplicate, with an initial PCR using standard primers, followed by a second amplification using longer primers containing adapters for cluster generation and indexes so that different patient samples could be multiplexed. Download figure Open in new tab Figure 2: Overview of the statistical analysis. A) The total cfDNA is modelled as a mixture of maternal ϕ i and foetal ψ i components with proportions determined by the foetal fraction w . B) Estimated maternal and foetal expected fractions of mutated reads for each amplicon in a random subject. C) The estimated probability density functions summarising the expected maternal and foetal fractions in B. The stronger the shift E [ P ] of the foetal component to the right of the maternal component, the more likely it is that the foetus has the disease. D) Overview of model training. We identified an optimal E [ P ] threshold equal to 0.62. The variance of the various performance metrics was estimated using the bootstrap. E) Overview of applying the calibrated model on both the training and test cohorts. F, G, H) Overview of model performance with decreasing foetal fraction. A false negative arises at a foetal fraction of 0.5%, while false positives arise at foetal fractions less than 4%. Four samples (in green) with E [ P ] scores close to the threshold are also characterised by small foetal fractions (≤0.9%) Download figure Open in new tab Figure 3: Overview of training and test cohorts used in this study. A total of 57 subjects were recruited, all of which had previously undergone prenatal invasive testing. 29 subjects were used for calibrating the method, i.e. for finding an optimal threshold for the score E [ P ] (training phase). The remaining 28 subjects were used for evaluating the predictive capacity of the calibrated model (test phase). Overall, the method returned 5 false positives and a single false negative at foetal fractions ≤0.5%. Foetal fractions were determined using RT-PCR 9 , followed by library preparation and NGS. Identifying low-level variants using NGS is challenging due to errors during library preparation and sequencing. Additionally, library preparation of highly fragmented cfDNA is more difficult. To overcome these issues, we employed a PCR-based method that generates small amplicons for triplicate analysis to improve noise reduction and error correction. For the library preparation, the first step was a short 17-cycle PCR with short primers specifically designed to preserve the wild-type-to-sickle allele ratio of the template DNA, thus minimizing amplification bias. A second PCR followed with considerably longer primers containing the indices, primer binding sites and adapters required for sequencing. Paired-end ultra-deep sequencing was carried out using the MiSeq 300 v2 chemistry. Levels of the HbS allele were derived from FASTQ files using a Perl script. Up to 8 patient samples were sequenced per run, alongside three controls – with genotypes AA (normal), AS (carrier) and SS (disease) – for normalizing sample data. NGS is complemented by a bespoke statistical model taking as input the estimated foetal fraction w , the number of reads harbouring the HbS allele and the total number of reads covering the locus of the mutation in each amplicon in both the mother and a non-pregnant control with the same genotype as the mother. For each amplicon, we model the expected fraction of reads harbouring the HbS allele in the total cfDNA isolated from the mother as a mixture of foetal ( ψ i ) and maternal ( ϕ i ) components in proportions determined by the foetal fraction w ( Figure 2A ). The control shares the same expected fraction of mutated reads ϕ i with the mother, which helps increase the precision of estimates. For each case, we estimate the foetal ( ψ 1 , ψ 2 , ψ 3 ) and maternal ( ϕ 1 , ϕ 2 , ϕ 3 ) expected fractions of mutated reads per amplicon ( Figure 2B ) and probability density functions summarising each triplet of expected fractions ( Figure 2C ). For a carrier mother (AS), the density function for the maternal mutated read fractions is concentrated close to 50%, while for a foetus with increased mutated read fractions, the corresponding function is shifted to the right ( Figure 2C ). The expected magnitude of this shift, E [ P ], ranges between 0 and 1 and it predicts the foetus as HbSS (homozygote), if sufficiently high. In order to determine an optimal threshold value for E [ P ], we recruited 29 subjects with known foetal disease status from IPT. We applied NIPT and calculated an E [ P ] value for each. Then, for E [ P ] threshold values between 0 and 1, we predicted the foetal disease status in each case and calculated the sensitivity (True Positive Rate or TPR), specificity (True Negative Rate or TNR), and Matthews Correlation Coefficient (MCC), a balanced metric for measuring the performance of binary classifiers ( Figure 2D ). We estimated the variance of all three metrics using the bootstrap (Online Methods). Increasing values of the E [ P ] threshold had opposite effect on sensitivity and specificity, reaching an optimal value of 62% at which MCC was maximised ( Figure 2D ). A second MCC maximum at 77% associated with lower sensitivity (higher false negative rate) was ignored. At the optimal E [ P ] threshold, our NIPT correctly identified 8 true positives, 17 true negatives, 1 false negative and 3 false positives ( Table 1 ; Figure 3 ) achieving 89% sensitivity, 85% specificity, 73% positive predictive value (PPV) and 94% negative predictive value (NPV). In order to assess the predictive capacity of our NIPT, we recruited 28 additional subjects, which had also undergone IPT, and we applied our method on each using the previously determined optimal E [ P ] threshold of 62%. The test returned 7 true positives, 19 true negatives, 2 false positives and no false negatives, achieving 100% sensitivity and NPV, 91% specificity and 78% PPV ( Table 1 ; Figure 3 ). On all 57 subjects, the proposed NIPT called 5 false positives and 1 false negative ( Figure 2E ) corresponding to 94% sensitivity, 88% specificity, 75% PPV and 98% NPV. The single false negative can be attributed to the very low foetal fraction (0.5%) of the corresponding case ( Figure 2F ). Furthermore, four samples close to the E [ P ] threshold (green dots in Figures 2E-G ) are also associated with very low (≤0.9%) foetal fractions ( Figure 2G ). Overall, the NIPT demonstrates 100% sensitivity and NPV for foetal fractions higher than 0.5%, and 100% specificity and PPV for foetal fractions higher than or equal to 4% ( Figure 2H ). These results suggest that our methodology can be used as a non-invasive screening tool to exclude the presence of an affected baby. We expect that our test will be implemented similarly to non-invasive aneuploidy diagnosis for confidently excluding an affected pregnancy, thereby considerably reducing the number of IPT performed to confirmatory testing of positive results only. This will particularly also apply to countries with high incidence of sickle cell disease and a demand for prenatal diagnosis (i.e. Nigeria) where access to invasive testing is limited due to high cost and relatively low numbers of trained obstetricians. By comparison, phlebotomy is relatively inexpensive and readily available. This methodology could potentially make NIPT available to a larger number of beneficiaries. We believe that the future adoption of this NIPT will also have utility in other clinical scenarios: Treatment options for children and young adults who develop early complications from SCA include bone marrow or cord blood transplantation. Several clinical trials using autologous gene-editing of stem cells also reported successful outcomes 10 , 11 . Knowledge of the foetal genotype before birth would allow advance preparation for umbilical cord blood stem cell sampling for future cellular therapy, without placing the foetus at unnecessary risk through invasive testing. Neonatal screening programmes around the world face logistical challenges because of loss to follow-up of a significant number of babies after birth due to delays in obtaining post-natal screening test results. Knowledge of the foetal genotype prior to birth would allow targeting follow-up efforts to affected pregnancies only. Importantly, our method does not require prior knowledge of the father’s genotype, which will allow its straightforward incorporation into routine antenatal care with a rapid turnaround time. The test could easily be adapted for detection of dominant monogenic disorders, de novo mutations and other autosomal recessive conditions. Online Methods Ethics statement Mothers undergoing invasive prenatal diagnosis for SCA were consented for blood sampling in accordance with the Helsinki Declaration for service evaluation. Samples A total of 57 samples were collected from 2012 to 2014, where the gestational ages were between 8 and 17 weeks. Blood was collected in either EDTA or Streck tubes, and processed within 6 or 24 hours of collection, respectively. Blood samples were initially centrifuged at 3000 or 3400 rpm for 10 min at room temperature to separate the plasma. The plasma supernatant was transferred to a fresh tube and micro-centrifuged at 7000 or 14000 rpm for 10 min at room temperature prior to storage at −80°C. cfDNA was extracted using the QIAamp Circulating Nucleic Acid Kit according to the manufacturer’s instructions and eluted into a final volume of 70µl. Foetal fraction determination and library preparation Initially, a Real-Time Polymerase Chain Reaction (RT-PCR) assay was carried out to determine the foetal fraction in the total cfDNA present. This was achieved by assessing the methylation status of the RASSF1A promoter, which is a universal foetal DNA marker, as described previously 9 . For the library preparations, an amplicon-based approach was used where cfDNA samples were amplified by 3 different primer sets in singleplex PCRs. Amplification was carried out in two stages. For the initial PCR, there was 10µl input DNA per amplicon and upon completion of this step, a 2.5µl aliquot was taken and a second PCR performed to attach the appropriate adaptors, primer binding sites and barcodes for sequencing. For both amplifications, primers at a final concentration of 240nmol/L and 0.1 U/µl of Pwo DNA polymerase (Roche) were used with 2X QIAGEN Multiplex PCR Master Mix. The cycling conditions for both PCRs were as follows: 94°C for 3 minutes, then 17 cycles (primary PCR); 25 cycles for amplicons 1 and 2 or 22 cycles for amplicon 3 at 94°C for 45 seconds, 56°C for 45 seconds and 72°C for 1 minute followed by 72°C for 10 minutes (secondary PCR). Following the second PCR, a clean-up was performed and the samples quantified for sequencing using the Qubit Fluorometer. Paired-end ultra-deep sequencing was carried out on a MiSeq, where 84 bases were sequenced from both ends of the DNA fragments. Bioinformatics From the generated FASTQ files, reads above Q30 were processed with a custom Perl script. Using a 2bp amplicon ID tag, the reads were divided into 3 amplicons, based on the following rules: the 84bp sequenced amplicon 1 starts with the bases AC; amplicon 2 ends with the bases TC; and amplicon 3 ends with the bases AC. Using the script, both reads were interrogated to find variants in the DNA fragment at the HbS mutation site (chr11:5,248,232). Statistical analysis Statistical analysis was performed in R 12 and Stan 13 using a bespoke statistical model. The model takes as input a) an estimate of the foetal fraction w in the total cfDNA, b) the total number of reads R i covering the HbS locus and the number of reads r i carrying the HbS allele for each amplicon i , as determined from sequencing the total cfDNA extracted from the mother’s blood and c) the total number of reads S i covering the HbS locus and the number of reads s i carrying the HbS allele for each amplicon i , as determined from sequencing the total cfDNA extracted from the blood of a non-pregnant control subject with the same genotype as the mother (usually a carrier). Inference in the model aims at estimating the expected fractions of foetal- and maternal-specific DNA harbouring the HbS allele from which we predict the disease status of the foetus. We assume that the read counts r i and s i follow Binomial distributions with parameters θ i and ϕ i , the expected fractions of reads harbouring the HbS mutation in the mother and control, respectively: The fundamental intuition in the above model is to express θ i as a mixture of maternal ϕ i and foetal ψ i fractions in proportions determined by the fraction w of foetal DNA in the mother’s blood, as follows: In the above mixture, the expected fraction ϕ i is shared with the control and it can be estimated with relatively high confidence. We complete the model by imposing Beta priors (parametrized by mean and variance) on ϕ i and ψ i , as shown below: where μ ϕ ~ Uniform (0, 1), υ ϕ ~ Uniform (0, 1), and . Analogous expressions hold for parameters μ ψ , υ ψ and . The above model was encoded in Stan and estimated using an adaptive Hamiltonian Monte Carlo algorithm with default parameters 14 . Having estimated the foetal and maternal fractions of reads for each amplicon i , we need to decide whether the foetus is homozygous mutant (SS) and, therefore, the disease is present, or not (not SS). Assuming the mother is a carrier (i.e. she has genotype AS), we expect the fractions ϕ i to be on average close to 50%. If the foetus is SS, then the foetal fractions ψ i should be on average higher than ϕ i . We can quantify the difference between foetal and maternal expected read fractions by calculating the overlap between the above Beta distributions, as follows: This measures the average probability that the foetal fraction is higher than the maternal fraction. Since the parameters of the above Beta distributions are random variables themselves, P is also a random variable with expectation , where π ≡ π ( μ ϕ , υ ϕ , μ ψ , υ ψ ) stands for the posterior distribution of its arguments. Given the output of the MCMC sampler, this expectation can be readily calculated as follows: a) for the t -th posterior sample of the parameter vector ( μ ϕ , υ ϕ , μ ψ , υ ψ ), sample one or more values of ϕ from its generative Beta distribution and, for each such value, calculate the inner integral in the above equation as a regularised incomplete beta function (this can be done efficiently using the pbeta routine in R), b) take the average of this integral over all samples of ϕ . This constitutes a posterior sample of P , corresponding to the t -th posterior sample of the parameter vector. c) Repeat until all T posterior parameter samples are exhausted. The required expectation can be calculated as an average over all T posterior samples of P . Subsequently, we can say that the foetus has the disease (i.e. its genotype is SS), if the expectation of P exceeds a threshold c , i.e. E [ P ] > c . In order to find an optimal value of c , we utilise a training dataset with n=29 subjects and we identify the value of c in the range from 0 to 1 (in steps of 0.01), which maximizes the classification performance of the statistical model, as measured by the Matthews Correlation Coefficient (MCC). The variance of MCC was determined using 1000 bootstrap samples, each with size equal to the size of the training dataset ( Figure 2D ). The trained model was subsequently tested on a validation cohort with n=28 subjects. Footnotes ↵ * joint first co-authors References 1. ↵ Lo , Y. M. et al. Presence of fetal DNA in maternal plasma and serum . Lancet 350 , 485 – 487 ( 1997 ). OpenUrl CrossRef PubMed Web of Science 2. ↵ Lo , Y. M. et al. Quantitative analysis of fetal DNA in maternal plasma and serum: implications for noninvasive prenatal diagnosis . Am J Hum Genet 62 , 768 – 775 ( 1998 ). OpenUrl CrossRef PubMed Web of Science 3. ↵ Norton , M. E. et al. Cell-free DNA analysis for noninvasive examination of trisomy . N Engl J Med 372 , 1589 – 1597 ( 2015 ). OpenUrl CrossRef PubMed 4. ↵ Zhang , J. et al. Non-invasive prenatal sequencing for multiple Mendelian monogenic disorders using circulating cell-free fetal DNA . Nat Med ( 2019 ). 5. ↵ Papasavva , T. , Kalikas , I. , Kyrri , A. & Kleanthous , M. Arrayed primer extension for the noninvasive prenatal diagnosis of beta-thalassemia based on detection of single nucleotide polymorphisms . Ann. N. Y. Acad. Sci . 1137 , 302 – 308 ( 2008 ). OpenUrl CrossRef PubMed 6. ↵ Phylipsen , M. et al. Non-invasive prenatal diagnosis of beta-thalassemia and sickle-cell disease using pyrophosphorolysis-activated polymerization and melting curve analysis . Prenat. Diagn . 32 , 578 – 587 ( 2012 ). OpenUrl PubMed 7. ↵ March of Dimes global report on birth defects . Available at: http://www.marchofdimes.org/mission/march-of-dimes-global-report-on-birth-defects.aspx . (Accessed: 7th May 2019 ) 8. ↵ Barrett , A. N. , McDonnell , T. C. R. , Chan , K. C. A . & Chitty , L. S. Digital PCR analysis of maternal plasma for noninvasive detection of sickle cell anemia . Clin. Chem . 58 , 1026 – 1032 ( 2012 ). OpenUrl Abstract / FREE Full Text 9. ↵ White , H. E. , Dent , C. L. , Hall , V. J. , Crolla , J. A. & Chitty , L. S. Evaluation of a Novel Assay for Detection of the Fetal Marker RASSF1A: Facilitating Improved Diagnostic Reliability of Noninvasive Prenatal Diagnosis . PLOS ONE 7 , e45073 ( 2012 ). OpenUrl CrossRef PubMed 10. ↵ Ribeil , J.-A. et al. Gene Therapy in a Patient with Sickle Cell Disease . N Engl J Med 376 , 848 – 855 ( 2017 ). OpenUrl CrossRef PubMed 11. ↵ Nowogrodzki , A. Gene therapy targets sickle-cell disease . Nature 564 , S12 – S13 ( 2018 ). OpenUrl 12. ↵ R Core Team . R: A Language and Environment for Statistical Computing . ( R Foundation for Statistical Computing , 2019 ). 13. ↵ Carpenter , B. et al. Stan: A Probabilistic Programming Language . J. Stat. Softw. Artic . 76 , 1 – 32 ( 2017 ). OpenUrl 14. ↵ Hoffman , M. D. & Gelman , A. The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo . J Mach Learn Res 15 , 1593 – 1623 ( 2014 ). OpenUrl Back to top Previous Next Posted May 10, 2019. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders Anthony Cutts , Dimitrios V. Vavoulis , Mary Petrou , Frances Smith , Barnaby Clark , Shirley Henderson , Anna Schuh bioRxiv 635342; doi: https://doi.org/10.1101/635342 Share This Article: Copy Citation Tools A method for non-invasive prenatal diagnosis of monogenic autosomal recessive disorders Anthony Cutts , Dimitrios V. Vavoulis , Mary Petrou , Frances Smith , Barnaby Clark , Shirley Henderson , Anna Schuh bioRxiv 635342; doi: https://doi.org/10.1101/635342 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Genomics Subject Areas All Articles Animal Behavior and Cognition (7975) Biochemistry (18653) Bioengineering (14781) Bioinformatics (44202) Biophysics (22484) Cancer Biology (19609) Cell Biology (26767) Clinical Trials (138) Developmental Biology (13907) Ecology (20900) Epidemiology (2067) Evolutionary Biology (25338) Genetics (16108) Genomics (23413) Immunology (18625) Microbiology (42282) Molecular Biology (17958) Neuroscience (92996) Paleontology (694) Pathology (2973) Pharmacology and Toxicology (5068) Physiology (8074) Plant Biology (15917) Scientific Communication and Education (2093) Synthetic Biology (4538) Systems Biology (10192) Zoology (2376) window.__CF$cv$params={r:'a3881491aec590cf',t:'MTc4ODk3NjY5OQ==',u:'01a087521fb07522b44ff266cb8d2a2b',ut:'s5F8hgrCnlmDUrfA9vzs717aqsP5h0MJTR.pK4bYHks-1788976701-1.2.1.1-QYVp2VYLRPNYiVJWQop2VyIH5e2GEdOUVwpzR.7qhIWfXzNjlbD3RuhKkVENbITNRzG0iLy8TPH_ebA9rW3qenm5UltUAidKmiPb1n4dd1M',i:60};(function(){if(!document.body)return;var s=document.createElement('script');s.src='/cdn-cgi/challenge-platform/scripts/precursor/main.js';document.head.appendChild(s);})();

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

⚙ Ask this paper AI returns verbatim quotes from the full text · source: preprint-html ⓘ

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. The paper's references may be in our DB but unresolved to ``paper_id`` (resolution happens at ingest when the cited DOI matches a row we already have). Run the cross-source citation reconcile pass to retry.

Source provenance

europepmc
last seen: 2026-05-19T01:45:01.086888+00:00
unpaywall
last seen: 2026-05-29T02:00:03.542394+00:00
License: CC-BY-NC-ND-4.0