On the study of microbial transcriptomes using second- and third-generation sequencing technologies

TPX2 upregulates MMP13 to promote the progression of lipopolysaccharide-induced osteoarthritis

Purpose

This study seeks to identify potential clinical biomarkers for osteoarthritis (OA) using bioinformatics and investigate OA mechanisms through cellular assays.

Methods

Differentially Expressed Genes (DEGs) from GSE52042 (four OA samples, four control samples) were screened and analyzed with protein-protein interaction (PPI) analysis. Overlapping genes in GSE52042 and GSE206848 (seven OA samples, and seven control samples) were identified and evaluated using Gene Set Enrichment Analysis (GSEA) and clinical diagnostic value analysis to determine the hub gene. Finally, whether and how the hub gene impacts LPS-induced OA progression was explored by in vitro experiments, including Western blotting (WB), co-immunoprecipitation (Co-IP), flow cytometry, etc.

Result

Bioinformatics analysis of DEGs (142 up-regulated and 171 down-regulated) in GSE52042 identified two overlapping genes (U2AF2, TPX2) that exhibit significant clinical diagnostic value. These genes are up-regulated in OA samples from both GSE52042 and GSE206848 datasets. Notably, TPX2, which AUC = 0.873 was identified as the hub gene. In vitro experiments have demonstrated that silencing TPX2 can alleviate damage to chondrocytes induced by lipopolysaccharide (LPS). Furthermore, there is a protein interaction between TPX2 and MMP13 in OA. Excessive MMP13 can attenuate the effects of TPX2 knockdown on LPS-induced changes in OA protein expression, cell growth, and apoptosis.

Conclusion

In conclusion, our findings shed light on the molecular mechanisms of OA and suggested TPX2 as a potential therapeutic target. TPX2 could promote the progression of LPS-induced OA by up-regulating the expression of MMP13, which provides some implications for clinical research.

Gut Microbiota in Children with Hand Foot and Mouth Disease on 16S rRNA Gene Sequencing

Hand foot and mouth disease (HFMD) is a contagious and seasonal viral disease in children. The gut microbiota of HFMD children is not clear now. The study aimed to explore the gut microbiota of HFMD children. The 16S rRNA gene of the gut microbiota of ten HFMD patients and ten healthy children were sequenced on the NovaSeq and PacBio platforms respectively. There were significant differences in gut microbiota between the patients and healthy children. The diversity and abundance of gut microbiota in HFMD patients were lower than that in healthy children. The species Roseburia inulinivorans and Romboutsia timonensis were more abundant in healthy children than those in HFMD patients, which suggests that the two species may be used as probiotics for adjusting the gut microbiota of HFMD patients. Meanwhile, the results of 16S rRNA gene sequences from the two platforms were different. The NovaSeq platform identified more microbiota and has the characteristics of high throughput, short time and low price. However, the NovaSeq platform has low resolution at the species level. The PacBio platform has high resolution based on its long reads length, which is more suitable for species-level analysis. But, the shortcomings of the high price and low throughput of PacBio still need to be overcome. With the development of sequencing technology, the reduction in sequencing price and the increase in throughput will promote the third-generation sequencing technology used in the study of gut microbes.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

Nanopore sequencing of RNA and cDNA molecules in Escherichia coli

High-throughput sequencing dramatically changed our view of transcriptome architectures and allowed for ground-breaking discoveries in RNA biology. Recently, sequencing of full-length transcripts based on the single-molecule sequencing platform from Oxford Nanopore Technologies (ONT) was introduced and is widely used to sequence eukaryotic and viral RNAs. However, experimental approaches implementing this technique for prokaryotic transcriptomes remain scarce. Here, we present an experimental and bioinformatic workflow for ONT RNA-seq in the bacterial model organism Escherichia coli, which can be applied to any microorganism. Our study highlights critical steps of library preparation and computational analysis and compares the results to gold standards in the field. Furthermore, we comprehensively evaluate the applicability and advantages of different ONT-based RNA sequencing protocols, including direct RNA, direct cDNA, and PCR-cDNA. We find that (PCR)-cDNA-seq offers improved yield and accuracy compared to direct RNA sequencing. Notably, (PCR)-cDNA-seq is suitable for quantitative measurements and can be readily used for simultaneous and accurate detection of transcript 5' and 3' boundaries, analysis of transcriptional units, and transcriptional heterogeneity. In summary, based on our comprehensive study, we show nanopore RNA-seq to be a ready-to-use tool allowing rapid, cost-effective, and accurate annotation of multiple transcriptomic features. Thereby nanopore RNA-seq holds the potential to become a valuable alternative method for RNA analysis in prokaryotes.

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

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

A chromosome-scale genome assembly of a diploid alfalfa, the progenitor of autotetraploid alfalfa

Alfalfa (Medicago sativa L.) is one of the most important and widely cultivated forage crops. It is commonly used as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Based on Illumina, Nanopore and Hi-C data, we assembled a chromosome-scale assembly of Medicago sativa spp. caerulea (voucher PI464715), the direct diploid progenitor of autotetraploid alfalfa. The assembled genome comprises 793.2 Mb of genomic sequence and 47,202 annotated protein-coding genes. The contig N50 length is 3.86 Mb. This genome is almost twofold larger and contains more annotated protein-coding genes than that of its close relative, Medicago truncatula (420 Mb and 44,623 genes). The more expanded gene families compared with those in M. truncatula and the expansion of repetitive elements rather than whole-genome duplication (i.e., the two species share the ancestral Papilionoideae whole-genome duplication event) may have contributed to the large genome size of M. sativa spp. caerulea. Comparative and evolutionary analyses revealed that M. sativa spp. caerulea diverged from M. truncatula ~5.2 million years ago, and the chromosomal fissions and fusions detected between the two genomes occurred during the divergence of the two species. In addition, we identified 489 resistance (R) genes and 82 and 85 candidate genes involved in the lignin and cellulose biosynthesis pathways, respectively. The near-complete and accurate diploid alfalfa reference genome obtained herein serves as an important complement to the recently assembled autotetraploid alfalfa genome and will provide valuable genomic resources for investigating the genomic architecture of autotetraploid alfalfa as well as for improving breeding strategies in alfalfa.

Next-generation technologies for studying host-pathogen interactions: a focus on dual transcriptomics, CRISPR/Cas9 screening and organs-on-chips

Pathogens constantly interact with their hosts and the environment, and therefore have evolved unique virulence mechanisms to target and breach host defense barriers and manipulate host immune response to establish an infection. Advances in technologies that allow genome mining, gene editing such as CRISPR/Cas9, genomic, epigenomic and transcriptomic studies such as dual RNA-seq, coupled with bioinformatics, have accelerated the field of host-pathogen interactions within a broad range of infection models. Underpinning of the molecular changes that accompany invasion of eukaryotic cells with pathogenic microorganisms at the intersection of host, pathogen and their local environment has provided a better understanding of infectious disease mechanisms and antimicrobial strategies. The recent evolution of physiologically relevant three-dimensional (3-D) tissue/organ models and microfluidic organ-on-chip devices also provided a window to a more predictive framework of infectious disease processes. These approaches combined hold the potential to highly impact discovery of novel drug targets and vaccine candidates of the future. Here, we review three of the available and emerging technologies-dual RNA-seq, CRISPR/Cas9 screening and organs-on-chips, applicable to the high throughput study and deciphering of interaction networks between pathogens and their hosts that are critical for the development of novel therapeutics.

Two's company: studying interspecies relationships with dual RNA-seq

Organisms do not exist isolated from each other, but constantly interact. Cells can sense the presence of interaction partners by a range of receptors and, via complex regulatory networks, specifically react by changing the expression of many of their genes. Technological advances in next-generation sequencing over the recent years now allow us to apply RNA sequencing to two species at the same time (dual RNA-seq), and thus to directly study the gene expression of two interacting species without the need to physically separate cells or RNA. In this review, we give an overview over the latest studies in interspecies interactions made possible by dual RNA-seq, ranging from pathogenic to symbiotic relationships. We summarize state-of-the-art experimental techniques, bioinformatic data analysis and data interpretation, while also highlighting potential problems and pitfalls starting from the selection of meaningful time points and number of reads to matters of rRNA depletion. A short outlook on new trends in the field of dual RNA-seq concludes this review, looking at sequencing of non-coding RNAs during host-pathogen interactions and the prediction of molecular interspecies interactions networks.

Taxonomy and systematics of plant probiotic bacteria in the genomic era

Recent decades have predicted significant changes within our concept of plant endophytes, from only a small number specific microorganisms being able to colonize plant tissues, to whole communities that live and interact with their hosts and each other. Many of these microorganisms are responsible for health status of the plant, and have become known in recent years as plant probiotics. Contrary to human probiotics, they belong to many different phyla and have usually had each genus analysed independently, which has resulted in lack of a complete taxonomic analysis as a group. This review scrutinizes the plant probiotic concept, and the taxonomic status of plant probiotic bacteria, based on both traditional and more recent approaches. Phylogenomic studies and genes with implications in plant-beneficial effects are discussed. This report covers some representative probiotic bacteria of the phylum Proteobacteria, Actinobacteria, Firmicutes and Bacteroidetes, but also includes minor representatives and less studied groups within these phyla which have been identified as plant probiotics.