Advantages of genome sequencing by long-read sequencer using SMRT technology in medical area

Detection of hemophilia a genetic variants using third-generation long-read sequencing

BACKGROUND

Hemophilia A (HA) is an X-linked recessive genetic disorder caused by pathogenic variations of the factor VIII -encoding gene, F8 gene. Due to the large size and diverse types of variations in the F8 gene, causative mutations in F8 cannot be simultaneously detected in one step by traditional molecular analysis, and genetic molecular diagnosis and prenatal screening of HA still face significant difficulties and challenges in clinical practice. Therefore, we aimed to develop and validate an efficient, accurate, and time-saving method for the genetic detection of HA.

METHODS

A comprehensive analysis of hemophilia A (CAHEA) method based on long-range PCR and long-read sequencing (LRS) was used to detect F8 gene mutations in 14 clinical HA samples. The LRS results were compared with those of the conventional methods to evaluate the accuracy and sensitivity of the proposed approach.

RESULTS

The CAHEA method successfully identified 14 F8 variants in all probands, including 3 small insertion deletions, 4 single nucleotide variants, and 7 intron 22 inversions in a "one-step" manner, of which 2 small deletions have not been reported previously. Moreover, this method provided an opportunity to analyze the mechanism of rearrangement and the pathogenicity of F8 variants. The LRS results were validated and found to be in 100% agreement with those obtained using the conventional method.

CONCLUSION

Our proposed LRS-based F8 gene detection method is an accurate and reproducible genetic screening and diagnostic method with significant clinical value. It provides efficient, comprehensive, and accurate genetic screening and diagnostic services for individuals at high risk of HA as well as for premarital and prenatal populations.

The Third-Generation Sequencing Challenge: Novel Insights for the Omic Sciences

The understanding of the human genome has been greatly improved by the advent of next-generation sequencing technologies (NGS). Despite the undeniable advantages responsible for their widespread diffusion, these methods have some constraints, mainly related to short read length and the need for PCR amplification. As a consequence, long-read sequencers, called third-generation sequencing (TGS), have been developed, promising to overcome NGS. Starting from the first prototype, TGS has progressively ameliorated its chemistries by improving both read length and base-calling accuracy, as well as simultaneously reducing the costs/base. Based on these premises, TGS is showing its potential in many fields, including the analysis of difficult-to-sequence genomic regions, structural variations detection, RNA expression profiling, DNA methylation study, and metagenomic analyses. Protocol standardization and the development of easy-to-use pipelines for data analysis will enhance TGS use, also opening the way for their routine applications in diagnostic contexts.

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.

Whole Genome Sequencing: Applications in Clinical Bacteriology

The success in determining the whole genome sequence of a bacterial pathogen was first achieved in 1995 by determining the complete nucleotide sequence of Haemophilus influenzae Rd using the chain-termination method established by Sanger et al. in 1977 and automated by Hood et al. in 1987. However, this technology was laborious, costly, and time-consuming. Since 2004, high-throughput next-generation sequencing technologies have been developed, which are highly efficient, require less time, and are cost-effective for whole genome sequencing (WGS) of all organisms, including bacterial pathogens. In recent years, the data obtained using WGS technologies coupled with bioinformatics analyses of the sequenced genomes have been projected to revolutionize clinical bacteriology. WGS technologies have been used in the identification of bacterial species, strains, and genotypes from cultured organisms and directly from clinical specimens. WGS has also helped in determining resistance to antibiotics by the detection of antimicrobial resistance genes and point mutations. Furthermore, WGS data have helped in the epidemiological tracking and surveillance of pathogenic bacteria in healthcare settings as well as in communities. This review focuses on the applications of WGS in clinical bacteriology.

When Livestock Genomes Meet Third-Generation Sequencing Technology: From Opportunities to Applications

Third-generation sequencing technology has found widespread application in the genomic, transcriptomic, and epigenetic research of both human and livestock genetics. This technology offers significant advantages in the sequencing of complex genomic regions, the identification of intricate structural variations, and the production of high-quality genomes. Its attributes, including long sequencing reads, obviation of PCR amplification, and direct determination of DNA/RNA, contribute to its efficacy. This review presents a comprehensive overview of third-generation sequencing technologies, exemplified by single-molecule real-time sequencing (SMRT) and Oxford Nanopore Technology (ONT). Emphasizing the research advancements in livestock genomics, the review delves into genome assembly, structural variation detection, transcriptome sequencing, and epigenetic investigations enabled by third-generation sequencing. A comprehensive analysis is conducted on the application and potential challenges of third-generation sequencing technology for genome detection in livestock. Beyond providing valuable insights into genome structure analysis and the identification of rare genes in livestock, the review ventures into an exploration of the genetic mechanisms underpinning exemplary traits. This review not only contributes to our understanding of the genomic landscape in livestock but also provides fresh perspectives for the advancement of research in this domain.

Unraveling metagenomics through long-read sequencing: a comprehensive review

The study of microbial communities has undergone significant advancements, starting from the initial use of 16S rRNA sequencing to the adoption of shotgun metagenomics. However, a new era has emerged with the advent of long-read sequencing (LRS), which offers substantial improvements over its predecessor, short-read sequencing (SRS). LRS produces reads that are several kilobases long, enabling researchers to obtain more complete and contiguous genomic information, characterize structural variations, and study epigenetic modifications. The current leaders in LRS technologies are Pacific Biotechnologies (PacBio) and Oxford Nanopore Technologies (ONT), each offering a distinct set of advantages. This review covers the workflow of long-read metagenomics sequencing, including sample preparation (sample collection, sample extraction, and library preparation), sequencing, processing (quality control, assembly, and binning), and analysis (taxonomic annotation and functional annotation). Each section provides a concise outline of the key concept of the methodology, presenting the original concept as well as how it is challenged or modified in the context of LRS. Additionally, the section introduces a range of tools that are compatible with LRS and can be utilized to execute the LRS process. This review aims to present the workflow of metagenomics, highlight the transformative impact of LRS, and provide researchers with a selection of tools suitable for this task.

Pathogenic bacteria and fungi in bioaerosols from specialized hospitals in Shandong province, East China

Bacteria and fungi are abundant and ubiquitous in bioaerosols in hospital environments. Understanding the distribution and diversity of microbial communities within bioaerosols is critical for mitigating their detrimental effects. Our knowledge on the composition of bacteria or fungi in bioaerosols is limited, especially the potential pathogens present in fine particulate matter (PM2.5) from specialized hospitals. Thirty p.m.2.5 filter samples were collected from five hospitals (i.e., oral, dermatology, chest, eye, and general hospitals) in Shandong Province, East China. The diversity of bacteria and fungi was analyzed at the species level using single-molecule real-time sequencing of the 16 S and internal transcribed spacer 1 (ITS) ribosomal genes, respectively. Significant differences were detected across sampling sites in terms of microbial diversity and community composition in PM2.5 as well as pollution concentrations. The range of PM2.5 concentrations observed in hospital halls was higher, ranging from 39.0 to 46.2 μg/m3, compared to the wards where the concentrations ranged from 10.7 to 25.2 μg/m3. Furthermore, microbial variations in PM2.5 bioaerosols were associated with hospital type. The most dominant pathogens identified were Vibrio metschnikovii, Staphylococcus epidermidis, Staphylococcus haemolyticus, Fusarium pseudensiforme, and Aspergillus ruber. Among these, A. ruber was identified as an opportunistic fungus in a hospital setting for the first time. Nine potentially novel strains of F. pseudensiforme, showing 84.5%-92.0% ITS sequence similarity to known Fusarium isolates, were identified in PM2.5 samples from all hospitals (excluding an eye hospital). This study highlights the importance of hospital environments in shaping microbial aerosol communities. To the best of our knowledge, this is the first study to provide insights into the bacterial and fungal biodiversity of PM2.5 in specialized hospitals, enriching research in healthcare environmental microbiology and carrying significant public health implications.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPSeq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPSeq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate Long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPseq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection (UMAP) analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen (HLA) molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPseq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

Comparative genomics and DNA methylation analysis of Pseudomonas aeruginosa clinical isolate PA3 by single-molecule real-time sequencing reveals new targets for antimicrobials

Introduction

Pseudomonas aeruginosa (P.aeruginosa) is an important opportunistic pathogen with broad environmental adaptability and complex drug resistance. Single-molecule real-time (SMRT) sequencing technique has longer read-length sequences, more accuracy, and the ability to identify epigenetic DNA alterations.

Methods

This study applied SMRT technology to sequence a clinical strain P. aeruginosa PA3 to obtain its genome sequence and methylation modification information. Genomic, comparative, pan-genomic, and epigenetic analyses of PA3 were conducted.

Results

General genome annotations of PA3 were discovered, as well as information about virulence factors, regulatory proteins (RPs), secreted proteins, type II toxin-antitoxin (TA) pairs, and genomic islands. A genome-wide comparison revealed that PA3 was comparable to other P. aeruginosa strains in terms of identity, but varied in areas of horizontal gene transfer (HGT). Phylogenetic analysis showed that PA3 was closely related to P. aeruginosa 60503 and P. aeruginosa 8380. P. aeruginosa's pan-genome consists of a core genome of roughly 4,300 genes and an accessory genome of at least 5,500 genes. The results of the epigenetic analysis identified one main methylation sites, N6-methyladenosine (m6A) and 1 motif (CATNNNNNNNTCCT/AGGANNNNNNNATG). 16 meaningful methylated sites were picked. Among these, purH, phaZ, and lexA are of great significance playing an important role in the drug resistance and biological environment adaptability of PA3, and the targeting of these genes may benefit further antibacterial studies.

Disucssion

This study provided a detailed visualization and DNA methylation information of the PA3 genome and set a foundation for subsequent research into the molecular mechanism of DNA methyltransferase-controlled P. aeruginosa pathogenicity.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

New algorithms for accurate and efficient de novo genome assembly from long DNA sequencing reads

Building de novo genome assemblies for complex genomes is possible thanks to long-read DNA sequencing technologies. However, maximizing the quality of assemblies based on long reads is a challenging task that requires the development of specialized data analysis techniques. We present new algorithms for assembling long DNA sequencing reads from haploid and diploid organisms. The assembly algorithm builds an undirected graph with two vertices for each read based on minimizers selected by a hash function derived from the k-mer distribution. Statistics collected during the graph construction are used as features to build layout paths by selecting edges, ranked by a likelihood function. For diploid samples, we integrated a reimplementation of the ReFHap algorithm to perform molecular phasing. We ran the implemented algorithms on PacBio HiFi and Nanopore sequencing data taken from haploid and diploid samples of different species. Our algorithms showed competitive accuracy and computational efficiency, compared with other currently used software. We expect that this new development will be useful for researchers building genome assemblies for different species.

Genome-Wide Pathway Exploration of the Epidermidibacterium keratini EPI-7T

Functional cosmetics industries using skin microbiome screening and beneficial materials isolated from key microorganisms are receiving increasing attention. Since Epidermidibacterium keratini EPI-7T was first discovered in human skin, previous studies have confirmed that it can produce a new pyrimidine compound, 1,1′-biuracil, having anti-aging effects on human skin. Therefore, we conducted genomic analyses to judge the use value of E. keratini EPI-7T and provide up-to-date information. Whole-genome sequencing analysis of E. keratini EPI-7T was performed to generate new complete genome and annotation information. E. keratini EPI-7T genome was subjected to comparative genomic analysis with a group of closely-related strains and skin flora strains through bioinformatic analysis. Furthermore, based on annotation information, we explored metabolic pathways for valuable substances that can be used in functional cosmetics. In this study, the whole-genome sequencing (WGS) and annotation results of E. keratini EPI-7T were improved, and through comparative analysis, it was confirmed that the E. keratini EPI-7T has more metabolite-related genes than comparison strains. In addition, we annotated the vital genes for biosynthesis of 20 amino acids, orotic acid, riboflavin (B2) and chorismate. In particular, we were able to prospect that orotic acid could accumulate inside E. keratini EPI-7T under uracil-enriched conditions. Therefore, through a genomics approach, this study aims to provide genetic information for the hidden potential of E. keratini EPI-7T and the strain development and biotechnology utilization to be conducted in further studies.

Full-Length Transcriptome of the Great Himalayan Leaf-Nosed Bats (Hipposideros armiger) Optimized Genome Annotation and Revealed the Expression of Novel Genes

The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats' echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.

Full-length transcriptome of in Medicago sativa L. roots in response to drought stress

Background: Alfalfa (Medicago sativa L.), serves as a legume with high drought tolerance, is a major forage crop with a high biomass of production. However, the molecular mechanism of Alfalfa in response to drought stress are still unclear. Results: We constructed the first full-length transcriptome for Alfalfa root. 21.53Gb clean data were obtained by further data filtering, in which incorporate 566,076 reads of Insert (ROI), and 409,291 full length reads non-Chimeric (FLNC) sequences. Combined with second-generation sequencing (SGS), there were 2615, 6011, and 4617 differentially expressed genes (DEGs) in three comparisons. KEGG pathway analysis showed enrichment of ribosome, glutathione metabolism, and biosynthesis of amino acids are among the DEGs. The majority of transcription factors (TFs) from DEGs were AP2/ERF-ERF (37), C2H2 (32), and bHLH (22) bZIP (22), followed by C3H (19), MYB (18), WRKY (18), GRAS (16), and NAC (15). 32 C2H2 genes were differentially expressed in three groups. In addition, TFs annotated as C3H (19), MYB (18), GRAS (16), and NAC (15) also changed significantly in expression in the three comparisons. We found 24 genes participate in the abscisic acid (ABA) and auxin hormone signaling pathway in response to drought stress, and monitored the expression patterns of these related genes. Conclusion: The present study enhanced our understanding of the genetic diversity and complexity, and provides greater insight into the fundamental transcriptome reprogramming of Alfalfa under drought.

Chromosome-level genome assembly of goose provides insight into the adaptation and growth of local goose breeds

BACKGROUND

Anatidae contains numerous waterfowl species with great economic value, but the genetic diversity basis remains insufficiently investigated. Here, we report a chromosome-level genome assembly of Lion-head goose (Anser cygnoides), a native breed in South China, through the combination of PacBio, Bionano, and Hi-C technologies.

FINDINGS

The assembly had a total genome size of 1.19 Gb, consisting of 1,859 contigs with an N50 length of 20.59 Mb, generating 40 pseudochromosomes, representing 97.27% of the assembled genome, and identifying 21,208 protein-coding genes. Comparative genomic analysis revealed that geese and ducks diverged approximately 28.42 million years ago, and geese have undergone massive gene family expansion and contraction. To identify genetic markers associated with body weight in different geese breeds, including Wuzong goose, Huangzong goose, Magang goose, and Lion-head goose, a genome-wide association study was performed, yielding an average of 1,520.6 Mb of raw data that detected 44,858 single-mucleotide polymorphisms (SNPs). Genome-wide association study showed that 6 SNPs were significantly associated with body weight and 25 were potentially associated. The significantly associated SNPs were annotated as LDLRAD4, GPR180, and OR, enriching in growth factor receptor regulation pathways.

CONCLUSIONS

We present the first chromosome-level assembly of the Lion-head goose genome, which will expand the genomic resources of the Anatidae family, providing a basis for adaptation and evolution. Candidate genes significantly associated with different goose breeds may serve to understand the underlying mechanisms of weight differences.

Quasispecies Fitness Partition to Characterize the Molecular Status of a Viral Population. Negative Effect of Early Ribavirin Discontinuation in a Chronically Infected HEV Patient

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

Next generation sequencing technologies to explore the diversity of germplasm resources: achievements and trends in tomato

Tomato is one of the major vegetable crops grown worldwide and a model species for genetic and biological research. Progress in genomic technologies made possible the development of forefront methods for high-scale sequencing, providing comprehensive insight into the genetic architecture of germplasm resources. This review revisits next-generation sequencing strategies and applications to investigate the diversity of tomato, describing the common platforms used for SNP genotyping of large collections, de novo sequencing, and whole genome resequencing. Significant findings in evolutionary history are outlined, thus discussing how genomics has provided new hints about the processes behind domestication. Finally, achievement and perspectives on pan-genome construction and graphical pan-genome development toward precise mining of the natural variation to be exploited for breeding purposes are presented.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

PacBio Full-Length and Illumina Transcriptomes of the Gill Reveal the Molecular Response of Corbicula fluminea under Aerial Exposure

Air exposure is a common stress for Corbicula fluminea, an economically important freshwater shellfish consumed in China, during aquaculture and transportation. However, little is known about its molecular responses to air exposure. Therefore, this study used a combination of PacBio full-length and Illumina transcriptomes to investigate its molecular responses to air exposure. A total of 36,772 transcripts were obtained using PacBio sequencing. Structural analysis identified 32,069 coding sequences, 1906 transcription factors, 8873 simple sequence repeats, and 17,815 long non-coding RNAs. Subcellular localization analysis showed that most transcripts were located in the cytoplasm and nucleus. After 96-h of air exposure, 210 differentially expressed genes (DEGs) in the gill were obtained via Illumina sequencing. Among these DEGs, most of the genes related to glycolysis, tricarboxylic acid cycle, lipid metabolism, and amino acid metabolism were upregulated. Additionally, many DEGs associated with immunity, cytoskeleton reorganization, autophagy, and ferroptosis were identified. These findings indicated that metabolic strategy change, immune response, cytoskeleton reconstruction, autophagy, and ferroptosis might be the important mechanisms that C. fluminea use to cope with air exposure. This study will enrich the gene resources of C. fluminea and provide valuable data for studying the molecular mechanisms coping with air exposure in C. fluminea and other freshwater mollusks.

Nanopore Electrochemistry for Pathogen Detection

Pathogen infections have seriously threatened human health, and there is an urgent demand for rapid and efficient pathogen identification to provide instructions in clinical diagnosis and therapeutic intervention. Recently, nanopore technology, a rapidly maturing technology which delivers ultrasensitive sensing and high throughput in real-time and at low cost, has achieved success in pathogen detection. Furthermore, the remarkable development of nanopore sequencing, for example, the MinION sequencer from Oxford Nanopore Technologies (ONT) as a competitive sequencing technology, has facilitated the rapid analysis of disease-related microbiomes at the whole-genome level and on a large scale. Here, we highlighted recent advances in nanopore approaches for pathogen detection at the single-molecule level. We also overviewed the applications of nanopore sequencing in pathogenic bacteria identification and diagnosis. In the end, we discussed the challenges and future developments of nanopore technology as promising tools for the management of infections, which may be helpful to aid understanding as well as decision-making.

Epigenetic tumor heterogeneity in the era of single-cell profiling with nanopore sequencing

Nanopore sequencing has brought the technology to the next generation in the science of sequencing. This is achieved through research advancing on: pore efficiency, creating mechanisms to control DNA translocation, enhancing signal-to-noise ratio, and expanding to long-read ranges. Heterogeneity regarding epigenetics would be broad as mutations in the epigenome are sensitive to cause new challenges in cancer research. Epigenetic enzymes which catalyze DNA methylation and histone modification are dysregulated in cancer cells and cause numerous heterogeneous clones to evolve. Detection of this heterogeneity in these clones plays an indispensable role in the treatment of various cancer types. With single-cell profiling, the nanopore sequencing technology could provide a simple sequence at long reads and is expected to be used soon at the bedside or doctor's office. Here, we review the advancements of nanopore sequencing and its use in the detection of epigenetic heterogeneity in cancer.

Mining of Thousands of Prokaryotic Genomes Reveals High Abundance of Prophages with a Strictly Narrow Host Range

Phages and prophages are one of the principal modulators of microbial populations. However, much of their diversity is still poorly understood. Here, we extracted 33,624 prophages from 13,713 complete prokaryotic genomes to explore the prophage diversity and their relationships with their host. Our results reveal that prophages were present in 75% of the genomes studied. In addition, Enterobacterales were significantly enriched in prophages. We also found that pathogens are a significant reservoir of prophages. Finally, we determined that the prophage relatedness and the range of genomic hosts were delimited by the evolutionary relationships of their hosts. On a broader level, we got insights into the prophage population, identified in thousands of publicly available prokaryotic genomes, by comparing the prophage distribution and relatedness between them and their hosts.

IMPORTANCE Phages and prophages play an essential role in controlling their host populations either by modulating the host abundance or providing them with genes that benefit the host. The constant growth in next-generation sequencing technology has caused the development of powerful computational tools to identify phages and prophages with high precision. Making it possible to explore the prophage populations integrated into host genomes on a large scale. However, it is still a new and under-explored area, and efforts are still required to identify prophage populations to understand their dynamics with their hosts.

Migrating to Long-Read Sequencing for Clinical Routine BCR-ABL1 TKI Resistance Mutation Screening

Objective

The aim of this project was to implement long-read sequencing for BCR-ABL1 TKI resistance mutation screening in a clinical setting for patients undergoing treatment for chronic myeloid leukemia.

Materials and Methods

Processes were established for registering and transferring samples from the clinic to an academic sequencing facility for long-read sequencing. An automated analysis pipeline for detecting mutations was established, and an information system was implemented comprising features for data management, analysis and visualization. Clinical validation was performed by identifying BCR-ABL1 TKI resistance mutations by Sanger and long-read sequencing in parallel. The developed software is available as open source via GitHub at https://github.com/pharmbio/clamp.

Results

The information system enabled traceable transfer of samples from the clinic to the sequencing facility, robust and automated analysis of the long-read sequence data, and communication of results from sequence analysis in a reporting format that could be easily interpreted and acted upon by clinical experts. In a validation study, all 17 resistance mutations found by Sanger sequencing were also detected by long-read sequencing. An additional 16 mutations were found only by long-read sequencing, all of them with frequencies below the limit of detection for Sanger sequencing. The clonal distributions of co-existing mutations were automatically resolved through the long-read data analysis. After the implementation and validation, the clinical laboratory switched their routine protocol from using Sanger to long-read sequencing for this application.

Conclusions

Long-read sequencing delivers results with higher sensitivity compared to Sanger sequencing and enables earlier detection of emerging TKI resistance mutations. The developed processes, analysis workflow, and software components lower barriers for adoption and could be extended to other applications.

Analysis of direct-acting antiviral-resistant hepatitis C virus haplotype diversity by single-molecule and long-read sequencing

The method of analyzing individual resistant hepatitis C virus (HCV) by a combination of haplotyping and resistance-associated substitution (RAS) has not been fully elucidated because conventional sequencing has only yielded short and fragmented viral genomes. We performed haplotype analysis of HCV mutations in 12 asunaprevir/daclatasvir treatment-failure cases using the Oxford Nanopore sequencer. This enabled single-molecule long-read sequencing using rolling circle amplification (RCA) for correction of the sequencing error. RCA of the circularized reverse-transcription polymerase chain reaction products successfully produced DNA longer than 30 kilobase pairs (kb) containing multiple tandem repeats of a target 3 kb HCV genome. The long-read sequencing of these RCA products could determine the original sequence of the target single molecule as the consensus nucleotide sequence of the tandem repeats and revealed the presence of multiple viral haplotypes with the combination of various mutations in each host. In addition to already known signature RASs, such as NS3-D168 and NS5A-L31/Y93, there were various RASs specific to a different haplotype after treatment failure. The distribution of viral haplotype changed over time; some haplotypes disappeared without acquiring resistant mutations, and other haplotypes, which were not observed before treatment, appeared after treatment. Conclusion: The combination of various mutations other than the known signature RAS was suggested to influence the kinetics of individual HCV quasispecies in the direct-acting antiviral treatment. HCV haplotype dynamic analysis will provide novel information on the role of HCV diversity within the host, which will be useful for elucidating the pathological mechanism of HCV-related diseases.

Integrated SMRT and Illumina Sequencing Provide New Insights into Crocin Biosynthesis of Gardenia jasminoides

Crocins are valuable bioactive components of gardenia fruit, and their biosynthesis and accumulation have attracted widespread interest. Studies have investigated the biosynthesis and accumulation of crocin based on Illumina sequencing, but there is a lack of reports based on full-length transcriptome sequencing. Utilising SMRT sequencing and high-performance liquid chromatography (HPLC), we explored crocin biosynthesis and accumulation in the fruit of Gardenia jasminoides. HPLC analysis showed that crocins specifically exist in fruit and that the content of crocins increases gradually during fruit development. SMRT sequencing generated 46,715 high-quality full-length isoforms, including 5230 novel isoforms that are not present in the G. jasminoides genome. Furthermore, a total of 46 genes and 91 lncRNAs were involved in the biosynthesis and accumulation of crocin. The qRT-PCR indicated that genes involved in crocin biosynthesis reached a peak in the NOV stage. These findings contributed to our understanding of crocin biosynthesis and accumulation.

Whole-Genome Sequence and Comparative Analysis of Trichoderma asperellum ND-1 Reveal Its Unique Enzymatic System for Efficient Biomass Degradation

The lignocellulosic enzymes of Trichoderma asperellum have been intensely investigated toward efficient conversion of biomass into high-value chemicals/industrial products. However, lack of genome data is a remarkable hurdle for hydrolase systems studies. The secretory enzymes of newly isolated T. asperellum ND-1 during lignocellulose degradation are currently poorly known. Herein, a high-quality genomic sequence of ND-1, obtained by both Illumina HiSeq 2000 sequencing platforms and PacBio single-molecule real-time, has an assembly size of 35.75 Mb comprising 10,541 predicted genes. Secretome analysis showed that 895 proteins were detected, with 211 proteins associated with carbohydrate-active enzymes (CAZymes) responsible for biomass hydrolysis. Additionally, T. asperellum ND-1, T. atroviride IMI 206040, and T. virens Gv-298 shared 801 orthologues that were not identified in T. reesei QM6a, indicating that ND-1 may play critical roles in biological-control. In-depth analysis suggested that, compared with QM6a, the genome of ND-1 encoded a unique enzymatic system, especially hemicellulases and chitinases. Moreover, after comparative analysis of lignocellulase activities of ND-1 and other fungi, we found that ND-1 displayed higher hemicellulases (particularly xylanases) and comparable cellulases activities. Our analysis, combined with the whole-genome sequence information, offers a platform for designing advanced T. asperellum ND-1 strains for industrial utilizations, such as bioenergy production.

Full-Length Transcriptome Analysis of the Halophyte Nitraria sibirica Pall

Background: Nitraria sibirica Pall. is one of the pioneer tree species in saline–alkali areas due to its extreme salt tolerance. However, the lack of information on its genome limits the further exploration of the molecular mechanisms in N. sibirica under salt stress. Methods: In this study, we used single-molecule real-time (SMRT) technology based on the PacBio Iso-Seq platform to obtain transcriptome data from N. sibirica under salt treatment for the first time, which is helpful for our in-depth analysis of the salt tolerance and molecular characteristics of N. sibirica. Results: Our results suggested that a total of 234,508 circular consensus sequences (CCSs) with a mean read length of 2121 bp were obtained from the 19.26 Gb raw data. Furthermore, based on transcript cluster analysis, 93,713 consensus isoforms were obtained, including 92,116 high-quality isoforms. After removing redundant sequences, 49,240 non-redundant transcripts were obtained from high-quality isoforms. A total of 37,261 SSRs, 1816 LncRNAs and 47,314 CDSs, of which 40,160 carried complete ORFs, were obtained. Based on our transcriptome data, we also analyzed the coding genes of H+-PPase, and the results of both bioinformatics and functional analyses indicated that the gene prediction via full-length transcripts obtained by SMRT technology is reliable and effective. In summary, our research data obtained by SMRT technology provides more reliable and accurate information for the further analysis of the regulatory network and molecular mechanism of N. sibirica under salt stress.

Identification of genes involved in drought tolerance in seedlings of the desert grass, Psammochloa villosa (Poaceae), based on full-length isoform sequencing and de novo assembly from short reads

Psammochloa villosa is a perennial herbaceous plant that is dominant within arid regions of the Inner Mongolian Plateau and the Qinghai-Tibet Plateau in China, where it is an endemic species and exhibits strong drought tolerance and wind resistance. To study drought tolerance in P. villosa and determine its molecular basis, we simulated high and moderate drought stress in a controlled environment and then analyzed transcriptome sequences by combining long-read sequences from a representative, wild-grown individual with short reads from the treatment groups. We obtained 184,076 high-quality isoforms as a reference and 168,650 genes (91.6%), which we were able to annotate according to public databases. Ultimately, we obtained 119,005 unigenes representing the transcriptome of P. villosa under drought stress and, among these, we identified 3089 differentially expressed genes and 1484 transcription factors. Physiologically, P. villosa that was exposed to high and moderate drought stress had reduced germination rates and shorter buds but generated more chlorophyll, which is atypical under drought stress and possibly reflects an adaptation of these plants to their arid environment. We inferred that significantly upregulated genes were annotated as 'Chlorophyll a-b binding protein' and 'Light-harvesting chlorophyll-protein' among drought and control groups. Broadly, our analyses revealed that drought stress triggered many genome-level responses, especially related to mitigation of radical oxygen species (ROS), which increase in concentration under drought stress. In particular, in the high drought stress group compared with the control, GO enrichment analysis revealed a significant enrichment of upregulated genes (n = 10) involved in mitigation of oxidative stress. Similarly, using KEGG we found significant enrichment of genes in the phenylpropanoid biosynthesis pathway (11 genes), which yields phenols that scavenge ROS. We also inferred that many genes involved in metabolism of arginine and proline, which may serve as both scavengers of ROS and osmoprotectants that interact with stress response genes based on our protein-protein interaction network analysis. We verified the relative expression levels of eight genes associated with mitigation of ROS, DNA repair, and transmembrane transporter activity using qRT-PCR, and the results were consistent with our inferences from transcriptomes. This study provides insights into the genomic and physiological basis of drought tolerance in P. villosa and represents a resource for development of the species as a forage crop via molecular breeding within arid lands.

Current Updates on Expanded Carrier Screening: New Insights in the Omics Era

Genetic carrier screening has been successfully used over the last decades to identify individuals at risk of transmitting specific DNA variants to their newborns, thus having an affected child. Traditional testing has been offered based on familial and/or ethnic backgrounds. The development of high-throughput technologies, such as next-generations sequencing, able to allow the study of large genomic regions in a time and cost-affordable way, has moved carrier screening toward a more comprehensive and extensive approach, i.e., expanded carrier screening (ECS). ECS simultaneously analyses several disease-related genes and better estimates individuals’ carrier status. Indeed, it is not influenced by ethnicity and is not limited to a subset of mutations that may arise from poor information in some populations. Moreover, if couples carry out ECS before conceiving a baby, it allows them to obtain a complete estimation of their genetic risk and the possibility to make an informed decision regarding their reproductive life. Despite these advantages, some weakness still exists regarding, for example, the number of genes and the kind of diseases to be analyzed and the interpretation and communication of the obtained results. Once these points are fixed, it is expectable that ECS will become an ever more frequent practice in clinical settings.

Long-Reads-Based Metagenomics in Clinical Diagnosis With a Special Focus on Fungal Infections

Identification of the causative infectious agent is essential in the management of infectious diseases, with the ideal diagnostic method being rapid, accurate, and informative, while remaining cost-effective. Traditional diagnostic techniques rely on culturing and cell propagation to isolate and identify the causative pathogen. These techniques are limited by the ability and the time required to grow or propagate an agent in vitro and the facts that identification based on morphological traits are non-specific, insensitive, and reliant on technical expertise. The evolution of next-generation sequencing has revolutionized genomic studies to generate more data at a cheaper cost. These are divided into short- and long-read sequencing technologies, depending on the length of reads generated during sequencing runs. Long-read sequencing also called third-generation sequencing emerged commercially through the instruments released by Pacific Biosciences and Oxford Nanopore Technologies, although relying on different sequencing chemistries, with the first one being more accurate both platforms can generate ultra-long sequence reads. Long-read sequencing is capable of entirely spanning previously established genomic identification regions or potentially small whole genomes, drastically improving the accuracy of the identification of pathogens directly from clinical samples. Long-read sequencing may also provide additional important clinical information, such as antimicrobial resistance profiles and epidemiological data from a single sequencing run. While initial applications of long-read sequencing in clinical diagnosis showed that it could be a promising diagnostic technique, it also has highlighted the need for further optimization. In this review, we show the potential long-read sequencing has in clinical diagnosis of fungal infections and discuss the pros and cons of its implementation.

JC Polyomavirus whole genome sequencing at the single molecule level reveals emerging neurotropic populations in Progressive Multifocal Leucoencephalopathy

BACKGROUND

JC polyomavirus (JCV) mostly causes asymptomatic persistent renal infections but may give rise in immunosuppressed patients to neurotropic variants which replicate in the brain causing progressive multifocal leukoencephalopathy (PML). Rearrangements in the JCV genome regulator non-coding control region (NCCR) and missense mutations in the viral capsid VP1 gene differentiate neurotropic variants from virus excreted in urine.

METHODS

To investigate intra-host emergence of JCV neurotropic populations in PML, we deep sequenced JCV whole genome recovered from cerebrospinal fluid (CSF) and urine samples from 32 HIV- and non HIV-infected PML patients at the single-molecule level.

RESULTS

JCV strains distributed among 6 out of 7 known genotypes. Common patterns of NCCR rearrangements included an initial deletion mostly located in a short 10-nucleotide sequence, followed by duplications/insertions. Multiple NCCR variants present in individual CSF samples shared at least one rearrangement suggesting they stemmed from a unique viral population. NCCR variants independently acquired single or double PML-specific adaptive VP1 mutations. NCCR variants recovered from urine and CSF displayed opposite deletion or duplication patterns in binding sites for transcription factors.

DISCUSSION

Long read deep sequencing shed light on emergence of neurotropic JCV populations in PML.

Composition and Dynamics of H1N1 and H7N9 Influenza A Virus Quasispecies in a Co-infected Patient Analyzed by Single Molecule Sequencing Technology

A human co-infected with H1N1 and H7N9 subtypes influenza A virus (IAV) causes a complex infectious disease. The identification of molecular-level variations in composition and dynamics of IAV quasispecies will help to understand the pathogenesis and provide guidance for precision medicine treatment. In this study, using single-molecule real-time sequencing (SMRT) technology, we successfully acquired full-length IAV genomic sequences and quantified their genotypes abundance in serial samples from an 81-year-old male co-infected with H1N1 and H7N9 subtypes IAV. A total of 26 high diversity nucleotide loci was detected, in which the A-G base transversion was the most abundant substitution type (67 and 64%, in H1N1 and H7N9, respectively). Seven significant amino acid variations were detected, such as NA:H275Y and HA: R222K in H1N1 as well as PB2:E627K and NA: K432E in H7N9, which are related to viral drug-resistance or mammalian adaptation. Furtherly, we retrieved 25 H1N1 and 22 H7N9 genomic segment haplotypes from the eight samples based on combining high-diversity nucleotide loci, which provided a more concise overview of viral quasispecies composition and dynamics. Our approach promotes the popularization of viral quasispecies analysis in a complex infectious disease, which will boost the understanding of viral infections, pathogenesis, evolution, and precision medicine.

A whole genome duplication drives the genome evolution of Phytophthora betacei, a closely related species to Phytophthora infestans

BACKGROUND

Pathogens of the genus Phytophthora are the etiological agents of many devastating diseases in several high-value crops and forestry species such as potato, tomato, cocoa, and oak, among many others. Phytophthora betacei is a recently described species that causes late blight almost exclusively in tree tomatoes, and it is closely related to Phytophthora infestans that causes the disease in potato crops and other Solanaceae. This study reports the assembly and annotation of the genomes of P. betacei P8084, the first of its species, and P. infestans RC1-10, a Colombian strain from the EC-1 lineage, using long-read SMRT sequencing technology.

RESULTS

Our results show that P. betacei has the largest sequenced genome size of the Phytophthora genus so far with 270 Mb. A moderate transposable element invasion and a whole genome duplication likely explain its genome size expansion when compared to P. infestans, whereas P. infestans RC1-10 has expanded its genome under the activity of transposable elements. The high diversity and abundance (in terms of copy number) of classified and unclassified transposable elements in P. infestans RC1-10 relative to P. betacei bears testimony of the power of long-read technologies to discover novel repetitive elements in the genomes of organisms. Our data also provides support for the phylogenetic placement of P. betacei as a standalone species and as a sister group of P. infestans. Finally, we found no evidence to support the idea that the genome of P. betacei P8084 follows the same gene-dense/gense-sparse architecture proposed for P. infestans and other filamentous plant pathogens.

CONCLUSIONS

This study provides the first genome-wide picture of P. betacei and expands the genomic resources available for P. infestans. This is a contribution towards the understanding of the genome biology and evolutionary history of Phytophthora species belonging to the subclade 1c.

Unraveling the complexity of faba bean (Vicia faba L.) transcriptome to reveal cold-stress-responsive genes using long-read isoform sequencing technology

Faba bean (Vicia faba L.), a globally important grain legume providing a stable source of dietary protein, was one of the earliest plant cytogenetic models. However, the lack of draft genome annotations and unclear structural information on mRNA transcripts have impeded its genetic improvement. To address this, we sequenced faba bean leaf transcriptome using the PacBio single-molecule long-read isoform sequencing platform. We identified 28,569 nonredundant unigenes, ranging from 108 to 9669 bp, with a total length of 94.5 Mb. Many unigenes (3597, 12.5%) had 2-20 isoforms, indicating a highly complex transcriptome. Approximately 96.5% of the unigenes matched sequences in public databases. The predicted proteins and transcription factors included NB-ARC, Myb_domain, C3H, bHLH, and heat shock proteins, implying that this genome has an abundance of stress resistance genes. To validate our results, we selected WCOR413-15785, DHN2-12403, DHN2-14197, DHN2-14797, COR15-14478, and HVA22-15 unigenes from the ICE-CBF-COR pathway to analyze their expression patterns in cold-treated samples via qRT-PCR. The expression of dehydrin-related genes was induced by cold stress. The assembled data provide the first insights into the deep sequencing of full-length RNA from faba bean at the single-molecule level. This study provides an important foundation to improve gene modeling and protein prediction.

Nematode Genome Announcement: Draft Genome of Meloidogyne chitwoodi, an Economically Important Pest of Potato in the Pacific Northwest

Meloidogyne chitwoodi is one of the most devastating pests of potato in the U.S. Pacific Northwest (PNW). Nematode-infected tubers develop external as well as internal defects, making the potato tubers unmarketable, and resulting in economic losses. Draft genome assemblies of three M. chitwoodi genotypes-race 1, race 2 and race 1 pathotype Roza-were generated using Illumina and PacBio Sequel RS II sequencing. The final assemblies consist of 30, 39, and 38 polished contigs for race 1, race 2 and race 1 pathotype Roza, respectively, with average N₅₀ of 2.37 Mb and average assembled genome size of approximately 47.41 Mb. On average, 10,508 genes were annotated for each genome. Benchmarking universal single-copy ortholog (BUSCO) analysis indicated that 69.80% of the BUSCOs were complete whereas 68.80, 0.93, and 12.67% were single copy, duplicated, and fragmented, respectively. These highly contiguous genomes will enrich resources to study potato-nematode interactions and enhance breeding efforts to develop nematode-resistant potato varieties for the PNW.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome sequence, transcriptome, and annotation of rodent malaria parasite Plasmodium yoelii nigeriensis N67

Background

Rodent malaria parasites are important models for studying host-malaria parasite interactions such as host immune response, mechanisms of parasite evasion of host killing, and vaccine development. One of the rodent malaria parasites is Plasmodium yoelii, and multiple P. yoelii strains or subspecies that cause different disease phenotypes have been widely employed in various studies. The genomes and transcriptomes of several P. yoelii strains have been analyzed and annotated, including the lethal strains of P. y. yoelii YM (or 17XL) and non-lethal strains of P. y. yoelii 17XNL/17X. Genomic DNA sequences and cDNA reads from another subspecies P. y. nigeriensis N67 have been reported for studies of genetic polymorphisms and parasite response to drugs, but its genome has not been assembled and annotated.

Results

We performed genome sequencing of the N67 parasite using the PacBio long-read sequencing technology, de novo assembled its genome and transcriptome, and predicted 5383 genes with high overall annotation quality. Comparison of the annotated genome of the N67 parasite with those of YM and 17X parasites revealed a set of genes with N67-specific orthology, expansion of gene families, particularly the homologs of the Plasmodium chabaudi erythrocyte membrane antigen, large numbers of SNPs and indels, and proteins predicted to interact with host immune responses based on their functional domains.

Conclusions

The genomes of N67 and 17X parasites are highly diverse, having approximately one polymorphic site per 50 base pairs of DNA. The annotated N67 genome and transcriptome provide searchable databases for fast retrieval of genes and proteins, which will greatly facilitate our efforts in studying the parasite biology and gene function and in developing effective control measures against malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07555-9.

NMD inhibition by 5-azacytidine augments presentation of immunogenic frameshift-derived neoepitopes

Frameshifted protein sequences elicit tumor-specific T cell-mediated immune responses in microsatellite-unstable (MSI) cancers if presented by HLA class I molecules. However, their expression and presentation are limited by nonsense-mediated RNA decay (NMD). We employed an unbiased immunopeptidomics workflow to analyze MSI HCT-116 cells and identified >10,000 HLA class I-presented peptides including five frameshift-derived InDel neoepitopes. Notably, pharmacological NMD inhibition with 5-azacytidine stabilizes frameshift-bearing transcripts and increases the HLA class I-mediated presentation of InDel neoepitopes. The frameshift mutation underlying one of the identified InDel neoepitopes is highly recurrent in MSI colorectal cancer cell lines and primary patient samples, and immunization with the corresponding neoepitope induces strong CD8+ T cell responses in an HLA-A∗02:01 transgenic mouse model. Our data show directly that pharmacological NMD inhibition augments HLA class I-mediated presentation of immunogenic frameshift-derived InDel neoepitopes thus highlighting the clinical potential of NMD inhibition in anti-cancer immunotherapy strategies.

16S rRNA long-read nanopore sequencing is feasible and reliable for endometrial microbiome analysis

RESEARCH QUESTION

Full-length 16S rRNA gene sequencing using nanopore technology is a fast alternative to conventional short-read 16S rRNA gene sequencing with low initial investment costs that has been used for various microbiome studies but has not yet been investigated as an alternative approach for endometrial microbiome analysis. Is in-situ 16S rRNA gene long-read sequencing using portable nanopore sequencing technology feasible and reliable for endometrial microbiome analysis?

DESIGN

A prospective experimental study based on 33 patients seeking infertility treatment between January and October 2019. A 16S rRNA gene long-read nanopore sequencing protocol for analysing endometrial microbiome samples was established, including negative controls for contamination evaluation and positive controls for bias evaluation. Contamination caused by kit and exterior sources was identified and excluded from the analysis. Endometrial samples from 33 infertile patients were sequenced using the optimized long-read nanopore sequencing protocol and compared with conventional short-read sequencing carried out by external laboratories.

RESULTS

Of the 33 endometrial patient samples, 23 successfully amplified (69.7%) and their microbiome was assessed using nanopore sequencing. Of those 23 samples, 14 (60.9%) were Lactobacillus-dominated (>80% of reads mapping to Lactobacillus), with 10 samples resulting in more than 90% Lactobacillus reads. Our long-read nanopore sequencing revealed results similar to two conventional short-read sequencing approaches and to long-read sequencing validation carried out in external laboratories.

CONCLUSION

In this pilot study, 16S rRNA gene long-read nanopore sequencing was established to analyse the endometrial microbiome in situ that could be widely applied owing to its cost efficiency and portable character.

Whole genome sequencing and comparative genomics analysis of Pectobacterium carotovorum identifies key pathogenic genes

Based on Single moleculereal time(SMRT)sequencing technology, the high-quality whole genome sequence of Pectobacterium carotovorum (PC1) was obtained by the PacBio RS II sequencer. The genome is a single circular chromosome of 5.3 Mb in size, containing three kinds of m6A methylation modification by SMRT Portal analysis. Genome annotation showed that 575 virulence factor genes, 304 drug resistance genes, 774 pathogen genes, 7 secretory systems and 22 pairs of two-component regulatory system could be relevant to bacterial pathogenicity. In addition, the average nucleotide identities (ANI) analysisshowed that the PC1 exhibited the highest homology with the Pectobacteriumcarotovorumsubsp.carotovorumstrain BP201601.1 (NZ_CP034236). There are 28 unique gene families to PC1 using cluster analysis of gene families. According to the analysis of key pathogenic genes, we have obtained three kinds of highly conserved genes related to cell wall degrading enzymes, including 19 pectinase genes, 25 cellulase genes and 22 protease genes. Our studies have provided a theoretical basis for investigation of bacterial soft rot and biological specific bactercides of PC1.

Comparison of Structural and Short Variants Detected by Linked-Read and Whole-Exome Sequencing in Multiple Myeloma

Linked-read sequencing was developed to aid the detection of large structural variants (SVs) from short-read sequencing efforts. We performed a systematic evaluation to determine if linked-read exome sequencing provides more comprehensive and clinically relevant information than whole-exome sequencing (WES) when applied to the same set of multiple myeloma patient samples. We report that linked-read sequencing detected a higher number of SVs (n = 18,455) than WES (n = 4065). However, linked-read predictions were dominated by inversions (92.4%), leading to poor detection of other types of SVs. In contrast, WES detected 56.3% deletions, 32.6% insertions, 6.7% translocations, 3.3% duplications and 1.2% inversions. Surprisingly, the quantitative performance assessment suggested a higher performance for WES (AUC = 0.791) compared to linked-read sequencing (AUC = 0.766) for detecting clinically validated cytogenetic alterations. We also found that linked-read sequencing detected more short variants (n = 704) compared to WES (n = 109). WES detected somatic mutations in all MM-related genes while linked-read sequencing failed to detect certain mutations. The comparison of somatic mutations detected using linked-read, WES and RNA-seq revealed that WES and RNA-seq detected more mutations than linked-read sequencing. These data indicate that WES outperforms and is more efficient than linked-read sequencing for detecting clinically relevant SVs and MM-specific short variants.

Verification of CRISPR editing and finding transgenic inserts by Xdrop indirect sequence capture followed by short- and long-read sequencing

Validation of CRISPR-Cas9 editing typically explores the immediate vicinity of the gene editing site and distal off-target sequences, which has led to the conclusion that CRISPR-Cas9 editing is very specific. However, an increasing number of studies suggest that on-target unintended editing events like deletions and insertions are relatively frequent but unfortunately often missed in the validation of CRISPR-Cas9 editing. The deletions may be several kilobases-long and only affect one allele. The gold standard in molecular validation of gene editing is direct sequencing of relatively short PCR amplicons. This approach allows the detection of small editing events but fails in detecting large rearrangements, in particular when only one allele is affected. Detection of large rearrangements requires that an extended region is analyzed and the characterization of events may benefit from long-read sequencing. Here we implemented Xdrop™, a new microfluidic technology that allows targeted enrichment of long regions (~100 kb) using just a single standard PCR primer set. Sequencing of the enriched CRISPR-Cas9 gene-edited region in four cell lines on long- and short-read sequencing platforms unravelled unknown and unintended genome editing events. The analysis revealed accidental kilobases-large insertions in three of the cell lines, which remained undetected using standard procedures. We also applied the targeted enrichment approach to identify the integration site of a transgene in a mouse line. The results demonstrate the potential of this technology in gene editing validation as well as in more classic transgenics.

Recent Advanced Technologies for the Characterization of Xenobiotic-Degrading Microorganisms and Microbial Communities

Global environmental contamination with a complex mixture of xenobiotics has become a major environmental issue worldwide. Many xenobiotic compounds severely impact the environment due to their high toxicity, prolonged persistence, and limited biodegradability. Microbial-assisted degradation of xenobiotic compounds is considered to be the most effective and beneficial approach. Microorganisms have remarkable catabolic potential, with genes, enzymes, and degradation pathways implicated in the process of biodegradation. A number of microbes, including Alcaligenes, Cellulosimicrobium, Microbacterium, Micrococcus, Methanospirillum, Aeromonas, Sphingobium, Flavobacterium, Rhodococcus, Aspergillus, Penecillium, Trichoderma, Streptomyces, Rhodotorula, Candida, and Aureobasidium, have been isolated and characterized, and have shown exceptional biodegradation potential for a variety of xenobiotic contaminants from soil/water environments. Microorganisms potentially utilize xenobiotic contaminants as carbon or nitrogen sources to sustain their growth and metabolic activities. Diverse microbial populations survive in harsh contaminated environments, exhibiting a significant biodegradation potential to degrade and transform pollutants. However, the study of such microbial populations requires a more advanced and multifaceted approach. Currently, multiple advanced approaches, including metagenomics, proteomics, transcriptomics, and metabolomics, are successfully employed for the characterization of pollutant-degrading microorganisms, their metabolic machinery, novel proteins, and catabolic genes involved in the degradation process. These technologies are highly sophisticated, and efficient for obtaining information about the genetic diversity and community structures of microorganisms. Advanced molecular technologies used for the characterization of complex microbial communities give an in-depth understanding of their structural and functional aspects, and help to resolve issues related to the biodegradation potential of microorganisms. This review article discusses the biodegradation potential of microorganisms and provides insights into recent advances and omics approaches employed for the specific characterization of xenobiotic-degrading microorganisms from contaminated environments.

Full-length transcriptome sequences of ridgetail white prawn Exopalaemon carinicauda provide insight into gene expression dynamics during thermal stress

Marine heat waves and extreme high temperature become more frequent and intense in these years, which affected the survival of aquaculture animals. The ridgetail white prawn Exopalaemon carinicauda is an important economic species in eastern China, which has remarkable thermal tolerance. However, there has been little study of its thermal-adaptation mechanisms due to the complex genetic structure and unknown genome. To better understand the molecular mechanisms of E. carinicauda to adapt to the changing temperature, a combination of Illumina-based short reads RNA-seq and single molecule real-time-based full-length transcriptome sequencing was used in this study. In total, 17,212 unigenes from high-quality transcripts of E. carinicauda were generated and 14,663 complete ORFs were detected with an average length of 1980 bp. In addition, the transcriptome profiles of E. carinicauda treated with 34 °C heat stress for 6 and 24 h were analyzed. These differentially expressed genes were primarily enriched in oxidation-reduction process (Gene Ontology enrichment, GO) and the pathways of starch and sucrose metabolism (Kyoto Encyclopedia of Genes and Genomes enrichment, KEGG) after 6 h thermal stress, which indicated that E. carinicauda was suffering the attack by reactive oxygen species. After 24 h thermal stress, these differentially expressed genes were enriched in the pathway of lysosome, glycine, serine and threonine metabolism, fatty acid metabolism (KEGG), which indicated the oxidative stress was decreased. Interestingly, 40 genes for hemocyanin were found to be downregulated after 6 h heat stress, which indicated that the immunocompetence of E. carinicauda decreased after short term thermal stress (6 h). After 24 h thermal stress, E. carinicauda showed transcriptional adaptation to high temperature by upregulating of 11 genes encoding molecular chaperones, including HSP40 and HSP90 which were firstly reported to be related to thermal stress in E. carinicauda. These results promote a better understanding of the thermal-adaptation mechanism of E. carinicauda.

Antarctic desert soil bacteria exhibit high novel natural product potential, evaluated through long-read genome sequencing and comparative genomics

Actinobacteria and Proteobacteria are important producers of bioactive natural products (NP), and these phyla dominate in the arid soils of Antarctica, where metabolic adaptations influence survival under harsh conditions. Biosynthetic gene clusters (BGCs) which encode NPs, are typically long and repetitious high G + C regions difficult to sequence with short-read technologies. We sequenced 17 Antarctic soil bacteria from multi-genome libraries, employing the long-read PacBio platform, to optimize capture of BGCs and to facilitate a comprehensive analysis of their NP capacity. We report 13 complete bacterial genomes of high quality and contiguity, representing 10 different cold-adapted genera including novel species. Antarctic BGCs exhibited low similarity to known compound BGCs (av. 31%), with an abundance of terpene, non-ribosomal peptide and polyketide-encoding clusters. Comparative genome analysis was used to map BGC variation between closely related strains from geographically distant environments. Results showed the greatest biosynthetic differences to be in a psychrotolerant Streptomyces strain, as well as a rare Actinobacteria genus, Kribbella, while two other Streptomyces spp. were surprisingly similar to known genomes. Streptomyces and Kribbella BGCs were predicted to encode antitumour, antifungal, antibacterial and biosurfactant-like compounds, and the synthesis of NPs with antibacterial, antifungal and surfactant properties was confirmed through bioactivity assays.

Differential structures and enterotype-like clusters of Bifidobacterium responses to probiotic fermented milk consumption across subjects using a Bifidobacterium-target procedure

The health-promoting attributes of bifidobacteria have piqued interest of researchers worldwide. However, scant published studies are available pertinent to bifidobacteria in microbiota/metagenomics datasets due to its intrinsic low abundance and limitations of detection methods. In this work, we designed a procedure to optimize the detection of the bifidobacterial population in complex biological samples with single-molecule real-time sequencing (SMRT) technology, including one primer pair designated as Bif-6 and a Bifidobacterium-specific database. The optimized procedure detected 14 bifidobacterial species/subspecies in ten human stool samples (2024 sequences per sample) and eight breast milk samples (3473 sequences per sample), respectively. Furthermore, by using the optimized procedure of SMRT, we investigated the effect of a 4-week-intervention of probiotic fermented milk (PFM; 200 g/day) on the gut bifidobacteria population of adults. The results showed that consuming PFM changed the structure and enterotype-like clusters of Bifidobacterium. After the consumption of PFM, the level of gut Bifidobacterium animalis increased significantly, replacing several originally dominating taxa in some subjects, including B. catenulatum, B. breve, and B. bifidum. On the other hand, B. adolescentis was, unaffectedly, the representative species in subjects having an original enterotype-like cluster of B. adolescentis. In conclusion, our work designed a procedure for detecting the bifidobacterial population in complex samples. By applying the currently designed procedure, we found that the PFM intervention changed the bifidobacterial enterotype-like cluster of some subjects, and such change was dependent on the basal bifidobacterial population.

The shared frameshift mutation landscape of microsatellite-unstable cancers suggests immunoediting during tumor evolution

The immune system can recognize and attack cancer cells, especially those with a high load of mutation-induced neoantigens. Such neoantigens are abundant in DNA mismatch repair (MMR)-deficient, microsatellite-unstable (MSI) cancers. MMR deficiency leads to insertion/deletion (indel) mutations at coding microsatellites (cMS) and to neoantigen-inducing translational frameshifts. Here, we develop a tool to quantify frameshift mutations in MSI colorectal and endometrial cancer. Our results show that frameshift mutation frequency is negatively correlated to the predicted immunogenicity of the resulting peptides, suggesting counterselection of cell clones with highly immunogenic frameshift peptides. This correlation is absent in tumors with Beta-2-microglobulin mutations, and HLA-A*02:01 status is related to cMS mutation patterns. Importantly, certain outlier mutations are common in MSI cancers despite being related to frameshift peptides with functionally confirmed immunogenicity, suggesting a possible driver role during MSI tumor evolution. Neoantigens resulting from shared mutations represent promising vaccine candidates for prevention of MSI cancers.