Quantitative and qualitative mutational impact of ionizing radiation on normal cells.

The comprehensive genomic impact of ionizing radiation (IR), a carcinogen, on healthy somatic cells remains unclear. Using large-scale whole-genome sequencing (WGS) of clones expanded from irradiated murine and human single cells, we revealed that IR induces a characteristic spectrum of short insertions or deletions (indels) and structural variations (SVs), including balanced inversions, translocations, composite SVs (deletion-insertion, deletion-inversion, and deletion-translocation composites), and complex genomic rearrangements (CGRs), including chromoplexy, chromothripsis, and SV by breakage-fusion-bridge cycles. Our findings suggest that 1 Gy IR exposure causes an average of 2.33 mutational events per Gb genome, comprising 2.15 indels, 0.17 SVs, and 0.01 CGRs, despite a high level of inter-cellular stochasticity. The mutational burden was dependent on total irradiation dose, regardless of dose rate or cell type. The findings were further validated in IR-induced secondary cancers and single cells without clonalization. Overall, our study highlights a comprehensive and clear picture of IR effects on normal mammalian genomes.

Comprehensive characterization of maternal, fetal, and neonatal microbiomes supports prenatal colonization of the gastrointestinal tract.

In this study, we aimed to comprehensively characterize the microbiomes of various samples from pregnant women and their neonates, and to explore the similarities and associations between mother-neonate pairs, sample collection sites, and obstetrical factors. We collected samples from vaginal discharge and amniotic fluid in pregnant women and umbilical cord blood, gastric liquid, and meconium from neonates. We identified 19,597,239 bacterial sequences from 641 samples of 141 pregnant women and 178 neonates. By applying rigorous filtering criteria to remove contaminants, we found evidence of microbial colonization in traditionally considered sterile intrauterine environments and the fetal gastrointestinal track. The microbiome distribution was strongly grouped by sample collection site, rather than the mother-neonate pairs. The distinct bacterial composition in meconium, the first stool passed by newborns, supports that microbial colonization occurs during normal pregnancy. The microbiome in neonatal gastric liquid was similar, but not identical, to that in maternal amnionic fluid, as expected since fetuses swallow amnionic fluid in utero and their urine returns to the fluid under normal physiological conditions. Establishing a microbiome library from various samples formed only during pregnancy is crucial for understanding human development and identifying microbiome modifications in obstetrical complications.

ClinPharmSeq: A targeted sequencing panel for clinical pharmacogenetics implementation.

The accurate identification of genetic variants contributing to therapeutic drug response or adverse effects is the first step in implementation of precision drug therapy. Targeted sequencing has recently become a common methodology for large-scale studies of genetic variation thanks to its favorable balance between low cost, high throughput, and deep coverage. Here, we present ClinPharmSeq, a targeted sequencing panel of 59 genes with associations to pharmacogenetic (PGx) phenotypes, as a platform to explore the relationship between drug response and genetic variation, both common and rare. For validation, we sequenced DNA from 64 ethnically diverse Coriell samples with ClinPharmSeq to call star alleles (haplotype patterns) in 27 genes using the bioinformatics tool PyPGx. These reference samples were extensively characterized by multiple laboratories using PGx testing assays and, more recently, whole genome sequencing. We found that ClinPharmSeq can consistently generate deep-coverage data (mean = 274x) with high uniformity (30x or above = 94.8%). Our genotype analysis identified a total of 185 unique star alleles from sequencing data, and showed that diplotype calls from ClinPharmSeq are highly concordant with that from previous publications (97.6%) and whole genome sequencing (97.9%). Notably, all 19 star alleles with complex structural variation including gene deletions, duplications, and hybrids were recalled with 100% accuracy. Altogether, these results demonstrate that the ClinPharmSeq platform offers a feasible path for broad implementation of PGx testing and optimization of individual drug treatments.

Indigenous Ancestry and Admixture in the Uruguayan Population.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in "Charrúas" and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

Genomic and transcriptomic analyses reveal a tandem amplification unit of 11 genes and mutations in mismatch repair genes in methotrexate-resistant HT-29 cells.

DHFR gene amplification is commonly present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. In this study, we proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real-time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). We identified an amplification unit spanning 11 genes, from the DHFR gene to the ATP6AP1L gene position, with high adjusted interaction frequencies on chromosome 5 (~2.2 Mbp) and a twenty-fold tandemly amplified region, and novel inversions at the start and end positions of the amplified region as well as frameshift insertions in most of the MSH and MLH genes were detected. These mutations might stimulate chromosomal breakage and cause the dysregulation of mismatch repair. Characterizing the tandem gene-amplified unit may be critical for identifying the mechanisms that trigger genomic rearrangements. These findings may provide new insight into the mechanisms underlying the amplification process and the evolution of drug resistance.

Classification of High-Grade Serous Ovarian Carcinoma by Epithelial-to-Mesenchymal Transition Signature and Homologous Recombination Repair Genes.

High-grade serous ovarian cancer (HGSOC) is one of the deadliest cancers that can occur in women. This study aimed to investigate the molecular characteristics of HGSOC through integrative analysis of multi-omics data. We used fresh-frozen, chemotherapy-naïve primary ovarian cancer tissues and matched blood samples of HGSOC patients and conducted next-generation whole-exome sequencing (WES) and RNA sequencing (RNA-seq). Genomic and transcriptomic profiles were comprehensively compared between patients with germline BRCA1/2 mutations and others with wild-type BRCA1/2. HGSOC samples initially divided into two groups by the presence of germline BRCA1/2 mutations showed mutually exclusive somatic mutation patterns, yet the implementation of high-dimensional analysis of RNA-seq and application of epithelial-to-mesenchymal (EMT) index onto the HGSOC samples revealed that they can be divided into two subtypes; homologous recombination repair (HRR)-activated type and mesenchymal type. Patients with mesenchymal HGSOC, characterized by the activation of the EMT transcriptional program, low genomic alteration and diverse cell-type compositions, exhibited significantly worse overall survival than did those with HRR-activated HGSOC (p = 0.002). In validation with The Cancer Genome Atlas (TCGA) HGSOC data, patients with a high EMT index (≥the median) showed significantly worse overall survival than did those with a low EMT index (

Evaluation of low-pass genome sequencing in polygenic risk score calculation for Parkinson's disease.

BACKGROUND: Low-pass sequencing (LPS) has been extensively investigated for applicability to various genetic studies due to its advantages over genotype array data including cost-effectiveness. Predicting the risk of complex diseases such as Parkinson's disease (PD) using polygenic risk score (PRS) based on the genetic variations has shown decent prediction accuracy. Although ultra-LPS has been shown to be effective in PRS calculation, array data has been favored to the majority of PRS analysis, especially for PD. RESULTS: Using eight high-coverage WGS, we assessed imputation approaches for downsampled LPS data ranging from 0.5 × to 7.0 × . We demonstrated that uncertain genotype calls of LPS diminished imputation accuracy, and an imputation approach using genotype likelihoods was plausible for LPS. Additionally, comparing imputation accuracies between LPS and simulated array illustrated that LPS had higher accuracies particularly at rare frequencies. To evaluate ultra-low coverage data in PRS calculation for PD, we prepared low-coverage WGS and genotype array of 87 PD cases and 101 controls. Genotype imputation of array and downsampled LPS were conducted using a population-specific reference panel, and we calculated risk scores based on the PD-associated SNPs from an East Asian meta-GWAS. The PRS models discriminated cases and controls as previously reported when both LPS and genotype array were used. Also strong correlations in PRS models for PD between LPS and genotype array were discovered. CONCLUSIONS: Overall, this study highlights the potentials of LPS under 1.0 × followed by genotype imputation in PRS calculation and suggests LPS as attractive alternatives to genotype array in the area of precision medicine for PD.

Whole genome sequencing reveals a frameshift mutation and a large deletion in YY1AP1 in a girl with a panvascular artery disease.

BACKGROUND: Rare diseases are pathologies that affect less than 1 in 2000 people. They are difficult to diagnose due to their low frequency and their often highly heterogeneous symptoms. Rare diseases have in general a high impact on the quality of life and life expectancy of patients, which are in general children or young people. The advent of high-throughput sequencing techniques has improved diagnosis in several different areas, from pediatrics, achieving a diagnostic rate of 41% with whole genome sequencing (WGS) and 36% with whole exome sequencing, to neurology, achieving a diagnostic rate between 47 and 48.5% with WGS. This evidence has encouraged our group to pursue a molecular diagnosis using WGS for this and several other patients with rare diseases. RESULTS: We used whole genome sequencing to achieve a molecular diagnosis of a 7-year-old girl with a severe panvascular artery disease that remained for several years undiagnosed. We found a frameshift variant in one copy and a large deletion involving two exons in the other copy of a gene called YY1AP1. This gene is related to Grange syndrome, a recessive rare disease, whose symptoms include stenosis or occlusion of multiple arteries, congenital heart defects, brachydactyly, syndactyly, bone fragility, and learning disabilities. Bioinformatic analyses propose these mutations as the most likely cause of the disease, according to its frequency, in silico predictors, conservation analyses, and effect on the protein product. Additionally, we confirmed one mutation in each parent, supporting a compound heterozygous status in the child. CONCLUSIONS: In general, we think that this finding can contribute to the use of whole genome sequencing as a diagnosis tool of rare diseases, and in particular, it can enhance the set of known mutations associated with different diseases.

Direct conversion of adult human fibroblasts into functional endothelial cells using defined factors.

We aimed to directly convert adult human dermal fibroblasts (aHDFs) into functional endothelial cells (ECs). Lentiviral vectors encoding endothelial transcription factors (TFs) were constructed. We examined whether five TFs (FOXO1, ER71, KLF2, TAL1, and LMO2) used for the generation of mouse induced ECs (iECs) could convert the aHDFs into human iECs. Twenty-eight days after transduction with lentiviral constructs, 32.1 ± 5.1% cells expressed vascular endothelial (VE)-cadherin. Factor screening revealed that only three factors (3F: ER71, KLF2, and TAL1) were necessary to induce VE-cadherin (+) cells (49.4 ± 3.5%). However, whole transcriptome sequencing showed that VE-cadherin (+) cells were not completely reprogrammed. Mature iECs double-positive for VE-cadherin/Pecam1 (DP cells) with a cobblestone appearance were obtained at a frequency of only 5.1 ± 0.6%. Using whole transcriptome analysis, the potential factors that could block the conversion were screened. Among candidates TWIST1-knockdown enhanced efficiency of conversion. Rosiglitazone, an inhibitor of epithelial-mesenchymal transition (EMT), also improved the conversion efficiency. Moreover, a 2nd second-stage conversion process, in which VE-cadherin (+) cells were incubated for additional two weeks, further enhanced the efficiency. The final protocol for 6 weeks yielded a conversion rate of 19.6 ± 3.0% iECs, defined by DP cells depicting the nature of mature ECs in various analyses. Further analyses revealed that the genetic and epigenetic profiles of iECs resembled those of functional ECs. Collectively, aHDFs can be converted into functional ECs through the transduction of ER71, KLF2, and TAL1, combined with two EMT inhibitors (siTWIST1 and rosiglitazone), followed by 2nd stage conversion.

Alterations of transcriptome signatures in head trauma-related neurodegenerative disorders.

Chronic traumatic encephalopathy (CTE) is a neurodegenerative disease that is associated with repetitive traumatic brain injury (TBI). CTE is known to share similar neuropathological features with Alzheimer's disease (AD), but little is known about the molecular properties in CTE. To better understand the neuropathological mechanism of TBI-related disorders, we conducted transcriptome sequencing analysis of CTE including AD and CTE with AD (CTE/AD) post-mortem human brain samples. Through weighted gene co-expression network analysis (WGCNA) and principal component analysis (PCA), we characterized common and unique transcriptome signatures among CTE, CTE/AD, and AD. Interestingly, synapse signaling-associated gene signatures (such as synaptotagmins) were commonly down-regulated in CTE, CTE/AD, and AD. Quantitative real-time PCR (qPCR) and Western blot analyses confirmed that the levels of synaptotagmin 1 (SYT1) were markedly decreased in CTE and AD compared to normal. In addition, calcium/calmodulin-dependent protein kinase II (CaMKII), protein kinase A (PKA), protein kinase C (PKC), and AMPA receptor genes that play a pivotal role in memory function, were down-regulated in head trauma-related disorders. On the other hand, up-regulation of cell adhesion molecules (CAMs) associated genes was only found in CTE. Our results indicate that dysregulation of synaptic transmission- and memory function-related genes are closely linked to the pathology of head injury-related disorder and AD. Alteration of CAMs-related genes may be specific pathological markers for the CTE pathology.

The tumor immune microenvironmental analysis of 2,033 transcriptomes across 7 cancer types.

Understanding the tumor microenvironment is important to efficiently identify appropriate patients for immunotherapies in a variety of cancers. Here, we presented the tumor microenvironmental analysis of 2,033 cancer samples across 7 cancer types: colon adenocarcinoma, skin cutaneous melanoma, kidney renal papillary cell carcinoma, sarcoma, pancreatic adenocarcinoma, glioblastoma multiforme, and pheochromocytoma / paraganglioma from The Cancer Genome Atlas cohort. Unsupervised hierarchical clustering based on the gene expression profiles separated the cancer samples into two distinct clusters, and characterized those into immune-competent and immune-deficient subtypes using the estimated abundances of infiltrated immune and stromal cells. We demonstrated differential tumor microenvironmental characteristics of immune-competent subtypes across 7 cancer types, particularly immunosuppressive tumor microenvironment features in kidney renal papillary cell carcinoma with significant poorer survival rates and immune-supportive features in sarcoma and skin cutaneous melanoma. Additionally, differential genomic instability patterns between the subtypes were found across the cancer types, and discovered that immune-competent subtypes in most of cancer types had significantly higher immune checkpoint gene expressions. Overall, this study suggests that our subtyping approach based on transcriptomic data could contribute to precise prediction of immune checkpoint inhibitor responses in a wide range of cancer types.

Epigenome signatures landscaped by histone H3K9me3 are associated with the synaptic dysfunction in Alzheimer's disease.

The pathogenesis of Alzheimer's disease (AD) and the commonest cause of dementia in the elderly remain incompletely understood. Recently, epigenetic modifications have been shown to play a potential role in neurodegeneration, but the specific involvement of epigenetic signatures landscaped by heterochromatin has not been studied in AD. Herein, we discovered that H3K9me3-mediated heterochromatin condensation is elevated in the cortex of sporadic AD postmortem brains. In order to identify which epigenomes are modulated by heterochromatin, we performed H3K9me3-chromatin immunoprecipitation (ChIP)-sequencing and mRNA-sequencing on postmortem brains from normal subjects and AD patients. The integrated analyses of genome-wide ChIP- and mRNA-sequencing data identified epigenomes that were highly occupied by H3K9me3 and inversely correlated with their mRNA expression levels in AD. Biological network analysis further revealed H3K9me3-landscaped epigenomes to be mainly involved in synaptic transmission, neuronal differentiation, and cell motility. Together, our data show that the abnormal heterochromatin remodeling by H3K9me3 leads to down-regulation of synaptic function-related genes, suggesting that the epigenetic alteration by H3K9me3 is associated with the synaptic pathology of sporadic AD.

NARD: whole-genome reference panel of 1779 Northeast Asians improves imputation accuracy of rare and low-frequency variants.

Here, we present the Northeast Asian Reference Database (NARD), including whole-genome sequencing data of 1779 individuals from Korea, Mongolia, Japan, China, and Hong Kong. NARD provides the genetic diversity of Korean (n = 850) and Mongolian (n = 384) ancestries that were not present in the 1000 Genomes Project Phase 3 (1KGP3). We combined and re-phased the genotypes from NARD and 1KGP3 to construct a union set of haplotypes. This approach established a robust imputation reference panel for Northeast Asians, which yields the greatest imputation accuracy of rare and low-frequency variants compared with the existing panels. NARD imputation panel is available at https://nard.macrogen.com/ .

Integrative analysis of genomic and transcriptomic characteristics associated with progression of aggressive thyroid cancer.

Anaplastic thyroid cancer (ATC) and advanced differentiated thyroid cancers (DTCs) show fatal outcomes, unlike DTCs. Here, we demonstrate mutational landscape of 27 ATCs and 86 advanced DTCs by massively-parallel DNA sequencing, and transcriptome of 13 ATCs and 12 advanced DTCs were profiled by RNA sequencing. TERT, AKT1, PIK3CA, and EIF1AX were frequently co-mutated with driver genes (BRAFV600E and RAS) in advanced DTCs as well as ATC, but tumor suppressors (e.g., TP53 and CDKN2A) were predominantly altered in ATC. CDKN2A loss was significantly associated with poor disease-specific survival in patients with ATC or advanced DTCs, and up-regulation of CD274 (PD-L1) and PDCD1LG2 (PD-L2). Transcriptome analysis revealed a fourth molecular subtype of thyroid cancer (TC), ATC-like, which hardly reflects the molecular signatures in DTC. Furthermore, the activation of JAK-STAT signaling pathway could be a potential druggable target in RAS-positive ATC. Our findings provide insights for precision medicine in patients with advanced TCs.

Peroxiredoxin5 Controls Vertebrate Ciliogenesis by Modulating Mitochondrial Reactive Oxygen Species.

AIMS: Peroxiredoxin5 (Prdx5), a thioredoxin peroxidase, is an antioxidant enzyme that is widely studied for its antioxidant properties and protective roles in neurological and cardiovascular disorders. This study is aimed at investigating the functional significance of Prdx5 in mitochondria and at analyzing its roles in ciliogenesis during the process of vertebrate development. RESULTS: We found that several Prdx genes were strongly expressed in multiciliated cells in developing Xenopus embryos, and their peroxidatic functions were crucial for normal cilia development. Depletion of Prdx5 increased levels of cellular reactive oxygen species (ROS), consequently leading to mitochondrial dysfunction and abnormal cilia formation. Proteomic and transcriptomic approaches revealed that excessive ROS accumulation on Prdx5 depletion subsequently reduced the expression level of pyruvate kinase (PK), a key metabolic enzyme in energy production. We further confirmed that the promotor activity of PK was significantly reduced on Prdx5 depletion and that the reduction in PK expression and its promoter activity led to ciliary defects observed in Prdx5-depleted cells. INNOVATION: Our data revealed the novel relationship between ROS and Prdx5 and the consequent effects of this interaction on vertebrate ciliogenesis. The normal process of ciliogenesis is interrupted by the Prdx5 depletion, resulting in excessive ROS levels and suggesting cilia as vulnerable targets of ROS. CONCLUSION: Prdx5 plays protective roles in mitochondria and is critical for normal cilia development by regulating the levels of ROS. The loss of Prdx5 is associated with excessive production of ROS, resulting in mitochondrial dysfunction and aberrant ciliogenesis.

Comprehensive analysis of the tumor immune micro-environment in non-small cell lung cancer for efficacy of checkpoint inhibitor.

Characterizing the molecular immune subtype and micro-environment of lung cancer is necessary to understand immunogenic interactions between infiltrating immune and stromal cells, and how tumor cells overcome immune checkpoint blockades. This study seeks to identify computational methodologies for subtyping gene expression-based tumor-immune micro-environment interactions, which differentiate non-small cell lung cancer (NSCLC) into immune-defective and immune-competent subtypes. Here, 101 lung squamous cell carcinomas (LUSCs) and 87 lung adenocarcinomas (LUADs) tumor samples have been analyzed. Several micro-environmental factors differentially induce LUAD or LUSC immune subtypes, as well as immune checkpoint expression. In particular, tumor-associated macrophages (TAMs) are key immune cells play a vital role in inflammation and cancer micro-environments of LUSCs; whereas, regulatory B cells are immunosuppressive and tumorigenic in LUADs. Additionally, cytolytic activity upon CD8+ T cell activation is decreased by the abundance of B cells and macrophages in immune-competent subtypes. Therefore, identifying immune subtypes in lung cancer and their impact on tumor micro-environment will lead to clinical tools for assessing LUADs and LUSCs in patients, as well as maximize the efficacy of immune checkpoint inhibitors.

Development of a common platform for the noninvasive prenatal diagnosis of X-linked diseases.

OBJECTIVE: The aim of this study was to develop a common targeted massively parallel sequencing platform for the noninvasive prenatal diagnosis (NIPD) of multiple X-linked diseases. METHOD: The custom capture probe was designed to target 33 genes and recombination hotspots. We tested the carrier mother and male proband pair of 6 families. Plasma DNA of the pregnant carrier mother was collected at different gestational weeks and sequenced. The fetal genotype of each family was determined by estimating the imbalance between the 2 maternal haplotypes constructed using a common custom-designed platform. RESULTS: The targeted sequencing of the maternal, proband, and fetal genomic DNAs and maternal plasma DNAs resulted in uniform coverage across the target region. Three to 5 recombination points were observed in each sample. However, these recombination points did not affect the haplotype dosage analysis for fetal genotype prediction. Consequently, all fetal genotypes in the 6 families obtained from haplotype dosage analysis of maternal plasma sequencing data were predicted correctly. CONCLUSIONS: Since a single platform that covers multiple diseases may prevent the need for disease-specific probes for the NIPD of individual disorders, this approach may provide a practical advantage for clinically implementing the NIPD of X-linked diseases.

Integrative analysis of oncogenic fusion genes and their functional impact in colorectal cancer.

BACKGROUND: Fusion genes are good candidates of molecular targets for cancer therapy. However, there is insufficient research on the clinical implications and functional characteristics of fusion genes in colorectal cancer (CRC). METHODS: In this study, we analysed RNA sequencing data of CRC patients (147 tumour and 47 matched normal tissues) to identify oncogenic fusion genes and evaluated their role in CRC. RESULTS: We validated 24 fusion genes, including novel fusions, by three algorithms and Sanger sequencing. Fusions from most patients were mutually exclusive CRC oncogenes and included tumour suppressor gene mutations. Eleven fusion genes from 13 patients (8.8%) were determined as oncogenic fusion genes by analysing their gene expression and function. To investigate their oncogenic impact, we performed proliferation and migration assays of CRC cell lines expressing fusion genes of GTF3A-CDK8, NAGLU- IKZF3, RNF121- FOLR2, and STRN-ALK. Overexpression of these fusion genes increased cell proliferation except GTF3A-CDK8. In addition, overexpression of NAGLU-IKZF3 enhanced migration of CRC cells. We demonstrated that NAGLU-IKZF3, RNF121-FOLR2, and STRN-ALK had tumourigenic effects in CRC. CONCLUSION: In summary, we identified and characterised oncogenic fusion genes and their function in CRC, and implicated NAGLU-IKZF3 and RNF121-FOLR2 as novel molecular targets for personalised medicine development.

Targeted linked-read sequencing for direct haplotype phasing of maternal DMD alleles: a practical and reliable method for noninvasive prenatal diagnosis.

For the noninvasive prenatal diagnosis (NIPD) of X-linked recessive diseases such as Duchenne muscular dystrophy (DMD), maternal haplotype phasing is a critical step for dosage analysis of the inherited allele. Until recently, the proband-based indirect haplotyping method has been preferred despite its limitations for use in clinical practice. Here, we describe a method for directly determining the maternal haplotype without requiring the proband's DNA in DMD families. We used targeted linked-read deep sequencing (mean coverage of 692×) of gDNA from 5 mothers to resolve their haplotypes and predict the mutation status of the fetus. The haplotype of DMD alleles in the carrier mother was successfully phased through a targeted linked-read sequencing platform. Compared with the proband-based phasing method, linked-read sequencing was more accurate in differentiating whether the recombination events occurred in the proband or in the fetus. The predicted inheritance of the DMD mutation was diagnosed correctly in all 5 families in which the mutation had been confirmed using amniocentesis or chorionic villus sampling. Direct haplotyping by this targeted linked-read sequencing method could be used as a phasing method for the NIPD of DMD, especially when the genomic DNA of the proband is unavailable.