Kagami Ogata syndrome: a small deletion refines critical region for imprinting.

Kagami-Ogata syndrome is a rare imprinting disorder and its phenotypic overlap with multiple different etiologies hampers diagnosis. Genetic etiologies include paternal uniparental isodisomy (upd(14)pat), maternal allele deletions of differentially methylated regions (DMR) in 14q32.2 or pure primary epimutations. We report a patient with Kagami-Ogata syndrome and an atypical diagnostic odyssey with several negative standard-of-care genetic tests followed by epigenetic testing using methylation microarray and a targeted analysis of whole-genome sequencing to reveal a 203 bp deletion involving the MEG3 transcript and MEG3:TSS-DMR. Long-read sequencing enabled the simultaneous detection of the deletion, phasing, and biallelic hypermethylation of the MEG3:TSS-DMR region in a single assay. This case highlights the challenges in the sequential genetic testing paradigm, the utility of long-read sequencing as a single comprehensive diagnostic assay, and the smallest reported deletion causing Kagami-Ogata syndrome allowing important insights into the mechanism of imprinting effects at this locus.

Comprehensive Gene Panel Testing for Hearing Loss in Children: Understanding Factors Influencing Diagnostic Yield.

OBJECTIVE: To evaluate factors influencing the diagnostic yield of comprehensive gene panel testing (CGPT) for hearing loss (HL) in children and to understand the characteristics of undiagnosed probands. STUDY DESIGN: This was a retrospective cohort study of 474 probands with childhood-onset HL who underwent CGPT between 2016 and 2020 at a single center. Main outcomes and measures included the association between clinical variables and diagnostic yield and the genetic and clinical characteristics of undiagnosed probands. RESULTS: The overall diagnostic yield was 44% (209/474) with causative variants involving 41 genes. While the diagnostic yield was high in the probands with congenital, bilateral, and severe HL, it was low in those with unilateral, noncongenital, or mild HL; cochlear nerve deficiency; preterm birth; neonatal intensive care unit admittance; certain ancestry; and developmental delay. Follow-up studies on 49 probands with initially inconclusive CGPT results changed the diagnostic status to likely positive or negative outcomes in 39 of them (80%). Reflex to exome sequencing on 128 undiagnosed probands by CGPT revealed diagnostic findings in 8 individuals, 5 of whom had developmental delays. The remaining 255 probands were undiagnosed, with 173 (173/255) having only a single variant in the gene(s) associated with autosomal recessive HL and 28% (48/173) having a matched phenotype. CONCLUSION: CGPT efficiently identifies the genetic etiologies of HL in children. CGPT-undiagnosed probands may benefit from follow-up studies or expanded testing.

Long-read sequencing for molecular diagnostics in constitutional genetic disorders.

Long-read sequencing (LRS) has been around for more than a decade, but widespread adoption of the technology has been slow due to the perceived high error rates and high sequencing cost. This is changing due to the recent advancements to produce highly accurate sequences and the reducing costs. LRS promises significant improvement over short read sequencing in four major areas: (1) better detection of structural variation (2) better resolution of highly repetitive or nonunique regions (3) accurate long-range haplotype phasing and (4) the detection of base modifications natively from the sequencing data. Several successful applications of LRS have demonstrated its ability to resolve molecular diagnoses where short-read sequencing fails to identify a cause. However, the argument for increased diagnostic yield from LRS remains to be validated. Larger cohort studies may be required to establish the realistic boundaries of LRS's clinical utility and analytical validity, as well as the development of standards for clinical applications. We discuss the limitations of the current standard of care, and contrast with the applications and advantages of two major LRS platforms, PacBio and Oxford Nanopore, for molecular diagnostics of constitutional disorders, and present a critical argument about the potential of LRS in diagnostic settings.

Genome sequencing increases diagnostic yield in clinically diagnosed Alagille syndrome patients with previously negative test results.

PURPOSE: Detection of all major classes of genomic variants in a single test would decrease cost and increase the efficiency of genomic diagnostics. Genome sequencing (GS) has the potential to provide this level of comprehensive detection. We sought to demonstrate the utility of GS in the molecular diagnosis of 18 patients with clinically defined Alagille syndrome (ALGS), who had a negative or inconclusive result by standard-of-care testing. METHODS: We performed GS on 16 pathogenic variant-negative probands and two probands with inconclusive results (of 406 ALGS probands) and analyzed the data for sequence, copy-number, and structural variants in JAG1 and NOTCH2. RESULTS: GS identified four novel pathogenic alterations including a copy-neutral inversion, a partial deletion, and a promoter variant in JAG1, and a partial NOTCH2 deletion, for an additional diagnostic yield of 0.9%. Furthermore, GS resolved two complex rearrangements, resulting in identification of a pathogenic variant in 97.5% (n = 396/406) of patients after GS. CONCLUSION: GS provided an increased diagnostic yield for individuals with clinically defined ALGS who had prior negative or incomplete genetic testing by other methods. Our results show that GS can detect all major classes of variants and has potential to become a single first-tier diagnostic test for Mendelian disorders.

Impaired Redox and Protein Homeostasis as Risk Factors and Therapeutic Targets in Toxin-Induced Biliary Atresia.

BACKGROUND & AIMS: Extrahepatic biliary atresia (BA) is a pediatric liver disease with no approved medical therapy. Recent studies using human samples and experimental modeling suggest that glutathione redox metabolism and heterogeneity play a role in disease pathogenesis. We sought to dissect the mechanistic basis of liver redox variation and explore how other stress responses affect cholangiocyte injury in BA. METHODS: We performed quantitative in situ hepatic glutathione redox mapping in zebrafish larvae carrying targeted mutations in glutathione metabolism genes and correlated these findings with sensitivity to the plant-derived BA-linked toxin biliatresone. We also determined whether genetic disruption of HSP90 protein quality control pathway genes implicated in human BA altered biliatresone toxicity in zebrafish and human cholangiocytes. An in vivo screening of a known drug library was performed to identify novel modifiers of cholangiocyte injury in the zebrafish experimental BA model, with subsequent validation. RESULTS: Glutathione metabolism gene mutations caused regionally distinct changes in the redox potential of cholangiocytes that differentially sensitized them to biliatresone. Disruption of human BA-implicated HSP90 pathway genes sensitized zebrafish and human cholangiocytes to biliatresone-induced injury independent of glutathione. Phosphodiesterase-5 inhibitors and other cyclic guanosine monophosphate signaling activators worked synergistically with the glutathione precursor N-acetylcysteine in preventing biliatresone-induced injury in zebrafish and human cholangiocytes. Phosphodiesterase-5 inhibitors enhanced proteasomal degradation and required intact HSP90 chaperone. CONCLUSION: Regional variation in glutathione metabolism underlies sensitivity to the biliary toxin biliatresone and may account for the reported association between BA transplant-free survival and glutathione metabolism gene expression. Human BA can be causatively linked to genetic modulation of protein quality control. Combined treatment with N-acetylcysteine and cyclic guanosine monophosphate signaling enhancers warrants further investigation as therapy for BA.

Thermally Stable Low-Loss Polymer Dielectrics Enabled by Attaching Cross-Linkable Antioxidant to Polypropylene.

Polymer dielectrics with low-loss and high-temperature tolerance are extremely desirable as electrical energy storage materials for advanced electronics and electrical power applications. They can allow fast switching rates during power conversion and therefore achieve high power densities without thermal issues. Here, we explore polypropylene (PP), the state of the art dielectric polymer, and present an innovative approach to substantially improve the thermal stability and concurrently reduce the dielectric loss of PP. In particular, cross-linkable antioxidant groups, hindered phenol (HP), are incorporated into PP via well-controlled chemical synthesis. The grafted HP can simultaneously serve as radical scavenger and cross-linker, thereby constraining thermally decomposed radicals and charge transport in the synthesized PP-HP copolymer. As a result, the upper-temperature limit of PP-HP is greatly extended to 190 °C and the electrical loss is even gradually reduced upon thermal annealing. The copolymer after heating under 190 °C exhibits better dielectric properties than the PP without any thermal treatment. The experimental results indicate that the PP-HP copolymers are promising materials for high-temperature, low-loss, and high-voltage dielectric applications.

Exome Sequencing in Individuals with Isolated Biliary Atresia.

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.

A highly sensitive and specific workflow for detecting rare copy-number variants from exome sequencing data.

BACKGROUND: Exome sequencing (ES) is a first-tier diagnostic test for many suspected Mendelian disorders. While it is routine to detect small sequence variants, it is not a standard practice in clinical settings to detect germline copy-number variants (CNVs) from ES data due to several reasons relating to performance. In this work, we comprehensively characterized one of the most sensitive ES-based CNV tools, ExomeDepth, against SNP array, a standard of care test in clinical settings to detect genome-wide CNVs. METHODS: We propose a modified ExomeDepth workflow by excluding exons with low mappability prior to variant calling to drastically reduce the false positives originating from the repetitive regions of the genome, and an iterative variant calling framework to assess the reproducibility. We used a cohort of 307 individuals with clinical ES data and clinical SNP array to estimate the sensitivity and false discovery rate of the CNV detection using exome sequencing. Further, we performed targeted testing of the STRC gene in 1972 individuals. To reduce the number of variants for downstream analysis, we performed a large-scale iterative variant calling process with random control cohorts to assess the reproducibility of the CNVs. RESULTS: The modified workflow presented in this paper reduced the number of total variants identified by one third while retaining a higher sensitivity of 97% and resulted in an improved false discovery rate of 11.4% compared to the default ExomeDepth pipeline. The exclusion of exons with low mappability removes 4.5% of the exons, including a subset of exons (0.6%) in disease-associated genes which are intractable by short-read next-generation sequencing (NGS). Results from the reproducibility analysis showed that the clinically reported variants were reproducible 100% of the time and that the modified workflow can be used to rank variants from high to low confidence. Targeted testing of 30 CNVs identified in STRC, a challenging gene to ascertain by NGS, showed a 100% validation rate. CONCLUSIONS: In summary, we introduced a modification to the default ExomeDepth workflow to reduce the false positives originating from the repetitive regions of the genome, created a large-scale iterative variant calling framework for reproducibility, and provided recommendations for implementation in clinical settings.

Protein-elongating mutations in MYH11 are implicated in a dominantly inherited smooth muscle dysmotility syndrome with severe esophageal, gastric, and intestinal disease.

Gastrointestinal motility disorders include a spectrum of mild to severe clinical phenotypes that are caused by smooth muscle dysfunction. We investigated the genetic etiology of severe esophageal, gastric, and colonic dysmotility in two unrelated families with autosomal dominant disease presentation. Using exome sequencing, we identified a 2 base pair insertion at the end of the myosin heavy chain 11 (MYH11) gene in all affected members of Family 1 [NM_001040113:c.5819_5820insCA(p.Gln1941Asnfs*91)] and a 1 base pair deletion at the same genetic locus in Proband 2 [NM_001040113:c.5819del(p.Pro1940Hisfs*91)]. Both variants are predicted to result in a similarly elongated protein product. Heterozygous dominant negative MYH11 pathogenic variants have been associated with thoracic aortic aneurysm and dissection while biallelic null alleles have been associated with megacystis microcolon intestinal hypoperistalsis syndrome. This report highlights heterozygous protein-elongating MYH11 variants affecting the SM2 isoforms of MYH11 as a cause for severe gastrointestinal dysmotility, and we hypothesize that the mechanistic pathogenesis of this disease, dominant hypercontractile loss-of-function, is distinct from those implicated in other diseases involving MYH11 dysfunction.

Association of Enterovirus D68 with Acute Flaccid Myelitis, Philadelphia, Pennsylvania, USA, 2009-2018.

Acute flaccid myelitis (AFM) is a polio-like disease that results in paralysis in previously healthy persons. Although the definitive cause of AFM remains unconfirmed, enterovirus D68 (EV-D68) is suspected based on 2014 data demonstrating an increase in AFM cases concomitant with an EV-D68 outbreak. We examined the prevalence in children and the molecular evolution of EV-D68 for 2009-2018 in Philadelphia, Pennsylvania, USA. We detected widespread EV-D68 circulation in 2009, rare detections in 2010 and 2011, and then biennial circulation, only in even years, during 2012-2018. Prevalence of EV-D68 significantly correlated with AFM cases during this period. Finally, whole-genome sequencing revealed early detection of the B1 clade in 2009 and continued evolution of the B3 clade from 2016 to 2018. These data reinforce the need to improve surveillance programs for nonpolio enterovirus to identify possible AFM triggers and predict disease prevalence to better prepare for future outbreaks.

Alagille syndrome mutation update: Comprehensive overview of JAG1 and NOTCH2 mutation frequencies and insight into missense variant classification.

Alagille syndrome is an autosomal dominant disease with a known molecular etiology of dysfunctional Notch signaling caused primarily by pathogenic variants in JAGGED1 (JAG1), but also by variants in NOTCH2. The majority of JAG1 variants result in loss of function, however disease has also been attributed to lesser understood missense variants. Conversely, the majority of NOTCH2 variants are missense, though fewer of these variants have been described. In addition, there is a small group of patients with a clear clinical phenotype in the absence of a pathogenic variant. Here, we catalog our single-center study, which includes 401 probands and 111 affected family members amassed over a 27-year period, to provide updated mutation frequencies in JAG1 and NOTCH2 as well as functional validation of nine missense variants. Combining our cohort of 86 novel JAG1 and three novel NOTCH2 variants with previously published data (totaling 713 variants), we present the most comprehensive pathogenic variant overview for Alagille syndrome. Using this data set, we developed new guidance to help with the classification of JAG1 missense variants. Finally, we report clinically consistent cases for which a molecular etiology has not been identified and discuss the potential for next generation sequencing methodologies in novel variant discovery.

Genome maps across 26 human populations reveal population-specific patterns of structural variation.

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

Identification of Polycystic Kidney Disease 1 Like 1 Gene Variants in Children With Biliary Atresia Splenic Malformation Syndrome.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

Utility and limitations of exome sequencing as a genetic diagnostic tool for children with hearing loss.

PURPOSE: Hearing loss (HL) is the most common sensory disorder in children. Prompt molecular diagnosis may guide screening and management, especially in syndromic cases when HL is the single presenting feature. Exome sequencing (ES) is an appealing diagnostic tool for HL as the genetic causes are highly heterogeneous. METHODS: ES was performed on a prospective cohort of 43 probands with HL. Sequence data were analyzed for primary and secondary findings. Capture and coverage analysis was performed for genes and variants associated with HL. RESULTS: The diagnostic rate using ES was 37.2%, compared with 15.8% for the clinical HL panel. Secondary findings were discovered in three patients. For 247 genes associated with HL, 94.7% of the exons were targeted for capture and 81.7% of these exons were covered at 20× or greater. Further analysis of 454 randomly selected HL-associated variants showed that 89% were targeted for capture and 75% were covered at a read depth of at least 20×. CONCLUSION: ES has an improved yield compared with clinical testing and may capture diagnoses not initially considered due to subtle clinical phenotypes. Technical challenges were identified, including inadequate capture and coverage of HL genes. Additional considerations of ES include secondary findings, cost, and turnaround time.

NIPBL+/- haploinsufficiency reveals a constellation of transcriptome disruptions in the pluripotent and cardiac states.

Cornelia de Lange syndrome (CdLS) is a complex disorder with multiple structural and developmental defects caused by mutations in structural and regulatory proteins involved in the cohesin complex. NIPBL, a cohesin regulatory protein, has been identified as a critical protein responsible for the orchestration of transcriptomic regulatory networks necessary for embryonic development. Mutations in NIPBL are responsible for the majority of cases of CdLS. Through RNA-sequencing of human induced pluripotent stem cells and in vitro-derived cardiomyocytes, we identified hundreds of mRNAs, pseudogenes, and non-coding RNAs with altered expression in NIPBL+/- patient-derived cells. We demonstrate that NIPBL haploinsufficiency leads to upregulation of gene sets identified in functions related to nucleosome, chromatin assembly, RNA modification and downregulation of Wnt signaling, cholesterol biosynthesis and vesicular transport in iPSC and cardiomyocytes. Mutations in NIPBL result in the dysregulation of many genes responsible for normal heart development likely resulting in the variety of structural cardiac defects observed in the CdLS population.

Utility and limitations of exome sequencing in the molecular diagnosis of pediatric inherited platelet disorders.

Inherited platelet disorders (IPD) are a heterogeneous group of rare disorders that affect platelet number and function and often predispose to other significant medical complications. In spite of the identification of over 50 IPD disease-associated genes, a molecular diagnosis is only identified in a minority (10%) of affected patients without a clinically suspected etiology. We studied a cohort of 21 pediatric patients with suspected IPDs by exome sequencing (ES) to: (1) examine the performance of the exome test for IPD genes, (2) determine if this exome-wide diagnostic test provided a higher diagnostic yield than has been previously reported, (3) to evaluate the frequency of variants of uncertain significance identified, and (4) to identify candidate variants for functional evaluation in patients with an uncertain or negative diagnosis. We established a high priority gene list of 53 genes, evaluated exome capture kit performance, and determined the coverage for these genes and disease-related variants. We identified likely disease causing variants in 5 of the 21 probands (23.8%) and variants of uncertain significance in 52% of patients studied. In conclusion, ES has the potential to molecularly diagnose causes of IPD, and to identify candidate genes for functional evaluation. Robust exome sequencing also requires that coverage of genes known to be associated with clinical findings of interest need to be carefully examined and supplemented if necessary. Clinicians who undertake ES should understand the limitations of the test and the full significance of results that may be returned.

Flexible robust binder-free carbon nanotube membranes for solid state and microcapacitor application.

We present a liquid phase post synthesis self-assemble protocol that transforms trillions of carbon nanotubes (CNTs) in powder form into densely packed flexible, robust and binder-free macroscopic membranes with a hierarchical pore structure. We employ charge transfer engineering to spontaneously disperse the CNTs in a liquid medium. The processing protocol has limited or no impact on the intrinsic properties of the CNTs. As the thickness of the CNT membrane is increased, we observed a gradual transition from high flexibility to buckling and brittleness in the flexural properties of the membranes. The binder-free CNT membranes have bulk mass density greater than that of water (1.0 g cm-3). We correlate the mass of the CNTs in the membrane to the thickness of the membrane and obtained a bulk mass density of ∼1.11 g cm-3 ± 0.03 g cm-3. We demonstrate the use of the CNT membranes as electrode in a pristine and oxidized single/stacked solid-state capacitor as well as pristine interdigitated microcapacitor that show time constant of ∼32 ms with no degradation in performance even after 10 000 cycles. The capacitors show very good temperature dependence over a wide range of temperatures with good cycling performance up to 90 °C. The specific capacitance of the pseudocapacitive CNT electrode at room temperature was 72 F g-1 and increased to 100 F g-1 at 70 °C. The leakage current of bipolar stacked solid state capacitor was ∼100 nA cm-2 at 2.5 V when held for 72 h.

Neurodevelopmental Disorders Caused by De Novo Variants in KCNB1 Genotypes and Phenotypes.

Importance: Knowing the range of symptoms seen in patients with a missense or loss-of-function variant in KCNB1 and how these symptoms correlate with the type of variant will help clinicians with diagnosis and prognosis when treating new patients. Objectives: To investigate the clinical spectrum associated with KCNB1 variants and the genotype-phenotype correlations. Design, Setting, and Participants: This study summarized the clinical and genetic information of patients with a presumed pathogenic variant in KCNB1. Patients were identified in research projects or during clinical testing. Information on patients from previously published articles was collected and authors contacted if feasible. All patients were seen at a clinic at one of the participating institutes because of presumed genetic disorder. They were tested in a clinical setting or included in a research project. Main Outcomes and Measures: The genetic variant and its inheritance and information on the patient's symptoms and characteristics in a predefined format. All variants were identified with massive parallel sequencing and confirmed with Sanger sequencing in the patient. Absence of the variant in the parents could be confirmed with Sanger sequencing in all families except one. Results: Of 26 patients (10 female, 15 male, 1 unknown; mean age at inclusion, 9.8 years; age range, 2-32 years) with developmental delay, 20 (77%) carried a missense variant in the ion channel domain of KCNB1, with a concentration of variants in region S5 to S6. Three variants that led to premature stops were located in the C-terminal and 3 in the ion channel domain. Twenty-one of 25 patients (84%) had seizures, with 9 patients (36%) starting with epileptic spasms between 3 and 18 months of age. All patients had developmental delay, with 17 (65%) experiencing severe developmental delay; 14 (82%) with severe delay had behavioral problems. The developmental delay was milder in 4 of 6 patients with stop variants and in a patient with a variant in the S2 transmembrane element rather than the S4 to S6 region. Conclusions and Relevance: De novo KCNB1 missense variants in the ion channel domain and loss-of-function variants in this domain and the C-terminal likely cause neurodevelopmental disorders with or without seizures. Patients with presumed pathogenic variants in KCNB1 have a variable phenotype. However, the type and position of the variants in the protein are (imperfectly) correlated with the severity of the disorder.

High-throughput single-molecule mapping links subtelomeric variants and long-range haplotypes with specific telomeres.

Accurate maps and DNA sequences for human subtelomere regions, along with detailed knowledge of subtelomere variation and long-range telomere-terminal haplotypes in individuals, are critical for understanding telomere function and its roles in human biology. Here, we use a highly automated whole genome mapping technology in nano-channel arrays to analyze large terminal human chromosome segments extending from chromosome-specific subtelomere sequences through subtelomeric repeat regions to terminal (TTAGGG)n repeat tracts. We establish detailed maps for subtelomere gap regions in the human reference sequence, detect many new large subtelomeric variants and demonstrate the feasibility of long-range haplotyping through segmentally duplicated subtelomere regions. These features make the method a uniquely valuable new tool for improving the quality of genome assemblies in complex DNA regions. Based on single molecule mapping of telomere-terminal DNA fragments, we provide proof of principle for a novel method to estimate telomere lengths linked to distinguishable telomeric haplotypes; this single-telomere genotyping method may ultimately enable delineation of human cis elements involved in telomere length regulation.