Short Read Alignment Using SOAP2

Single-cell transcriptomics: background, technologies, applications, and challenges

Single-cell sequencing was developed as a high-throughput tool to elucidate unusual and transient cell states that are barely visible in the bulk. This technology reveals the evolutionary status of cells and differences between populations, helps to identify unique cell subtypes and states, reveals regulatory relationships between genes, targets and molecular mechanisms in disease processes, tumor heterogeneity, the state of the immune environment, etc. However, the high cost and technical limitations of single-cell sequencing initially prevented its widespread application, but with advances in research, numerous new single-cell sequencing techniques have been discovered, lowering the cost barrier. Many single-cell sequencing platforms and bioinformatics methods have recently become commercially available, allowing researchers to make fascinating observations. They are now increasingly being used in various industries. Several protocols have been discovered in this context and each technique has unique characteristics, capabilities and challenges. This review presents the latest advancements in single-cell transcriptomics technologies. This includes single-cell transcriptomics approaches, workflows and statistical approaches to data processing, as well as the potential advances, applications, opportunities and challenges of single-cell transcriptomics technology. You will also get an overview of the entry points for spatial transcriptomics and multi-omics.

Variation in Bombyx mori immune response against fungal pathogen Beauveria bassiana with variability in cell wall β-1,3-glucan

Entomopathogenic fungi are protected by a cell wall with dynamic structure for adapting to various environmental conditions. β-1,3-Glucan recognition proteins activate the innate immune system of insects by recognizing surface molecules of fungi. However, the associations between pathogenicity and the different components of entomopathogenic fungal cell walls remain unclear. Three Beauveria bassiana strains were selected that have significantly differing virulence against Bombyx mori. The molecular mechanisms underlying the immune response in B. mori were investigated using RNA sequencing, which revealed differences in the immune response to different B. bassiana strains at 12 h post-infection. Immunofluorescence assays revealed that β-1,3-glucan content had an opposite trend to that of fungal virulence. β-1,3-Glucan injection upregulated BmβGRP4 expression and significantly reduced the virulence of the high-virulence strain but not that of the medium-virulence or low-virulence strains. BmβGRP4 silencing in B. mori with RNA interference resulted in the opposite virulence pattern, indicating that the virulence of B. bassiana was affected by the cell walls' content of β-1,3-glucan, which could be recognized by BmβGRP4. Furthermore, interference with the gene CnA (calcineurin catalytic A subunit) involved in β-1,3-glucan synthesis eliminated differences in virulence between B. bassiana strains. These results indicate that strains of a single species of pathogenic fungi that have differing cell wall components are recognized differently by the innate immune system of B. mori.

A new and rare type of hepatocellular carcinoma: Survival and gene analysis of portal vein tumour thrombus-type hepatocellular carcinoma

BACKGROUND

Portal vein tumour thrombus (PVTT) in patients with hepatocellular carcinoma (HCC) is known as a major complication associated with poor survival. We clinically defined a new and rare type of HCC, PVTT-type HCC (PVTT-HCC), in a small group of HCC patients with HCC presenting only as PVTT without a demonstrable parenchyma tumour. The clinicopathological and biological features of PVTT-HCC are not clear.

METHODS

The data for patients who had PVTT-HCC with histologically confirmed HCC from January 2004 to December 2012 at the Eastern Hepatobiliary Surgery Hospital were retrospectively analysed. The survival outcomes of patients with PVTT-HCC were compared with those of HCC patients with PVTT (HCC-PVTT). Propensity score matching (PSM) analysis was performed to match patients at a ratio of 1:3. Then, we performed RNA-Seq analysis of liver samples from PVTT-HCC and HCC-PVTT patients to identify and compare differentially expressed genes and biological pathways between the two groups.

RESULTS

We observed and collected 10 rare cases of PVTT-HCC and performed a prospective cohort study to compare overall survival (OS) between PVTT-HCC and HCC-PVTT. PVTT invaded the main portal vein in 10 PVTT-HCC patients. Univariate and multivariate analyses demonstrated that ChildPugh (A/B), different treatments (LR/non-LR), and different groups were independent risk factors for OS. The median OS was 10.3 months (95 % CI = 6.7-13.8) in the HCC-PVTT group and 7.5 months (95 % CI = 2.8-12.1) in the PVTT-HCC group (P = 0.042). From RNA-Seq, 1630 differentially expressed genes were obtained, of which 731 were upregulated and 899 downregulated in PVTT-HCC compared with HCC-PVTT.

CONCLUSIONS

The survival outcomes of patients with PVTT-HCC were worse than those of patients with HCC-PVTT. RNA-Seq demonstrated differential gene expression between PVTT-HCC and HCC-PVTT, indicating that the former may have distinguishing biological characteristics and be a new and rare type of HCC.

Directed Evolution Detects Supernumerary Centric Chromosomes Conferring Resistance to Azoles in Candida auris

Candida auris exhibits resistance to multiple antifungal drug classes and sterilization agents, posing threats to the immunocompromised worldwide. Among the four major geographical clades, the East Asian clade 2 isolates of C. auris are mostly drug susceptible. In this study, we experimentally evolved one such drug-susceptible isolate for multiple generations in the presence of the antifungal compound fluconazole and analyzed changes in the karyotype, DNA sequence, and gene expression profiles in three evolved drug-resistant isolates. Next-generation sequencing and electrophoretic karyotyping confirm the presence of segmental aneuploidy as supernumerary chromosomes originating from centromere-inclusive chromosomal duplication events in two such cases. A 638-kb region and a 675-kb region, both of which originated from chromosome 5 and contained its centromere region, are instances of supernumerary chromosome formation identified in two evolved fluconazole-resistant isolates. Loss of the supernumerary chromosomes from the drug-resistant isolates results in a complete reversal of fluconazole susceptibility. Transcriptome analysis of the third isolate identified overexpression of drug efflux pumps as a possible non-aneuploidy-driven mechanism of drug resistance. Together, this study reveals how both aneuploidy-driven and aneuploidy-independent mechanisms may operate in parallel in an evolving population of C. auris in the presence of an antifungal drug, in spite of starting from the same strain grown under similar conditions, to attain various levels of fluconazole resistance. IMPORTANCE Fungal pathogens develop drug resistance through multiple pathways by acquiring gene mutations, increasing the copy number of genes, or altering gene expression. In this study, we attempt to understand the mechanisms of drug resistance in the recently emerged superbug, C. auris. One approach to studying this aspect is identifying various mechanisms operating in drug-resistant clinical isolates. An alternative approach is to evolve a drug-susceptible isolate in the presence of an antifungal compound and trace the changes that result in drug resistance. Here, we evolve a drug-susceptible isolate of C. auris in the laboratory in the presence of a widely used antifungal compound, fluconazole. In addition to the already known changes like overexpression of drug efflux pumps, this study identifies a novel mechanism of azole resistance by the emergence of additional chromosomes through segmental duplication of chromosomal regions, including centromeres. The centric supernumerary chromosome helps stable amplification of a set of genes with an extra copy to confer fluconazole resistance.

Draft Genome Sequence of an Epiphytic Strain, Bacillus sp. Strain WL1, Isolated from the Surface of Nostoc flagelliforme Colonies in Yinchuan, Ningxia, China

Bacillus species are Gram-positive, aerobic, spore-forming bacteria that are widely spread in soil, dust, and water. One strain, Bacillus sp. strain WL1, was isolated from the surface of the cyanobacterium Nostoc flagelliforme in Yinchuan, Ningxia, China. The draft genome sequence of this strain was 5.62 Mbp.

Performance evaluation of pipelines for mapping, variant calling and interval padding, for the analysis of NGS germline panels

Background

Next-generation sequencing (NGS) represents a significant advancement in clinical genetics. However, its use creates several technical, data interpretation and management challenges. It is essential to follow a consistent data analysis pipeline to achieve the highest possible accuracy and avoid false variant calls. Herein, we aimed to compare the performance of twenty-eight combinations of NGS data analysis pipeline compartments, including short-read mapping (BWA-MEM, Bowtie2, Stampy), variant calling (GATK-HaplotypeCaller, GATK-UnifiedGenotyper, SAMtools) and interval padding (null, 50 bp, 100 bp) methods, along with a commercially available pipeline (BWA Enrichment, Illumina®). Fourteen germline DNA samples from breast cancer patients were sequenced using a targeted NGS panel approach and subjected to data analysis.

Results

We highlight that interval padding is required for the accurate detection of intronic variants including spliceogenic pathogenic variants (PVs). In addition, using nearly default parameters, the BWA Enrichment algorithm, failed to detect these spliceogenic PVs and a missense PV in the TP53 gene. We also recommend the BWA-MEM algorithm for sequence alignment, whereas variant calling should be performed using a combination of variant calling algorithms; GATK-HaplotypeCaller and SAMtools for the accurate detection of insertions/deletions and GATK-UnifiedGenotyper for the efficient detection of single nucleotide variant calls.

Conclusions

These findings have important implications towards the identification of clinically actionable variants through panel testing in a clinical laboratory setting, when dedicated bioinformatics personnel might not always be available. The results also reveal the necessity of improving the existing tools and/or at the same time developing new pipelines to generate more reliable and more consistent data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04144-1.

Complete Genome Sequence of Curtobacterium sp. Strain YC1, Isolated from the Surface of Nostoc flagelliforme Colonies in Yinchuan, Ningxia, China

The genus Curtobacterium belongs to the family Microbacteriaceae, within the phylum Actinobacteria. This genus includes a wide range of Gram-positive species associated with plants and soils. Here, we report the genome sequence of a new strain, Curtobacterium sp. strain YC1, which was isolated from the surface of Nostoc flagelliforme colonies. The genome of this strain contains one chromosome and one plasmid, and its size is 3.4 Mb.

RNA sequencing: new technologies and applications in cancer research

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.

Deoxynivalenol globally affects the selection of 3' splice sites in human cells by suppressing the splicing factors, U2AF1 and SF1

Deoxynivalenol (DON) is one of the most abundant mycotoxins and has adverse effects on several biological processes, posing risks of protein synthesis-disrupting effects and ribotoxic response. Therefore, chronic exposure to DON would fundamentally reshape the global expression pattern. Whether DON causes toxic effects on mRNA splicing, a fundamental biological process, remains unclear. In this study, we found that administration of the relative low dosage of DON dramatically changed the alternative splicing of pre-mRNA in HepG2 cells. The overall number of transcripts with aberrant selection of 3' splice sites was significantly increased in DON-exposed HepG2 cells. This effect was further confirmed in two other human cell lines, HEK293 and Caco-2, suggesting that this DON-induced alteration in splicing patterns was universal in human cells. Among these DON-induced changes in alternative splicing, the expression levels of two related splicing factors, SF1 and U2AF1, which are essential for 3' splice site recognitions, were strongly suppressed. Overexpression of either of the two splicing factors strongly alleviated the DON-induced aberrant selection of 3' splice sites. Moreover, SF1 was required for human cell proliferation in DON exposure, and the restoration of SF1 expression partially reinstated the proliferation potential for DON-treated cells. In conclusion, our study suggests that DON, even at a low dosage, has great potential to change gene expression globally by affecting not only protein synthesis but also mRNA processing in human cells.

Accuracy assessment of fusion transcript detection via read-mapping and de novo fusion transcript assembly-based methods

BACKGROUND

Accurate fusion transcript detection is essential for comprehensive characterization of cancer transcriptomes. Over the last decade, multiple bioinformatic tools have been developed to predict fusions from RNA-seq, based on either read mapping or de novo fusion transcript assembly.

RESULTS

We benchmark 23 different methods including applications we develop, STAR-Fusion and TrinityFusion, leveraging both simulated and real RNA-seq. Overall, STAR-Fusion, Arriba, and STAR-SEQR are the most accurate and fastest for fusion detection on cancer transcriptomes.

CONCLUSION

The lower accuracy of de novo assembly-based methods notwithstanding, they are useful for reconstructing fusion isoforms and tumor viruses, both of which are important in cancer research.

Genetic Diversity Analysis of Olive Germplasm (Olea europaea L.) With Genotyping-by-Sequencing Technology

Olive (Olea europaea L.) is a very important edible oil crop and has been cultivated for about 4,000 years in the Mediterranean area. Due to its nutritional and economic importance, researches on germplasm characterization received extensive attention. In this study, using the genotyping-by-sequencing (GBS) technology, we carried out genetic diversity analysis on 57 olive cultivars with different geographical origins. In total, 73,482 high-quality single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) > 5%, call rate > 50%, and heterozygosity rate < 10% were obtained at the whole genome level. Genetic structure and phylogenetic analysis showed that the 57 olive cultivars could be classified into two groups (Group I and Group II). No clear geographical distributions of cultivars were observed generally between the two groups. The average nucleotide diversities (π) specific for Group I and Group II were 0.317 and 0.305. The fixation index (F_ST) between Group I and Group II was 0.033. In Group II, cultivars could be further divided into two subgroups (Group IIa and Group IIb), which seem to be associated with their fruit sizes. The five Chinese-bred cultivars were all clustered in Group II, showing a closer genetic relationship with those from the central Mediterranean region and limited genetic background. It is therefore necessary for Chinese olive breeding programs to incorporate other genetic basis by utilizing germplasm from the other regions particularly from the east Mediterranean region as breeding parents. The results showed that GBS is an effective marker choice for cultivar characterization and genetic diversity analysis in olive and will help us better understand the genetic backgrounds of the crop.

Genomic Features, Comparative Genomic Analysis, and Antimicrobial Susceptibility Patterns of Chryseobacterium arthrosphaerae Strain ED882-96 Isolated in Taiwan

Bacteria belonging to the genus Chryseobacterium are ubiquitously distributed in natural environments, plants, and animals. Except C. indologenes and C. gleum, other Chryseobacterium species rarely cause human diseases. This study reported the whole-genome features, comparative genomic analysis, and antimicrobial susceptibility patterns of C. arthrosphaerae ED882-96 isolated in Taiwan. Strain ED882-96 was collected from the blood of a patient who had alcoholic liver cirrhosis and was an intravenous drug abuser. This isolate was initially identified as C. indologenes by using matrix-assisted laser desorption ionization–time of flight mass spectrometry. The analysis of 16S ribosomal RNA gene sequence revealed that ED882-96 shared 100% sequence identity with C. arthrosphaerae type strain CC-VM-7^T. The results of whole-genome sequencing of ED882-96 showed two chromosome contigs and one plasmid. The total lengths of the draft genomes of chromosome and plasmid were 4,249,864 bp and 435,667 bp, respectively. The findings of both in silico DNA–DNA hybridization and average nucleotide identity analyses clearly demonstrated that strain ED882-96 was a species of C. arthrosphaerae. A total of 83 potential virulence factor homologs were predicted in the whole-genome sequencing of strain ED882-96. This isolate was resistant to all tested antibiotics, including β-lactams, β-lactam/β-lactamase inhibitor combinations, aminoglycosides, fluoroquinolones, tetracycline, glycylcycline, and trimethoprim-sulfamethoxazole. Only one antibiotic resistance gene was recognized in the plasmid. By contrast, many antibiotic resistance genes were identified in the chromosome. The findings of this study suggest that strain ED882-96 is a highly virulent and multidrug-resistant pathogen. Knowledge regarding genomic characteristics and antimicrobial susceptibility patterns provides valuable insights into this uncommon species.

Comparative Transcriptome Reveals Benzenoid Biosynthesis Regulation as Inducer of Floral Scent in the Woody Plant Prunus mume

Mei (Prunus mume) is a peculiar woody ornamental plant famous for its inviting fragrance in winter. However, in this valuable plant, the mechanism behind floral volatile development remains poorly defined. Therefore, to explore the floral scent formation, a comparative transcriptome was conducted in order to identify the global transcripts specifying flower buds and blooming flowers of P. mume. Differentially expressed genes were identified between the two different stages showing great discrepancy in floral volatile production. Moreover, according to the expression specificity among the organs (stem, root, fruit, leaf), we summarized one gene cluster regulating the benzenoid floral scent. Significant gene changes were observed in accordance with the formation of benzenoid, thus pointing the pivotal roles of genes as well as cytochrome-P450s and short chain dehydrogenases in the benzenoid biosynthetic process. Further, transcription factors like EMISSION OF BENZENOID I and ODORANT I performed the same expression pattern suggesting key roles in the management of the downstream genes. Taken together, these data provide potential novel anchors for the benzenoid pathway, and the insight for the floral scent induction and regulation mechanism in woody plants.