High accuracy of single-molecule real-time sequencing in detecting a rare α-globin fusion gene in carrier screening population.

INTRODUCTION: The α-globin fusion gene between the HBA2 and HBAP1 genes becomes clinically important in thalassemia screening because this fusion gene can cause severe hemoglobin (Hb) H disease when combining with α0 -thalassemia (α0 -thal). Due to its uncommon rearrangement in the α gene cluster without dosage changes, this fusion gene is undetectable by common molecular testing approaches used for α-thal diagnosis. METHODS: In this study, we used the single-molecule real-time (SMRT) sequencing technique to detect this fusion gene in 23 carriers identified by next-generation sequencing (NGS) among 16,504 screened individuals. Five primers for α and ß thalassemia were utilized. RESULTS: According to the NGS results, the 23 carriers include 14 pure heterozygotes, eight compound heterozygotes with common α-thal alleles, and one homozygote. By using SMRT, the fusion mutant was successfully detected in all 23 carriers. Furthermore, SMRT corrected the diagnosis in two "pure" heterozygotes: one was compound heterozygote with anti-3.7 triplication, and the other was homozygote. CONCLUSION: Our results indicate that SMRT is a superior method compared to NGS in detecting the α fusion gene, attributing to its efficient, accurate, and one-step properties.

Long-read HiFi sequencing of NUDT15: Phased full-gene haplotyping and pharmacogenomic allele discovery.

To determine the phase of NUDT15 sequence variants for more comprehensive star (*) allele diplotyping, we developed a novel long-read single-molecule real-time HiFi amplicon sequencing method. A 10.5 kb NUDT15 amplicon assay was validated using reference material positive controls and additional samples for specimen type and blinded accuracy assessment. Triplicate NUDT15 HiFi sequencing of two reference material samples had nonreference genotype concordances of >99.9%, indicating that the assay is robust. Notably, short-read genome sequencing of a subset of samples was unable to determine the phase of star (*) allele-defining NUDT15 variants, resulting in ambiguous diplotype results. In contrast, long-read HiFi sequencing phased all variants across the NUDT15 amplicons, including a *2/*9 diplotype that previously was characterized as *1/*2 in the 1000 Genomes Project v3 data set. Assay throughput was also tested using 8.5 kb amplicons from 100 Ashkenazi Jewish individuals, which identified a novel NUDT15 *1 suballele (c.-121G>A) and a rare likely deleterious coding variant (p.Pro129Arg). Both novel alleles were Sanger confirmed and assigned as *1.007 and *20, respectively, by the PharmVar Consortium. Taken together, NUDT15 HiFi amplicon sequencing is an innovative method for phased full-gene characterization and novel allele discovery, which could improve NUDT15 pharmacogenomic testing and subsequent phenotype prediction.

Context-dependent DNA polymerization effects can masquerade as DNA modification signals.

BACKGROUND: Single molecule measurements of DNA polymerization kinetics provide a sensitive means to detect both secondary structures in DNA and deviations from primary chemical structure as a result of modified bases. In one approach to such analysis, deviations can be inferred by monitoring the behavior of DNA polymerase using single-molecule, real-time sequencing with zero-mode waveguide. This approach uses a Single Molecule Real Time (SMRT)-sequencing measurement of time between fluorescence pulse signals from consecutive nucleosides incorporated during DNA replication, called the interpulse duration (IPD). RESULTS: In this paper we present an analysis of loci with high IPDs in two genomes, a bacterial genome (E. coli) and a eukaryotic genome (C. elegans). To distinguish the potential effects of DNA modification on DNA polymerization speed, we paired an analysis of native genomic DNA with whole-genome amplified (WGA) material in which DNA modifications were effectively removed. Adenine modification sites for E. coli are known and we observed the expected IPD shifts at these sites in the native but not WGA samples. For C. elegans, such differences were not observed. Instead, we found a number of novel sequence contexts where IPDs were raised relative to the average IPDs for each of the four nucleotides, but for which the raised IPD was present in both native and WGA samples. CONCLUSION: The latter results argue strongly against DNA modification as the underlying driver for high IPD segments for C. elegans, and provide a framework for separating effects of DNA modification from context-dependent DNA polymerase kinetic patterns inherent in underlying DNA sequence for a complex eukaryotic genome.

Genome Resource of a Hypervirulent Strain C9-3 of Xanthomonas oryzae pv. oryzae Causing Bacterial Blight of Rice.

Xanthomonas oryzae pv. oryzae is the causal agent of bacterial blight, one of the most devastating diseases of rice. Here, a hypervirulent strain, C9-3, defeating Xa1, Xa10, xa13, and Xa23 resistance genes, was used to extract genomic DNA for single molecule real-time (SMRT) sequencing. After assembly, the genome consists of a single-circular chromosome with the size of 4,924,298 bp with G+C content of 63.7% and contains 4,715 genes. Annotation and analysis of the TALE genes using a suite of applications named AnnoTALE suggested that 17 transcription activator-like effectors, including 15 typical TALEs and 2 iTALEs/truncTALEs, were encoded in the genome. The approach and genome resource will contribute to the discovery of new virulence effectors and understanding on rice-X. oryzae pv. oryzae interactions.

Pathogenic 12-kb copy-neutral inversion in syndromic intellectual disability identified by high-fidelity long-read sequencing.

We report monozygotic twin girls with syndromic intellectual disability who underwent exome sequencing but with negative pathogenic variants. To search for variants that are unrecognized by exome sequencing, high-fidelity long-read genome sequencing (HiFi LR-GS) was applied. A 12-kb copy-neutral inversion was precisely identified by HiFi LR-GS after trio-based variant filtering. This inversion directly disrupted two genes, CPNE9 and BRPF1, the latter of which attracted our attention because pathogenic BRPF1 variants have been identified in autosomal dominant intellectual developmental disorder with dysmorphic facies and ptosis (IDDDFP), which later turned out to be clinically found in the twins. Trio-based HiFi LR-GS together with haplotype phasing revealed that the 12-kb inversion occurred de novo on the maternally transmitted chromosome. This study clearly indicates that submicroscopic copy-neutral inversions are important but often uncharacterized culprits in monogenic disorders and that long-read sequencing is highly advantageous for detecting such inversions involved in genetic diseases.

A Novel 5 kb Deletion in the ß-Globin Gene Cluster Identified in a Chinese Patient.

ß-Thalassemia (ß-thal), a highly prevalent disease in tropical and subtropical regions of Southern China, is caused mainly by point mutations in the ß-globin gene cluster. However, large deletions have also been found to contribute to some types of ß-thal. We identified a novel 5 kb deletion in the ß-globin cluster in a Chinese patient using multiplex ligation-dependent probe amplification (MLPA), and characterized it with single molecule real-time (SMRT) sequencing, gap-polymerase chain reaction (gap-PCR) and Sanger sequencing. The deletion was located between positions 5226189 and 5231091 on chromosome 11 (GRCh38), extending from 4 kb upstream of the 5' untranslated region (5'UTR) to the second intron of the ß-globin gene. The patient with this deletion presented with microcytosis and hypochromic red cells, as well as relatively high Hb F and Hb A2 levels. Our research indicated that SMRT sequencing is a useful tool for accurate detection of large deletions. Our study broadens the spectrum of deletional ß-thalassemias and provides a perspective for further study of the function of the ß-globin cluster.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds.

BACKGROUND: The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map. RESULTS: Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features. CONCLUSIONS: The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Investigating the bacterial microbiota of traditional fermented dairy products using propidium monoazide with single-molecule real-time sequencing.

Traditional fermented dairy foods have been the major components of the Mongolian diet for millennia. In this study, we used propidium monoazide (PMA; binds to DNA of nonviable cells so that only viable cells are enumerated) and single-molecule real-time sequencing (SMRT) technology to investigate the total and viable bacterial compositions of 19 traditional fermented dairy foods, including koumiss from Inner Mongolia (KIM), koumiss from Mongolia (KM), and fermented cow milk from Mongolia (CM); sample groups treated with PMA were designated PKIM, PKM, and PCM. Full-length 16S rRNA sequencing identified 195 bacterial species in 121 genera and 13 phyla in PMA-treated and untreated samples. The PMA-treated and untreated samples differed significantly in their bacterial community composition and α-diversity values. The predominant species in KM, KIM, and CM were Lactobacillus helveticus, Streptococcus parauberis, and Lactobacillus delbrueckii, whereas the predominant species in PKM, PKIM, and PCM were Enterobacter xiangfangensis, Lactobacillus helveticus, and E. xiangfangensis, respectively. Weighted and unweighted principal coordinate analyses showed a clear clustering pattern with good separation and only minor overlapping. In addition, a pure culture method was performed to obtain lactic acid bacteria resources in dairy samples according to the results of SMRT sequencing. A total of 102 LAB strains were identified and Lb. helveticus (68.63%) was the most abundant, in agreement with SMRT sequencing results. Our results revealed that the bacterial communities of traditional dairy foods are complex and vary by type of fermented dairy product. The PMA treatment induced significant changes in bacterial community structure.

Adaptive Strategies of Bacillus thuringiensis Isolated from Acid Mine Drainage Site in Sabah, Malaysia.

The adaptive process in bacteria is driven by specific genetic elements which regulate phenotypic characteristics such as tolerance to high metal ion concentrations and the secretion of protective biofilms. Extreme environments such as those associated with heavy metal pollution and extremes of acidity offer opportunities to study the adaptive mechanisms of microorganisms. This study focused on the genome analysis of Bacillus thuringiensis (Bt MCMY1), a gram positive rod shaped bacterium isolated from an acid mine drainage site in Sabah, Malaysia by using a combination of Single Molecule Real Time DNA Sequencing, Scanning Electron Microscopy (SEM) and Fourier Transform Infrared Spectroscopy (FTIR). The genome size of Bt MCMY1 was determined to be 5,458,152 bases which was encoded on a single chromosome. Analysis of the genome revealed genes associated with resistance to Copper, Mercury, Arsenic, Cobalt, Zinc, Cadmium and Aluminum. Evidence from SEM and FTIR indicated that the bacterial colonies form distinct films which bear the signature of polyhydroxyalkanoates (PHA) and this finding was supported by the genome data indicating the presence of a genetic pathway associated with the biosynthesis of PHAs. This is the first report of a Bacillus sp. isolated from an acid mine drainage site in Sabah, Malaysia and the genome sequence will provide insights into the manner in which B. thuringiensis adapts to acid mine drainage.

Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6.

The cytochrome P450-2D6 (CYP2D6) enzyme metabolizes ∼25% of common medications, yet homologous pseudogenes and copy number variants (CNVs) make interrogating the polymorphic CYP2D6 gene with short-read sequencing challenging. Therefore, we developed a novel long-read, full gene CYP2D6 single molecule real-time (SMRT) sequencing method using the Pacific Biosciences platform. Long-range PCR and CYP2D6 SMRT sequencing of 10 previously genotyped controls identified expected star (*) alleles, but also enabled suballele resolution, diplotype refinement, and discovery of novel alleles. Coupled with an optimized variant-calling pipeline, CYP2D6 SMRT sequencing was highly reproducible as triplicate intra- and inter-run nonreference genotype results were completely concordant. Importantly, targeted SMRT sequencing of upstream and downstream CYP2D6 gene copies characterized the duplicated allele in 15 control samples with CYP2D6 CNVs. The utility of CYP2D6 SMRT sequencing was further underscored by identifying the diplotypes of 14 samples with discordant or unclear CYP2D6 configurations from previous targeted genotyping, which again included suballele resolution, duplicated allele characterization, and discovery of a novel allele and tandem arrangement. Taken together, long-read CYP2D6 SMRT sequencing is an innovative, reproducible, and validated method for full-gene characterization, duplication allele-specific analysis, and novel allele discovery, which will likely improve CYP2D6 metabolizer phenotype prediction for both research and clinical testing applications.

High-resolution community profiling of arbuscular mycorrhizal fungi.

Community analyses of arbuscular mycorrhizal fungi (AMF) using ribosomal small subunit (SSU) or internal transcribed spacer (ITS) DNA sequences often suffer from low resolution or coverage. We developed a novel sequencing based approach for a highly resolving and specific profiling of AMF communities. We took advantage of previously established AMF-specific PCR primers that amplify a c. 1.5-kb long fragment covering parts of SSU, ITS and parts of the large ribosomal subunit (LSU), and we sequenced the resulting amplicons with single molecule real-time (SMRT) sequencing. The method was applicable to soil and root samples, detected all major AMF families and successfully discriminated closely related AMF species, which would not be discernible using SSU sequences. In inoculation tests we could trace the introduced AMF inoculum at the molecular level. One of the introduced strains almost replaced the local strain(s), revealing that AMF inoculation can have a profound impact on the native community. The methodology presented offers researchers a powerful new tool for AMF community analysis because it unifies improved specificity and enhanced resolution, whereas the drawback of medium sequencing throughput appears of lesser importance for low-diversity groups such as AMF.

Unveiling genomic features linked to traits of plant growth-promoting bacterial communities from sugarcane.

Microorganisms are ubiquitous, and those inhabiting plants have been the subject of several studies. Plant-associated bacteria exhibit various biological mechanisms that enable them to colonize host plants and, in some cases, enhance their fitness. In this study, we describe the genomic features predicted to be associated with plant growth-promoting traits in six bacterial communities isolated from sugarcane. The use of highly accurate single-molecule real-time sequencing technology for metagenomic samples from these bacterial communities allowed us to recover 17 genomes. The taxonomic assignments for the binned genomes were performed, revealing taxa distributed across three main phyla: Bacillota, Bacteroidota, and Pseudomonadota, with the latter being the most representative. Subsequently, we functionally annotated the metagenome-assembled genomes (MAGs) to characterize their metabolic pathways related to plant growth-promoting traits. Our study successfully identified the enrichment of important functions related to phosphate and potassium acquisition, modulation of phytohormones, and mechanisms for coping with abiotic stress. These findings could be linked to the robust colonization of these sugarcane endophytes.

Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations-Application to HIV-1 quasispecies.

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing, which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence polymerase-chain reaction (PCR) amplicons derived from cDNA templates tagged with unique molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR. The use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Production of highly accurate sequences from the large datasets produced from SMRT-UMI sequencing is facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline). PORPIDpipeline automatically filters and parses circular consensus reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination, heteroduplex formation, or early cycle PCR errors. The optimized SMRT-UMI sequencing and PORPIDpipeline methods presented here represent a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus quasispecies in a virus transmitter-recipient pair of individuals.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features.

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Distribution characteristics and clinical phenotype analyses of hemoglobin variants in the Z region of Central Guangxi, Southern China.

OBJECTIVE: To investigate the molecular diagnosis of hemoglobin variants in Z region by Capillary electrophoresis in Central Guangxi, Southern China, and analyze their distribution and phenotypic characteristics, to provide a reference for clinical consultation and prenatal diagnosis for couples. METHODS: A total of 23,709 subjects were collected for blood routine analysis, hemoglobin analysis, and common α- and ß-globin gene loci in Chinese population. The hemoglobin electrophoresis components were divided into Zone 1-Zone 15 (Z1-Z15) by Capillary zone electrophoresis (CE). For samples not clearly detected by the conventional technology, Sanger sequencing, and multiplex ligation-dependent probe amplification (MLPA) were used. Single-molecule real-time (SMRT) sequencing technology was used to analyze rare-type genes in a sample with a structural variation. RESULTS: Ten rare hemoglobin variants distributed in Z region were detected in 23,709 samples, including Hb Cibeles, which was reported for the first time in Asia; Hb J-Broussais, Hb G-Honolulu and J-Wenchang-Wuming, they were first reported in Guangxi; 1 case of Hb Anti-Lepore Liuzhou, which was a newly discovered hemoglobin variant; hemoglobin variants Hb G-Siriraj, Hb Handsworth, Hb Q-Thailand, Hb Ube-2, Hb NewYork were also detected. CONCLUSION: There are a few studies on rare hemoglobin variants in Z region in Southern China. Ten rare hemoglobin variants were found in this study. The hematological phenotype and component content of hemoglobin variants are related to the occurrence of thalassemia. This study enriched the data of rare hemoglobin variants in Southern China and provided a comprehensive data basis for prenatal diagnosis of hemoglobin variants in this area.

The Long and Short of Genome Sequencing: Using a Hybrid Sequencing Strategy to Sequence Oral Microbial Genomes.

Since our chapter on genome sequencing using the GS-FLX pyrosequencer in the First Edition of this book, significant advances have been made in next-generation DNA sequencing (NGS) technology. Not only has the GS-FLX become extinct, but the more recent introduction and establishment of the so-called third-generation DNA sequencers by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) has revolutionized genomics yet again by generating ultra-long (>100,000 basepair) sequence reads concomitant with an incredible reduction in cost per sequenced basepair. Unfortunately, the ultra-high sequence yields of third-generation sequencers are compromised by their inherent sequencing error rates, prompting an alternative sequencing strategy, i.e., a hybrid sequencing strategy, which combines PacBio/ONT primary datasets with complementary datasets generated by mainstream short-read NGS platforms, e.g., Illumina or Ion Torrent. Although the concept of a hybrid sequencing strategy is not new, existing yields and accuracy of ultra-long and short-read sequencing technologies makes such a strategy achievable, resulting in complete genome sequences in one hit. In this chapter, we describe our updated laboratory and bioinformatic protocols that will allow the average research group to obtain complete oral microbial genome sequences assembled from a combination of DNA sequence data generated by NGS and third-generation platforms.

Hb Q-Thailand heterozygosity unlinked with the (-α4.2/) α+-thalassemia deletion allele identified by long-read SMRT sequencing: hematological and molecular analyses.

OBJECTIVE: In the present study, two unrelated cases of Hb Q-Thailand heterozygosity unlinked with the (-α4.2/) α+-thalassemia deletion allele were identified by long-read single molecule real-time (SMRT) sequencing in southern China. The aim of this study was to report the hematological and molecular features as well as diagnostic aspects of the rare manifestation. METHODS: Hematological parameters and hemoglobin analysis results were recorded. A suspension array system for routine thalassemia genetic analysis and long-read SMRT sequencing were applied in parallel for thalassemia genotyping. Traditional methods, including Sanger sequencing, multiplex gap-polymerase chain reaction (gap-PCR) and multiplex ligation-dependent probe amplification (MLPA), were used together to confirm the thalassemia variants. RESULTS: Long-read SMRT sequencing was used to diagnose two Hb Q-Thailand heterozygous patients for whom the hemoglobin variant was unlinked to the (-α4.2/) allele for the first time. The hitherto undescribed genotypes were verified by traditional methods. Hematological parameters were compared with those of Hb Q-Thailand heterozygosity linked with the (-α4.2/) deletion allele in our study. For the positive control samples, long-read SMRT sequencing revealed a linkage relationship between the Hb Q-Thailand allele and the (-α4.2/) deletion allele. CONCLUSIONS: Identification of the two patients confirms that the linkage relationship between the Hb Q-Thailand allele and the (-α4.2/) deletion allele is a common possibility but not a certainty. Remarkably, as it is superior to traditional methods, SMRT technology may eventually serve as a more comprehensive and precise method that holds promising prospects in clinical practice, especially for rare variants.

A novel 15.8 kb deletion α-thalassemia confirmed by long-read single-molecule real-time sequencing: Hematological phenotypes and molecular characterization.

BACKGROUND: Thalassemia is the most frequent recessive Mendelian inherited monogenic disease worldwide, and is characterized by the impaired synthesis of globin chains due to disease-causing variants in α- or ß-globin genes. There are many conventional methods to diagnose thalassemia but all of them have limitations. CASE REPORT: We present the case of a 37-year-old female with abnormal values of routine hematological indices who was admitted for genetic screening of thalassemia. Genomic DNA was extracted and used for genetic assays covering the known and potential novel genotypes in HBA and HBB genes using a suspension-array system, gap-polymerase chain reaction (Gap-PCR), PCR-reverse dot blot (PCR-RDB) and multiplex ligation-dependent probe amplification (MLPA). Finally, using long-read single-molecule real-time (SMRT) sequencing, we first confirmed the case with a novel 15.8 kb deletion located in the HBA gene (Chr16:163886-179768, GRch38/hg38). CONCLUSIONS: Our results showed that long-read SMRT sequencing has great advantages in the detection of rare α-globin gene variants. This study may provide a reference protocol for the use of long-read SMRT sequencing for the detection of known and potential novel genotypes of thalassemia in the population and improve the accuracy of genetic counseling and prenatal diagnosis.

Detection of rare thalassemia mutations using long-read single-molecule real-time sequencing.

Gap- polymerase chain reaction (PCR), reverse dot-blot assay (RDB), real-time PCR based multicolor melting curve analysis (MMCA assay), multiplex ligation-dependent probe amplification (MLPA) and Sanger sequencing are conventional methods to diagnose thalassemia but all of them have limitations. In this study, we applied single-molecule real-time (SMRT) sequencing following multiplex long-range PCR to uncover rare mutations in nine patients and their family members. The patients with different results between Gap-PCR and MMCA assay or with phenotype not matching genotype were included. Using SMRT sequencing, we first identified the carriers with αααanti3.7/HKαα, -α762bpα/αα (chr16:172,648-173,409), ααfusion/αQSα (in a trans configuration), two cases with novel gene rearrangements and another case with a novel 341 bp insertion in α-globin gene cluster, respectively. One carrier with --SEA/αααanti4.2, and two carriers with the coexistence of globin variant and an α-globin gene duplication were also found. Most importantly, we could determine two defects in α-globin gene cluster being a cis or trans configuration in a single test. Our results showed that SMRT has great advantages in detection of α-globin gene triplications, rare deletions and determination of a cis or trans configuration. SMRT is a comprehensive and one-step method for thalassemia screening and diagnosis, especially for detection of rare thalassemia mutations.

Dynamic Interplay between Structural Variations and 3D Genome Organization in Pancreatic Cancer.

Structural variations (SVs) are the greatest source of variations in the genome and can lead to oncogenesis. However, the identification and interpretation of SVs in human cancer remain technologically challenging. Here, long-read sequencing is first employed to depict the signatures of structural variations in carcinogenesis of human pancreatic ductal epithelium. Then widespread reprogramming of the 3D chromatin architecture is revealed by an in situ Hi-C technique. Integrative analyses indicate that the distribution pattern of SVs among the 3D genome is highly cell-type specific and the bulk remodeling effects of SVs in the chromatin organization partly depend on intercellular genomic heterogeneity. Meanwhile, contact domains tend to minimize these disrupting effects of SVs within local adjacent genomic regions to maintain overall stability. Notably, complex genomic rearrangements involving two key driver genes CDKN2A and SMAD4 are identified, and their influence on the expression of oncogenes MIR31HG, MYO5B, etc., are further elucidated from both a linear view and 3D perspective. Overall, this work provides a genome-wide resource and highlights the impact, complexity, and dynamicity of the interplay between structural variations and high-order chromatin organization, which expands the current understanding of the pathogenesis of SVs in human cancer.