A Pilot Study of Assessing Whole Genome Sequencing in Newborn Screening in Unselected Children in China
Min Jian,Xiaohong Wang,Yuanyuan Sui,Mingyan Fang,Chenchen Feng,Yingping Huang,Chunhua Liu,Ruidong Guo,Yuanning Guan,Yuxiao Gao,Zhiwei Wang,Shuli Li,Bochen Cheng,Lina Sun,Fenghua Cui,Jia Guo,Ying Zhan,Guohong Zhang,Ling Zheng,Fengxia Su,Wei Xue,Puyi Qian,Shaobo Gao,Jiayu Chen,Lingyao Guan,Haorong Lu,Karsten Kristiansen,Xin Jin,Fang Chen,Yuhuan Zhao,Lennart Hammarstrom,Xiaojing Jiang,Junnian Liu,Ya Gao
DOI: https://doi.org/10.1002/ctm2.843
IF: 8.554
2022-01-01
Clinical and Translational Medicine
Abstract:We investigated whether screening by whole genome sequencing (WGS) in unselected newborns provides more information of potentially curable or treatable medical conditions than routine newborn screening (NBS). We demonstrated that compared with routine NBS, WGS produced fewer false positive results and identified more actionable pathogenic or likely pathogenic variants in the selective 246 genes. Previously, WGS has been used to identify mutated genes in newborn children with a suspected disease.1 However, sequencing of apparently healthy newborns has remained controversial due to technical concerns and ethical issues.2 In this study, 321 non-pre-selected newborns from a cohort of pregnant women in Qingdao, China were recruited (Table 1). DNA from 303 umbilical cord blood samples and 18 umbilical cords was extracted for 40X WGS. For data interpretation, we selected 251 genes associated with 59 Mendelian disorders, 164 primary immunodeficiency diseases (PIDs) and five pharmacogenetic (PGx) genes, following the guidelines by the Recommended Uniform Screening Panel (RUSP), the International Union of Immunologic Societies (IUIS) Expert Committee for Primary Immunodeficiency, the Dutch Pharmacogenetics Working Group (DPWG), and the Clinical Pharmacogenetics Implementation Consortium (CPIC).3-5 Sequencing protocol, data analysis pipeline, and criteria for sequence variants interpretation following the ACMG/AMP guidelines are described in the Supporting Information. The WGS results were compared with NBS results, including the mandatory checks of hearing impairment and four metabolic diseases, the metabolic testing of 48 inherited metabolic diseases (IMDs), and the genetic screening for 20 hearing loss loci incorporated into the local NBS program in China.6, 7 Among the analysed DNA samples of 321 newborns, the average sequencing depth was 47.42X (28.84X–82.90X) and the average coverage was 99.48% (99.01%–99.89%) (Figure 1A). For the 59 Mendelian disorders, a total of 131 pathogenic or likely pathogenic (P/LP) mutations and 5 pathogenic copy number variations were detected in 107 of the 321 newborns (33.33%), corresponding to 106 carriers of 28 diseases and 1 patient with phenylketonuria (PKU) (Figure 1B and Table 2). The 25.23% of newborns (n = 81) carried one P/LP mutations, and 7.17% and 0.93% of newborns (n = 23 and n = 3) carried two or three P/LP mutations, respectively. Hearing loss, methylmalonic acidemia (MMA), primary congenital hypothyroidism (CH), and PKU were diseases with the most carriers, while GJB2 (28/321, 8.72%), MMACHC (11/321, 3.43%), DUOX2 (10/321, 3.12%), PAH (8/321, 2.49%) and SLC26A4 (8/321, 2.49%) were the top five genes with the highest carrier frequencies of P/LP mutations (Table 2). For the 164 PIDs, 9 heterozygous P/LP variants in 6 genes were identified in 9 newborns (2.80%), all in a heterozygous state (Table 2). Four newborns were shown to carry heterozygous variant of unknown significance (VUS) in the gene SLC25A13 (c.2T>C, p.M1T), which was predicted as start loss and likely affecting the initiator methionine of the SLC25A13 mRNA. Two newborns carried a VUS in ASS1(c.-4C>T, p.?). Although these VUSs were not included in the final report to the participants, follow-up of the children with VUSs will be conducted till 3 years of age. Sanger sequencing confirmed 143 out of 145 mutations identified by WGS, resulting in an accuracy of 98.62%. Carriers of SMN1 mutations were validated by multiplex ligation-dependent probe amplification and real-time quantitative PCR, showing that five out of the six predicted carriers were true. Of the 321 newborns, 312 (97.20%) had the results of 48 IMDs screening and genetic hearing loss screening on 20 loci, which identified one newborn with PKU and one infant with increased blood level of isovalerylcarnitine (Table 1). In addition, 18 carriers harbouring 20 pathogenic mutations causing hearing impairment were detected by genetic hearting loss screening, albeit all 321 children passed the physical hearing screening at hospital (Supporting Information). The newborn WGS also identified the PKU case and 18 hearing loss carriers (Figure 2A). However, the child with increased level of isovalerylcarnitine was confirmed to be a carrier of 3-methylcrotonyl-CoA carboxylase deficiency by WGS. In addition, WGS identified two infants carrying compound heterozygous P/LP variants in GJB2 (Figure 3) and four children carrying pathogenic mutations in MT-RNR1 (c1095T > C), suggesting an increased risk of late-onset deafness or drug-induced hearing loss, respectively. Although currently non-symptomatic, the two newborns with GJB2 variants were scheduled to undergo hearing tests every 6 months, and the four newborns with the m.1095T mutation in MT-RNR1 were advised to avoid using aminoglycosides. Interestingly, we observed that 313 newborns (97.51%) carried at least one actionable PGx variant (Figure 1C). This result is in line with a European 44 000 biobank participants study, where 99.8% of the participants had a genotype associated with increased risks to at least one medication.8 Furthermore, we found three common PGx variants in the Qingdao cohort, CYP2D6*10 (48.60%), NUDT15*3 (13.08%) and UGT1A1*6 (21.18%) (Figure 2B and Table 3), that showed significant frequency differences as compared to East Asian populations (p < 0.05). An important aspect when screening for disorders in a given population is the use of a matched control database as variants can be highly specific for a given ethnic group.9 Most databases published to date are based on individuals of European descent and many populations have limited or poor representation. Limitations of the current study are the small sample size and restricted metabolic tests. A large-scale NBS effort is needed to validate our findings and fully investigate the treatable or curable medical conditions in newborns. The technical challenge of newborn WGS is to screen genes with high homology due to the misalignment of short-read sequencing. Therefore, the customized pipeline is needed to improve the accuracy and sensitivity of SNVs at genes with high-level homology. Albeit the present cost and turnaround time of WGS is several times more than the present NBS methods, in the forseeable future the pitfalls of WGS cost and turnaround time will likely facilitate the application of newborn WGS in NBS programs. In our study, selective identification of genomic data, where therapeutic options are available, did not violate the Wilson–Jungner criteria10. Our work provides a basis for future research on expanding screening genes and diseases in newborn screening program. Given adequate cost-effectiveness, WGS should be considered in future newborn screening programs. Further discussion of the interpretation accuracy and ethical use of genomic information needs to take place on a global scale. We appreciate the participation of the volunteers and their families. Without their support, this work would not have been possible. This work was also supported by the China National GeneBank (CNGB). This study was funded by the National Natural Science Foundation of China (No.31800765), the Shenzhen Municipal Government of China (JCY20170817145047361) and the Guangdong Provincial Key Laboratory of Genome Read and Write (No. 2017B030301011). The authors declare no conflict of interest. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.