LoRDEC: accurate and efficient long read error correction

Transcriptome Level Reveals the Triterpenoid Saponin Biosynthesis Pathway of Bupleurum falcatum L

Bupleurum falcatum L. is frequently used in traditional herbal medicine in Asia. Saikosaponins (SSs) are the main bioactive ingredients of B. falcatum, but the biosynthetic pathway of SSs is unclear, and the biosynthesis of species-specific phytometabolites is little known. Here we resolved the transcriptome profiles of B. falcatum to identify candidate genes that might be involved in the biosynthesis of SSs. By isoform sequencing (Iso-Seq) analyses of the whole plant, a total of 26.98 Gb of nucleotides were obtained and 124,188 unigenes were identified, and 81,594 unigenes were successfully annotated. A total of 1033 unigenes of 20 families related to the mevalonate (MVA) pathway and methylerythritol phosphate (MEP) pathway of the SS biosynthetic pathway were identified. The WGCNA (weighted gene co-expression network analysis) of these unigenes revealed that only the co-expression module of MEmagenta, which contained 343 unigenes, was highly correlated with the biosynthesis of SSs. Comparing differentially expressed gene analysis and the WGCNA indicated that 130 out of 343 genes of the MEmagenta module exhibited differential expression levels, and genes with the most "hubness" within this module were predicted. Manipulation of these genes might improve the biosynthesis of SSs.

Conjugative IncHI2/HI2A plasmids harbouring mcr-9 in colistin-susceptible Escherichia coli isolated from diseased pigs in Japan

Colistin is a last resort antimicrobial used for the treatment of gram-negative bacterial infections. Plasmid-mediated colistin resistance (mcr) genes are a cause of global concern, and, thus far, mcr-1-10 have been identified. In a previous study, we screened mcr-1-5 in Escherichia coli derived from diseased pigs in Japan and reported a high prevalence of mcr-1, -3 and -5. However, the previous report on the prevalence of mcr genes was inaccurate. In the present study, we aimed to clarify the prevalence of all reported variants of mcr in E. coli derived from the diseased pigs, which were previously screened for mcr-1-5. Additionally, we also characterized the mcr-9-positive E. coli , which was detected in this study. We screened mcr in 120 E. coli strains from diseased pigs and mcr-positive E. coli and an mcr-carrying plasmid were also characterized. One mcr-9-positive colistin-susceptible E. coli strain was detected (0.8 %). Plasmid-mediated mcr-9 was transferred to E. coli ML4909 as the recipient strain, and it was located on IncHI2/HI2A plasmid p387_L with other antimicrobial resistance genes (ARGs). The region harbouring ARGs including mcr-9, was similar to that on the Klebsiella pneumoniae chromosome harbouring mcr-9 isolated in Japan. mcr-3, -5 and -9 were detected (4.2 %) in colistin-susceptible strains. mcr-9 was found to be disseminated via the plasmid IncHI2/HI2A p387_L and transferred and inserted into chromosomes via a transposon. Our results suggest that mcr genes should be monitored regularly, regardless of their susceptibility to colistin.

The combination of SMRT sequencing and Illumina sequencing highlights organ-specific and age-specific expression patterns of miRNAs in Sika Deer

Due to the lack of high-quality Sika Deer (Cervus nippon) transcriptome and sRNAome across multiple organs or development stages, it is impossible to comprehensively analyze the mRNA and miRNA regulatory networks related to growth, development and immunity response. In this study, we used single molecule-real time sequencing (SMRT-seq) and Illumina sequencing methods to generate transcriptome and sRNAome from ten tissues and four age groups of Sika Deer to help us understand molecular characteristics and global miRNA expression profiles. The results showed that a total of 240,846 consensus transcripts were generated with an average length of 2,784 bp. 4,329 Transcription factors (TFs), 109,000 Simple Sequence Repeats (SSRs) and 18,987 Long non-coding RNAs (LncRNAs) were identified. Meanwhile, 306 known miRNAs and 143 novel miRNAs were obtained. A large number of miRNAs showed organ-specific and age-specific differential expression patterns. In particular, we found that the organ-specific miRNAs were enriched in the brain, some of which shared only between the brain and adrenal. These miRNAs were involved in maintaining specific functions within the brain and adrenal. By constructing miRNA96mRNA interaction networks associated with Sika Deer immunity, we found that miRNAs (miR-148a, miR-26a, miR-214, let-7b, etc.) and mRNAs (CD6, TRIM38, C3, CD163, etc.) might play an important role in the immune response of Sika Deer spleen. Together, our study generated an improved transcript annotation for Sika Deer by SMRT-seq and revealed the role of miRNA in regulating the growth, development and immunity response of Sika Deer.

Transcriptome Profiling of Stem-Differentiating Xylem in Response to Abiotic Stresses Based on Hybrid Sequencing in Cunninghamia lanceolata

Cunninghamia lanceolata (C. lanceolata) belongs to Gymnospermae, which are fast-growing and have desirable wood properties. However, C. lanceolata's stress resistance is little understood. To unravel the physiological and molecular regulation mechanisms under environmental stresses in the typical gymnosperm species of C. lanceolata, three-year-old plants were exposed to simulated drought stress (polyethylene glycol 8000), salicylic acid, and cold treatment at 4 °C for 8 h, 32 h, and 56 h, respectively. Regarding the physiological traits, we observed a decreased protein content and increased peroxidase upon salicylic acid and polyethylene glycol treatment. Superoxide dismutase activity either decreased or increased at first and then returned to normal under the stresses. Regarding the molecular regulation, we used both nanopore direct RNA sequencing and short-read sequencing to reveal a total of 5646 differentially expressed genes in response to different stresses, of which most had functions in lignin catabolism, pectin catabolism, and xylan metabolism, indicating that the development of stem-differentiating xylem was affected upon stress treatment. Finally, we identified a total of 51 AP2/ERF, 29 NAC, and 37 WRKY transcript factors in C. lanceolata. The expression of most of the NAC TFs increased under cold stress, and the expression of most of the WRKY TFs increased under cold and SA stress. These results revealed the transcriptomics responses in C. lanceolata to short-term stresses under this study's experimental conditions and provide preliminary clues about stem-differentiating xylem changes associated with different stresses.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

VeChat: correcting errors in long reads using variation graphs

Error correction is the canonical first step in long-read sequencing data analysis. Current self-correction methods, however, are affected by consensus sequence induced biases that mask true variants in haplotypes of lower frequency showing in mixed samples. Unlike consensus sequence templates, graph-based reference systems are not affected by such biases, so do not mistakenly mask true variants as errors. We present VeChat, as an approach to implement this idea: VeChat is based on variation graphs, as a popular type of data structure for pangenome reference systems. Extensive benchmarking experiments demonstrate that long reads corrected by VeChat contain 4 to 15 (Pacific Biosciences) and 1 to 10 times (Oxford Nanopore Technologies) less errors than when being corrected by state of the art approaches. Further, using VeChat prior to long-read assembly significantly improves the haplotype awareness of the assemblies. VeChat is an easy-to-use open-source tool and publicly available at https://github.com/HaploKit/vechat .

The gastrodin biosynthetic pathway in Pholidota chinensis Lindl. revealed by transcriptome and metabolome profiling

Pholidota chinensis Lindl. is an epiphytic or lithophytic perennial herb of Orchidaceae family used as a garden flower or medicinal plant to treat high blood pressure, dizziness and headache in traditional Chinese medicine. Gastrodin (GAS) is considered as a main bioactive ingredient of this herb but the biosynthetic pathway remains unclear in P. chinensis. To elucidate the GAS biosynthesis and identify the related genes in P. chinensis, a comprehensive analysis of transcriptome and metabolome of roots, rhizomes, pseudobulbs and leaves were performed by using PacBio SMART, Illumina Hiseq and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS/MS). A total of 1,156 metabolites were identified by UPLC-MS/MS, of which 345 differential metabolites were mainly enriched in phenylpropanoid/phenylalanine, flavone and flavonol biosynthesis. The pseudobulbs make up nearly half of the fresh weight of the whole plant, and the GAS content in the pseudobulbs was also the highest in four tissues. Up to 23,105 Unigenes were obtained and 22,029 transcripts were annotated in the transcriptome analysis. Compared to roots, 7,787, 8,376 and 9,146 differentially expressed genes (DEGs) were identified in rhizomes, pseudobulbs and leaves, respectively. And in total, 80 Unigenes encoding eight key enzymes for GAS biosynthesis, were identified. Particularly, glycosyltransferase, the key enzyme of the last step in the GAS biosynthetic pathway had 39 Unigenes candidates, of which, transcript28360/f2p0/1592, was putatively identified as the most likely candidate based on analysis of co-expression, phylogenetic analysis, and homologous searching. The metabolomics and transcriptomics of pseudobulbs versus roots showed that 8,376 DEGs and 345 DEMs had a substantial association based on the Pearson's correlation. This study notably enriched the metabolomic and transcriptomic data of P. chinensis, and it provides valuable information for GAS biosynthesis in the plant.

Hybrid assembly using long reads resolves repeats and completes the genome sequence of a laboratory strain of Staphylococcus aureus subsp. Aureus RN4220

Staphylococcus aureus RN4220 has been extensively used by staphylococcal researchers as an intermediate strain for genetic manipulation due to its ability to accept foreign DNA. Despite its wide use in laboratories, its complete genome is not available. In this study, we used a hybrid genome assembly approach using minION long reads and Illumina short reads to sequence the complete genome of S. aureus RN4220. The comparative analysis of the annotated complete genome showed the presence of 39 genes fragmented in the previous assembly, many of which were located near the repeat regions. Using RNA-Seq reads, we showed that a higher number of reads could be mapped to the complete genome than the draft genome and the gene expression profile obtained using the complete genome also differs from that obtained from the draft genome. Furthermore, by comparative transcriptomic analysis, we showed the correlation between expression levels of staphyloxanthin biosynthetic genes and the production of yellow pigment. This study highlighted the importance of long reads in completing microbial genomes, especially those possessing repetitive elements.

•

S. aureus RN4220 is used as an intermediate strain for genetic manipulation.

•

We completed its genome and found 39 fragmented genes in previous genome assembly.

•

RNA-Seq analysis improved mapping of the reads with the use of complete genome.

•

Expression of staphyloxanthin biosynthetic genes was correlated with its production.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

RNA-Seq and Iso-Seq Reveal the Important Role of COMT and CCoAOMT Genes in Accumulation of Scopoletin in Noni (Morinda citrifolia)

Scopoletin, the main component of clinical drugs and the functional component of health products, is highly abundant in noni fruit (Morinda citrifolia). Multiple enzyme genes regulate scopoletin accumulation. In the present study, differentially expressed genes of noni were analyzed by RNA sequencing (RNA-Seq) and the full-length genes by isoform-sequencing (Iso-Seq) to find the critical genes in the scopoletin accumulation mechanism pathway. A total of 32,682 full-length nonchimeric reads (FLNC) were obtained, out of which 16,620 non-redundant transcripts were validated. Based on KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation and differential expression analysis, two differentially expressed genes, caffeic acid 3-O-methyltransferase (COMT) and caffeoyl-CoA O-methyltransferase (CCoAOMT), were found in the scopoletin accumulation pathway of noni. Real-time quantitative polymerase chain reaction (q-PCR), phylogenetic tree analysis, gene expression analysis, and the change in scopoletin content confirmed that these two proteins are important in this pathway. Based on these results, the current study supposed that COMT and CCoAOMT play a significant role in the accumulation of scopoletin in noni fruit, and COMT (gene number: gene 7446, gene 8422, and gene 6794) and CCoAOMT (gene number: gene 12,084) were more significant. These results provide the importance of COMT and CCoAOMT and a basis for further understanding the accumulation mechanism of scopoletin in noni.

HaploDMF: viral haplotype reconstruction from long reads via deep matrix factorization

MOTIVATION

Lacking strict proofreading mechanisms, many RNA viruses can generate progeny with slightly changed genomes. Being able to characterize highly similar genomes (i.e. haplotypes) in one virus population helps study the viruses' evolution and their interactions with the host/other microbes. High-throughput sequencing data has become the major source for characterizing viral populations. However, the inherent limitation on read length by next-generation sequencing makes complete haplotype reconstruction difficult.

RESULTS

In this work, we present a new tool named HaploDMF that can construct complete haplotypes using third-generation sequencing (TGS) data. HaploDMF utilizes a deep matrix factorization model with an adapted loss function to learn latent features from aligned reads automatically. The latent features are then used to cluster reads of the same haplotype. Unlike existing tools whose performance can be affected by the overlap size between reads, HaploDMF is able to achieve highly robust performance on data with different coverage, haplotype number and error rates. In particular, it can generate more complete haplotypes even when the sequencing coverage drops in the middle. We benchmark HaploDMF against the state-of-the-art tools on simulated and real sequencing TGS data on different viruses. The results show that HaploDMF competes favorably against all others.

AVAILABILITY AND IMPLEMENTATION

The source code and the documentation of HaploDMF are available at https://github.com/dhcai21/HaploDMF.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

Chromosomal coharboring of blaIMP-60 and mcr-9 in Enterobacter asburiae isolated from a Japanese woman with empyema: a case report

BACKGROUND

Polymyxin E (colistin) is a last-resort antibiotic to treat infections caused by carbapenemase-producing Enterobacteriaceae (CPE). However, reports of CPEs resistant to colistin have been increasing, and the mcr genes are emerging as resistance mechanisms. Among them, plasmid-mediate mcr-9 is known to be associated with colistin resistance, whereas reports on chromosomal mcr-9 and its association with colistin resistance in humans are few.

CASE PRESENTATION

We identified Enterobacter asburiae harboring mcr-9 and bla_IMP-60 in the pleural fluid of a patient with empyema. The long-read sequencing technique revealed that these genes were located on its chromosome. Despite the lack of exposure to colistin, the organism showed microcolonies in the inhibition circle in the E-test and disk diffusion test. Antibiotic susceptibility testing by broth microdilution confirmed its resistance to colistin.

CONCLUSION

Our case report showed that mcr-9 can be present not only on plasmids but also on the chromosome in E. asburiae, and that the presence of mcr-9 on its chromosome may influence its susceptibility to colistin.

Fec: a fast error correction method based on two-rounds overlapping and caching

The third-generation sequencing technology has advanced genome analysis with long-read length, but the reads need error correction due to the high error rate. Error correction is a time-consuming process especially when the sequencing coverage is high. Generally, for a pair of overlapping reads A and B, the existing error correction methods perform a base-level alignment from B to A when correcting the read A. And another base-level alignment from A to B is performed when correcting the read B. However, based on our observation, the base-level alignment information can be reused. In this article, we present a fast error correction tool Fec, using two-rounds overlapping and caching. Fec can be used independently or as an error correction step in an assembly pipeline. In the first round, Fec uses a large window size (20) to quickly find enough overlaps to correct most of the reads. In the second round, a small window size (5) is used to find more overlaps for the reads with insufficient overlaps in the first round. When performing base-level alignment, Fec searches the cache first. If the alignment exists in the cache, Fec takes this alignment out and deduces the second alignment from it. Otherwise, Fec performs base-level alignment and stores the alignment in the cache. We test Fec on nine datasets, and the results show that Fec has 1.24-38.56 times speed-up compared to MECAT, CANU and MINICNS on five PacBio datasets and 1.16-27.8 times speed-up compared to NECAT and CANU on four nanopore datasets.

AVAILABILITY AND IMPLEMENTATION

Fec is available at https://github.com/zhangjuncsu/Fec.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Camellia nitidissima Petals

Flavonoids are the main medicinal ingredients in Camellia nitidissima, but the regulatory mechanism of flavonoid biosynthesis in flowers is unclear; therefore, the flavonoids in C. nitidissima have not been effectively used. The present study performed full-length transcriptome sequencing of C. nitidissima flower. Furthermore, the reported RNA-sequencing data of C. nitidissima petals were reanalyzed using the full-length transcriptome as a reference, and the regulatory mechanism of flavonoid synthesis in petals was elucidated. The analysis identified 43,350 isoforms annotated in non-redundant protein (Nr), Kyoto Encyclopedia of Genes and Genomes (KEGG), EuKaryotic Orthologous Groups (KOG), and Swiss-Prot databases, among which 34,602 aligned to Camellia sinensis genes. A total of 11,857 differentially expressed genes (DEGs), including 112 related to flavonoid synthesis, were identified by pairwise comparison. Subsequently, analysis of the phylogeny and the conserved motifs of R2R3-MYB using the proteins sequences identified three R2R3-MYB transcription factors that regulated flavonoid biosynthesis. Weighted gene co-expression network analysis (WGCNA) identified phenylalanine ammonia-lyase (PAL) and 4-coumarate: CoA ligase(4CL) as the hub genes and showed that bHLH79 interacted with PAL. Finally, validated the expression of seven DEGs involved in flavonoid biosynthesis using real-time quantitative PCR (qRT-PCR). Thus, the present study generated and used the full-length transcriptome as the reference to analyze the transcriptome of petals and proposed a possible regulatory mechanism of flavonoid synthesis in C. nitidissima. The study's findings unravel the genetic mechanisms underlying flavonoid synthesis and suggest candidate genes for the genetic improvement of C. nitidissima.

Arun P, Cui G, Bay LK, van Oppen MJH, Webster NS, Aranda M. The coral Acropora loripes genome reveals an alternative pathway for cysteine biosynthesis in animals

The metabolic capabilities of animals have been derived from well-studied model organisms and are generally considered to be well understood. In animals, cysteine is an important amino acid thought to be exclusively synthesized through the transsulfuration pathway. Corals of the genus Acropora have lost cystathionine β-synthase, a key enzyme of the transsulfuration pathway, and it was proposed that Acropora relies on the symbiosis with dinoflagellates of the family Symbiodiniaceae for the acquisition of cysteine. Here, we identify the existence of an alternative pathway for cysteine biosynthesis in animals through the analysis of the genome of the coral Acropora loripes. We demonstrate that these coral proteins are functional and synthesize cysteine in vivo, exhibiting previously unrecognized metabolic capabilities of animals. This pathway is also present in most animals but absent in mammals, arthropods, and nematodes, precisely the groups where most of the animal model organisms belong to, highlighting the risks of generalizing findings from model organisms.

Full-Length Transcriptome Maps of Reef-Building Coral Illuminate the Molecular Basis of Calcification, Symbiosis, and Circa-Dian Genes

Coral transcriptomic data largely rely on short-read sequencing, which severely limits the understanding of coral molecular mechanisms and leaves many important biological questions unresolved. Here, we sequence the full-length transcriptomes of four common and frequently dominant reef-building corals using the PacBio Sequel II platform. We obtain information on reported gene functions, structures, and expression profiles. Among them, a comparative analysis of biomineralization-related genes provides insights into the molecular basis of coral skeletal density. The gene expression profiles of the symbiont Symbiodiniaceae are also isolated and annotated from the holobiont sequence data. Finally, a phylogenetic analysis of key circadian clock genes among 40 evolutionarily representative species indicates that there are four key members in early metazoans, including cry genes; Clock or Npas2; cyc or Arntl; and tim, while per, as the fifth member, occurs in Bilateria. In summary, this work provides a foundation for further work on the manipulation of skeleton production or symbiosis to promote the survival of these important organisms.

Identification of candidate sex-specific genomic regions in male and female Asian arowana genomes

BACKGROUND

Asian arowana, Scleropages formosus, is one of the most expensive aquarium fish species worldwide. Its sex, however, cannot be distinguished clearly at any development stage, which impedes captive breeding and species protection for this endangered aquarium fish.

RESULTS

To discover molecular clues to the sex of Asian arowana, we sequenced 26.5 Gb of PacBio HiFi reads and 179.2 Gb of Hi-C reads for 1 male fish and also sequenced 106.5 Gb of Illumina reads, 36.0 Gb of PacBio Sequel reads, and 80.7 Gb of Hi-C reads for 1 female individual. The final male and female genome assemblies were approximately 756.8 Mb and 781.5 Mb in length and contained 25,262 and 25,328 protein-coding genes, respectively. We also resequenced the genomes of 15 male and 15 female individuals with approximately 722.1 Gb of Illumina reads. A genome-wide association study identified several potentially divergent regions between male and female individuals. In these regions, cd48 and cfap52 could be candidate genes for sex determination of Asian arowana. We also found some structural variations in few chromosomes between male and female individuals.

CONCLUSION

We provided an improved reference genome assembly of female arowana and generated the first sequenced genome of 1 male individual. These valuable genetic resources and resequencing data may improve global aquarium fish research.

Scalable, ultra-fast, and low-memory construction of compacted de Bruijn graphs with Cuttlefish 2

The de Bruijn graph is a key data structure in modern computational genomics, and construction of its compacted variant resides upstream of many genomic analyses. As the quantity of genomic data grows rapidly, this often forms a computational bottleneck. We present Cuttlefish 2, significantly advancing the state-of-the-art for this problem. On a commodity server, it reduces the graph construction time for 661K bacterial genomes, of size 2.58Tbp, from 4.5 days to 17-23 h; and it constructs the graph for 1.52Tbp white spruce reads in approximately 10 h, while the closest competitor requires 54-58 h, using considerably more memory.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Transcriptome innovations in primates revealed by single-molecule long-read sequencing

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Characterization of the gene expression profile response to drought stress in Haloxylon using PacBio single-molecule real-time and Illumina sequencing

Haloxylon ammodendron and Haloxylon persicum are important drought-tolerant plants in northwest China. The whole-genome sequencing of H. ammodendron and H. persicum grown in their natural environment is incomplete, and their transcriptional regulatory network in response to drought environment remains unclear. To reveal the transcriptional responses of H. ammodendron and H. persicum to an arid environment, we performed single-molecule real-time (SMRT) and Illumina RNA sequencing. In total, 20,246,576 and 908,053 subreads and 435,938 and 210,334 circular consensus sequencing (CCS) reads were identified by SMRT sequencing of H. ammodendron and H. persicum, and 15,238 and 10,135 unigenes, respectively, were successfully obtained. In addition, 9,794 and 7,330 simple sequence repeats (SSRs) and 838 and 71 long non-coding RNAs were identified. In an arid environment, the growth of H. ammodendron was restricted; plant height decreased significantly; basal and branch diameters became thinner and hydrogen peroxide (H₂O₂) content and peroxidase (POD) activity were increased. Under dry and wet conditions, 11,803 and 15,217 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum, respectively. There were 319 and 415 DEGs in the signal transduction pathways related to drought stress signal perception and transmission, including the Ca²⁺ signal pathway, the ABA signal pathway, and the MAPK signal cascade. In addition, 217 transcription factors (TFs) and 398 TFs of H. ammodendron and H. persicum were differentially expressed, including FAR1, MYB, and AP2/ERF. Bioinformatic analysis showed that under drought stress, the expression patterns of genes related to active oxygen [reactive oxygen species (ROS)] scavenging, functional proteins, lignin biosynthesis, and glucose metabolism pathways were altered. Thisis the first full-length transcriptome report concerning the responses of H. ammodendron and H. persicum to drought stress. The results provide a foundation for further study of the adaptation to drought stress. The full-length transcriptome can be used in genetic engineering research.

Long-Read Sequencing Annotation of the Transcriptome in DNA-PK Inactivated Cells

The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) with a Ku70/Ku80 heterodimer constitutes the intact DNA-PK kinase, which is an upstream component of the DNA repair machinery that signals the DNA damage, orchestrates the DNA repair, and serves to maintain genome integrity. Beyond its role in DNA damage repair, the DNA-PK kinase is also implicated in transcriptional regulation and RNA metabolism, with an illuminated impact on tumor progression and therapeutic responses. However, the efforts to identify DNA-PK regulated transcriptomes are limited by short-read sequencing to resolve the full complexity of the transcriptome. Therefore, we leveraged the PacBio Single Molecule, Real-Time (SMRT) Sequencing platform to study the transcriptome after DNA-PK inactivation to further underscore the importance of its role in diseases. Our analysis revealed additional novel transcriptome and complex gene structures in the DNA-PK inactivated cells, identifying 8,355 high-confidence new isoforms from 3,197 annotated genes and 523 novel genes. Among them, 380 lncRNAs were identified. We validated these findings using computational approaches and confirmatory transcript quantification with short-read sequencing. Several novel isoforms representing distinct splicing events have been validated through PCR experiments. Our analyses provide novel insights into DNA-PK function in transcriptome regulation and RNA metabolism.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

A chromosome-level genome assembly of the jade perch (Scortum barcoo)

Endemic to Australia, jade perch (Scortum barcoo) is a highly profitable freshwater bass species. It has extraordinarily high levels of omega-3 polyunsaturated fatty acids (PUFAs), which detailed genes involved in are largely unclear. Meanwhile, there were four chromosome-level bass species have been previous sequenced, while the bass ancestor genome karyotypes have not been estimated. Therefore, we sequenced, assembled and annotated a genome of jade perch to characterize the detailed genes for biosynthesis of omega-3 PUFAs and to deduce the bass ancestor genome karyotypes. We constructed a chromosome-level genome assembly with 24 pairs of chromosomes, 657.7 Mb in total length, and the contig and the scaffold N50 of 4.8 Mb and 28.6 Mb respectively. We also identified repetitive elements (accounting for 19.7% of the genome assembly) and predicted 26,905 protein-coding genes. Meanwhile, we performed genome-wide localization and characterization of several important genes encoding some key enzymes in the biosynthesis pathway of PUFAs. These genes may contribute to the high concentration of omega-3 in jade perch. Moreover, we conducted a series of comparative genomic analyses among four representative bass species at a chromosome level, resulting in a series of sequences of a deductive bass ancestor genome.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina Sequencing • Oxford Nanopore Sequencing

Transcriptomic comparison sheds new light on regulatory networks for dimorphic flower development in response to photoperiod in Viola prionantha

BACKGROUND

Chasmogamous (CH)-cleistogamous (CL) dimorphic flowers are developed in Viola prionantha. However, the environmental and genetic factors necessary for the CH-CL transition are unknown.

RESULTS

In the present work, short-day (SD) conditions induced CH flowers, whereas long days (LDs) triggered CL flowers in V. prionantha. Compared to fully developed CH flowers, CL flowers had less mature stamens, no nectar glands, and immature petals. Comparative transcriptomics revealed differentially expressed genes (DEGs) during CL and CH development. Core genes in the photoperiod pathway, such as V. prionantha orthologs of GIGANTEA (GI), CONSTANS (CO), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), which promote floral induction, were highly expressed in CL flowers, whereas UNUSUAL FLORAL ORGANS (UFO) and B-class MADS-box genes for floral organ identity and development showed an opposite alteration. Moreover, genes in the glycolytic process, sucrose metabolic process, and fatty acid biosynthetic process were all highly expressed in CH flowers. Interestingly, V. prionantha orthologs of the B-class MADS-box genes APETALA3 (AP3) and PISTILLATA (PI) might relate to these sugar-fatty acid processes and were co-expressed with GAIP-B-like and YABBY5 (YAB5), which regulate the development of the petal, stamen, and nectary. Compared to CH flowers, DEGs and hub genes in the most significantly correlated modules of the gene co-expression network, which are involved in abiotic and biotic responses, were upregulated in CL flowers.

CONCLUSIONS

We proposed an integrative model for transcription regulation of genes in the photoperiod pathway, floral organ development, stress response, and sugar-fatty acid processes to determine CH-CL flower development in V. prionantha. Particularly, under LDs, activated GI may induce genes involved in the stress-response pathways, and then downregulated AP3 and PI or UFO to inhibit the sugar-fatty acid metabolic processes, together forming CL flowers. In contrast, CH flowers were produced under SDs. This work provides novel insights into the developmental evolution of dimorphic flowers in Viola.

Whole-Genome Sequencing and Comparative Genome Analysis of Fusarium solani-melongenae Causing Fusarium Root and Stem Rot in Sweetpotatoes

Sweetpotato (Ipomoea batatas) is the eighth most important crop globally. However, the production and quality of sweetpotatoes are threatened by Fusarium diseases that are prevalent around the world. In this study, a Fusarium species that causes root and stem rot in sweetpotatoes was studied. The pathogenic fungus CRI 24-3 was isolated and sequenced using third- and next-generation sequencing techniques and a 49.6 Mb chromosome-level draft genome containing 15,374 putative coding genes were obtained. Molecular phylogenetic analysis showed that CRI 24-3 was an F. solani-melongenae strain within clade 3 of the F. solani species complex (FSSC). CRI 24-3 showed a relatively high number of virulence factors, such as carbohydrate-active enzymes (CAZymes), pathogen-host interaction (PHI) proteins, and terpene synthases (TSs), compared with the number of those identified in other sequenced FSSC members. Comparative genome analysis revealed considerable conservation and unique characteristics between CRI 24-3 and other FSSC species. In conclusion, the findings in the current study provide important genetic information about F. solani-melongenae and should be useful in the exploration of pathogenicity mechanisms and the development of Fusarium disease management strategies.

IMPORTANCE Fusarium root and stem rot in sweetpotato are prevalent in the main sweetpotato-growing areas in China, and fungal isolation, morphological characteristics, and molecular phylogenetic analysis of the disease causal agent (F. solani-melongenae isolate CRI 24-3) were systematically studied. The genome sequence of F. solani-melongenae isolates CRI 24-3 was first reported, which should provide a basis for genome assembly of other closely related Fusarium species. Carbohydrate-active enzymes predicted in CRI 24-3 may be important to convert the substantial polysaccharides to sustainable and renewable energy. Moreover, other virulence factors facilitating Fusarium diseases, including effectors and toxic secondary metabolites, are ideal objects for pathogenicity mechanism research and molecular targets for fungicide development. The findings of comparative genome analysis of CRI 24-3 and 15 sequenced members of the F. solani species complex help promote an integral understanding of genomic features and evolutionary relationships in Fusarium.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

A survey of transcriptome complexity using full-length isoform sequencing in the tea plant Camellia sinensis

Tea is one of the most popular beverages and its leaves are rich in catechins, contributing to the diverse flavor as well as beneficial for human health. However, the study of the post-transcriptional regulatory mechanism affecting the synthesis of catechins remains insufficient. Here, we sequenced the transcriptome using PacBio sequencing technology and obtained 63,111 full-length high-quality isoforms, including 1302 potential novel genes and 583 highly reliable fusion transcripts. We also identified 1204 lncRNAs with high quality, containing 188 known and 1016 novel lncRNAs. In addition, 311 mis-annotated genes were corrected based on the high-quality Isoseq reads. A large number of alternative splicing (AS) events (3784) and alternative polyadenylation (APA) genes (18,714) were analyzed, accounting for 8.84% and 43.7% of the total annotated genes, respectively. We also found that 2884 genes containing AS and APA features exhibited higher expression levels than other genes. These genes are mainly involved in amino acid biosynthesis, carbon fixation in photosynthetic organisms, phenylalanine, tyrosine, tryptophan biosynthesis, and pyruvate metabolism, suggesting that they play an essential role in the catechins content of tea polyphenols. Our results further improved the level of genome annotation and indicated that post-transcriptional regulation plays a crucial part in synthesizing catechins.

Transcript Complexity and New Insights of Restorer Line in CMS-D8 Cotton Through Full-Length Transcriptomic Analysis

Hybrid utilization has proficiently increased crop production worldwide. The cytoplasmic male sterility (CMS) system has emerged as an efficient tool for commercial hybrid cotton seed production. The restorer line with dominant Rf₂ gene can restore the fertility of the CMS-D8 sterile line. However, the molecular mechanism of fertility restoration remains unclear in CMS-D8 cotton that limits wider utilization of three-line hybrid breeding. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to understand fertility restoration mechanism of CMS-D8 cotton. In total, 228,106 full-length non-chimeric transcriptome sequences were obtained from anthers of developing flowering buds. The analysis results identified 3,174 novel isoforms, 2,597 novel gene loci, 652 long non-coding RNAs predicted from novel isoforms, 7,234 alternative splicing events, 114 fusion transcripts, and 1,667 genes with alternative polyadenylation. Specially, two novel genes associated with restoration function, Ghir_D05.742.1 and m64033_190821_201011/21103726/ccs were identified and showed significant higher levels of expression in restorer line than sterile and maintainer lines. Our comparative full-length transcriptome analysis provides new insights into the molecular function of Rf₂ fertility restorer gene. The results of this study offer a platform for fertility restoration candidate gene discovery in CMS-D8 cotton.

Genome-Wide Identification of Genes Related to Biosynthesis of Phenolic Acid Derivatives in Bletilla striata at Different Suspension Culture Stages

To screen the genes regulating the biosynthesis of phenolic acid derivatives from the genome of Bletilla striata, we designed a suspension culture system to sample the cells for the following experiments. The contents of four phenolic acid derivatives were determined by high-performance liquid chromatography, and several full-length transcriptome sequencings of RNA samples at 10 time points were performed for bioinformatics analysis. The correlation analysis was used to identify and verify the key DEGs involved in the biosynthesis of the four phenolic acid derivatives. The results showed that the contents of p-hydroxybenzylalcohol (HBA), Dactylorhin A, Militarine, and Coelonin peaked at 33 days postinoculation (Dpi), 18 Dpi, 39 Dpi, and 39 Dpi of the culture system, respectively. Based on transcriptome data, 80 DEGs involved in the biosynthesis of phenolic acid derivatives were obtained. The KEGG pathway enrichment analysis classified them mostly into five metabolic pathways: phenylpropane biosynthesis, starch and sucrose metabolic, cyanoamino acid metabolism, gluconeogenesis and glycolysis, and phenylalanine metabolism. qPCR analysis revealed that the relative gene expression levels were consistent with the overall trend of transcriptome sequencing results. Among them, 14, 18, 23, and 41 unigenes were found to be involved in the synthesis of HBA, Dactylorhin A, Coelonin, and Militarine, respectively. These unigenes laid a solid foundation for elucidating the biosynthesis mechanism of phenolic acid derivatives in suspension cells of B. striata.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

Symbiodiniaceae microRNAs and their targeting sites in coral holobionts: A transcriptomics-based exploration

Corals should make excellent models for cross-kingdom research because of their natural animal-photobiont holobiont composition, yet a lack of studies and experimental data restricts their use. Here we integrate new full-length transcriptomes and small RNAs of four common reef-building corals with the published Cladocopium genomes to gain deeper insight into gene regulation in coral-Symbiodiniaceae holobionts. Eleven novel Symbiodiniaceae miRNAs get identified, and enrichment results of their target genes show that they might play a role in downregulating rejection from host coral cells, protecting symbiont from autophagy and apoptosis in parallel. This work provides evidence for the early origin of cross-kingdom regulation as a mechanism of self-defense autotrophs can use against heterotrophs, sheds more light on coral-Symbiodiniaceae holobionts, and contributes valuable data for further coral research.

Transcriptomic and proteomic analyses of Cucurbita ficifolia Bouché (Cucurbitaceae) response to Fusarium oxysporum f.sp. cucumerium

Background

Fusarium oxysporum f. sp. cucumerinum (FOC) is the causal agent of cucumber Fusarium wilt, which can cause extensive damages and productivity losses. Cucurbita ficifolia Bouché (Cucurbitaceae) is usually used as rootstock for cucumber because of its excellent resistance to Fusarium wilt. Our previous study found that C.ficifolia has high FOC resistance, the underlying mechanism of which is unclear.

Results

Transcriptome and proteome profiling was performed on the basis of RNA-Seq and isobaric tag for relative and absolute quantitation technology to explore the molecular mechanisms of the response of Cucurbita ficifolia Bouché to Fusarium oxysporum f. sp. cucumerium infection. Comparative analyses revealed that 1850 genes and 356 protein species were differentially regulated at 2d and 4d after FOC inoculation. However, correlation analysis revealed that only 11 and 39 genes were differentially regulated at both the transcriptome and proteome levels after FOC inoculation at 2d and 4d, respectively. After FOC inoculation, plant hormones signal transduction, transcription factors were stimulated, whereas wax biosynthesis and photosynthesis were suppressed. Increased synthesis of oxidative-redox proteins is involved in resistance to FOC.

Conclusions

This study is the first to reveal the response of C. ficifolia leaf to FOC infection at the transcriptome and proteome levels, and to show that FOC infection activates plant hormone signaling and transcription factors while suppressing wax biosynthesis and photosynthesis. The accumulation of oxidative-redox proteins also plays an important role in the resistance of C. ficifolia to FOC. Results provide new information regarding the processes of C. ficifolia leaf resistance to FOC and will contribute to the breeding of cucumber rootstock with FOC resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08674-7.

Comparative analysis of mitochondrial genomes provides insights into the mechanisms underlying an S-type cytoplasmic male sterility (CMS) system in wheat (Triticum aestivum L.)

Cytoplasmic male sterility (CMS) has been widely used in crop cross breeding. There has been much research on wheat CMS. However, the correlation between S-type CMS and mitochondrial genome remains elusive. Herein, we sequenced the mitochondrial genome of wheat CMS line and compared it with the maintainer line. The results showed that the mitochondrial genome of CMS line encoded 26 tRNAs, 8 rRNAs, and 35 protein-coding genes, and the cob encoding complex III in which the protein coding gene is mutated. This protein is known to affect reactive oxygen (ROS) production. The analysis of ROS metabolism in developing anthers showed that the deficiency of antioxidants and antioxidant enzymes in the sterile system aggravated membrane lipid oxidation, resulting in ROS accumulation, and influencing the anther development. Herein, cob is considered as a candidate causative gene sequence for CMS.

Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L

Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca²⁺ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.

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.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus)

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

A rust fungus effector directly binds plant pre-mRNA splice site to reprogram alternative splicing and suppress host immunity

Alternative splicing (AS) is a crucial post-transcriptional regulatory mechanism in plant resistance. However, whether and how plant pathogens target splicing in their host remains mostly unknown. For example, although infection by Puccinia striiformis f. sp. tritici (Pst), a pathogenic fungus that severely affects the yield of wheat worldwide, has been shown to significantly influence the levels of alternatively spliced transcripts in the host, the mechanisms that govern this process, and its functional consequence have not been examined. Here, we identified Pst_A23 as a new Pst arginine-rich effector that localizes to host nuclear speckles, nuclear regions enriched in splicing factors. We demonstrated that transient expression of Pst_A23 suppresses plant basal defence dependent on the Pst_A23 nuclear speckle localization and that this protein plays an important role in virulence, stable silencing of which improves wheat stripe rust resistance. Remarkably, RNA-Seq data revealed that AS patterns of 588 wheat genes are altered in Pst_A23-overexpressing lines compared to control plants. To further examine the direct relationship between Pst_A23 and AS, we confirmed direct binding between two RNA motifs predicted from these altered splicing sites and Pst_A23 in vitro. The two RNA motifs we chose occur in the cis-element of TaXa21-H and TaWRKY53, and we validated that Pst_A23 overexpression results in decreased functional transcripts of TaXa21-H and TaWRKY53 while silencing of TaXa21-H and TaWRKY53 impairs wheat resistance to Pst. Overall, this represents formal evidence that plant pathogens produce 'splicing' effectors, which regulate host pre-mRNA splicing by direct engagement of the splicing sites, thereby interfering with host immunity.

RegScaf: a regression approach to scaffolding

MOTIVATION

Crucial to the correctness of a genome assembly is the accuracy of the underlying scaffolds that specify the orders and orientations of contigs together with the gap distances between contigs. The current methods construct scaffolds based on the alignments of 'linking' reads against contigs. We found that some 'optimal' alignments are mistaken due to factors such as the contig boundary effect, particularly in the presence of repeats. Occasionally, the incorrect alignments can even overwhelm the correct ones. The detection of the incorrect linking information is challenging in any existing methods.

RESULTS

In this study, we present a novel scaffolding method RegScaf. It first examines the distribution of distances between contigs from read alignment by the kernel density. When multiple modes are shown in a density, orientation-supported links are grouped into clusters, each of which defines a linking distance corresponding to a mode. The linear model parameterizes contigs by their positions on the genome; then each linking distance between a pair of contigs is taken as an observation on the difference of their positions. The parameters are estimated by minimizing a global loss function, which is a version of trimmed sum of squares. The least trimmed squares estimate has such a high breakdown value that it can automatically remove the mistaken linking distances. The results on both synthetic and real datasets demonstrate that RegScaf outperforms some popular scaffolders, especially in the accuracy of gap estimates by substantially reducing extremely abnormal errors. Its strength in resolving repeat regions is exemplified by a real case. Its adaptability to large genomes and TGS long reads is validated as well.

AVAILABILITY AND IMPLEMENTATION

RegScaf is publicly available at https://github.com/lemontealala/RegScaf.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-molecule real-time sequencing of the full-length transcriptome of purple garlic (Allium sativum L. cv. Leduzipi) and identification of serine O-acetyltransferase family proteins involved in cysteine biosynthesis

BACKGROUND

Garlic (Allium sativum L.), whose bioactive components are mainly organosulfur compounds (OSCs), is a herbaceous perennial widely consumed as a green vegetable and a condiment. Yet, the metabolic enzymes involved in the biosynthesis of OSCs are not identified in garlic.

RESULTS

Here, a full-length transcriptome of purple garlic was generated via PacBio and Illumina sequencing, to characterize the garlic transcriptome and identify key proteins mediating the biosynthesis of OSCs. Overall, 22.56 Gb of clean data were generated, resulting in 454 698 circular consensus sequence (CCS) reads, of which 83.4% (379 206) were identified as being full-length non-chimeric reads - their further transcript clustering facilitated identification of 36 571 high-quality consensus reads. Once corrected, their genome-wide mapping revealed that 6140 reads were novel isoforms of known genes, and 2186 reads were novel isoforms from novel genes. We detected 1677 alternative splicing events, finding 2902 genes possessing either two or more poly(A) sites. Given the importance of serine O-acetyltransferase (SERAT) in cysteine biosynthesis, we investigated the five SERAT homologs in garlic. Phylogenetic analysis revealed a three-tier classification of SERAT proteins, each featuring a serine acetyltransferase domain (N-terminal) and one or two hexapeptide transferase motifs. Template-based modeling showed that garlic SERATs shared a common homo-trimeric structure with homologs from bacteria and other plants. The residues responsible for substrate recognition and catalysis were highly conserved, implying a similar reaction mechanism. In profiling the five SERAT genes' transcript levels, their expression pattern varied significantly among different tissues.

CONCLUSION

This study's findings deepen our knowledge of SERAT proteins, and provide timely genetic resources that could advance future exploration into garlic's genetic improvement and breeding. © 2021 Society of Chemical Industry.

Transcriptome Analysis of Otodectes cynotis in Different Developmental Stages

The mite Otodectes cynotis is distributed worldwide and parasitism the ear canals of cats and dogs, causing otitis externa. Molecular biology of O. cynotis is poorly understood, with only a few genes being deposited in public databases. In the present study, we aimed to perform transcriptome analysis of O. cynotis using SMRT and Illumina sequencing of RNA from different development stages. SMRT-Seq of O. cynotis demonstrated 5,431 final transcripts, including 406 long non-coding RNAs and 2,698 differentially expressed genes (DEGs), including 1,357 up-regulated genes and 1,341 down-regulated genes between adult mites and nymph/larva. A total of 397 putative allergen genes were detected, 231 of which were DEGs. Among them, 77 were homologous of known mite allergens. The expression level of allergen genes hints at the pathogenicity of mites in different life stages, and the protein interaction network analysis could identify possible key genes in the pathogenic mechanism. Intriguingly, Gene Ontology analysis showed that most of the (DEGs) were associated with the terms hydrolase activity and proteolysis. Kyoto Encyclopedia of genes and genomes (KEGG) analysis identified drug metabolism-cytochrome P450 signal pathway as one of the top pathways. SMRT-Seq of the full-length transcriptome of O. cynotis was performed first, and a valuable resource was acquired through the combination analysis with the Illumina sequencing data. The results of our analyses provide new information for further research into Otodectes cynotis.

Full-Length Transcriptome Sequencing Reveals Tissue-Specific Gene Expression Profile of Mangrove Clam Geloina erosa

Mollusca is the second largest animal phylum and represents one of the most evolutionarily successful animal groups. Geloina erosa, a species of Corbiculidae, plays an important role in mangrove ecology. It is highly adaptable and can withstand environmental pollution and microbial infections. However, there is no reference genome or full-length transcriptome available for G. erosa. This impedes the study of the biological functions of its different tissues because transcriptome research requires reference genome or full-length transcriptome as a reference to improve accuracy. In this study, we applied a combination of Illumina and PacBio single-molecule real-time sequencing technologies to sequence the full-length transcriptomes of G. erosa tissues. Transcriptomes of nine samples obtained from three tissues (hepatopancreas, gill, and muscle) were sequenced using Illumina. Furthermore, we obtained 87,310 full-length reads non-chimeric sequences. After removing redundancy, 22,749 transcripts were obtained. The average Q score of 30 was 94.48%. In total, 271 alternative splicing events were predicted. There were 14,496 complete regions and 3,870 lncRNAs. Differential expression analysis revealed tissue-specific physiological functions. The gills mainly express functions related to filtration, metabolism, identifying pathogens and activating immunity, and neural activity. The hepatopancreas is the main tissue related to metabolism, it also involved in the immune response. The muscle mainly express functions related to muscle movement and control, it contains more energy metabolites that gill and hepatopancreas. Our research provides an important reference for studying the gene expression of G. erosa under various environmental stresses. Moreover, we present a reliable sequence that will provide an excellent foundation for further research on G. erosa.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense)

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Reconstructing viral haplotypes using long reads

MOTIVATION

Most RNA viruses lack strict proofreading during replication. Coupled with a high replication rate, some RNA viruses can form a virus population containing a group of genetically related but different haplotypes. Characterizing the haplotype composition in a virus population is thus important to understand viruses' evolution. Many attempts have been made to reconstruct viral haplotypes using next-generation sequencing (NGS) reads. However, the short length of NGS reads cannot cover distant single-nucleotide variants, making it difficult to reconstruct complete or near-complete haplotypes. Given the fast developments of third-generation sequencing technologies, a new opportunity has arisen for reconstructing full-length haplotypes with long reads.

RESULTS

In this work, we developed a new tool, RVHaplo to reconstruct haplotypes for known viruses from long reads. We tested it rigorously on both simulated and real viral sequencing data and compared it against other popular haplotype reconstruction tools. The results demonstrated that RVHaplo outperforms the state-of-the-art tools for viral haplotype reconstruction from long reads. Especially, RVHaplo can reconstruct the rare (1% abundance) haplotypes that other tools usually missed.

AVAILABILITY AND IMPLEMENTATION

The source code and the documentation of RVHaplo are available at https://github.com/dhcai21/RVHaplo.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data

Tenderness is an important indicator of meat quality. Novel isoforms associated with meat tenderness and the role of the CCCTC-binding factor (CTCF) in regulating alternative splicing to produce isoforms in sheep are largely unknown. The current project studied six sheep from two crossbred populations (Dorper × Hu × Hu, DHH and Dorper × Dorper × Hu, DDH) with divergent meat tenderness. Pooled Iso-seq data were used to annotate the sheep genomes. Then, the updated genome annotation and six RNA-seq data were combined to identify differentially expressed isoforms (DEIs) in muscles between DHH and DDH. These data were also combined with peaks detected from CTCF ChIP-seq data to investigate the regulatory role of CTCF for the alternative splicing. As a result, a total of 624 DEIs were identified between DDH and DHH. For example, isoform 7.524.18 transcribed from CAPN3 may be associated with meat tenderness. In addition, a total of 86 genes were overlapped between genes with transcribed DEIs and genes in differential peaks identified by CTCF ChIP-seq. Among these overlapped genes, ANKRD23 produces different isoforms which may be regulated by CTCF via methylation. As preliminary research, our results identified novel isoforms associated with meat tenderness and revealed the possible regulating mechanisms of alternative splicing to produce isoforms.

A chromosomal-level reference genome of the widely utilized Coccidioides posadasii laboratory strain "Silveira"

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

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.