Molecular mechanism of Serratia marcescens Bizio infection in Reticulitermes chinensis Snyder based on full-length SMRT transcriptome sequencing

Reticulitermes chinensis Snyder is an important pest in forestry and construction and is widely distributed in China. We found that Serratia marcescens Bizio strain SM1 has insecticidal activity to R. chinensis, but the pathogenic mechanism of SM1 to R. chinensis is not clear. Therefore, full-length transcriptome sequencing was performed on R. chinensis infected with SM1 and the control group. A total of 230 differentially expressed genes were identified by comparing SM1 infection group and the control group, among which 103 were downregulated and 127 were upregulated. We found downregulated genes in nine metabolic pathway categories, among which carbohydrate metabolism had the most downregulated genes, followed by energy metabolism and amino acid metabolism. We also found that some downregulated genes were related to pattern recognition receptors, cellular immunity, and humoral immunity, indicating that R. chinensis immunity was negatively affected by SM1 infection. In addition, some genes in signal transduction and genetic information processing pathways were downregulated. In this study, high-throughput full-length transcriptome analysis was used to analyse the pathogenic mechanism of SM1 to R. chinensis. The results of this study provide useful information for exploring the relationship between SM1 and R. chinensis, and provide theoretical support for the future application of SM1 and the prevention and treatment of R. chinensis.

The novel alleles, HLA-G*01:04:09 and HLA-DPB1*04:01:01:136, detected in a hematopoietic cell donor

Two novel alleles, HLA-G*01:04:09 and HLA-DPB1*04:01:01:136, were identified in a single healthy individual.

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

The novel HLA-F*01:16 and HLA-F*01:17 alleles identified in hematopoietic cell donors

Two novel non-classical HLA class I alleles have been characterized, HLA-F*01:16 and -F*01:17.

Fifty novel HLA-DPB1 alleles identified in a UK cohort of unrelated hematopoietic cell donors and recipients

HLA-DPB1 is the classical HLA class II genes with the least recorded variation on the IPD-IMGT/HLA Database, suggesting the full extent of its diversity is perhaps yet to be characterized. Here, a full-gene typing strategy was employed to genotype a UK cohort of 1470 HCT recipients (n = 744) and donors (n = 726). In total, 2940 full-length HLA-DPB1 sequences were generated, comprising 193 distinct alleles. Of these, 107 sequences contained novel variation, totaling 49 unique intronic HLA-DPB1 alleles, and one coding variant (HLA-DPB1*1188:01). Full-gene sequencing resulted in zygosity changes for 129 individuals by identifying two distinct intronic variants of the same coding allele. We verified the existence of nine unconfirmed alleles and extended the sequence of two existing alleles on the IPD-IMGT/HLA Database.

Mycobacterium tuberculosis complex molecular networks and their regulation: Implications of strain heterogeneity on epigenetic diversity and transcriptome regulation

Tuberculosis has been a public health crisis since the 1900, which has caused the highest mortalities due to a single bacterial infection worldwide, that was recently further complicated by the Coronavirus disease 2019 pandemic. The causative agent of Tuberculosis, Mycobacterium tuberculosis, belongs to a genetically well-characterized family of strains known as the Mycobacterium tuberculosis complex, which has complicated progress made towards eradicating Tuberculosis due to pathogen-specific phenotypic differences in the members of this complex. Mycobacterium tuberculosis complex strains are genetically diverse human- and animal-adapted pathogens belonging to 7 lineages (Indo-Oceanic, East-Asian, East-African Indian, Euro-American, M. africanum West Africa 1, M. africanum West Africa 2 and Ethopia), respectively and the recently identified Lineage 8 and M. africanum Lineage 9. Genomic studies have revealed that Mycobacterium tuberculosis complex members are ∼99 % similar, however, due to selective pressure and adaptation to human host, they are prone to mutations that have resulted in development of drug resistance and phenotypic heterogeneity that impact strain virulence. Furthermore, members of the Mycobacterium tuberculosis complex have preferred geographic locations and possess unique phenotypic characteristics that is linked to their pathogenicity. Due to the recent advances in development next generation sequencing platforms, several studies have revealed epigenetic changes in genomic regions combined with "unique" gene regulatory mechanisms through non-coding RNAs that are responsible for strain-specific behaviour on in vitro and in vivo infection models. The current review provides up to date epigenetic patterns, gene regulation through non-coding RNAs, together with implications of these mechanisms in down-stream proteome and metabolome, which may be responsible for "unique" responses to infection by members of the Mycobacterium tuberculosis complex. Understanding lineage-specific molecular mechanisms during infection may provide novel drug targets and disease control measures towards World Health organization END-TB strategy.

Construction of heat stress regulation networks based on Illumina and SMRT sequencing data in potato

Potato (Solanum tuberosum L.) is one of the most important tuber food crops in the world; however, the cultivated potatoes are susceptible to high temperature, by which potato production is adversely affected. Understanding the coping mechanism of potato to heat stress is essential to secure yield and expand adaptability under environmental conditions with rising temperature. However, the lack of heat-related information has significantly limited the identification and application of core genes. To gain deeper insights into heat tolerance genes, next-generation sequencing and single-molecule real-time sequencing were used to learn the transcriptional response of potato to heat stress and 13,159 differentially expressed genes (DEGs) were identified in this study. All DEGs were grouped into 12 clusters using the K-means clustering algorithm. Gene Ontology enrichment analysis revealed that they were involved in temperature signaling, phytohormone, and protein modification. Among them, there were 950 differentially expressed transcription factors (DETFs). According to the network analysis of DETFs at the sixth hour under heat stress, we found some genes that were previously reported to be associated with photoperiodic tuberization, StCO (CONSTANS), tuber formation, StBEL11 (BEL1-LIKE 11), and earliness in potato, StCDF1 (CYCLING DOF FACTOR 1) responding to temperature. Furthermore, we verified the relative expression levels using quantitative real-time polymerase chain reaction, and the results were consistent with the inferences from transcriptomes. In addition, there were 22,125 alternative splicing events and 2,048 long non-coding RNAs. The database and network established in this study will extend our understanding of potato response to heat stress. It ultimately provided valuable resources for molecular analysis of heat stress response in potato and cultivation of potato varieties with heat tolerance.

A Survey of Archaeal Restriction-Modification Systems

When compared with bacteria, relatively little is known about the restriction-modification (RM) systems of archaea, particularly those in taxa outside of the haloarchaea. To improve our understanding of archaeal RM systems, we surveyed REBASE, the restriction enzyme database, to catalog what is known about the genes and activities present in the 519 completely sequenced archaeal genomes currently deposited there. For 49 (9.4%) of these genomes, we also have methylome data from Single-Molecule Real-Time (SMRT) sequencing that reveal the target recognition sites of the active m6A and m4C DNA methyltransferases (MTases). The gene-finding pipeline employed by REBASE is trained primarily on bacterial examples and so will look for similar genes in archaea. Nonetheless, the organizational structure and protein sequence of RM systems from archaea are highly similar to those of bacteria, with both groups acquiring systems from a shared genetic pool through horizontal gene transfer. As in bacteria, we observe numerous examples of "persistent" DNA MTases conserved within archaeal taxa at different levels. We experimentally validated two homologous members of one of the largest "persistent" MTase groups, revealing that methylation of C(m5C)WGG sites may play a key epigenetic role in Crenarchaea. Throughout the archaea, genes encoding m6A, m4C, and m5C DNA MTases, respectively, occur in approximately the ratio 4:2:1.

SMRT Sequencing Enables High-Throughput Identification of Novel AAVs from Capsid Shuffling and Directed Evolution

The use of AAV capsid libraries coupled with various selection strategies has proven to be a remarkable approach for generating novel AAVs with enhanced and desired features. The inability to reliably sequence the complete capsid gene in a high-throughput manner has been the bottleneck of capsid engineering. As a result, many library strategies are confined to localized and modest alterations in the capsid, such as peptide insertions or single variable region (VR) alterations. The caveat of short reads by means of next-generation sequencing (NGS) hinders the diversity of capsid library construction, shifting the field away from whole-capsid modifications. We generated AAV capsid shuffled libraries of naturally occurring AAVs and applied directed evolution in both mice and non-human primates (NHPs), with the goal of yielding AAVs that are compatible across both species for translational applications. We recovered DNA from the tissues of injected animal and used single molecule real-time (SMRT) sequencing to identify variants enriched in the central nervous system (CNS). We provide insights and considerations for variant identification by comparing bulk tissue sequencing to that of isolated nuclei. Our work highlights the potential advantages of whole-capsid engineering, as well as indispensable methodological improvements for the analysis of recovered capsids, including the nuclei-enrichment step and SMRT sequencing.

Utilization of multiple genetic methods for prenatal diagnosis of rare thalassemia variants

Background: Thalassemia is the most prevalent monogenic disorder caused by an imbalance between the α- and β-globin chains as a result of pathogenic variants in the α- or β-globin genes. Novel or complex structural changes in globin genes are major hurdles for genetic consulting and prenatal diagnosis. Methods: From 2020 to 2022, genetic analysis was performed on 1,316 families suspected of having children with thalassemia major, including 42 pregnant couples suspected of being thalassemia carriers with rare variants. Multiple techniques including multiplex ligation-dependent probe amplification (MLPA), Sanger sequencing, targeted next-generation sequencing, and single-molecule real-time (SMRT) sequencing were used to diagnose rare thalassemia. Results: The rate of prenatal diagnosis for rare thalassemia variants was 3.19% (42/1,316). The most prevalent alleles of α- and β-thalassemia are Chinese ^Gγ(^Aγδβ)⁰and -- ^THAI deletion. In addition, ten rare complex genotypes include one Chinese ^Gγ(^Aγδβ)⁰ deletion combined with HBG1-HBG2 fusion, two rare deletions at HBB gene (hg38, Chr11: 5224211-5232470, hg38, Chr11: 5224303-5227790), one complete 7,412 bp fusion gene for anti-Lepore Hong Kong, two complex rearrangements of the α-globin gene cluster, two novel duplications, and two rare large deletions in the α-globin gene cluster. Conclusion: Accurate gene diagnosis for probands with combined molecular biology techniques is the key to prenatal diagnosis of rare thalassemia.

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.

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.

Combined CRISPR toolkits reveal the domestication landscape and function of the ultra-long and highly repetitive silk genes

Highly repetitive sequences play a major structural and function role in the genome. In the present study, we developed Cas9-assisted cloning and SMRT sequencing of long repetitive sequences (CACS) to sequence and manipulate highly repetitive genes from eukaryotic genomes. CACS combined Cas9-mediated cleavage of a target segment from an intact genome, Gibson assembly cloning, and PacBio SMRT sequencing. Applying CACS, we directly cloned and sequenced the complete sequences of fibroin heavy chain (FibH) genes from 17 domesticated (Bombyx mori) and 7 wild (Bombyx mandarina) silkworms. Our analysis revealed the unique fine structure organization, genetic variations, and domestication dynamics of FibH. We also demonstrated that the length of the repetitive regions determined the mechanical properties of silk fiber, which was further confirmed by Cas9 editing of FibH. CACS is a simple, robust, and efficient approach, providing affordable accessibility to highly repetitive regions of a genome. STATEMENT OF SIGNIFICANCE: Silkworm silk is the earliest and most widely used animal fiber, and its excellent performance mainly depends on the fibroin heavy chain (FibH) protein. The FibH gene is the main breakthrough in understanding the formation mechanism and improvement of silk fiber. In the study, we developed a CACS method for characterizing the fine structure and domestication landscape of 24 silkworm FibH genes. We used CRISPR/Cas9 to edit the repetitive sequence of FibH genes, revealing the relationship between FibH genes and mechanical properties of silkworm silk. Our study is helpful in modifying silk genes to manipulate other valuable highly repetitive sequences, and provides insight for silkworm breeding.

HLA/MHC and KIR characterization in humans and non-human primates using Oxford Nanopore Technologies and Pacific Biosciences sequencing platforms

The gene products of the HLA/MHC and KIR multigene families are important modulators of the immune system and are associated with health and disease. Characterization of the genes encoding these receptors has been integrated into different biomedical applications, including transplantation and reproduction biology, immune therapies and in fundamental research into disease susceptibility or resistance. Conventional short-read sequencing strategies have shown their value in high throughput typing, but are insufficient to uncover the entire complexity of the highly polymorphic HLA/MHC and KIR gene systems. The implementation of single-molecule and real-time sequencing platforms, offered by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), revolutionized the fields of genomics and transcriptomics. Using fundamentally distinct principles, these platforms generate long-read data that can unwire the plasticity of the HLA/MHC and KIR genes, including high-resolution characterization of genes, alleles, phased haplotypes, transcription levels and epigenetics modification patterns. These insights might have profound clinical relevance, such as improved matching of donors and patients in clinical transplantation, but could also lift disease association studies to a higher level. Even more, a comprehensive characterization may refine animal models in preclinical studies. In this review, the different HLA/MHC and KIR characterization approaches using PacBio and ONT platforms are described and discussed.

Full-Length Transcriptome Analysis of the Ichthyotoxic Harmful Alga Heterosigma akashiwo (Raphidophyceae) Using Single-Molecule Real-Time Sequencing

The raphidophyte Heterosigma akashiwo is a harmful algal species. The bloom of this organism has been associated with the massive mortality of fish in many coastal waters. To investigate the molecular mechanism of H. akashiwo blooms, having a reliable reference transcriptome of this species is essential. Therefore, in this study, a full-length transcriptome of H. akashiwo was obtained by single-molecule real-time sequencing. In total, 45.44 Gb subread bases were generated, and 16,668 unigenes were obtained after the sequencing data processing. A total of 8666 (52.00%) unigenes were successfully annotated using seven public databases. Among them, mostly phosphorus and nitrogen metabolism genes were detected. Moreover, there were 300 putative transcription factors, 4392 putative long non-coding RNAs, and 7851 simple sequence repeats predicted. This study provides a valuable reference transcriptome for understanding how H. akashiwo blooms at a molecular level.

Whole-genome single molecule real-time sequencing of SARS-CoV-2 Omicron

New variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome can only be identified using accurate sequencing methods. Single molecule real-time (SMRT) sequencing has been used to characterize Alpha and Delta variants, but not Omicron variants harboring numerous mutations in the SARS-CoV-2 genome. This study assesses the performance of a target capture SMRT sequencing protocol for whole genome sequencing (WGS) of SARS-CoV-2 Omicron variants and compared it to that of an amplicon SMRT sequencing protocol optimized for Omicron variants. The failure rate of the target capture protocol (6%) was lower than that of the amplicon protocol (34%, p < 0.001) on our data set, and the median genome coverage with the target capture protocol (98.6% [interquartile range (IQR): 86-99.4]) was greater than that with the amplicon protocol (76.6% [IQR: 66-89.6], [p < 0.001]). The percentages of samples with >95% whole genome coverage were 64% with the target capture protocol and 19% with the amplicon protocol (p < 0.05). The clades of 96 samples determined with both protocols were 93% concordant and the lineages of 59 samples were 100% concordant. Thus, target capture SMRT sequencing appears to be an efficient method for WGS, genotyping and detecting mutations of SARS-CoV-2 Omicron variants.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

86 novel HLA-E alleles discovered through full-gene sequencing of 6227 hematopoietic cell transplant patients and unrelated donors

Until recently the number of alleles of the nonclassical HLA class I gene HLA-E documented in the IPD-IMGT/HLA Database was small and as a result, the gene was often not considered to be notably polymorphic. Here, we describe our work in identifying and submitting 86 novel HLA-E alleles after full-gene single-molecule real-time (SMRT) DNA sequencing of 6227 DNA samples. These samples were comprised of 2468 patients undergoing hematopoietic cell transplantation and 3759 unrelated potential donors. A total of 111 unique HLA-E alleles were detected in this cohort. The majority of novel alleles (79.1%) contained polymorphisms in intronic regions, highlighting the significant undiscovered variation present in the noncoding regions of the HLA-E gene.

The transcriptomic response of Hyphantria cunea (Drury) to the infection of Serratia marcescens Bizio based on full-length SMRT transcriptome sequencing

Hyphantria cunea (Drury) is a globally important forest pest. We found that the Serratia marcescens Bizio strain SM1 had insecticidal activity against H. cunea, but the transcriptomic response of H. cunea to SM1 were not clear. Therefore, we performed full-length sequencing of the transcriptomes of H. cunea larvae infected with SM1 and the control group. A total of 1,183 differentially expressed genes (DEGs) were identified by comparing the group infected with SM1 and the control group, including 554 downregulated genes and 629 upregulated genes. We found many downregulated genes in metabolic pathways. Furthermore, some of these downregulated genes were involved in cellular immunity, melanization, and detoxification enzymes, which showed that SM1 weakened H. cunea immunity. In addition, genes in the juvenile hormone synthesis pathway were upregulated, which was detrimental to the survival of H. cunea. This research analyzed the transcriptomic response of H. cunea to SM1 by high-throughput full-length transcriptome sequencing. The results provide useful information to explore the relationship between S. marcescens and H. cunea, and theoretical support for the application of S. marcescens and the control of H. cunea in the future.

Identification of thyroid hormone response genes in the remodeling of dorsal muscle during Microhyla fissipes metamorphosis

INTRODUCTION

Extensive morphological, biochemical, and cellular changes occur during anuran metamorphosis, which is triggered by a single hormone, thyroid hormone (TH). The function of TH is mainly mediated through thyroid receptor (TR) by binding to the specific thyroid response elements (TREs) of direct response genes, in turn regulating the downstream genes in the cascade. The remodeling of dorsal skeletal muscle during anuran metamorphosis provides the perfect model to identify the immediate early and direct response genes that are important during apoptosis, proliferation, and differentiation of the muscle.

METHODS

In our current study, we performed Illumina sequencing combined with single-molecule real-time (SMRT) sequencing in the dorsal muscle of Microhyla fissipes after TH, cycloheximide (CHX), and TH_CHX treatment.

RESULTS AND DISCUSSION

We first identified 1,245 differentially expressed transcripts (DETs) after TH exposure, many of which were involved in DNA replication, protein processing in the endoplasmic reticulum, cell cycle, apoptosis, p53 signaling pathway, and protein digestion and absorption. In the comparison of the TH group vs. control group and TH_CHX group vs. CHX group overlapping gene, 39 upregulated and 6 downregulated genes were identified as the TH directly induced genes. Further analysis indicated that AGGTCAnnTnAGGTCA is the optimal target sequence of target genes for TR/RXR heterodimers in M. fissipes. Future investigations on the function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian dorsal muscle remodeling. These full-length and high-quality transcriptomes in this study also provide an important foundation for future studies in M. fissipes metamorphosis.

Full-Length Transcriptome of Myotis pilosus as a Reference Resource and Mining of Auditory and Immune Related Genes

Rickett's big-footed bat, Myotis pilosus, which belongs to the family Vespertilionida, is the only known piscivorous bat in East Asia. Accurate whole genome and transcriptome annotations are essential for the study of bat biological evolution. The lack of a whole genome for M. pilosus has limited our understanding of the molecular mechanisms underlying the species' evolution, echolocation, and immune response. In the present work, we sequenced the entire transcriptome using error-corrected PacBio single-molecule real-time (SMRT) data. Then, a total of 40 GB of subreads were generated, including 29,991 full-length non-chimeric (FLNC) sequences. After correction by Illumina short reads and de-redundancy, we obtained 26,717 error-corrected isoforms with an average length of 3018.91 bp and an N50 length of 3447 bp. A total of 1528 alternative splicing (AS) events were detected by transcriptome structural analysis. Furthermore, 1032 putative transcription factors (TFs) were identified, with additional identification of several long non-coding RNAs (lncRNAs) with high confidence. Moreover, several key genes, including PRL-2, DPP4, Glul, and ND1 were also identified as being associated with metabolism, immunity, nervous system processes, and auditory perception. A multitude of pattern recognition receptors was identified, including NLR, RLR, SRCR, the antiviral molecule IRF3, and the IFN receptor subunit IFNAR1. High-quality reference genomes at the transcriptome level may be used to quantify gene or transcript expression, evaluate alternative splicing levels, identify novel transcripts, and enhance genome annotation in bats.

Single-Molecule Real-Time Sequencing of Full-Length Transcriptome and Identification of Genes Related to Male Development in Cannabis sativa

Female Cannabis sativa plants have important therapeutic properties. The sex ratio of the dioecious cannabis is approximately 1:1. Cultivating homozygous female plants by inducing female plants to produce male flowers is of great practical significance. However, the mechanism underlying cannabis male development remains unclear. In this study, single-molecule real-time (SMRT) sequencing was performed using a mixed sample of female and induced male flowers from the ZYZM1 cannabis variety. A total of 15,241 consensus reads were identified, and 13,657 transcripts were annotated across seven public databases. A total of 48 lncRNAs with an average length of 986.54 bp were identified. In total, 8202 transcripts were annotated as transcription factors, the most common of which were bHLH transcription factors. Moreover, tissue-specific expression pattern analysis showed that 13 MADS transcription factors were highly expressed in male flowers. Furthermore, 232 reads of novel genes were predicted and enriched in lipid metabolism, and qRT-PCR results showed that CER1 may be involved in the development of cannabis male flowers. In addition, 1170 AS events were detected, and two AS events were further validated. Taken together, these results may improve our understanding of the complexity of full-length cannabis transcripts and provide a basis for understanding the molecular mechanism of cannabis male development.

Detection and Quantitation of DNA Damage on a Genome-wide Scale Using RADAR-seq

The formation and persistence of DNA damage can impact biological processes such as DNA replication and transcription. To maintain genome stability and integrity, organisms rely on robust DNA damage repair pathways. Techniques to detect and locate DNA damage sites across a genome enable an understanding of the consequences of DNA damage as well as how damage is repaired, which can have key diagnostic and therapeutic implications. Importantly, advancements in technology have enabled the development of high-throughput sequencing-based DNA damage detection methods. These methods require DNA enrichment or amplification steps that limit the ability to quantitate the DNA damage sites. Further, each of these methods is typically tailored to detect only a specific type of damage. RAre DAmage and Repair (RADAR) sequencing is a DNA sequencing workflow that overcomes these limitations and enables detection and quantitation of DNA damage sites in any organism on a genome-wide scale. RADAR-seq works by replacing DNA damage sites with a patch of modified bases that can be directly detected by Pacific Biosciences Single-Molecule Real Time sequencing. Here, we present three protocols that enable detection of thymine dimers and ribonucleotides in bacterial and archaeal genomes. Basic Protocol 1 enables construction of a reference genome required for RADAR-seq analyses. Basic Protocol 2 describes how to locate, quantitate, and compare thymine dimer levels in Escherichia coli exposed to varying amounts of UV light. Basic Protocol 3 describes how to locate, quantitate, and compare ribonucleotide levels in wild-type and ΔRNaseH2 Thermococcus kodakarensis. Importantly, all three protocols provide in-depth steps for data analysis. Together they serve as proof-of-principle experiments that will allow users to adapt the protocols to locate and quantitate a wide variety of DNA damage sites in any organism. © 2022 New England Biolabs. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Constructing a reference genome utilizing SMRT sequencing Basic Protocol 2: Mapping and quantitating genomic thymine dimer formation in untreated versus UV-irradiated E. coli using RADAR-seq Basic Protocol 3: Mapping and quantitating genomic ribonucleotide incorporation in wildtype versus ΔRNaseH2 T. kodakarensis using RADAR-seq.

Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.

Contrasting Responses of Arbuscular Mycorrhizal Fungal Families to Simulated Climate Warming and Drying in a Semiarid Shrubland

We carried out a 4-year manipulative field experiment in a semiarid shrubland in southeastern Spain to assess the impacts of experimental warming (W), rainfall reduction (RR), and their combination (W + RR) on the composition and diversity of arbuscular mycorrhizal fungal (AMF) communities in rhizosphere soil of H. syriacum and G. struthium shrubs using single-molecule real-time (SMRT) DNA sequencing. Across climate treatments, we encountered 109 AMF operational taxonomic units (OTUs) that were assigned to four families: Glomeraceae (93.94%), Gigasporaceae (2.19%), Claroideoglomeraceae (1.95%), and Diversisporaceae (1.92%). AMF community composition and diversity at OTU level were unaffected by the climate manipulation treatments, except for a significant decrease in AMF OTU richness in the W treatment relative to the control. However, we found a significant decrease of AMF family richness in all climate manipulation treatments relative to the control treatment. Members of the Gigasporaceae and Diversisporaceae families appeared to be highly vulnerable to intensification of heat and drought stress, as their abundances decreased by 67% and 77%, respectively, in the W + RR treatment relative to current ambient conditions. In contrast, the relative abundance and dominance of the Glomeraceae family within the AMF community increased significantly under the W + RR treatment, with Glomeraceae being the indicator family for the W + RR treatment. The interaction between warming and rainfall reduction had a significant effect on AMF community structure at family level. These findings provide new insights to help in the conservation of the soil biodiversity facing climate change in dryland ecosystems.

Whole-genome sequencing of SARS-CoV-2: Comparison of target capture and amplicon single molecule real-time sequencing protocols

Fast, accurate sequencing methods are needed to identify new variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome. Single-molecule real-time (SMRT) Pacific Biosciences (PacBio) provides long, highly accurate sequences by circular consensus reads. This study compares the performance of a target capture SMRT PacBio protocol for whole-genome sequencing (WGS) of SARS-CoV-2 to that of an amplicon PacBio SMRT sequencing protocol. The median genome coverage was higher (p < 0.05) with the target capture protocol (99.3% [interquartile range, IQR: 96.3-99.5]) than with the amplicon protocol (99.3% [IQR: 69.9-99.3]). The clades of 65 samples determined with both protocols were 100% concordant. After adjusting for C_t values, S gene coverage was higher with the target capture protocol than with the amplicon protocol. After stratification on C_t values, higher S gene coverage with the target capture protocol was observed only for samples with C_t  > 17 (p < 0.01). PacBio SMRT sequencing protocols appear to be suitable for WGS, genotyping, and detecting mutations of SARS-CoV-2.

Integration analysis of PacBio SMRT- and Illumina RNA-seq reveals P450 genes involved in thiamethoxam detoxification in Bradysia odoriphaga

The sciarid fly Bradysia odoriphaga is a serious pest of Chinese chive (Liliaceae). Neonicotinoid insecticides including thiamethoxam have been used for B. odoriphaga control. However, thiamethoxam resistance in B. odoriphaga has developed in recent years. To identify potential genes involved in detoxification metabolism of thiamethoxam in B. odoriphaga, a PacBio single-molecule real-time (SMRT) transcriptome sequencing and Illumina RNA-seq analysis on thiamethoxam treated B. odoriphaga were performed to explore differentially expressed genes in B. odoriphaga. After SMRT sequencing, analysis of Illumina RNA-Seq data showed a total of 172 differentially expressed genes (DEGs) after thiamethoxam treatment, among which eight upregulated DEGs were P450 genes that may be related to thiamethoxam metabolism. The qRT-PCR results of the eight up-regulated P450 unigenes after thiamethoxam treatment were consistent with RNA-Seq data. Furthermore, oral delivery mediated RNA interference of the eight upregulated P450 transcripts followed by insecticide bioassay was conducted, and three P450 unigenes were verified to be related to thiamethoxam detoxification in B. odoriphaga. This study provides new information about the P450 genes involved in thiamethoxam detoxification in B. odoriphaga.

Full-length transcriptome sequencing analysis and characterization, development and validation of microsatellite markers in Kengyilia melanthera

As a typical psammophyte of the Triticeae, Kengyilia melanthera possesses high feeding potential and great utilization values in desertification control in the Qinghai-Tibet Plateau. However, few gene function and genetic studies have been performed in K. melanthera. In this study, single-molecule real-time sequencing technology was used to obtain the full-length transcriptome sequence of K. melanthera, following the functional annotation of transcripts and prediction of coding sequences (CDSs), transcription factors (TFs), and long noncoding RNA (lncRNA) sequences. Meanwhile, a total of 42,433 SSR loci were detected, with 5'-UTRs having the most SSR loci and trinucleotide being the most abundant type. In total, 108,399 SSR markers were designed, and 300 SSR markers were randomly selected for diversity verification of K. melanthera. A total of 49 polymorphic SSR markers were used to construct the genetic relationships of 56 K. melanthera accessions, among which 21 SSR markers showed good cross-species transferability among the related species. In conclusion, the full-length transcriptome sequence of the K. melanthera will assist gene prediction and promote molecular biology and genomics research, and the polymorphic SSR markers will promote molecular-assisted breeding and related research of K. melanthera and its relatives.

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.

Identification of the novel HLA-E*01:03:02:25 allele in an acute lymphoblastic leukemia patient

HLA-E*01:03:02:25 sequenced from blood and buccal samples of an ALL patient.

A chromosome scale tomato genome built from complementary PacBio and Nanopore sequences alone reveals extensive linkage drag during breeding

The assembly and scaffolding of plant crop genomes facilitate the characterization of genetically diverse cultivated and wild germplasm. The cultivated tomato (Solanum lycopersicum) has been improved through the introgression of genetic material from related wild species, including resistance to pandemic strains of tobacco mosaic virus (TMV) from Solanum peruvianum. Here we applied PacBio HiFi and ONT Nanopore sequencing to develop independent, highly contiguous and complementary assemblies of an inbred TMV-resistant tomato variety. We show specific examples of how HiFi and ONT datasets can complement one another to improve assembly contiguity. We merged the HiFi and ONT assemblies to generate a long-read-only assembly where all 12 chromosomes were represented as 12 contiguous sequences (N50 = 68.5 Mbp). This chromosome scale assembly did not require scaffolding using an orthogonal data type. The merged assembly was validated by chromosome conformation capture data and is highly consistent with previous tomato genome assemblies that made use of genetic maps and Hi-C for scaffolding. Our long-read-only assembly reveals that a complex series of structural variants linked to the TMV resistance gene likely contributed to linkage drag of a 64.1-Mbp region of the S. peruvianum genome during tomato breeding. Through marker studies and ONT-based comprehensive haplotyping we show that this minimal introgression region is present in six cultivated tomato hybrid varieties developed in three commercial breeding programs. Our results suggest that complementary long read technologies can facilitate the rapid generation of near-complete genome sequences.

A high-quality reference genome for the fish pathogen Streptococcus iniae

Fish mortality caused by Streptococcus iniae is a major economic problem in aquaculture in warm and temperate regions globally. There is also risk of zoonotic infection by S. iniae through handling of contaminated fish. In this study, we present the complete genome sequence of S. iniae strain QMA0248, isolated from farmed barramundi in South Australia. The 2.12 Mb genome of S. iniae QMA0248 carries a 32 kb prophage, a 12 kb genomic island and 92 discrete insertion sequence (IS) elements. These include nine novel IS types that belong mostly to the IS3 family. Comparative and phylogenetic analysis between S. iniae QMA0248 and publicly available complete S. iniae genomes revealed discrepancies that are probably due to misassembly in the genomes of isolates ISET0901 and ISNO. Long-range PCR confirmed five rRNA loci in the PacBio assembly of QMA0248, and, unlike S. iniae 89353, no tandemly repeated rRNA loci in the consensus genome. However, we found sequence read evidence that the tandem rRNA repeat existed within a subpopulation of the original QMA0248 culture. Subsequent nanopore sequencing revealed that the tandem rRNA repeat was the most prevalent genotype, suggesting that there is selective pressure to maintain fewer rRNA copies under uncertain laboratory conditions. Our study not only highlights assembly problems in existing genomes, but provides a high-quality reference genome for S. iniae QMA0248, including manually curated mobile genetic elements, that will assist future S. iniae comparative genomic and evolutionary studies.

Genome-wide lone strand adenine methylation in Deinococcus radiodurans R1: Regulation of gene expression through DR0643-dependent adenine methylation

DNA methylation is a covalent modification of adenine or cytosine in the genome of an organism and is found in diverse microbes including the radiation resistant bacterium Deinococcus radiodurans R1. Although earlier findings have confirmed repression or de-repression of certain genes in adenine methyltransferase (DR_0643/Dam1_DR) deficient D. radiodurans mutant however, the overall regulatory aspects of Dam1_DR-mediated adenine methylation remain mostly unexplored. In the present study, we compared the genome-wide methylome and the corresponding transcriptome of D. radiodurans WT and Δdam1 mutant to explore the correlation between methylation and gene expression. In D. radiodurans, deletion of DR_0643 ORF (Δdam1) led to hypomethylation of 512 genes resulting in differential expression of 168 genes (99 genes are upregulated and 69 genes are downregulated). The modification patterns deduced for Dam1_DR (DR_0643) and Dam2_DR (DR_2267) were non-palindromic and atypical. Moreover, we observed methylation at opportunistic sites that show adenine methylation only in D. radiodurans Δdam1 and not in D. radiodurans WT. Correlation between the methylome and transcriptome suggests that hypomethylation at Dam1_DR specific sites had both negative as well as a positive effects on gene expression. Pathways such as amino acid metabolism, transport, oxidative phosphorylation, quorum sensing, signal transduction, two-component system, glycolysis/gluconeogenesis, TCA cycle, glyoxylate and dicarboxylate metabolism were modulated by Dam1_DR-mediated adenine methylation in D. radiodurans. Processes such as DNA repair, recombination, ATPase and transmembrane transporter activity were enriched when Dam1_DR mutant was subjected to radiation stress. We further evaluated the molecular interactions and mode of binding between Dam1_DR protein and S-adenosyl methionine using molecular docking followed by MD simulation. To get a better insight into the methylation mechanism, the Dam1_DR-SAM complex was also docked with a DNA molecule to elucidate DNA-Dam1_DR structural interaction during methyl-group transfer reaction. In summary, our work presents comprehensive and integrative approaches to investigate both functional and structural aspects of DNA adenine methyltransferase (Dam1_DR) in D. radiodurans biology.

Long range PCR-based deep sequencing for haplotype determination in mixed HCMV infections

Background

Short read sequencing has been used extensively to decipher the genome diversity of human cytomegalovirus (HCMV) strains, but falls short to reveal individual genomes in mixed HCMV strain populations. Novel third-generation sequencing platforms offer an extended read length and promise to resolve how distant polymorphic sites along individual genomes are linked. In the present study, we established a long amplicon PacBio sequencing workflow to identify the absolute and relative quantities of unique HCMV haplotypes spanning over multiple hypervariable sites in mixtures. Initial validation of this approach was performed with defined HCMV DNA templates derived from cell-culture enriched viruses and was further tested for its suitability on patient samples carrying mixed HCMV infections.

Results

Total substitution and indel error rate of mapped reads ranged from 0.17 to 0.43% depending on the stringency of quality trimming. Artificial HCMV DNA mixtures were correctly determined down to 1% abundance of the minor DNA source when the total HCMV DNA input was 4 × 10⁴ copies/ml. PCR products of up to 7.7 kb and a GC content < 55% were efficiently generated when DNA was directly isolated from patient samples. In a single sample, up to three distinct haplotypes were identified showing varying relative frequencies. Alignments of distinct haplotype sequences within patient samples showed uneven distribution of sequence diversity, interspersed by long identical stretches. Moreover, diversity estimation at single polymorphic regions as assessed by short amplicon sequencing may markedly underestimate the overall diversity of mixed haplotype populations.

Conclusions

Quantitative haplotype determination by long amplicon sequencing provides a novel approach for HCMV strain characterisation in mixed infected samples which can be scaled up to cover the majority of the genome by multi-amplicon panels. This will substantially improve our understanding of intra-host HCMV strain diversity and its dynamic behaviour.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08272-z.

The gene space of European mistletoe (Viscum album)

European mistletoe (Viscum album) is a hemiparasitic flowering plant that is known for its very special life cycle and extraordinary biochemical properties. Particularly, V. album has an unusual mode of cellular respiration that takes place in the absence of mitochondrial complex I. However, insights into the molecular biology of V. album so far are very limited. Since the genome of V. album is extremely large (estimated 600 times larger than the genome of the model plant Arabidopsis thaliana) it has not been sequenced up to now. We here report sequencing of the V. album gene space (defined as the space including and surrounding genic regions, encompassing coding as well as 5' and 3' non-coding regions). mRNA fractions were isolated from different V. album organs harvested in summer or winter and were analyzed via single-molecule real-time sequencing. We determined sequences of 39 092 distinct open reading frames encoding 32 064 V. album proteins (designated V. album protein space). Our data give new insights into the metabolism and molecular biology of V. album, including the biosynthesis of lectins and viscotoxins. The benefits of the V. album gene space information are demonstrated by re-evaluating mass spectrometry-based data of the V. album mitochondrial proteome, which previously had been evaluated using the A. thaliana genome sequence. Our re-examination allowed the additional identification of nearly 200 mitochondrial proteins, including four proteins related to complex I, which all have a secondary function not related to respiratory electron transport. The V. album gene space sequences are available at the NCBI.

Complete Genome Sequence of the Plant-Pathogenic Fungus Colletotrichum lupini

Colletotrichum is a fungal genus (Ascomycota, Sordariomycetes, Glomerellaceae) that includes many economically important plant pathogens that cause devastating diseases of a wide range of plants. In this work, using a combination of long- and short-read sequencing technologies, we sequenced the genome of Colletotrichum lupini RB221, isolated from white lupin (Lupinus albus) in France during a survey in 2014. The genome was assembled into 11 nuclear chromosomes and a mitochondrial genome with a total assembly size of 63.41 Mb and 36.55 kb, respectively. In total, 18,324 protein-encoding genes have been predicted, of which only 39 are specific to C. lupini. This resource will provide insight into pathogenicity factors and will help provide a better understanding of the evolution and genome structure of this important plant pathogen.[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.

Transcriptome profiling reveals gene regulation programs underlying tail development in the Ornamented Pygmy frog Microhyla fissipes

Introduction: Tadpole tail develops from the tailbud, an apparently homogenous mass of cells at the posterior of the embryo. While much progress has been made in understanding the origin and the induction of the tailbud, the subsequent outgrowth and differentiation have received much less attention, particularly with regard to global gene expression changes. Methods: By using RNA-seq with SMRT and further analyses, we report the transcriptome profiles at four key stages of tail development, from a small tailbud to the onset of feeding (S18, S19, S21 and S28) in Microhyla fissipes, an anuran with a number of advantages for developmental and genetic studies. Results: We obtained 48,826 transcripts and discovered 8807 differentially expressed transcripts (DETs, q < 0.05) among these four developmental stages. We functionally classified these DETs by using GO and KEGG analyses and revealed 110 significantly enriched GO categories and 6 highly enriched KEGG pathways (Protein digestion and absorption; ECM-receptor interaction; Pyruvate metabolism; Fatty acid degradation; Valine, leucine and isoleucine degradation; and Glyoxylate and dicarboxylate metabolism) that are likely critically involved in developmental changes in the tail. In addition, analyses of DETs between any two individual stages demonstrated the involvement of distinct biological pathways/GO terms at different stages of tail development. Furthermore, the most dramatic changes in gene expression profile are those between S28 and any of the other three stages. The upregulated DETs at S28 are highly enriched in "myosin complex" and "potassium channel activity", which are important for muscle contraction, a critical function of the tail that the animal needs by the end of embryogenesis. Additionally, many DETs and enriched pathways discovered here during tail development, such as HDAC1, Hes1 and Hippo signaling pathway, have also been reported to be vital for the tissue/organ regeneration, suggesting conserved functions between development and regeneration. Conclusion: The present staudy provides a golbal overview of gene expression patterns and new insights into the mechanism involved in anuran tail development and regeneration.

Single-molecule real-time transcript sequencing of developing cotton anthers facilitates genome annotation and fertility restoration candidate gene discovery

Heterosis refers to the superior phenotypes observed in hybrids. Cytoplasmic male sterility (CMS) system plays an important role in cotton heterosis utilization. However, the global gene expression patterns of CMS-D2 and its interaction with the restorer gene Rf1 remain unclear. Here, the full-length transcript sequencing was performed in anthers of the CMS-D2 restorer line using PacBio single-molecule real-time sequencing technology. Combining PacBio SMRT long-read isoforms and Illumina RNA-seq data, 107,066 isoforms from 44,338 loci were obtained, including 10,086 novel isoforms of novel genes and 66,419 new isoforms of known genes. Totally 56,572 alternative splicing (AS) events, 1146 lncRNAs, 61 fusion transcripts and 10,466 genes exhibited alternative polyadenylation (APA), and 60,995 novel isoforms with predicted open reading frames (ORFs) were further identified. Furthermore, the specifically expressed genes in restorer line were selected and confirmed by qRT-PCR. These findings provide a basis for upland cotton genome annotation and transcriptome research, and will help to reveal the molecular mechanism of interaction between Rf1 and CMS-D2 cytoplasm.

SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra antiquata

Coelomactra antiquata is an important aquatic economic shellfish with high medicinal value. However, because C. antiquata has no reference genome, a lot of molecular biology research cannot be carried out, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. antiquata. In the present study, we conducted the first full-length transcriptome analysis of C. antiquata by using PacBio single-molecule real-time (SMRT) sequencing technology. The results identified a total of 39,209 unigenes with an average length of 2,732 bp, 23,338 CDSs, 251 AS events, 9,881 lncRNAs, 20,106 SSRs, and 2,316 TFs. Subsequently, 59.22% (23,220) of the unigenes were successfully annotated, of which 23,164, 18,711, 15,840, 13,534, and 13,474 unigenes could be annotated using NR, Swiss-prot, KOG, GO, and KEGG databases, respectively. This study lays the foundation for the follow-up research of molecular biology and provides a reference for studying the more medicinal value of C. antiquata.

The diversity and composition of the human gut lactic acid bacteria and bifidobacterial microbiota vary depending on age

Aging is associated with gut microbiota alterations, characterized by changes in intestinal microbial diversity and composition. However, no study has yet focused on investigating age-related changes in the low-abundant but potentially beneficial subpopulations of gut lactic acid bacteria (LAB) and Bifidobacterium. Our study found that the subjects' age correlated negatively with the alpha diversity of the gut bifidobacterial microbiota, and such correlation was not observed in the gut LAB subpopulation. Principal coordinate analysis (PCoA) and analysis of distribution of operational taxonomic units (OTUs) revealed that the structure and composition of the gut bifidobacterial subpopulation of the longevous elderly group were rather different from that of the other three age groups. The same analyses were applied to identify age-dependent characteristics of the gut LAB subpopulation, and the results revealed that the gut LAB subpopulation of young adults was significantly different from that of all three elderly groups. Our study identified several potentially beneficial bacteria (e.g., Bifidobacterium breve and Bifidobacterium longum) that were enriched in the longevous elderly group (P < 0.05), and the relative abundance of Bifidobacterium adolescentis decreased significantly with the increase in age (P < 0.05). Although both bifidobacteria and LAB are generally considered as health-promoting taxa, their age-dependent distribution varied from each other, suggesting their different life stage changes and potentially different functional roles. This study provided novel species-level gut bifidobacterial and LAB microbiota profiles of a large cohort of subjects and identified several age-or longevity-associated features and biomarkers. KEY POINTS: • The alpha diversity of the gut bifidobacterial microbiota decreased with age, while LAB did not change. • The structure and composition of the gut bifidobacterial subpopulation of the longevous elderly group were rather different from that of the other three age groups. • Several potentially beneficial bacteria (e.g., Bifidobacterium breve and Bifidobacterium longum) that were enriched in the longevous elderly group.

Single-Molecule Real-Time and Illumina Sequencing to Analyze Transcriptional Regulation of Flavonoid Synthesis in Blueberry

Blueberries (Vaccinium corymbosum) contain large amounts of flavonoids, which play important roles in the plant's ability to resist stress and can also have beneficial effects on human health when the fruits are eaten. However, the molecular mechanisms that regulate flavonoid synthesis in blueberries are still unclear. In this study, we combined two different transcriptome sequencing platforms, single-molecule real-time (SMRT) and Illumina sequencing, to elucidate the flavonoid synthetic pathways in blueberries. We analyzed transcript quantity, length, and the number of annotated genes. We mined genes associated with flavonoid synthesis (such as anthocyanins, flavonols, and proanthocyanidins) and employed fluorescence quantitative PCR to analyze the expression of these genes and their correlation with flavonoid synthesis. We discovered one R2R3 MYB transcription factor from the sequencing library, VcMYB1, that can positively regulate anthocyanin synthesis in blueberries. VcMYB1 is mainly expressed in colored (mature) fruits. Experiments showed that overexpression and transient expression of VcMYB1 promoted anthocyanin synthesis in Arabidopsis, tobacco (Nicotiana benthamiana) plants and green blueberry fruits. Yeast one-hybrid (Y1H) assay, electrophoretic mobility shift assay, and transient expression experiments showed that VcMYB1 binds to the MYB binding site on the promoter of the structural gene for anthocyanin synthesis, VcMYB1 to positively regulate the transcription of VcDFR, thereby promoting anthocyanin synthesis. We also performed an in-depth investigation of transcriptional regulation of anthocyanin synthesis. This study provides background information and data for studying the synthetic pathways of flavonoids and other secondary metabolites in blueberries.

A Type I Restriction Modification System Influences Genomic Evolution Driven by Horizontal Gene Transfer in Paenibacillus polymyxa

Considered a “Generally Recognized As Safe” (GRAS) bacterium, the plant growth–promoting rhizobacterium Paenibacillus polymyxa has been widely applied in agriculture and animal husbandry. It also produces valuable compounds that are used in medicine and industry. Our previous work showed the presence of restriction modification (RM) system in P. polymyxa ATCC 842. Here, we further analyzed its genome and methylome by using SMRT sequencing, which revealed the presence of a larger number of genes, as well as a plasmid documented as a genomic region in a previous report. A number of mobile genetic elements (MGEs), including 78 insertion sequences, six genomic islands, and six prophages, were identified in the genome. A putative lysozyme-encoding gene from prophage P6 was shown to express lysin which caused cell lysis. Analysis of the methylome and genome uncovered a pair of reverse-complementary DNA methylation motifs which were widespread in the genome, as well as genes potentially encoding their cognate type I restriction-modification system PpoAI. Further genetic analysis confirmed the function of PpoAI as a RM system in modifying and restricting DNA. The average frequency of the DNA methylation motifs in MGEs was lower than that in the genome, implicating a role of PpoAI in restricting MGEs during genomic evolution of P. polymyxa. Finally, comparative analysis of R, M, and S subunits of PpoAI showed that homologs of the PpoAI system were widely distributed in species belonging to other classes of Firmicute, implicating a role of the ancestor of PpoAI in the genomic evolution of species beyond Paenibacillus.

Transcriptome analysis identifies putative genes involved in triterpenoid biosynthesis in Platycodon grandiflorus

Comprehensive transcriptome analysis of different Platycodon grandiflorus tissues discovered genes related to triterpenoid saponin biosynthesis. Platycodon grandiflorus (Jacq.) A. DC. (P. grandiflorus), a traditional Chinese medicine, contains considerable triterpenoid saponins with broad pharmacological activities. Triterpenoid saponins are the major components of P. grandiflorus. Here, single-molecule real-time and next-generation sequencing technologies were combined to comprehensively analyse the transcriptome and identify genes involved in triterpenoid saponin biosynthesis in P. grandiflorus. We quantified four saponins in P. grandiflorus and found that their total content was highest in the roots and lowest in the stems and leaves. A total of 173,354 non-redundant transcripts were generated from the PacBio platform, and three full-length transcripts of β-amyrin synthase, the key synthase of β-amyrin, were identified. A total of 132,610 clean reads obtained from the DNBSEQ platform were utilised to explore key genes related to the triterpenoid saponin biosynthetic pathway in P. grandiflorus, and 96 differentially expressed genes were selected as candidates. The expression levels of these genes were verified by quantitative real-time PCR. Our reliable transcriptome data provide valuable information on the related biosynthesis pathway and may provide insights into the molecular mechanisms of triterpenoid saponin biosynthesis in P. grandiflorus.

Single-Molecule Real-Time and Illumina-Based RNA Sequencing Data Identified Vernalization-Responsive Candidate Genes in Faba Bean (Vicia faba L.)

Faba bean (Vicia faba L.) is one of the most widely grown cool season legume crops in the world. Winter faba bean normally has a vernalization requirement, which promotes an earlier flowering and pod setting than unvernalized plants. However, the molecular mechanisms of vernalization in faba bean are largely unknown. Discovering vernalization-related candidate genes is of great importance for faba bean breeding. In this study, the whole transcriptome of faba bean buds was profiled by using next-generation sequencing (NGS) and single-molecule, real-time (SMRT) full-length transcriptome sequencing technology. A total of 29,203 high-quality non-redundant transcripts, 21,098 complete coding sequences (CDS), 1,045 long non-coding RNAs (lncRNAs), and 12,939 simple sequence repeats (SSRs) were identified. Furthermore, 4,044 differentially expressed genes (DEGs) were identified through pairwise comparisons. By Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, these differentially expressed transcripts were found to be enriched in binding and transcription factor activity, electron carrier activity, rhythmic process, and receptor activity. Finally, 50 putative vernalization-related genes that played important roles in the vernalization of faba bean were identified; we also found that the levels of vernalization-responsive transcripts showed significantly higher expression levels in cold-treated buds. The expression of VfSOC1, one of the candidate genes, was sensitive to vernalization. Ectopic expression of VfSOC1 in Arabidopsis brought earlier flowering. In conclusion, the abundant vernalization-related transcripts identified in this study will provide a basis for future researches on the vernalization and faba bean breeding and established a reference full-length transcriptome for future studies on faba bean.

Growth trait gene analysis of kuruma shrimp (Marsupenaeus japonicus) by transcriptome study

Growth traits are a vital standard for the animal culture industry. The molecular mechanism of growth traits remains poorly understood, especially in aquaculture, which hinders the development of the selective breeding industry. Genomic resources discovered by next-generation sequencing (NGS) have been widely applied in certain species. However, accurate assembly and downstream analysis by NGS are still major challenges for species without reference genomes. In this study, a comparative transcriptome analysis of an economic crustacean species (Marsupenaeus japonicus) between a fast growth group and slow growth group at different stages was performed by SMRT (single molecule real time) and NGS. A high-quality full-length transcriptome (e.g., mean length of unigenes was longer than those unigenes assembled by Illumina clean reads from previous reports, and annotation rate was higher than Illumina sequencing in the same studies) was generated and analyzed. Several differentially expressed genes (DEGs) related to growth were identified and validated by quantitative real-time PCR (qPCR). The results showed that compared with the late stage, more DEGs were identified at the early stage, indicating that the growth-related physiological activity differences between different individuals at the early stage were higher than at the late stage. Moreover, 215 DEGs were shared between the early stage and late stage, and 109 had divergent functions during development. These 109 genes may play an important role in regulating the specific growth rate (SGR) of kuruma shrimp. In addition, twelve growth-related pathways were shared between the two comparative groups. Among these pathways, the fly Hippo signaling pathway and its key gene Mj14-3-3-like were identified for the first time to be involved in growth traits in crustaceans. Further analysis showed that Mj14-3-3-like was significantly downregulated in the fast growth group at the early stage and late stage; its expression level was reduced to its lowest level at the intermolt stage (C), the most important growth stage in shrimp, suggesting that Mj14-3-3-like may inhibit the growth of kuruma shrimp. Our study helps to elucidate the genes involved in the molecular mechanisms governing growth traits in kuruma shrimp, which is valuable for future shrimp developmental research.

Full-Length Transcriptome Analysis Reveals Candidate Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is an important medicinal plant widely utilized for moxibustion heat therapy in China. The terpenoid biosynthesis process in A. argyi is speculated to play a key role in conferring its medicinal value. However, the molecular mechanism underlying terpenoid biosynthesis remains unclear, in part because the reference genome of A. argyi is unavailable. Moreover, the full-length transcriptome of A. argyi has not yet been sequenced. Therefore, in this study, de novo transcriptome sequencing of A. argyi's root, stem, and leaf tissues was performed to obtain those candidate genes related to terpenoid biosynthesis, by combining the PacBio single-molecule real-time (SMRT) and Illumina sequencing NGS platforms. And more than 55.4 Gb of sequencing data and 108,846 full-length reads (non-chimeric) were generated by the Illumina and PacBio platform, respectively. Then, 53,043 consensus isoforms were clustered and used to represent 36,820 non-redundant transcripts, of which 34,839 (94.62%) were annotated in public databases. In the comparison sets of leaves vs roots, and leaves vs stems, 13,850 (7,566 up-regulated, 6,284 down-regulated) and 9,502 (5,284 up-regulated, 4,218 down-regulated) differentially expressed transcripts (DETs) were obtained, respectively. Specifically, the expression profile and KEGG functional enrichment analysis of these DETs indicated that they were significantly enriched in the biosynthesis of amino acids, carotenoids, diterpenoids and flavonoids, as well as the metabolism processes of glycine, serine and threonine. Moreover, multiple genes encoding significant enzymes or transcription factors related to diterpenoid biosynthesis were highly expressed in the A. argyi leaves. Additionally, several transcription factor families, such as RLK-Pelle_LRR-L-1 and RLK-Pelle_DLSV, were also identified. In conclusion, this study offers a valuable resource for transcriptome information, and provides a functional genomic foundation for further research on molecular mechanisms underlying the medicinal use of A. argyi leaves.

SMRT and Illumina RNA Sequencing and Characterization of a Key NAC Gene LoNAC29 during the Flower Senescence in Lilium oriental 'Siberia'

Lily (Lilium spp.) is an important cut flower around the world. Flower senescence in lilies is characterized by the wilting and abscission of tepals, which results in a decrease in flower quality and huge economic loss. However, the mechanism underlying flower senescence in lilies is largely unknown. In this study, single-molecule, real-time (SMRT) and Illumina sequencing were carried out in L. oriental ‘Siberia’. Sequencing yielded 73,218 non-redundant transcripts, with an N50 of 3792 bp. These data were further integrated with three published transcriptomes through cogent analysis, which yielded 62,960 transcripts, with an increase in N50 of 3935 bp. Analysis of differentially expressed genes showed that 319 transcription factors were highly upregulated during flower senescence. The expression of twelve NAC genes and eleven senescence-associated genes (SAGs) showed that LoNAC29 and LoSAG39 were highly expressed in senescent flowers. Transient overexpression of LoNAC29 and LoSAG39 in tepals of lily notably accelerated flower senescence, and the promoter activity of LoSAG39 was strongly induced by LoNAC29. This work supported new evidence for the molecular mechanism of flower senescence and provided better sequence data for further study in lilies.

Molecular classification of a complex structural rearrangement of the RB1 locus in an infant with sporadic, isolated, intracranial, sellar region retinoblastoma

Retinoblastoma is a childhood cancer of the retina involving germline or somatic alterations of the RB Transcriptional Corepressor 1 gene, RB1. Rare cases of sellar-suprasellar region retinoblastoma without evidence of ocular or pineal tumors have been described. A nine-month-old male presented with a sellar-suprasellar region mass. Histopathology showed an embryonal tumor with focal Flexner-Wintersteiner-like rosettes and loss of retinoblastoma protein (RB1) expression by immunohistochemistry. DNA array-based methylation profiling confidently classified the tumor as pineoblastoma group A/intracranial retinoblastoma. The patient was subsequently enrolled on an institutional translational cancer research protocol and underwent comprehensive molecular profiling, including paired tumor/normal exome and genome sequencing and RNA-sequencing of the tumor. Additionally, Pacific Biosciences (PacBio) Single Molecule Real Time (SMRT) sequencing was performed from comparator normal and disease-involved tissue to resolve complex structural variations. RNA-sequencing revealed multiple fusions clustered within 13q14.1-q21.3, including a novel in-frame fusion of RB1-SIAH3 predicted to prematurely truncate the RB1 protein. SMRT sequencing revealed a complex structural rearrangement spanning 13q14.11-q31.3, including two somatic structural variants within intron 17 of RB1. These events corresponded to the RB1-SIAH3 fusion and a novel RB1 rearrangement expected to correlate with the complete absence of RB1 protein expression. Comprehensive molecular analysis, including DNA array-based methylation profiling and sequencing-based methodologies, were critical for classification and understanding the complex mechanism of RB1 inactivation in this diagnostically challenging tumor.

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.

Bifidobacterium Genome Assembly and Methylome Analysis Using Pacbio SMRT Sequencing

In this chapter, we present a generic method to achieve de novo genome assembly and methylome analysis of Bifidobacterium genomes using the Pacbio SMRT sequencing and SMRT Link pipeline. The methods described here cover the de novo Pacbio or hybrid Illumina-Pacbio assembly ideally intended to obtain a complete genome sequence, followed by the detection of base modifications and methylome analysis.The identified DNA motifs obtained by methylome analysis can be used to predict active restriction-modification (RM) systems in Bifidobacterium strains. The presence of active RM systems and knowledge on their target motifs may guide selection of suitable cloning or mutagenesis vectors for bifidobacteria.