High-throughput single-molecule telomere characterization

Techniques for assessing telomere length: A methodological review

Telomeres are located at the ends of chromosomes and have specific sequences with a distinctive structure that safeguards genes. They possess capping structures that protect chromosome ends from fusion events and ensure chromosome stability. Telomeres shorten in length during each cycle of cell division. When this length reaches a certain threshold, it can lead to genomic instability, thus being implicated in various diseases, including cancer and neurodegenerative disorders. The possibility of telomeres serving as a biomarker for aging and age-related disease is being explored, and their significance is still under study. This is because post-mitotic cells, which are mature cells that do not undergo mitosis, do not experience telomere shortening due to age. Instead, other causes, for example, exposure to oxidative stress, can directly damage the telomeres, causing genomic instability. Nonetheless, a general agreement has been established that measuring telomere length offers valuable insights and forms a crucial foundation for analyzing gene expression and epigenetic data. Numerous approaches have been developed to accurately measure telomere lengths. In this review, we summarize various methods and their advantages and limitations for assessing telomere length.

Cholinergic Mechanisms in Gastrointestinal Neoplasia

Acetylcholine-activated receptors are divided broadly into two major structurally distinct classes: ligand-gated ion channel nicotinic and G-protein-coupled muscarinic receptors. Each class encompasses several structurally related receptor subtypes with distinct patterns of tissue expression and post-receptor signal transduction mechanisms. The activation of both nicotinic and muscarinic cholinergic receptors has been associated with the induction and progression of gastrointestinal neoplasia. Herein, after briefly reviewing the classification of acetylcholine-activated receptors and the role that nicotinic and muscarinic cholinergic signaling plays in normal digestive function, we consider the mechanics of acetylcholine synthesis and release by neuronal and non-neuronal cells in the gastrointestinal microenvironment, and current methodology and challenges in measuring serum and tissue acetylcholine levels accurately. Then, we critically evaluate the evidence that constitutive and ligand-induced activation of acetylcholine-activated receptors plays a role in promoting gastrointestinal neoplasia. We focus primarily on adenocarcinomas of the stomach, pancreas, and colon, because these cancers are particularly common worldwide and, when diagnosed at an advanced stage, are associated with very high rates of morbidity and mortality. Throughout this comprehensive review, we concentrate on identifying novel ways to leverage these observations for prognostic and therapeutic purposes.

Elucidation of the molecular mechanism of the breakage-fusion-bridge (BFB) cycle using a CRISPR-dCas9 cellular model

Chromosome instability (CIN) is frequently observed in many tumors. The breakage-fusion-bridge (BFB) cycle has been proposed to be one of the main drivers of CIN during tumorigenesis and tumor evolution. However, the detailed mechanisms for the individual steps of the BFB cycle warrants further investigation. Here, we demonstrated that a nuclease-dead Cas9 (dCas9) coupled with a telomere-specific single-guide RNA (sgTelo) can be used to model the BFB cycle. First, we showed that targeting dCas9 to telomeres using sgTelo impeded DNA replication at telomeres and induced a pronounced increase of replication stress and DNA damage. Using Single-Molecule Telomere Assay via Optical Mapping (SMTA-OM), we investigated the genome-wide features of telomeres in the dCas9/sgTelo cells and observed a dramatic increase of chromosome end fusions, including fusion/ITS+ and fusion/ITS-.Consistently, we also observed an increase in the formation of dicentric chromosomes, anaphase bridges, and intercellular telomeric chromosome bridges (ITCBs). Utilizing the dCas9/sgTelo system, we uncovered many novel molecular and structural features of the ITCB and demonstrated that multiple DNA repair pathways are implicated in the formation of ITCBs. Our studies shed new light on the molecular mechanisms of the BFB cycle, which will advance our understanding of tumorigenesis, tumor evolution, and drug resistance.

A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Fluorogenic CRISPR for genomic DNA imaging

Genomic DNA exhibits high heterogeneity in terms of its dynamic within the nucleus, its structure and functional roles. CRISPR-based imaging approaches can image genomic loci in living cells. However, conventional CRISPR-based tools involve expressing constitutively fluorescent proteins, resulting in high background and nonspecific nucleolar signal. Here, we construct fluorogenic CRISPR (fCRISPR) to overcome these issues. fCRISPR is designed with dCas9, an engineered sgRNA, and a fluorogenic protein. Fluorogenic proteins are degraded unless they are bound to specific RNA hairpins. These hairpins are inserted into sgRNA, resulting in dCas9: sgRNA: fluorogenic protein ternary complexes that enable fluorogenic DNA imaging. With fCRISPR, we image various genomic DNA in different human cells with high signal-to-noise ratio and sensitivity. Furthermore, fCRISPR tracks chromosomes dynamics and length. fCRISPR also allows DNA double-strand breaks (DSBs) and repair to be tracked in real time. Taken together, fCRISPR offers a high-contrast and sensitive platform for imaging genomic loci.

Single-Molecule Telomere Assay via Optical Mapping (SMTA-OM) Can Potentially Define the ALT Positivity of Cancer

Telomeres play an essential role in protecting the ends of linear chromosomes and maintaining the integrity of the human genome. One of the key hallmarks of cancers is their replicative immortality. As many as 85-90% of cancers activate the expression of telomerase (TEL+) as the telomere maintenance mechanism (TMM), and 10-15% of cancers utilize the homology-dependent repair (HDR)-based Alternative Lengthening of Telomere (ALT+) pathway. Here, we performed statistical analysis of our previously reported telomere profiling results from Single Molecule Telomere Assay via Optical Mapping (SMTA-OM), which is capable of quantifying individual telomeres from single molecules across all chromosomes. By comparing the telomeric features from SMTA-OM in TEL+ and ALT+ cancer cells, we demonstrated that ALT+ cancer cells display certain unique telomeric profiles, including increased fusions/internal telomere-like sequence (ITS+), fusions/internal telomere-like sequence loss (ITS-), telomere-free ends (TFE), super-long telomeres, and telomere length heterogeneity, compared to TEL+ cancer cells. Therefore, we propose that ALT+ cancer cells can be differentiated from TEL+ cancer cells using the SMTA-OM readouts as biomarkers. In addition, we observed variations in SMTA-OM readouts between different ALT+ cell lines that may potentially be used as biomarkers for discerning subtypes of ALT+ cancer and monitoring the response to cancer therapy.

Detection of Transgene Location in the CYP2A13/2B6/2F1-transgenic Mouse Model using Optical Genome Mapping Technology

Most transgenic mouse models are generated through random integration of the transgene. The location of the transgene provides valuable information for assessing potential effects of the transgenesis on the host and for designing genotyping protocols that can amplify across the integration site, but it is challenging to identify. Here, we report the successful utility of optical genome mapping technology to identify the transgene insertion site in a CYP2A13/2B6/2F1-transgenic mouse model, which produces three human cytochrome P450 (P450) enzymes (CYP2A13, CYP2B6, and CYP2F1) that are encoded by neighboring genes on human chromosome 19. These enzymes metabolize many drugs, respiratory toxicants, and chemical carcinogens. Initial efforts to identify candidate insertion sites by whole genome sequencing was unsuccessful, apparently because the transgene is located in a region of the mouse genome that contains highly repetitive sequences. Subsequent utility of the optical genome mapping approach, which compares genome-wide marker distribution between the transgenic mouse genome and a reference mouse (GRCm38) or human (GRCh38) genome, localized the insertion site to mouse chromosome 14, between two marker positions at 4451324 base pair and 4485032 base pair. A transgene-mouse genome junction sequence was further identified through long-polymerase chain reaction amplification and DNA sequencing at GRCm38 Chr.14:4484726. The transgene insertion (∼2.4 megabase pair) contained 5-7 copies of the human transgenes, which replaced a 26.9-33.4 kilobase pair mouse genomic region, including exons 1-4 of Gm3182, a predicted and highly redundant gene. Finally, the sequencing results enabled the design of a new genotyping protocol that can distinguish between hemizygous and homozygous CYP2A13/2B6/2F1-transgenic mice. SIGNIFICANCE STATEMENT: This study characterizes the genomic structure of, and provides a new genotyping method for, a transgenic mouse model that expresses three human P450 enzymes, CYP2A13, CYP2B6, and CYP2F1, that are important in xenobiotic metabolism and toxicity. The demonstrated success in applying the optical genome mapping technology for identification of transgene insertion sites should encourage others to do the same for other transgenic models generated through random integration, including most of the currently available human P450 transgenic mouse models.

NPGREAT: assembly of human subtelomere regions with the use of ultralong nanopore reads and linked-reads

BACKGROUND

Human subtelomeric DNA regulates the length and stability of adjacent telomeres that are critical for cellular function, and contains many gene/pseudogene families. Large evolutionarily recent segmental duplications and associated structural variation in human subtelomeres has made complete sequencing and assembly of these regions difficult to impossible for many loci, complicating or precluding a wide range of genetic analyses to investigate their function.

RESULTS

We present a hybrid assembly method, NanoPore Guided REgional Assembly Tool (NPGREAT), which combines Linked-Read data with mapped ultralong nanopore reads spanning subtelomeric segmental duplications to potentially overcome these difficulties. Linked-Read sets of DNA sequences identified by matches with 1-copy subtelomere sequence adjacent to segmental duplications are assembled and extended into the segmental duplication regions using Regional Extension of Assemblies using Linked-Reads (REXTAL). Mapped telomere-containing ultralong nanopore reads are then used to provide contiguity and correct orientation for matching REXTAL sequence contigs as well as identification/correction of any misassemblies. Our method was tested for a subset of representative subtelomeres with ultralong nanopore read coverage in the haploid human cell line CHM13. A 10X Linked-Read dataset from CHM13 was combined with ultralong nanopore reads from the same genome to provide improved subtelomere assemblies. Comparison of Nanopore-only assemblies using SHASTA with our NPGREAT assemblies in the distal-most subtelomere regions showed that NPGREAT produced higher-quality and more complete assemblies than SHASTA alone when these regions had low ultralong nanopore coverage (such as cases where large segmental duplications were immediately adjacent to (TTAGGG) tracts).

CONCLUSION

In genomic regions with large segmental duplications adjacent to telomeres, NPGREAT offers an alternative economical approach to improving assembly accuracy and coverage using linked-read datasets when more expensive HiFi datasets of 10-20 kb reads are unavailable.

Multiplex structural variant detection by whole-genome mapping and nanopore sequencing

Identification of structural variants (SVs) breakpoints is important in studying mutations, mutagenic causes, and functional impacts. Next-generation sequencing and whole-genome optical mapping are extensively used in SV discovery and characterization. However, multiple platforms and computational approaches are needed for comprehensive analysis, making it resource-intensive and expensive. Here, we propose a strategy combining optical mapping and cas9-assisted targeted nanopore sequencing to analyze SVs. Optical mapping can economically and quickly detect SVs across a whole genome but does not provide sequence-level information or precisely resolve breakpoints. Furthermore, since only a subset of all SVs is known to affect biology, we attempted to type a subset of all SVs using targeted nanopore sequencing. Using our approach, we resolved the breakpoints of five deletions, five insertions, and an inversion, in a single experiment.

Characterization of full-length LINE-1 insertions in 154 genomes

Long interspersed nuclear elements (LINEs) are retrotransposons that contribute to genetic variation in the human genome. LINE-1 elements in larger-scale studies are challenging to identify using sequencing technologies due to cost and scalability. We developed an approach using optical mapping for detection of full-length LINE-1 insertions and 10× sequencing for confirmation. We found 51 true positive full-length LINE-1 insertions, of which 4 are novel insertions, in NA12878. Repeating our analysis on a larger sample set representing 26 populations, we identified 329 full-length LINE-1 elements, of which 123 are novel. 24.8% of these 329 LINE-1 insertions were shared amongst all 5 superpopulations (AFR, AMR, EUR, EAS, SAS). The African superpopulation has a higher percentage of population-specific LINE-1 insertions than any other superpopulation. These data indicate that our approach can provide high-speed, cost-effective, and increased accuracy for LINE-1 detection. These data also provide an insight into variations of LINE-1 elements between different populations.

Multicolor Whole-Genome Mapping in Nanochannels for Genetic Analysis

Analysis of structural variations (SVs) is important to understand mutations underlying genetic disorders and pathogenic conditions. However, characterizing SVs using short-read, high-throughput sequencing technology is difficult. Although long-read sequencing technologies are being increasingly employed in characterizing SVs, their low throughput and high costs discourage widespread adoption. Sequence motif-based optical mapping in nanochannels is useful in whole-genome mapping and SV detection, but it is not possible to precisely locate the breakpoints or estimate the copy numbers. We present here a universal multicolor mapping strategy in nanochannels combining conventional sequence-motif labeling system with Cas9-mediated target-specific labeling of any 20-base sequences (20mers) to create custom labels and detect new features. The sequence motifs are labeled with green fluorophores and the 20mers are labeled with red fluorophores. Using this strategy, it is possible to not only detect the SVs but also utilize custom labels to interrogate the features not accessible to motif-labeling, locate breakpoints, and precisely estimate copy numbers of genomic repeats. We validated our approach by quantifying the D4Z4 copy numbers, a known biomarker for facioscapulohumeral muscular dystrophy (FSHD) and estimating the telomere length, a clinical biomarker for assessing disease risk factors in aging-related diseases and malignant cancers. We also demonstrate the application of our methodology in discovering transposable long non-interspersed Elements 1 (LINE-1) insertions across the whole genome.

Single-molecule telomere length characterization by optical mapping in nano-channel array: Perspective and review on telomere length measurement

In humans, the telomere consists of tandem 5'TTAGGG3' DNA repeats on both ends of all 46 chromosomes. Telomere shortening has been linked to aging and age-related diseases. Similarly, telomere length changes have been associated with chemical exposure, molecular-level DNA damage, and tumor development. Telomere elongation has been associated to tumor development, caused due to chemical exposure and molecular-level DNA damage. The methods used to study these effects mostly rely on average telomere length as a biomarker. The mechanisms regulating subtelomere-specific and haplotype-specific telomere lengths in humans remain understudied and poorly understood, primarily because of technical limitations in obtaining these data for all chromosomes. Recent studies have shown that it is the short telomeres that are crucial in preserving chromosome stability. The identity and frequency of specific critically short telomeres potentially is a useful biomarker for studying aging, age-related diseases, and cancer. Here, we will briefly review the role of telomere length, its measurement, and our recent single-molecule telomere length measurement assay. With this assay, one can measure individual telomere lengths as well as identify their physically linked subtelomeric DNA. This assay can also positively detect telomere loss, characterize novel subtelomeric variants, haplotypes, and previously uncharacterized recombined subtelomeres. We will also discuss its applications in aging cells and cancer cells, highlighting the utility of the single molecule telomere length assay.

Customized optical mapping by CRISPR-Cas9 mediated DNA labeling with multiple sgRNAs

Whole-genome mapping technologies have been developed as a complementary tool to provide scaffolds for genome assembly and structural variation analysis (1,2). We recently introduced a novel DNA labeling strategy based on a CRISPR-Cas9 genome editing system, which can target any 20bp sequences. The labeling strategy is specifically useful in targeting repetitive sequences, and sequences not accessible to other labeling methods. In this report, we present customized mapping strategies that extend the applications of CRISPR-Cas9 DNA labeling. We first design a CRISPR-Cas9 labeling strategy to interrogate and differentiate the single allele differences in NGG protospacer adjacent motifs (PAM sequence). Combined with sequence motif labeling, we can pinpoint the single-base differences in highly conserved sequences. In the second strategy, we design mapping patterns across a genome by selecting sets of specific single-guide RNAs (sgRNAs) for labeling multiple loci of a genomic region or a whole genome. By developing and optimizing a single tube synthesis of multiple sgRNAs, we demonstrate the utility of CRISPR-Cas9 mapping with 162 sgRNAs targeting the 2Mb Haemophilus influenzae chromosome. These CRISPR-Cas9 mapping approaches could be particularly useful for applications in defining long-distance haplotypes and pinpointing the breakpoints in large structural variants in complex genomes and microbial mixtures.

Method comparison studies of telomere length measurement using qPCR approaches: A critical appraisal of the literature

Use of telomere length (TL) as a biomarker for various environmental exposures and diseases has increased in recent years. Various methods have been developed to measure telomere length. Polymerase chain reaction (PCR)-based methods remain wide-spread for population-based studies due to the high-throughput capability. While several studies have evaluated the repeatability and reproducibility of different TL measurement methods, the results have been variable. We conducted a literature review of TL measurement cross-method comparison studies that included a PCR-based method published between January 1, 2002 and May 25, 2020. A total of 25 articles were found that matched the inclusion criteria. Papers were reviewed for quality of methodologic reporting of sample and DNA quality, PCR assay characteristics, sample blinding, and analytic approaches to determine final TL. Overall, methodologic reporting was low as assessed by two different reporting guidelines for qPCR-based TL measurement. There was a wide range in the reported correlation between methods (as assessed by Pearson’s r) and few studies utilized the recommended intra-class correlation coefficient (ICC) for assessment of assay repeatability and methodologic comparisons. The sample size for nearly all studies was less than 100, raising concerns about statistical power. Overall, this review found that the current literature on the relation between TL measurement methods is lacking in validity and scientific rigor. In light of these findings, we present reporting guidelines for PCR-based TL measurement methods and results of analyses of the effect of assay repeatability (ICC) on statistical power of cross-sectional and longitudinal studies. Additional cross-laboratory studies with rigorous methodologic and statistical reporting, adequate sample size, and blinding are essential to accurately determine assay repeatability and replicability as well as the relation between TL measurement methods.

Advances in optical mapping for genomic research

Recent advances in optical mapping have allowed the construction of improved genome assemblies with greater contiguity. Optical mapping also enables genome comparison and identification of large-scale structural variations. Association of these large-scale genomic features with biological functions is an important goal in plant and animal breeding and in medical research. Optical mapping has also been used in microbiology and still plays an important role in strain typing and epidemiological studies. Here, we review the development of optical mapping in recent decades to illustrate its importance in genomic research. We detail its applications and algorithms to show its specific advantages. Finally, we discuss the challenges required to facilitate the optimization of optical mapping and improve its future development and application.

Single-molecule analysis of subtelomeres and telomeres in Alternative Lengthening of Telomeres (ALT) cells

BACKGROUND

Telomeric DNA is typically comprised of G-rich tandem repeat motifs and maintained by telomerase (Greider CW, Blackburn EH; Cell 51:887-898; 1987). In eukaryotes lacking telomerase, a variety of DNA repair and DNA recombination based pathways for telomere maintenance have evolved in organisms normally dependent upon telomerase for telomere elongation (Webb CJ, Wu Y, Zakian VA; Cold Spring Harb Perspect Biol 5:a012666; 2013); collectively called Alternative Lengthening of Telomeres (ALT) pathways. By measuring (TTAGGG) n tract lengths from the same large DNA molecules that were optically mapped, we simultaneously analyzed telomere length dynamics and subtelomere-linked structural changes at a large number of specific subtelomeric loci in the ALT-positive cell lines U2OS, SK-MEL-2 and Saos-2.

RESULTS

Our results revealed loci-specific ALT telomere features. For example, while each subtelomere included examples of single molecules with terminal (TTAGGG) n tracts as well as examples of recombinant telomeric single molecules, the ratio of these molecules was subtelomere-specific, ranging from 33:1 (19p) to 1:25 (19q) in U2OS. The Saos-2 cell line shows a similar percentage of recombinant telomeres. The frequency of recombinant subtelomeres of SK-MEL-2 (11%) is about half that of U2OS and Saos-2 (24 and 19% respectively). Terminal (TTAGGG) n tract lengths and heterogeneity levels, the frequencies of telomere signal-free ends, and the frequency and size of retained internal telomere-like sequences (ITSs) at recombinant telomere fusion junctions all varied according to the specific subtelomere involved in a particular cell line. Very large linear extrachromosomal telomere repeat (ECTR) DNA molecules were found in all three cell lines; these are in principle capable of templating synthesis of new long telomere tracts via break-induced repair (BIR) long-tract DNA synthesis mechanisms and contributing to the very long telomere tract length and heterogeneity characteristic of ALT cells. Many of longest telomere tracts (both end-telomeres and linear ECTRs) displayed punctate CRISPR/Cas9-dependent (TTAGGG) n labeling patterns indicative of interspersion of stretches of non-canonical telomere repeats.

CONCLUSION

Identifying individual subtelomeres and characterizing linked telomere (TTAGGG) n tract lengths and structural changes using our new single-molecule methodologies reveals the structural consequences of telomere damage, repair and recombination mechanisms in human ALT cells in unprecedented molecular detail and significant differences in different ALT-positive cell lines.

Telomere Length Measurement by Molecular Combing

Telomeres are repetitive regions of DNA bound by specialized proteins at the termini of linear chromosomes that prevent the natural chromosome ends from being recognized as DNA double strand breaks. Telomeric DNA is gradually eroded with each round of cell division, resulting in the accumulation of critically short or dysfunctional telomeres that eventually trigger cellular senescence. Consequently, telomere length is indicative of the proliferative capacity of a cell. Multiple methods exist to measure telomere length and telomere content, but a simple and reliable technique to accurately measure individual telomere lengths is currently lacking. We have developed the Telomere length Combing Assay (TCA) to measure telomere length on stretched DNA fibers. We used TCA to measure telomere erosion in primary human fibroblasts, and to detect telomere lengthening in response to activation of telomere maintenance pathways. TCA was also used to accurately measure telomere length in healthy individuals, and to identify critically short telomeres in patients with telomere biology disorders. TCA is performed on isolated DNA, negating the need for cycling cells. TCA is amenable to semi-automated image analysis, and can be fully automated using the Genomic Vision molecular combing platform. This not only precludes sampling bias, but also provides the potential for high-throughput applications and clinical development. TCA is a simple and versatile technique to measure the distribution of individual telomere lengths in a cell population, offering improved accuracy, and more detailed biological insight for telomere length measurement applications.

Human subtelomeric DNA methylation: regulation and roles in telomere function

Subtelomeres are the regions at chromosome ends, immediately adjacent to the terminal telomeric repeats. The majority of human subtelomeres are CpG-rich in their distal two kilobases, and are methylated during early embryonic development by the de novo DNA methyltransferase DNMT3B. The biological relevance of subtelomeric DNA methylation is highlighted by the presence of promoters for the long non-coding TERRA transcripts in these CpG-rich regions. Indeed, deviant subtelomeric methylation has been linked with abnormal telomeric phenotypes, as most strikingly found in ICF syndrome. Here we review recent studies that explore new aspects of subtelomeric methylation regulation and demonstrate the significance of maintaining proper DNA methylation at the extreme distal human subtelomeric regions.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Comprehensive Analysis of Human Subtelomeres by Whole Genome Mapping

Detailed comprehensive knowledge of the structures of individual long-range telomere-terminal haplotypes are needed to understand their impact on telomere function, and to delineate the population structure and evolution of subtelomere regions. However, the abundance of large evolutionarily recent segmental duplications and high levels of large structural variations have complicated both the mapping and sequence characterization of human subtelomere regions. Here, we use high throughput optical mapping of large single DNA molecules in nanochannel arrays for 154 human genomes from 26 populations to present a comprehensive look at human subtelomere structure and variation. The results catalog many novel long-range subtelomere haplotypes and determine the frequencies and contexts of specific subtelomeric duplicons on each chromosome arm, helping to clarify the currently ambiguous nature of many specific subtelomere structures as represented in the current reference sequence (HG38). The organization and content of some duplicons in subtelomeres appear to show both chromosome arm and population-specific trends. Based upon these trends we estimate a timeline for the spread of these duplication blocks.

A micropatterned substrate for on-surface enzymatic labelling of linearized long DNA molecules

Optical mapping of linearized DNA molecules is a promising new technology for sequence assembly and scaffolding, large structural variant detection, and diagnostics. This is currently achieved either using nanochannel confinement or by stretching single DNA molecules on a solid surface. While the first method necessitates DNA labelling before linearization, the latter allows for modification post-linearization, thereby affording increased process flexibility. Each method is constrained by various physical and chemical limitations. One of the most common techniques for linearization of DNA uses a hydrophobic surface and a receding meniscus, termed molecular combing. Here, we report the development of a microfabricated surface that can not only comb the DNA molecules efficiently but also provides for sequence-specific enzymatic fluorescent DNA labelling. By modifying a glass surface with two contrasting functionalities, such that DNA binds selectively to one of the two regions, we can control DNA extension, which is known to be critical for sequence-recognition by an enzyme. Moreover, the surface modification provides enzymatic access to the DNA backbone, as well as minimizing non-specific fluorescent dye adsorption. These enhancements make the designed surface suitable for large-scale and high-resolution single DNA molecule studies.

Nanopore Guided Assembly of Segmental Duplications near Telomeres

Human subtelomere regions are highly enriched in large segmental duplications and structural variants, leading to many gaps and misassemblies in these regions. We develop a novel method, NPGREAT (NanoPore Guided REgional Assembly Tool), which combines Nanopore ultralong read datasets and short-read assemblies derived from 10x linked-reads to efficiently assemble these subtelomere regions into a single continuous sequence. We show that with the use of ultralong Nanopore reads as a guide, the highly accurate shorter linked-read sequence contigs are correctly oriented, ordered, spaced and extended. In the rare cases where a linked-read sequence contig contains inaccurately assembled segments, the use of Nanopore reads allows for detection and correction of this error. We tested NPGREAT on four representative subtelomeres of the NA12878 human genome (10p, 16p, 19q and 20p). The results demonstrate that the final computed assembly of each subtelomere is accurate and complete.

OMMA enables population-scale analysis of complex genomic features and phylogenomic relationships from nanochannel-based optical maps

BACKGROUND

Optical mapping is an emerging technology that complements sequencing-based methods in genome analysis. It is widely used in improving genome assemblies and detecting structural variations by providing information over much longer (up to 1 Mb) reads. Current standards in optical mapping analysis involve assembling optical maps into contigs and aligning them to a reference, which is limited to pairwise comparison and becomes bias-prone when analyzing multiple samples.

FINDINGS

We present a new method, OMMA, that extends optical mapping to the study of complex genomic features by simultaneously interrogating optical maps across many samples in a reference-independent manner. OMMA captures and characterizes complex genomic features, e.g., multiple haplotypes, copy number variations, and subtelomeric structures when applied to 154 human samples across the 26 populations sequenced in the 1000 Genomes Project. For small genomes such as pathogenic bacteria, OMMA accurately reconstructs the phylogenomic relationships and identifies functional elements across 21 Acinetobacter baumannii strains.

CONCLUSIONS

With the increasing data throughput of optical mapping system, the use of this technology in comparative genome analysis across many samples will become feasible. OMMA is a timely solution that can address such computational need. The OMMA software is available at https://github.com/TF-Chan-Lab/OMTools.

Recompleting the Caenorhabditis elegans genome

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted ≥53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology.

Deciphering the complexity of the 4q and 10q subtelomeres by molecular combing in healthy individuals and patients with facioscapulohumeral dystrophy

BACKGROUND

Subtelomeres are variable regions between telomeres and chromosomal-specific regions. One of the most studied pathologies linked to subtelomeric imbalance is facioscapulohumeral dystrophy (FSHD). In most cases, this disease involves shortening of an array of D4Z4 macrosatellite elements at the 4q35 locus. The disease also segregates with a specific A-type haplotype containing a degenerated polyadenylation signal distal to the last repeat followed by a repetitive array of β-satellite elements. This classification applies to most patients with FSHD. A subset of patients called FSHD2 escapes this definition and carries a mutation in the SMCHD1 gene. We also recently described patients carrying a complex rearrangement consisting of a cis-duplication of the distal 4q35 locus identified by molecular combing.

METHODS

Using this high-resolution technology, we further investigated the organisation of the 4q35 region linked to the disease and the 10q26 locus presenting with 98% of homology in controls and patients.

RESULTS

Our analyses reveal a broad variability in size of the different elements composing these loci highlighting the complexity of these subtelomeres and the difficulty for genomic assembly. Out of the 1029 DNA samples analysed in our centre in the last 7 years, we also identified 54 cases clinically diagnosed with FSHD carrying complex genotypes. This includes mosaic patients, patients with deletions of the proximal 4q region and 23 cases with an atypical chromosome 10 pattern, infrequently found in the control population and never reported before.

CONCLUSION

Overall, this work underlines the complexity of these loci challenging the diagnosis and genetic counselling for this disease.

Translating genomics to the clinical diagnosis of disorders/differences of sex development

The medical and psychosocial challenges faced by patients living with Disorders/Differences of Sex Development (DSD) and their families can be alleviated by a rapid and accurate diagnostic process. Clinical diagnosis of DSD is limited by a lack of standardization of anatomical and endocrine phenotyping and genetic testing, as well as poor genotype/phenotype correlation. Historically, DSD genes have been identified through positional cloning of disease-associated variants segregating in families and validation of candidates in animal and in vitro modeling of variant pathogenicity. Owing to the complexity of conditions grouped under DSD, genome-wide scanning methods are better suited for identifying disease causing gene variant(s) and providing a clinical diagnosis. Here, we review a number of established genomic tools (karyotyping, chromosomal microarrays and exome sequencing) used in clinic for DSD diagnosis, as well as emerging genomic technologies such as whole-genome (short-read) sequencing, long-read sequencing, and optical mapping used for novel DSD gene discovery. These, together with gene expression and epigenetic studies can potentiate the clinical diagnosis of DSD diagnostic rates and enhance the outcomes for patients and families.

Cascaded Electrochemiluminescence Signal Amplifier for the Detection of Telomerase Activity from Tumor Cells and Tissues

Telomerase is closely linked to the physiological transformation of tumor cells and is commonly overexpressed in most types of tumor cells. Therefore, telomerase has become a potential biomarker for the process of tumorigenesis, progression, prognosis and metastasis. Thus, it is important to develop a simple, accurate and reliable method for detecting telomerase activity. As a high signal-to-noise ratio mode, electrochemiluminescence (ECL) has been widely applied in the field of biomedical analysis. Here, our objective was to construct an improved ECL signal amplifier for the detection of telomerase activity. Methods: A cascaded ECL signal amplifier was constructed to detect telomerase activity with high selectivity via controllable construction of a lysine-based dendric Ru(bpy)₃ ²⁺ polymer (DRP). The sensitivity, specificity and performance index were simultaneously evaluated by standard substance and cell and tissue samples. Results: With this cascaded ECL signal amplifier, high sensitivities of 100, 50, and 100 cells for three tumor cell lines (A549, MCF7 and HepG2 cell lines) were simultaneously achieved, and desirable specificity was also obtained. Furthermore, the excellent performance of this platform was also demonstrated in the detection of telomerase in tumor cells and tissues. Conclusion: This cascaded ECL signal amplifier has the potential to be a technological innovation in the field of telomerase activity detection.

Evolutionary rates of mammalian telomere-stability genes correlate with karyotype features and female germline expression

Telomeres protect the ends of eukaryotic chromosomes and are essential for cell viability. In mammals, telomere dynamics vary with life history traits (e.g. body mass and longevity), suggesting differential selection depending on physiological characteristics. Telomeres, in analogy to centromeric regions, also represent candidate meiotic drivers and subtelomeric DNA evolves rapidly. We analyzed the evolutionary history of mammalian genes implicated in telomere homeostasis (TEL genes). We detected widespread positive selection and we tested two alternative hypotheses: (i) fast evolution is driven by changes in life history traits; (ii) a conflict with selfish DNA elements at the female meiosis represents the underlying selective pressure. By accounting for the phylogenetic relationships among mammalian species, we show that life history traits do not contribute to shape diversity of TEL genes. Conversely, the evolutionary rate of TEL genes correlates with expression levels during meiosis and episodes of positive selection across mammalian species are associated with karyotype features (number of chromosome arms). We thus propose a telomere drive hypothesis, whereby (sub)telomeres and telomere-binding proteins are engaged in an intra-genomic conflict similar to the one described for centromeres.

CRISPR-based strategies for studying regulatory elements and chromatin structure in mammalian gene control

The development of high-throughput methods has enabled the genome-wide identification of putative regulatory elements in a wide variety of mammalian cells at an unprecedented resolution. Extensive genomic studies have revealed the important role of regulatory elements and genetic variation therein in disease formation and risk. In most cases, there is only correlative evidence for the roles of these elements and non-coding changes within these elements in pathogenesis. With the advent of genome- and epigenome-editing tools based on the CRISPR technology, it is now possible to test the functional relevance of the regulatory elements and alterations on a genomic scale. Here, we review the various CRISPR-based strategies that have been developed to functionally validate the candidate regulatory elements in mammals as well as the non-coding genetic variants found to be associated with human disease. We also discuss how these synthetic biology tools have helped to elucidate the role of three-dimensional nuclear architecture and higher-order chromatin organization in shaping functional genome and controlling gene expression.