Blog

An efficient molecular genetic testing strategy for incontinentia pigmenti based on single-tube long fragment read sequencing | npj Genomic Medicine

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

npj Genomic Medicine volume  9, Article number: 32 (2024 ) Cite this article Cancer Molecular Diagnostics Laboratory

An efficient molecular genetic testing strategy for incontinentia pigmenti based on single-tube long fragment read sequencing | npj Genomic Medicine

Incontinentia pigmenti (IP) is a rare X-linked dominant neuroectodermal dysplasia that primarily affects females. The only known causative gene is IKBKG, and the most common genetic cause is the recurrent IKBKG△ 4–10 deletion resulting from recombination between two MER67B repeats. Detection of variants in IKBKG is challenging due to the presence of a highly homologous non-pathogenic pseudogene IKBKGP1. In this study, we successfully identified four pathogenic variants in four IP patients using a strategy based on single-tube long fragment read (stLFR) sequencing with a specialized analysis pipeline. Three frameshift variants (c.519-3_519dupCAGG, c.1167dupC, and c.700dupT) were identified and subsequently validated by Sanger sequencing. Notably, c.519-3_519dupCAGG was found in both IKBKG and IKBKGP1, whereas the other two variants were only detected in the functional gene. The IKBKG△ 4–10 deletion was identified and confirmed in one patient. These results demonstrate that the proposed strategy can identify potential pathogenic variants and distinguish whether they are derived from IKBKG or its pseudogene. Thus, this strategy can be an efficient genetic testing method for IKBKG. By providing a comprehensive understanding of the whole genome, it may also enable the exploration of other genes potentially associated with IP. Furthermore, the strategy may also provide insights into other diseases with detection challenges due to pseudogenes.

Incontinentia pigmenti (IP, OMIM 308300), also known as Bloch-Sulzberger syndrome, is a rare neuroectodermal dysplasia characterized by various abnormalities of the skin, hair, teeth, eyes, and central nervous system1. It occurs primarily in females and can cause in utero lethality in males. Skin abnormalities in affected females evolve through four stages from infancy to adulthood: bullous stage, verrucous stage, hyperpigmentation stage, and atretic stage2. The hyperpigmented lesions typically fade after puberty. Other clinical features include alopecia, hypodontia, retinal hypervascularization, seizures, and central nervous system anomalies3,4. Occasionally, some affected males survive, which may be due to the presence of an additional X chromosome (47,XXY), somatic mosaicism, or hypomorphic variants5.

IP is caused by variants in the IKBKG (also known as NEMO) gene on Xq28. The gene is ~23 kb in length and consists of 10 exons. The detection rate of pathogenic variants is ~80%6,7. Deletions of exons 4–10, IKBKG△ 4–10, account for the majority of identified variants8,9. The complexity of variant detection is due to a highly homologous non-pathogenic pseudogene, IKBKGP1, which is located 31 kb distal to IKBKG in the opposite orientation10. In contrast to IKBKG, deletion of exons 4–10 in IKBKGP1 does not cause symptoms11. In addition, the existence of two 879 bp repeats, termed MER67B, one in intron 3 (MER67B1st) and one downstream of exon 10 (MER67B2nd), further complicates the diagnosis. Despite the advances in sequencing technology, molecular testing for IP is still mainly based on long-range polymerase chain reaction (PCR) and multiplex ligation-dependent probe amplification (MLPA)3,12. Therefore, alternative efficient approaches are needed to identify IKBKG-specific variations. Single-tube long fragment read (stLFR, MGI Tech) has been reported as an efficient technology that can utilize next-generation sequencing (NGS) platforms to sequence long DNA molecules13. By co-barcoding, short reads from the same long DNA fragment can be assembled according to the barcodes. The stLFR can retain long-range genomic information of ~20–300 kb in length. It may thus provide a solution for addressing complex genetic issues, such as identifying large structural variations (SVs) and distinguishing functional genes from pseudogenes. However, its related clinical applications have not yet been fully demonstrated.

Here, we present an efficient strategy for genetic testing of IP based on stLFR sequencing. We applied this strategy to four IP families and successfully identified four pathogenic variants in the IKBKG gene. The new strategy enables us to detect putative pathogenic variants throughout the genome and can differentiate specific variants of IKBKG from IKBKGP1. Thus, the strategy proposed here is able to overcome pseudogene-related complications in IP detection.

Four patients (Fig. 1) were detected using a strategy based on stLFR sequencing with a specialized analysis pipeline (Fig. 2). The quality of stLFR sequencing data was summarized in Table 1. The average genome sequencing depth was ~20×, ranging from 19.98× to 22.47×. The average long fragment length of the four samples was 39.29 kb, 22.41 kb, 35.4 kb, and 48.48 kb, respectively. Based on the routine analysis of stLFR, >3 million SNPs, >800 thousand insertion-deletion mutations (indels), and hundreds of copy number variants (CNVs) were identified for each sample, and none of them was located in the genomic region related to IP.

Affected individuals are shown in black, while arrows indicate the probands in each family (a–d). The genotype of each evaluated individual is displayed below his or her symbol. A plus sign denotes the normal allele, and a minus sign denotes the mutant allele.

a Schematic representation of the IKBKG gene and its downstream pseudogene IKBKGP1. The IKBKG gene contains 10 exons (exons 1–10), while IKBKGP1 contains only 7 exons (exons 3–10). The gray square arrows represent the highly homologous regions, and the green boxes represent the non-homologous regions. The MER67B repeated sequences are shown as yellow boxes. b Diagram of the analysis pipelines. The left side shows the routine bioinformatic pipeline of stLFR, and the right side shows a specialized analysis pipeline for IP.

After filtering with a specialized analysis pipeline, three intragenic variants were found in the IKBKG gene (NM_001099857.5: c.519-3_519dupCAGG, c.1167dupC and c.700dupT) (Fig. 3 and Table 2). Considering the inheritance pattern, allele frequency, predicted impact, and the annotation results of VarSome, these three variants were considered potential pathogenic variants for patients P1-P3. Notably, variant c.519-3_519dupCAGG was found in both IKBKG and IKBKGP1 in P1, whereas the other two variants were only detected in the functional gene. Subsequent variant testing was performed in the probands and their family members by Sanger sequencing (Figs. 1 and 4). In IP families 1 and 2, both the proband and the mother carried the variant, suggesting maternal transmission. In family 3, the variant proved to be de novo.

The c.519-3_519dupCAGG variant is found in both IKBKG (a) and IKBKGP1 (b) in P1. The c.1167dupC and c.700dupT variants are present only in the IKBKG gene in P2 (c) and P3 (d), respectively.

Variants are indicated by arrows. a c.519-3_519dupCAGG in P1; b c.1167dupC in P2; c c.700dupT in P3.

These frameshift variants were predicted to cause the loss of normal protein function either through protein truncation or nonsense-mediated mRNA decay. Two of them (c.519-3_519dupCAGG and c.1167dupC) have been previously reported7,14,15, and published functional studies have shown that c.1167dupC would cause a damaging effect on protein function16,17. None of these variants were observed in the gnomAD database. Therefore, they were considered pathogenic or likely pathogenic.

All assembly results of Supernova and SPAdes were demonstrated in Supplementary Figure 1. Two scaffolds (scaffold 21 in P1 and scaffold 51 in P4) with breakpoints near MER67B in IKBKG were detected by Supernova. They were subsequently excluded because the breakpoints were caused by N-gaps (Fig. 5). With SPAdes, two discontinuous or misassembled assemblers between two MER67B repeats were detected (NODE_3 in P2 and NODE_1 in P4) (Supplementary Figure 1 and Fig. 5). Compared to Supernova, the assemblers of SPAdes were longer, and no N was present. To further determine the authenticity of these two assemblers, reads passing through the breakpoint were extracted for long DNA fragment analysis. In P2, one of six reads was aligned to the MER67B1st sequence (Ref MER67B1st), and the remaining five were aligned to either the IKBKG△ 4–10 mutant sequence (Mut) or the MER67B2nd sequence (Ref MER67B2nd) (Supplementary Figure 2a). Because IKBKG△ 4–10 is caused by recombination between two MER67B repeat sequences, the reads supporting assemblers at the breakpoint could be a true deletion signal or simply mapped to Ref MER67B2nd. In P4, five extracted reads were mapped to Ref MER67B1st, and the other five were mapped to either Mut or Ref MER67B2nd (Supplementary Figure 3a). To verify the validity of these reads, we further used the co-barcoding information to trace all reads of these long DNA fragments. The distribution of all co-barcoded reads was illustrated in Supplementary Fig. 2b and 3b. In P2, multiple reads were mapped to the deletion region in three long fragments, demonstrating that there was indeed no deletion in these fragments. The remaining two fragments (1439_459_1531 and 851_274_1146) were considered as potential IKBKG△ 4–10, due to the absence of reads in the deletion region. In P4, three fragments turned out to be aligned to Ref MER67B2nd due to the presence of co-barcoded reads in the deletion region, while the other two fragments (531_1256_485 and 221_57_1365) were considered as potential IKBKG△ 4–10 (Supplementary Figure 3b). No CNVs were found in IKBKGP1 (data not shown).

Scaffolds detected by Supernova or SPAdes in P1 (a), P2 (b) and P4 (c and d). The horizontal lines represent the assembly sequences, and the red vertical lines represent the breakpoints. The gray rectangles are the N-gap regions, and the yellow rectangles are the recombined MER67B region. The genomic position (hg19) of the breakpoint is shown below each contig.

The IKBKG△ 4–10 deletions identified in P2 and P4 were then validated by MLPA (Fig. 6). The deletion in P2 turned out to be a false positive result, while the deletion in P4 was confirmed to be true. Since MLPA could not determine whether the deletion was in IKBKG or IKBKGP1, long-range PCR followed by Sanger sequencing was performed. The results showed that the deletion was in the IKBKG gene (data not shown). Combined with the results of stLFR analysis, IKBKG△ 4–10 in P4 was an IKBKG-specific deletion. MLPA confirmed that P4’s mother and sister also carried IKBKG△ 4–10 (data not shown).

a, b are MLPA results of P2 and P4, respectively.

Currently, the main challenge in IP detection is to differentiate the true IKBKG gene from its pseudogene IKBKGP1. To identify SNVs in the IKBKG gene, a common approach is to test all exons, exon-surrounding intronic regions, and promoter regions separately3,10. To eliminate the IKBKGP1 gene, a first round of long-range PCR would be performed18,19,20 (Supplementary Table 1). An alternative approach to identify variants in the functional IKBKG gene is direct sequencing of the cDNA16,19,21. However, sequencing each fragment can be a tedious task, and variants in regulatory, deep intronic, or non-coding regions remain undetected. Although NGS is widely used in clinical molecular testing, its application in IP is limited. Due to a 35.5 kb homologous region shared by IKBKG and IKBKGP110, it is extremely difficult for short-read technologies to identify specific variants. Attempts have been made using long-range PCR followed by NGS6,22. For CNVs, nested long-range PCR remains the gold standard method10,23,24. Another commonly used method is MLPA, which, however, cannot differentiate between IKBKG and its pseudogene3. Therefore, there is still a need to develop efficient testing approaches that can identify variants in both exons and introns. To fill the gaps, we proposed a strategy using stLFR in combination with unique barcodes to retrieve long fragment information. Routine analysis pipelines of stLFR can provide a comprehensive picture of the whole genome. Besides SNVs and CNVs, indels and other SVs can also be detected. However, routine analysis cannot identify variants in the homologous region of IKBKG due to mapping quality values of 0. Therefore, we established a specialized analysis pipeline for IKBKG. All putative pathogenic variants in the functional IKBKG gene can be detected and distinguished from IKBKGP1. Our data further showed that three of these variants (c.1167dupC, c.700dupT, and IKBKG△ 4–10) were present only in the IKBKG gene, and one (c.519-3_519dupCAGG) was present in both IKBKG and IKBKGP1. Additionally, our work may facilitate the discovery of novel putative variants or genes. Approximately one-fifth of all cases still have no known cause7,25. Individuals may have low-level somatic mosaicism6, pathogenic variants in other regions of the IKBKG gene7, or even in other genes. Our strategy can detect variants in these regions more effectively than conventional methods. Given that IKBKG is the only gene associated with IP and that some patients have a milder or different clinical phenotype6,26, investigation of other potential genetic causes should be considered. Although little data was used to analyze the IKBKG and IKBKGP1 genes, the remaining data could be used to identify other SVs and other causative genes. One drawback of stLFR is the cost, but it could provide a more comprehensive understanding of genome structural features and overcome the interference of highly homologous pseudogenes in the molecular testing of IPs. This, in turn, may improve IP detection and is worth considering.

In our proposed strategy, two assembly software were used to detect CNVs. Supernova is one of the most commonly used software for de novo assembly of linked reads and is capable of generating diploid assemblies of the human genome27,28. However, our results showed that Supernova would easily introduce Ns into the assembly results, making it difficult to detect the target CNVs. This may be due to the limited amount of data for assembly contigs. The sequencing depth of genomes in previous studies is generally more than 60X29,30,31, while it was about 20X in our study. Increasing the amount of sequencing data may improve the ability of Supernova to detect IKBKG△ 4–10. SPAdes, on the other hand, is a flexible assembler that is suitable for various data formats generated by different sequencing platforms and is capable of generating long and accurate assembly results from cross-species raw sequencing reads32,33. Although the amount of data was limited, our study showed that SPAdes successfully identified IKBKG△ 4–10. According to the MLPA verification results, the deletion in P4 was confirmed, while the one in P2 was a false positive. Although it may be frustrating to require MLPA and long-range PCR for verification, the strategy itself is logical. Insufficient data and low-coverage reads of long DNA fragments are possible causes. With the continuous reduction of sequencing costs, the specificity and positive predictive value of the pipeline for detecting assembly CNVs could be improved by increasing the amount of sequencing data.

In conclusion, we have proposed a feasible and promising strategy for IP testing using stLFR with a specialized analysis pipeline. It could provide a comprehensive understanding of the whole genome, and all putative pathogenic variants of IP could be detected. More importantly, specific variants of IKBKG could be distinguished from IKBKGP1. In addition, our strategy has the potential to uncover additional genes that may be associated with IP. The performance of our pipeline for IP testing could be further improved with the decreasing cost of sequencing. Furthermore, the strategy proposed here can address pseudogene-related issues in IP testing and provide insights into other diseases with detection challenges due to pseudogenes.

Five patients from four families were clinically diagnosed with IP based on typical skin manifestations (Fig. 1). P1 was further confirmed by skin biopsy. Both mothers of P1 and P2 had very mild hyperpigmentation. P3 had no family history. The mother of P4 had mild blisters at birth. Unilateral amblyopia and agenesis of permanent teeth occurred in one patient (P1), while one patient developed blindness in the right eye (P4). None of them exhibited any neurological disorders.

Written informed consent was obtained to participate and publish from both parents and patients. This study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee at Women’s Hospital, School of Medicine, Zhejiang University (20190038). The study is compliant with the Guidance of the Ministry of Science and Technology (MOST) for the Review and Approval of Human Genetic Resources.

The genomic DNA was extracted from peripheral blood leukocytes using the QIAGEN MagAttract HMW DNA Kit (QIAGEN, Germany) following the manufacturer’s protocol. An stLFR library was constructed using the MGIEasy stLFR Library Prep Kit (MGI Tech, China) according to the manufacturer’s instructions. Briefly, transposons were inserted into long DNA molecules. Subsequently, these transposon-inserted DNA sequences were hybridized with clonal barcoded beads and then ligated with barcoded oligo and adapters via splint oligo. After adding the library adapters, PCR amplification and circularization were performed to generate DNA nanoballs (DNBs). The prepared library was then sequenced on the MGISEQ-2000 platform (MGI Tech) with a 100 bp paired-end strategy. The expected raw data of each sample is 100GB or more.

The routine bioinformatic pipelines of stLFR sequencing were summarized in Fig. 2. Briefly, raw sequencing reads were first demultiplexed with the associated barcodes of long DNA molecules using the barcode split tool (GitHub). Routine pre-processing of the raw reads was then conducted, including filtering of low-quality reads, alignment (software: BWA, reference genome: GRCh37/HG19), and elimination of PCR duplicates. Variants, including SNVs, indels, CNVs, and SVs, were detected and annotated based on WGS pipelines using the MegaBOLT system (MGI Tech).

New reference genomes were constructed according to the genomic characteristics (Fig. 2a). Non-homologous regions and an additional 50 kb of extended regions were used to extract co-barcoded sequencing reads of long DNA fragments, which were designated as “non-homologous region 1” and “non-homologous region 2” for IKBKG and IKBKGP1, respectively. Sequence related to IKBKG was designated as reconstructed reference 1, while sequence related to IKBKGP1 was designated as reconstructed reference 2. Based on the primary data analysis, a specialized analysis pipeline was developed to distinguish the functional IKBKG gene from its pseudogene IKBKGP1 (Fig. 2b). Variants were detected by GATK HaplotypeCaller and filtered with hard filter patterns (Supplementary methods). The reconstructed reference positions were then converted to hg19 for subsequent annotation analysis. Finally, variants were annotated using unpublished in-house software. Annotated information included basic variant information, protein functional hazard predictions, nucleotide conservation predictions, ClinVar significance records, HGMD records, OMIM records, and pathogenicity predictions. Non-benign homozygous or heterozygous variants matching the inheritance pattern were selected as candidates. The VarSome website was then used to predict the pathogenicity of these variants, including SNVs and indels. Likely pathogenic and pathogenic variants of IKBKG were identified as potential IP-causing variants.

As demonstrated in Fig. 2b, barcodes of reads aligned to non-homologous reference regions were collected. Reads sharing these barcodes were then retrieved and assembled. The CNV detection strategy based on local de novo assembly was implemented using Supernova (version 2.1.1) and SPAdes. The analysis based on Supernova was conducted with an open-source software named “stlfr2supernova_pepeline” (https://github.com/BGI-Qingdao/stlfr2supernova_pipeline) according to the user’s manual. Meanwhile, the main parameters of SPAdes were ‘--careful --only-assembler -t 16’. The assembly results were then aligned separately to the human genome (hg19) using the QUAST software (version 5.1.0rc1) with the parameters ‘-t 20 --min-identity 95 --min-contig 500’. Based on the alignment results, discontinuous or misassembled assemblers were selected as potential CNVs. Subsequently, assemblers with breakpoints caused by Ns were excluded. For the remaining assemblers, reads passing through the breakpoint were extracted, and their co-barcoded long DNA fragments were traced. The distribution patterns of all reads belonging to these long DNA fragments were then analyzed. Finally, assemblers with no other reads distributed in the deletion region were considered candidate CNVs.

SNVs and indels were validated by PCR and Sanger sequencing, while CNVs were confirmed by MPLA using MLPA probemix P073-A1 (MRC Holland, Holland) and long-range PCR. Capillary electrophoresis was performed on ABI3500 (Life Technologies), and analyzed using Coffalyser (MRC Holland) or Gene Mapper software (Life Technologies).

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

The data that support the findings of this study have been deposited in the Genome Sequence Archive for Human with the accession code HRA006663.

The basic analysis pipelines of the stLFR sequencing data are publicly available on GitHub (https://github.com/MGI-tech-bioinformatics/stLFR_V1.3). The detection analysis of SNVs, indels, and CNVs that could distinguish the IKBKG gene from its pseudogene was based on open source software and is detailed in Methods and Supplementary Methods.

Scheuerle, AE & Ursini, M. v.Incontinence of pigment.in GeneReviews((R)) (eds. Adam, MP et al.) (Seattle, WA, 1993).

Bodemer, C. et al. Multidisciplinary consensus recommendations from a European network for the diagnosis and practical management of patients with incontinentia pigmenti. J. Eur. Acad. Dermatol. Venereol. 34, 1415–1424 (2020).

Article  CAS  PubMed  Google Scholar 

Kim, H. Y. et al. Importance of extracutaneous organ involvement in determining the clinical severity and prognosis of incontinentia pigmenti caused by mutations in the IKBKG gene. Exp. Dermatol. 30, 676–683 (2021).

Article  CAS  PubMed  Google Scholar 

Cammarata-Scalisi, F., Fusco, F. & Ursini, MV Pigment incontinence.Actas Dermosifiliogr.(Eng. Ed.) 110, 273–278 (2019).

Article  CAS  PubMed  Google Scholar 

Kenwrick, S. et al. Survival of male patients with incontinentia pigmenti carrying a lethal mutation can be explained by somatic mosaicism or Klinefelter syndrome. Am. J. Hum. Genet. 69, 1210–1217 (2001).

Article  CAS  PubMed  Google Scholar 

Kawai, M. et al. Molecular analysis of low-level mosaicism of the IKBKG mutation using the X chromosome inactivation pattern in Incontinentia Pigmenti. Mol. Genet. Genom. Med. 8, e1531 (2020).

Conte, M. I. et al. Insight into IKBKG/NEMO locus: report of new mutations and complex genomic rearrangements leading to incontinentia pigmenti disease. Hum. Mutat. 35, 165–177 (2014).

Article  CAS  PubMed  Google Scholar 

Aradhya, S. et al.A recurrent deletion in the ubiquitously expressed NEMO (IKK-gamma) gene accounts for the vast majority of incontinent pigment mutations.Hum.Mol.Genet.10, 2171–2179 (2001).

Article  CAS  PubMed  Google Scholar 

Fusco, F. et al. Molecular analysis of the genetic defect in a large cohort of IP patients and identification of novel NEMO mutations interfering with NF-kappaB activation. Hum. Mol. Genet. 13, 1763–1773 (2004).

Article  CAS  PubMed  Google Scholar 

Haque, M. N. et al. Analysis of IKBKG/NEMO gene in five Japanese cases of incontinentia pigmenti with retinopathy: fine genomic assay of a rare male case with mosaicism. J. Hum. Genet. 66, 205–214 (2021).

Article  CAS  PubMed  Google Scholar 

Francesca, F. et al. Microdeletion/duplication at the Xq28 IP locus causes a de novo IKBKG/NEMO/IKKgamma exon4_10 deletion in families with incontinentia pigmenti. Hum. Mutat. 30, 1284–1291 (2009).

Pizzamiglio, M. R. et al. Incontinentia pigmenti: learning disabilities are a fundamental hallmark of the disease. PLoS One 9, e87771 (2014).

Article  PubMed  PubMed Central  Google Scholar 

Wang, O. et al. Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly. Genome Res. 29, 798–808 (2019).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Zonana, J. et al. A novel X-linked disorder of immune deficiency and hypohidrotic ectodermal dysplasia is allelic to incontinentia pigmenti and due to mutations in IKK-gamma (NEMO). Am. J. Hum. Genet. 67, 1555–1562 (2000).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Jiang, J. et al. NEMO gene mutations in two chinese females with incontinentia pigmenti. Clin. Cosmet. Investig. Dermatol. 15, 815–821 (2022).

Article  PubMed  PubMed Central  Google Scholar 

Ohnishi, H. et al. Immunodeficiency in two female patients with incontinentia pigmenti with heterozygous NEMO mutation diagnosed by LPS unresponsiveness. J. Clin. Immunol. 37, 529–538 (2017).

Article  CAS  PubMed  Google Scholar 

Kawai, T. et al. Frequent somatic mosaicism of NEMO in T cells of patients with X-linked anhidrotic ectodermal dysplasia with immunodeficiency. Blood 119, 5458–5466 (2012).

Article  CAS  PubMed  Google Scholar 

Dangouloff-Ros, V. et al. Severe neuroimaging anomalies are usually associated with random X inactivation in leucocytes circulating DNA in X-linked dominant Incontinentia Pigmenti. Mol. Genet. Metab. 122, 140–144 (2017).

Article  CAS  PubMed  Google Scholar 

Sun, S. et al. A novel inhibitor of nuclear factor kappa-B kinase subunit gamma mutation identified in an incontinentia pigmenti patient with syndromic tooth agenesis. Arch. Oral. Biol. 101, 100–107 (2019).

Article  CAS  PubMed  Google Scholar 

Fryssira, H. et al.Incontinence pigment revisited.A novel nonsense mutation of the IKBKG gene.Acta Paediatr.100, 128–133 (2011).

Article  CAS  PubMed  Google Scholar 

Ramírez-Alejo, N. et al. Novel hypomorphic mutation in IKBKG impairs NEMO-ubiquitylation causing ectodermal dysplasia, immunodeficiency, incontinentia pigmenti, and immune thrombocytopenic purpura. Clin. Immunol. 160, 163–171 (2015).

Hull, S. et al. Somatic mosaicism of a novel IKBKG mutation in a male patient with incontinentia pigmenti. Am. J. Med. Genet. A 167, 1601–1604 (2015).

Article  CAS  PubMed  Google Scholar 

Bardaro, T. et al. Two cases of misinterpretation of molecular results in incontinentia pigmenti, and a PCR-based method to discriminate NEMO/IKKgamma dene deletion. Hum. Mutat. 21, 8–11 (2003).

Article  CAS  PubMed  Google Scholar 

Steffann, J. et al. A novel PCR approach for prenatal detection of the common NEMO rearrangement in incontinentia pigmenti. Prenat. Diagn. 24, 384–388 (2004).

Article  CAS  PubMed  Google Scholar 

Article  PubMed  PubMed Central  Google Scholar 

Wang, R., Lara-Corrales, I., Kannu, P. & Pope, E. Unraveling incontinentia pigmenti: a comparison of phenotype and genotype variants. J. Am. Acad. Dermatol 81, 1142–1149 (2019).

Article  CAS  PubMed  Google Scholar 

Weisenfeld, N. I., Kumar, V., Shah, P., Church, D. M. & Jaffe, D. B. Direct determination of diploid genome sequences. Genome Res. 27, 757–767 (2017).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Marks, P. et al. Resolving the full spectrum of human genome variation using linked-reads. Genome Res. 29, 635–645 (2019).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Peters, B. A., Liu, J. & Drmanac, R. Co-barcoded sequence reads from long DNA fragments: a cost-effective solution for ‘perfect genome’ sequencing. Front. Genet. 5, 466 (2014).

Zhang, L., Zhou, X., Weng, Z. & Sidow, A. Assessment of human diploid genome assembly with 10x linked-reads data. Gigascience 8, giz141 (2019).

Article  PubMed  PubMed Central  Google Scholar 

Mao, Q. et al. The whole genome sequences and experimentally phased haplotypes of over 100 personal genomes. Gigascience 5, 42 (2016).

Article  PubMed  PubMed Central  Google Scholar 

Prjibelski, A., Antipov, D., Meleshko, D., Lapidus, A. & Korobeynikov, A. Using SPAdes de novo assembler.Curr.Protoc.Bioinformatics 70, e102 (2020).

Article  CAS  PubMed  Google Scholar 

Ishengoma, E. & Rhode, C. Using SPAdes, AUGUSTUS, and BLAST in an automated pipeline for clustering homologous exome sequences. Curr. Protoc. 2, e449 (2022).

Article  CAS  PubMed  Google Scholar 

The work was supported by the National Natural Science Foundation of China (81901382 and 82171848), the Key Projects Jointly Constructed by the Ministry and the Province of Zhejiang Medical and Health Science and Technology Project (WKJ-ZJ-2127), and Zhejiang Provincial Natural Science Foundation of China (LY22H110004).

These authors contributed equally: Min Chen, Mei-Hua Tan.

Women’s Hospital, School of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310006, P. R. China

Min Chen, Yan-Mei Yang, Jia-Ling Yu, Ying-Zhi Huang, Yi-Xi Sun, Ye-Qing Qian, Kai Yan & Min-Yue Dong

Key Laboratory of Reproductive Genetics (Zhejiang University), Ministry of Education, Hangzhou, Zhejiang, 310006, P. R. China

Min Chen, Yan-Mei Yang, Jia-Ling Yu, Ying-Zhi Huang, Yi-Xi Sun, Ye-Qing Qian, Kai Yan & Min-Yue Dong

Key Laboratory of Women’s Reproductive Health of Zhejiang Province, Hangzhou, Zhejiang, 310006, P. R. China

Min Chen, Yan-Mei Yang, Jia-Ling Yu, Ying-Zhi Huang, Yi-Xi Sun, Ye-Qing Qian, Kai Yan & Min-Yue Dong

BGI Genomics, Shenzhen, Guangdong, 518083, P. R. China

Mei-Hua Tan & Li-Juan He

Lishui Maternity and Child Health Care Hospital, Lishui, Zhejiang, 323000, P. R. China

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

M.C. and M.H.T. contributed equally to this article. M.Y.D. conceived and coordinated the project and experimental design, recruited the patients, and revised the manuscript. M.C. contributed to the conception and design of the study, performed the experiments, analyzed data, and drafted the manuscript. M.H.T. performed the experiments, analyzed data, and drafted the manuscript. J.L. contributed to patient recruitment. Y.M.Y. and J.L.Y. participated in the experiment’s conduction. L.J.H. assisted with data analysis. Y.Z.H., Y.Q.Q., Y.X.S., and K.Y. contributed to variants interpretation and confirmation. All authors have read, revised, and approved the final manuscript.

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Chen, M., Tan, MH., Liu, J. et al. An efficient molecular genetic testing strategy for incontinentia pigmenti based on single-tube long fragment read sequencing. npj Genom. Med. 9, 32 (2024). https://doi.org/10.1038/s41525-024-00421-z

DOI: https://doi.org/10.1038/s41525-024-00421-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

npj Genomic Medicine (npj Genom. Med.) ISSN 2056-7944 (online)

An efficient molecular genetic testing strategy for incontinentia pigmenti based on single-tube long fragment read sequencing | npj Genomic Medicine

Ce Ivd Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.