Next-generation sequencing as a tool for breakpoint analysis in rearrangements of the globin gene clusters

Third-generation sequencing identified a novel complex variant in a patient with rare alpha-thalassemia

BACKGROUND

Thalassemias represent some of the most common monogenic diseases worldwide and are caused by variations in human hemoglobin genes which disrupt the balance of synthesis between the alpha and beta globin chains. Thalassemia gene detection technology is the gold standard to achieve accurate detection of thalassemia, but in clinical practice, most of the tests are only for common genotypes, which can easily lead to missing or misdiagnosis of rare thalassemia genotypes.

CASE PRESENTATION

We present the case of an 18-year-old Chinese female with abnormal values of routine hematological indices who was admitted for genetic screening for thalassemia. Genomic DNA was extracted and used for the genetic assays. Gap polymerase chain reaction and agarose gel electrophoresis were performed to detect HBA gene deletions, while PCR-reverse dot blot hybridization was used to detect point mutations in the HBA and HBB genes. Next-generation sequencing and third-generation sequencing (TGS) were used to identify known and potentially novel genotypes of thalassemia. We identified a novel complex variant αHb WestmeadαHb Westmeadαanti3.7/-α3.7 in a patient with rare alpha-thalassemia.

CONCLUSIONS

Our study identified a novel complex variant that expands the thalassemia gene variants spectrum. Meanwhile, the study suggests that TGS could effectively improve the specificity of thalassemia gene detection, and has promising potential for the discovery of novel thalassemia genotypes, which could also improve the accuracy of genetic counseling. Couples who are thalassemia carriers have the opportunity to reduce their risk of having a child with thalassemia.

Next-Generation Sequencing and Emerging Technologies

Genetic sequencing technologies are evolving at a rapid pace with major implications for research and clinical practice. In this review, the authors provide an updated overview of next-generation sequencing (NGS) and emerging methodologies. NGS has tremendously improved sequencing output while being more time and cost-efficient in comparison to Sanger sequencing. The authors describe short-read sequencing approaches, such as sequencing by synthesis, ion semiconductor sequencing, and nanoball sequencing. Third-generation long-read sequencing now promises to overcome many of the limitations of short-read sequencing, such as the ability to reliably resolve repeat sequences and large genomic rearrangements. By combining complementary methods with massively parallel DNA sequencing, a greater insight into the biological context of disease mechanisms is now possible. Emerging methodologies, such as advances in nanopore technology, in situ nucleic acid sequencing, and microscopy-based sequencing, will continue the rapid evolution of this area. These new technologies hold many potential applications for hematological disorders, with the promise of precision and personalized medical care in the future.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding

BACKGROUND

Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy.

RESULT

CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population.

CONCLUSION

We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.

Identification of small-sized intrachromosomal segments at the ends of INV-DUP-DEL patterns

The mechanism of chromosomal rearrangement associated with inverted-duplication-deletion (INV-DUP-DEL) pattern formation has been investigated by many researchers, and several possible mechanisms have been proposed. Currently, fold-back and subsequent dicentric chromosome formation has been established as non-recurrent INV-DUP-DEL pattern formation mechanisms. In the present study, we analyzed the breakpoint junctions of INV-DUP-DEL patterns in five patients using long-read whole-genome sequencing and detected 2.2-6.1 kb copy-neutral regions in all five patients. At the end of the INV-DUP-DEL, two patients exhibited chromosomal translocations, which are recognized as telomere capture, and one patient showed direct telomere healing. The remaining two patients had additional small-sized intrachromosomal segments at the end of the derivative chromosomes. These findings have not been previously reported but they may only be explained by the presence of telomere capture breakage. Further investigations are required to better understand the mechanisms underlying this finding.

Application of Targeted Next-Generation Sequencing for the Investigation of Thalassemia in a Developing Country: A Single Center Experience

Thalassemia is identified as a prevalent disease in Malaysia, known to be one of the developing countries. Fourteen patients with confirmed cases of thalassemia were recruited from the Hematology Laboratory. The molecular genotypes of these patients were tested using the multiplex-ARMS and GAP-PCR methods. The samples were repeatedly investigated using the Devyser Thalassemia kit (Devyser, Sweden), a targeted NGS panel targeting the coding regions of hemoglobin genes, namely the HBA1, HBA2, and HBB genes, which were used in this study. There were many different genetic variants found in 14 unrelated cases. Out of all fourteen cases, NGS was able to determine an additional -50 G>A (HBB:c.-100G>A) that were not identified by the multiplex-ARMS method, including HBA2 mutations, namely CD 79 (HBA2:c.239C>G). Other than that, CD 142 (HBA2:c.427T>C) and another non-deletional alpha thalassemia and alpha triplication were also not picked up by the GAP-PCR methods. We illustrated a broad, targeted NGS-based test that proposes benefits rather than using traditional screening or basic molecular methods. The results of this study should be heeded, as this is the first report on the practicality of targeted NGS concerning the biological and phenotypic features of thalassemia, especially in a developing population. Discovering rare pathogenic thalassemia variants and additional secondary modifiers may facilitate precise diagnosis and better disease prevention.

Molecular Capture of Mycobacterium tuberculosis Genomes Directly from Clinical Samples: A Potential Backup Approach for Epidemiological and Drug Susceptibility Inferences

The application of whole genome sequencing of Mycobacterium tuberculosis directly on clinical samples has been investigated as a means to avoid the time-consuming need for culture isolation that can lead to a potential prolonged suboptimal antibiotic treatment. We aimed to provide a proof-of-concept regarding the application of the molecular capture of M. tuberculosis genomes directly from positive sputum samples as an approach for epidemiological and drug susceptibility predictions. Smear-positive sputum samples (n = 100) were subjected to the SureSelectXT HS Target Enrichment protocol (Agilent Technologies, Santa Clara, CA, USA) and whole-genome sequencing analysis. A higher number of reads on target were obtained for higher smear grades samples (i.e., 3+ followed by 2+). Moreover, 37 out of 100 samples showed ≥90% of the reference genome covered with at least 10-fold depth of coverage (27, 9, and 1 samples were 3+, 2+, and 1+, respectively). Regarding drug-resistance/susceptibility prediction, for 42 samples, ≥90% of the >9000 hits that are surveyed by TB-profiler were detected. Our results demonstrated that M. tuberculosis genome capture and sequencing directly from clinical samples constitute a potential valid backup approach for phylogenetic inferences and resistance prediction, essentially in settings when culture is not routinely performed or for samples that fail to grow.

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.

Applications of next generation sequencing in the screening and diagnosis of thalassemia: A mini-review

Thalassemias are a group of inherited blood disorders that affects 5-7% of the world population. Comprehensive screening strategies are essential for the management and prevention of this disorder. Today, many clinical and research laboratories have widely utilized next-generation sequencing (NGS) technologies to identify diseases, from germline and somatic disorders to infectious diseases. Yet, NGS application in thalassemia is limited and has just recently surfaced due to current demands in seeking alternative DNA screening tools that are more efficient, versatile, and cost-effective. This review aims to understand the several aspects of NGS technology, including its most current and expanding uses, advantages, and limitations, along with the issues and solutions related to its integration into routine screening and diagnosis of thalassemias. Hitherto, NGS has been a groundbreaking technology that offers tremendous improvements as a diagnostic tool for thalassemia in terms of its higher throughput, accuracy, and adaptability. The superiority of NGS in detecting rare variants, solving complex hematological problems, and providing non-invasive alternatives to neonatal diagnosis cannot be overlooked. However, several pitfalls still preclude its use as a stand-alone technique over conventional methods.

Improved Characterization of Complex β-Globin Gene Cluster Structural Variants Using Long-Read Sequencing

Complex insertion-deletion (indel) events in the globin genes manifest in widely variable clinical phenotypes. Many are incompletely characterized because of a historic lack of efficient methods. A more complete assessment enables improved prediction of clinical impact, which guides emerging therapeutic choices. Current methods have limited capacity for breakpoint assignment and accurate assessment of mutation extent, especially in cases containing duplications or multiple deletions and insertions. Technology, such as long-read sequencing, holds promise for significant impact in the characterization of indel events because of read lengths that span large regions, resulting in improved resolution. Four known complex β-globin gene cluster indel types were assessed using single-molecule, real-time sequencing technology and showed high correlation with previous reports, including the Caribbean locus control deletion (g.5,305,478_5,310,336del), a large β-gene duplication containing the Hb S mutation (g.4,640,335_5,290,171dup with g.5,248,232T>A, c.20A>T; variant allele fraction, 64%), and two nested variants (double deletions with intervening inversion): the Indian ^Gγ(^Aγδβ)⁰-thalassemia (g.5,246,804-5,254,275del, g.5,254,276_5,269,600inv, and g.5,269,601_5,270,442del) and the Turkish/Macedonian (δβ)⁰ thalassemia (g.5,235,064_5,236,652del, g.5,236,653_5,244,280inv, and g.5,244,281_5,255,766del). Our data confirm long-read sequencing as an efficient and accurate method to identify these clinically significant complex events. Limitations include high-complexity sample preparation requirements, which hinder routine use in clinical laboratories. Continued improvements in sample and data workflow processes are needed to accommodate volumes in a tertiary clinical laboratory.

Identification of thalassemia gene cluster deletion by long-read whole-genome sequencing (LR-WGS)

INTRODUCTION

At present, a variety of molecular detection methods are obtained to diagnose thalassemia accurately. Although exome sequencing or specific panels have been widely used in clinical diagnosis of genetic diseases, the positive rate is about 25%-30%. Because the detection range is limited to exons and splice sites, and the read length is usually 100-150 bp, there are limitations in the detection of globin gene clusters with pseudogenes.

METHODS

In this study, seven thalassemia patients were selected to perform whole-genome sequencing (WGS) with long read at 400 bp to make accurate detection for thalassemia deletions. And we used PCR and Sanger sequencing to confirm the gene deletions in the patients.

RESULTS

WGS analysis detected a rare 172 kb deletion on the α-globin gene cluster at chr16: 57 009-330 001, 19 kb deletion at chr16: 215 396-234 699, 11 kb deletion at chr16:220 861-231 981; and 27 kb deletion on the β-globin gene deletion at chr11: 5 222 878-5 250 288, 21.4 kb deletion at chr11: 5 236 361-5 257 771, 78.9 kb deletion at chr11: 5 191 121-5 270 050. All the seven patients carried heterozygous deletions, including three in α-gene cluster, three in β-gene cluster, and one in both globin clusters.

CONCLUSION

Our results indicate that long-read WGS will be beneficial to the diagnosis of genetic diseases with pseudogenes or highly duplicated sequences and will enable clinical geneticists to inform high-risk couples and provide prenatal diagnosis.

Breakpoint junction analysis for complex genomic rearrangements with the caldera volcano-like pattern

Chromosomal triplications can be classified into recurrent and nonrecurrent triplications. Most of the nonrecurrent triplications are embedded in duplicated segments, and duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) has been established as one of the mechanisms of triplication. This study aimed to reveal the underlying mechanism of the TRP-DUP-TRP pattern of chromosomal aberrations, in which the appearance of moving averages obtained through array-based comparative genomic hybridization analysis is similar to the shadows of the caldera volcano-like pattern, which were first identified in two patients with neurodevelopmental disabilities. For this purpose, whole-genome sequencing using long-read Nanopore sequencing was carried out to confirm breakpoint junctions. Custom array analysis and Sanger sequencing were also used to detect all breakpoint junctions. As a result, the TRP-DUP-TRP pattern consisted of only two patterns of breakpoint junctions in both patients. In patient 1, microhomologies were identified in breakpoint junctions. In patient 2, more complex architectures with insertional segments were identified. Thus, replication-based mechanisms were considered as a mechanism of the TRP-DUP-TRP pattern.

Complex chromosomal rearrangements of human chromosome 21 in a patient manifesting clinical features partially overlapped with that of Down syndrome

The chromosomal region critical in Down syndrome has long been analyzed through genotype-phenotype correlation studies using data from many patients with partial trisomy 21. Owing to that, a relatively small region of human chromosome 21 (35.9 ~ 38.0 Mb) has been considered as Down syndrome critical region (DSCR). In this study, microarray-based comparative genomic hybridization analysis identified complex rearrangements of chromosome 21 in a patient manifesting clinical features partially overlapped with that of Down syndrome. Although the patient did not show up-slanting palpebral fissures and single transverse palmar creases, other symptoms were consistent with Down syndrome. Rearrangements were analyzed by whole-genome sequencing using Nanopore long-read sequencing. The analysis revealed that chromosome 21 was fragmented into seven segments and reassembled by six connected points. Among 12 breakpoints, 5 are located within the short region and overlapped with repeated segments. The rearrangement resulted in a maximum gain of five copies, but no region showed loss of genomic copy numbers. Breakpoint-junctions showed no homologous region. Based on these findings, chromoanasynthesis was considered as the mechanism. Although the distal 21q22.13 region was not included in the aberrant regions, some of the genes located on the duplicated regions, SOD1, SON, ITSN1, RCAN1, and RUNX1, were considered as possible candidate genes for clinical features of the patient. We discussed the critical region for Down syndrome, with the literature review.

Diagnosis of the accurate genotype of HKαα carriers in patients with thalassemia using multiplex ligation-dependent probe amplification combined with nested polymerase chain reaction

BACKGROUND

Patients carrying the HongKongαα (HKαα) allele and -α/ααα could be misdiagnosed as -α/αα by the current conventional thalassemia detection methods, leading to inaccurate genetic counseling and an incorrect prenatal diagnosis. This study was aimed to accurately analyze the genotypes of HKαα carriers and -α/ααα.

METHODS

Samples were collected in our hospital from July 2017 to October 2019. Twenty-four common types of Chinese thalassemia were screened by gap-polymerase chain reaction (Gap-PCR) and reverse dot blot (RDB). Anti-4.2 multiplex-PCR was used to confirm carriers of the ααα duplication with -α deletion. Two-round nested PCR and multiplex ligation-dependent probe amplification (MLPA) were applied to accurately identify and confirm their genotypes. For data analysis, we used descriptive statistics and Fisher's exact tests.

RESULTS

Two thousand five hundred and forty-four cases were identified as thalassemia in 5488 peripheral blood samples. The results showed that α, β, and αβ compound thalassemia were identified in 1190 (46.78%), 1286 (50.55%), and 68 (2.67%) cases, respectively. A total of 227 samples from thalassemia patients were identified as -α/αα by Gap-PCR, and the genotypes of two samples were uncertain. There was a difference between Gap-PCR and combined groups (Gap-PCR combined with nested PCR and MLPA) in detecting HKαα (P < 0.05). Among the 229 patients, 20 patients were identified as HKαα carriers and one was identified as -α/ααα by two-round nested PCR and MLPA, including 15 patients with HKαα/αα, three with HKαα/αα and β-thalassemia coinheritance, one with HKαα/--, one with HKαα/-α and β-thalassemia coinheritance, and one with -α/ααα and β-thalassemia coinheritance.

CONCLUSIONS

ααα and HKαα genotypes of patients carrying -α need to be detected to reduce the misdiagnosis rate of patients carrying HKαα and -α3.7/ααα alleles. More accurate genetic counseling can be provided in the clinic using nested PCR combined with MLPA.

Visual genotyping of thalassemia by using pyrrolidinyl peptide nucleic acid probes immobilized on carboxymethylcellulose-modified paper and enzyme-induced pigmentation

A simple probe pair was designed for the detection of hemoglobin E (HbE) genotype, a single-point mutation that leads to abnormal red blood cells commonly found in South East Asia. The key to differentiation is the use of a conformationally constrained peptide nucleic acid (PNA) that was immobilized on carboxymethylcellulose-modified paper. This was then used for target DNA binding and visualization by an enzyme-catalyzed pigmentation. The biotinylated target DNA bound to the immobilized probe was visually detected via alkaline phosphatase-linked streptavidin. This enzyme conjugate catalyzed the dephosphorylation of the substrate 5-bromo-4-chloro-3-indolyl phosphate, leading to a series of reactions that generate an intense, dark blue pigment. The test was validated with 100 DNA samples, which shows good discrimination among different genotypes (normal, HbE, and heterozygous) with 100% accuracy when optimal conditions of analysis were applied. The method does not require temperature control and can be performed at ambient temperature. This is an attractive feature for diagnosis in primary care, which accounts for a large part of affected population. Graphical abstract Schematic representation of a paper-based sensor for the detection of the gene Hemoglobin E. The interaction between an immobilized peptide nucleic acid and a DNA target leads to enzymatic pigmentation, allowing simple visual readout with up to 100% accuracy.

Diagnosis of Haemoglobinopathies: New Scientific Advances

The molecular defects underlying haemoglobinopathies are both deletions and point mutations in the alpha- or beta-globin genes or gene-clusters. To detect point mutations causing alpha- or beta-thalassaemia, direct sequencing is the method of choice to detect the widest spectrum of molecular defects. The most established approach in DNA diagnostics to screen for the most common deletion defects causing alpha-thalassaemia or beta-thalassaemia is gap- PCR, Multiplex Ligation-dependent Probe Amplification (MLPA) and Sanger Sequencing technology to detect breakpoint sequences of previously uncharacterized deletions/duplications. We demonstrate the recent advances in the determination of duplications and deletions causing alpha- or beta-thalassemia, using Next Generation Sequencing, array Comparative Genome Hybridization and Target Locus Amplification. We present three cases in which the use of advanced technologies allow the diagnosis of unexpected disease genotypes.

Molecular basis of α-thalassemia

α-Thalassemia is an inherited, autosomal recessive, disorder characterized by a microcytic hypochromic anemia. It is one of the most common monogenic gene disorders in the world population. The clinical severity varies from almost asymptomatic, to mild microcytic hypochromic, and to a lethal hemolytic condition, called Hb Bart's Hydrops Foetalis Syndrome. The molecular basis are usually deletions and less frequently, point mutations affecting the expression of one or more of the duplicated α-genes. The clinical variation and increase in disease severity is directly related to the decreased expression of one, two, three or four copies of the α-globin genes. Deletions and point mutations in the α-globin genes and their regulatory elements have been studied extensively in carriers and patients and these studies have given insight into the α-globin genes are regulated. By looking at naturally occurring deletions and point mutations, our knowledge of globin-gene regulation and expression will continue to increase and will lead to new targets of therapy.