Genome-wide identification and functional analysis of circRNAs in Trichophyton rubrum conidial and mycelial stages.

BACKGROUND: Circular RNAs (circRNAs) are a group of noncoding RNAs that participate in gene expression regulation in various pathways. The essential roles of circRNAs have been revealed in many species. However, knowledge of circRNAs in fungi is still not comprehensive. RESULTS: Trichophyton rubrum (T. rubrum) is considered a model organism of human pathogenic filamentous fungi and dermatophytes. In this study, we performed a genome-wide investigation of circRNAs in T. rubrum based on high-throughput sequencing and ultimately identified 4254 circRNAs. Most of these circRNAs were specific to the conidial or mycelial stage, revealing a developmental stage-specific expression pattern. In addition, 940 circRNAs were significantly differentially expressed between the conidial and mycelial stages. PCR experiments conducted on seven randomly selected differentially expressed (DE-) circRNAs confirmed the circularized structures and relative expression levels of these circRNAs. Based on their genome locations, most circRNAs originated from intergenic regions, unlike those in plants and animals. Furthermore, we constructed circRNA-miRNA-mRNA regulatory networks that included 661 DE-circRNAs targeting 140 miRNAs and further regulating 2753 mRNAs. The relative expression levels of two randomly selected circRNA-miRNA-mRNA axes were investigated by qRT-PCR, and the competing endogenous RNA (ceRNA) network theory was validated. Functional enrichment analysis of the target genes suggested that they were significantly involved in posttranscriptional processes and protein synthesis as well as some small-molecule metabolism processes. CircRNAs are relatively more conserved in closely related dermatophytes but rarely conserved in distantly related species. Tru_circ07138_001 is a highly conserved circRNA that was conserved in all ten dermatophytes analyzed in our study and three distantly related species. Its host gene TERG_07138 was also highly conserved in two of these distantly related species Gallus gallus and Caenorhabditis elegans. The specific role of this circRNA deserves further exploration. CONCLUSIONS: Our study is the first to provide a global profile of circRNAs in T. rubrum as well as dermatophytes. These results could serve as valuable resources for research on circRNA regulatory mechanisms in fungi and reveal new insights for further investigation of the physical characteristics of these significant human fungal pathogens.

Proteogenomic Analysis and Discovery of Immune Antigens in Mycobacterium vaccae.

Tuberculosis (TB) is one of the leading causes of death worldwide, especially in developing countries. Neonatal BCG vaccination occurs in various regions, but the level of protection varies in different populations. Recently, Mycobacterium vaccae is found to be an immunomodulating therapeutic agent that could confer a significant level of protection against TB. It is the only vaccine in a phase III trial from WHO to assess its efficacy and safety in preventing TB disease in people with latent TB infection. However, the mechanism of immunotherapy of M. vaccae remains poorly understood. In this study, the full genome of M. vaccae was obtained by next-generation sequencing technology, and a proteogenomic approach was successfully applied to further perform genome annotation using high resolution and high accuracy MS data. A total of 3,387 proteins (22,508 unique peptides) were identified, and 581 proteins annotated as hypothetical proteins in the genome database were confirmed. Furthermore, 38 novel protein products not annotated at the genome level were detected and validated. Additionally, the translational start sites of 445 proteins were confirmed, and 98 proteins were validated through extension of their translational start sites based on N terminus-derived peptides. The physicochemical characteristics of the identified proteins were determined. Thirty-five immunogenic proteins of M. vaccae were identified by immunoproteomic analysis, and 20 of them were selected to be expressed and validated by Western blot for immunoreactivity to serum from patients infected with M. tuberculosis The results revealed that eight of them showed strong specific reactive signals on the immunoblots. Furthermore, cellular immune response was further examined and one protein displayed a higher cellular immune level in pulmonary TB patients. Twelve identified immunogenic proteins have orthologous in H37Rv and BCG. This is the first study to obtain the full genome and annotation of M. vaccae using a proteogenomic approach, and some immunogenic proteins that were validated by immunoproteomic analysis could contribute to the understanding of the mechanism of M. vaccae immunotherapy.

Proteogenomic Analysis of Trichophyton rubrum Aided by RNA Sequencing.

Infections caused by dermatophytes, Trichophyton rubrum in particular, are among the most common diseases in humans. In this study, we present a proteogenomic analysis of T. rubrum based on whole-genome proteomics and RNA-Seq studies. We confirmed 4291 expressed proteins in T. rubrum and validated their annotated gene structures based on 35 874 supporting peptides. In addition, we identified 323 novel peptides (not present in the current annotated protein database of T. rubrum) that can be used to enhance current T. rubrum annotations. A total of 104 predicted genes supported by novel peptides were identified, and 127 gene models suggested by the novel peptides that conflicted with existing annotations were manually assigned based on transcriptomic evidence. RNA-Seq confirmed the validity of 95% of the total peptides. Our study provides evidence that confirms and improves the genome annotation of T. rubrum and represents the first survey of T. rubrum genome annotations based on experimental evidence. Additionally, our integrated proteomics and multisourced transcriptomics approach provides stronger evidence for annotation refinement than proteomic data alone, which helps to address the dilemma of one-hit wonders (uncertainties supported by only one peptide).

Mutations of novel influenza A(H10N8) virus in chicken eggs and MDCK cells.

The recent emergence of human infection with influenza A(H10N8) virus is an urgent public health concern. Genomic analysis showed that the virus was conserved in chicken eggs but presented substantial adaptive mutations in MDCK cells. Our results provide additional evidence for the avian origin of this influenza virus.

VFDB 2012 update: toward the genetic diversity and molecular evolution of bacterial virulence factors.

The virulence factor database (VFDB, http://www.mgc.ac.cn/VFs/) has served as a comprehensive repository of bacterial virulence factors (VFs) for >7 years. Bacterial virulence is an exciting and dynamic field, due to the availability of complete sequences of bacterial genomes and increasing sophisticated technologies for manipulating bacteria and bacterial genomes. The intricacy of virulence mechanisms offers a challenge, and there exists a clear need to decipher the 'language' used by VFs more effectively. In this article, we present the recent major updates of VFDB in an attempt to summarize some of the most important virulence mechanisms by comparing different compositions and organizations of VFs from various bacterial pathogens, identifying core components and phylogenetic clades and shedding new light on the forces that shape the evolutionary history of bacterial pathogenesis. In addition, the 2012 release of VFDB provides an improved user interface.

Identification and characterisation of non-coding small RNAs in the pathogenic filamentous fungus Trichophyton rubrum.

BACKGROUND: Accumulating evidence demonstrates that non-coding RNAs (ncRNAs) are indispensable components of many organisms and play important roles in cellular events, regulation, and development. RESULTS: Here, we analysed the small non-coding RNA (ncRNA) transcriptome of Trichophyton rubrum by constructing and sequencing a cDNA library from conidia and mycelia. We identified 352 ncRNAs and their corresponding genomic loci. These ncRNA candidates included 198 entirely novel ncRNAs and 154 known ncRNAs classified as snRNAs, snoRNAs and other known ncRNAs. Further bioinformatic analysis detected 96 snoRNAs, including 56 snoRNAs that had been annotated in other organisms and 40 novel snoRNAs. All snoRNAs belonged to two major classes--C/D box snoRNAs and H/ACA snoRNAs--and their potential target sites in rRNAs and snRNAs were predicted. To analyse the evolutionary conservation of the ncRNAs in T. rubrum, we aligned all 352 ncRNAs to the genomes of six dermatophytes and to the NCBI non-redundant nucleotide database (NT). The results showed that most of the identified snRNAs were conserved in dermatophytes. Of the 352 ncRNAs, 102 also had genomic loci in other dermatophytes, and 27 were dermatophyte-specific. CONCLUSIONS: Our systematic analysis may provide important clues to the function and evolution of ncRNAs in T. rubrum. These results also provide important information to complement the current annotation of the T. rubrum genome, which primarily comprises protein-coding genes.

Novel SARS-like betacoronaviruses in bats, China, 2011.

To clarify the evolutionary relationships among betavoronaviruses that infect bats, we analyzed samples collected during 2010-2011 from 14 insectivorous bat species in China. We identified complete genomes of 2 novel betacoronaviruses in Rhinolophus pusillus and Chaerephon plicata bats, which showed close genetic relationships with severe acute respiratory syndrome coronaviruses.

Full genome of influenza A (H7N9) virus derived by direct sequencing without culture.

An epidemic caused by influenza A (H7N9) virus was recently reported in China. Deep sequencing revealed the full genome of the virus obtained directly from a patient's sputum without virus culture. The full genome showed substantial sequence heterogeneity and large differences compared with that from embryonated chicken eggs.

Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces.

Bats are natural hosts for a large variety of zoonotic viruses. This study aimed to describe the range of bat viromes, including viruses from mammals, insects, fungi, plants, and phages, in 11 insectivorous bat species (216 bats in total) common in six provinces of China. To analyze viromes, we used sequence-independent PCR amplification and next-generation sequencing technology (Solexa Genome Analyzer II; Illumina). The viromes were identified by sequence similarity comparisons to known viruses. The mammalian viruses included those of the Adenoviridae, Herpesviridae, Papillomaviridae, Retroviridae, Circoviridae, Rhabdoviridae, Astroviridae, Flaviridae, Coronaviridae, Picornaviridae, and Parvovirinae; insect viruses included those of the Baculoviridae, Iflaviridae, Dicistroviridae, Tetraviridae, and Densovirinae; fungal viruses included those of the Chrysoviridae, Hypoviridae, Partitiviridae, and Totiviridae; and phages included those of the Caudovirales, Inoviridae, and Microviridae and unclassified phages. In addition to the viruses and phages associated with the insects, plants, and bacterial flora related to the diet and habitation of bats, we identified the complete or partial genome sequences of 13 novel mammalian viruses. These included herpesviruses, papillomaviruses, a circovirus, a bocavirus, picornaviruses, a pestivirus, and a foamy virus. Pairwise alignments and phylogenetic analyses indicated that these novel viruses showed little genetic similarity with previously reported viruses. This study also revealed a high prevalence and diversity of bat astroviruses and coronaviruses in some provinces. These findings have expanded our understanding of the viromes of bats in China and hinted at the presence of a large variety of unknown mammalian viruses in many common bat species of mainland China.

A comprehensive survey of bat sarbecoviruses across China in relation to the origins of SARS-CoV and SARS-CoV-2.

SARS-CoV and SARS-CoV-2 have been thought to originate from bats. In this study, we screened pharyngeal and anal swabs from 13 064 bats collected between 2016 and 2021 at 703 locations across China for sarbecoviruses, covering almost all known southern hotspots, and found 146 new bat sarbecoviruses. Phylogenetic analyses of all available sarbecoviruses show that there are three different lineages-L1 as SARS-CoV-related CoVs (SARSr-CoVs), L2 as SARS-CoV-2-related CoVs (SC2r-CoVs) and novel L-R (recombinants of L1 and L2)-present in Rhinolophus pusillus bats, in the mainland of China. Among the 146 sequences, only four are L-Rs. Importantly, none belong in the L2 lineage, indicating that circulation of SC2r-CoVs in China might be very limited. All remaining 142 sequences belong in the L1 lineage, of which YN2020B-G shares the highest overall sequence identity with SARS-CoV (95.8%). The observation suggests endemic circulations of SARSr-CoVs, but not SC2r-CoVs, in bats in China. Geographic analysis of the collection sites in this study, together with all published reports, indicates that SC2r-CoVs may be mainly present in bats of Southeast Asia, including the southern border of Yunnan province, but absent in all other regions within China. In contrast, SARSr-CoVs appear to have broader geographic distribution, with the highest genetic diversity and sequence identity to human sarbecoviruses along the southwest border of China. Our data provide the rationale for further extensive surveys in broader geographical regions within, and beyond, Southeast Asia in order to find the most recent ancestors of human sarbecoviruses.

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China.

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014-2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10-4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Complete genome sequence of the nitrogen-fixing and rhizosphere-associated bacterium Pseudomonas stutzeri strain DSM4166.

We present here the analysis of the whole-genome sequence of Pseudomonas stutzeri strain DSM4166, a diazotrophic isolate from the rhizosphere of a Sorghum nutans cultivar. To our knowledge, this is the second genome to be sequenced for P. stutzeri. The availability and analysis of the genome provide insight into the evolution of the nitrogen fixation property and identification of rhizosphere competence traits required in interactions with host plants.

Unbiased parallel detection of viral pathogens in clinical samples by use of a metagenomic approach.

Viral infectious diseases represent a major threat to public health and are among the greatest disease burdens worldwide. Rapid and accurate identification of viral agents is crucial for both outbreak control and estimating regional disease burdens. Recently developed metagenomic methods have proven to be powerful tools for simultaneous pathogen detection. Here, we performed a systematic study of the capability of the short-read-based metagenomic approach in the molecular detection of viral pathogens in nasopharyngeal aspirate samples from patients with acute lower respiratory tract infections (n = 16). Using the high-throughput capacity of ultradeep sequencing and a dedicated data interpretation method, we successfully identified seven species of known respiratory viral agents from 15 samples, a result that was consistent with results of conventional PCR assays. We also detected a coinfected case that was missed by regular PCR testing. Using the metagenomic data, 11 draft genomes of the abundantly detected viruses in the samples were reconstructed with 21.84% to 98.53% coverage. Our results show the power of the short-read-based metagenomic approach for accurate and parallel screening of viral pathogens. Although there are some inherent difficulties in applying this approach to clinical samples, including a lack of controls, limited specimen quantity, and high contamination rate, our work will facilitate further application of this unprecedented high-throughput method to clinical samples.

Genetic Characteristics of Human Parainfluenza Virus Types 1-4 From Patients With Clinical Respiratory Tract Infection in China.

Human parainfluenza viruses (HPIV1-4) cause acute respiratory tract infections, thereby impacting human health worldwide. However, there are no current effective antivirals or licensed vaccines for infection prevention. Moreover, sequence information for human parainfluenza viruses (HPIVs) circulating in China is inadequate. Therefore, to shed light on viral genetic diversity and evolution, we collected samples from patients infected with HPIV1-4 in China from 2012 to 2018 to sequence the viruses. We obtained 24 consensus sequences, comprising 1 for HPIV1, 2 for HPIV2, 19 for HPIV3, and 2 for HPIV4A. Phylogenetic analyses classified the 1 HPIV1 into clade 2, and the 2 HPIV4 sequences into cluster 4A. Based on the hemagglutinin-neuraminidase (HN) gene, a new sub-cluster was identified in one of the HPIV2, namely G1c, and the 19 HPIV3 sequences were classified into the genetic lineages of C3f and C3a. The results indicated that HPIV1-4 were co-circulated in China. Further, the lineages of sub-cluster C3 of HPIV3 were co-circulated in China. A recombination analysis indicated that a putative recombination event may have occurred in the HN gene of HPIV3. In the obtained sequences of HPIV3, we found that two amino acid substitution sites (R73K in the F protein of PUMCH14028/2014 and A281V in the HN protein of PUMCH13961/2014) and a negative selection site (amino acid position 398 in the F protein) corresponded to the previously reported neutralization-related sites. Moreover, amino acid substitution site (K108E) corresponded to the negative selection site (amino acid position 108) in the 10 F proteins of HPIV3. However, no amino acid substitution site corresponded to the glycosylation site in the obtained HPIV3 sequences. These results might help in studying virus evolution, developing vaccines, and monitoring HPIV-related respiratory diseases.

VFDB 2008 release: an enhanced web-based resource for comparative pathogenomics.

Virulence factor database (VFDB) was set up in 2004 dedicated for providing current knowledge of virulence factors (VFs) from various medical significant bacterial pathogens to facilitate pathogenomic research. Nowadays, complete genome sequences of almost all the major pathogenic microbes have been determined, which makes comparative genomics a powerful approach for uncovering novel virulence determinants and hidden aspects of pathogenesis. VFDB was therefore upgraded to present the enormous diversity of bacterial genomes in terms of virulence genes and their organization. The VFDB 2008 release includes the following new features; (i) detailed tabular comparison of virulence composition of a given genome with other genomes of the same genus, (ii) multiple alignments and statistical analysis of homologous VFs and (iii) graphical comparison of genomic organizations of virulence genes. Comparative analysis of the numerous VFs will improve our understanding of the nature and evolution of virulence, as well as the development of new therapeutic and preventive strategies. VFDB 2008 release offers more user-friendly tools for comparative pathogenomics and it is publicly accessible at http://www.mgc.ac.cn/VFs/.

Metagenomic analysis of the lung microbiome in pulmonary tuberculosis - a pilot study.

The lung microbiome plays an important role in the pathophysiological processes associated with pulmonary tuberculosis (PTB). However, only a few studies using 16S rDNA amplicon sequencing have been reported, and the interactions between Mycobacterium tuberculosis (MTB) and the lung microbiome remain poorly understood. Patients with respiratory symptoms and imaging abnormalities compatible with tuberculosis (TB) were enrolled. We analyzed the lung microbiome in bronchoalveolar lavage (BAL) samples from 30 MTB-positive (MTB+) subjects and 30 MTB negative (MTB-) subjects by shotgun metagenomic sequencing. Alpha diversity tended to be lower in the MTB+ group than in the MTB- group. There was a significant difference in beta diversity between the MTB+ and MTB- subjects. MTB+ lung samples were dominated by MTB, while MTB- samples were enriched with Streptococcus, Prevotella, Nesseria, Selenomonas and Bifidobacterium, which more closely resemble the microbial composition of a healthy lung. Network analysis suggested that MTB could greatly impact the microbial community structure. MTB+ and MTB- communities showed distinct functional signatures. Fungal communities were also found to be associated with the presence or absence of MTB. Furthermore, it was confirmed that 16S rDNA amplicon sequencing underrepresents Mycobacterium. This pilot study is the first to explore the interplay between MTB and the host microbiome by using metagenomic sequencing. MTB dominates the lung microbiome of MTB+ subjects, while MTB- subjects have a Streptococcus-enriched microbiome. The 16S approach underrepresents Mycobacterium and is not the best way to study the TB-associated microbiome.

An Optimized Metagenomic Approach for Virome Detection of Clinical Pharyngeal Samples With Respiratory Infection.

Respiratory virus infections are one of the major causes of acute respiratory disease or exacerbation of chronic obstructive pulmonary disease (COPD). However, next-generation sequencing has not been used for routine viral detection in clinical respiratory samples owing to its sophisticated technology. Here, several pharyngeal samples with COPD were collected to enrich viral particles using an optimized method (M3), which involved M1 with centrifugation, filtration, and concentration, M2 (magnetic beads) combined with mixed nuclease digestion, and M4 with no pretreatment as a control. Metagenomic sequencing and bioinformatics analyses showed that the M3 method for viral enrichment was superior in both viral sequencing composition and viral taxa when compared to M1, M2, and M4. M3 acquired the most viral reads and more complete sequences within 15-h performance, indicating that it might be feasible for viral detection in multiple respiratory samples in clinical practice. Based on sequence similarity analysis, 12 human viruses, including nine Anelloviruses and three coronaviruses, were characterized. Coronavirus OC43 with the largest number of viral reads accounted for nearly complete (99.8%) genome sequences, indicating that it may be a major viral pathogen involved in exacerbation of COPD.

Distinct lung microbial community states in patients with pulmonary tuberculosis.

An improved understanding of the lung microbiome may lead to better strategies to diagnose, treat, and prevent pulmonary tuberculosis (PTB). However, the characteristics of the lung microbiomes of patients with TB remain largely undefined. In this study, 163 bronchoalveolar lavage (BAL) samples were collected from 163 sputum-negative suspected PTB patients. Furthermore, 12 paired BAL samples were obtained from 12 Mycobacterium tuberculosis-positive (MTB+) patients before and after negative conversion following a two-month anti-TB treatment. The V3-V4 region of the 16S ribosomal RNA (rRNA) gene was used to characterize the microbial composition of the lungs. The results showed that the prevalence of MTB in the BAL samples was 42.9% (70/163) among the sputum-negative patients. The α-diversity of lung microbiota was significantly less diverse in MTB+ patients compared with Mycobacterium tuberculosis-negative (MTB-) patients. There was a significant difference in ß-diversity between MTB+ and MTB- patients. MTB+ patients were enriched with Anoxybacillus, while MTB- patients were enriched with Prevotella, Alloprevotella, Veillonella, and Gemella. There was no significant difference between the Anoxybacillus detection rates of MTB+ and MTB- patients. The paired comparison between the BAL samples from MTB+ patients and their negative conversion showed that BAL negative-conversion microbiota had a higher α-diversity. In conclusion, distinct features of airway microbiota could be identified between samples from patients with and without MTB. Our results imply links between lung microbiota and different clinical groups of active PTB.

Neoadjuvant PD-1 Immune Checkpoint Blockade Reverses Functional Immunodominance among Tumor Antigen-Specific T Cells.

PURPOSE: Surgical resection of primary tumor with regional lymphadenectomy remains the treatment of choice for patients with advanced human papillomavirus-negative head and neck squamous cell carcinoma. However, even when pathologic disease-free margins can be achieved, locoregional and/or distant disease relapse remains high. Perioperative immunotherapy may improve outcomes, but mechanistic data supporting the use of neoadjuvant or adjuvant treatment clinically are sparse. EXPERIMENTAL DESIGN: Two syngeneic models of oral cavity carcinoma with defined T-cell antigens were treated with programmed death receptor 1 (PD-1) mAb before or after surgical resection of primary tumors, and antigen-specific T-cell responses were explored with functional and in vivo challenge assays. RESULTS: We demonstrated that functional immunodominance developed among T cells targeting multiple independent tumor antigens. T cells specific for subdominant antigens expressed greater levels of PD-1. Neoadjuvant, but not adjuvant, PD-1 immune checkpoint blockade broke immunodominance and induced T-cell responses to dominant and subdominant antigens. Using tumors lacking the immunodominant antigen as a model of antigen escape, neoadjuvant PD-1 immune checkpoint blockade induced effector T-cell immunity against tumor cells lacking immunodominant but retaining subdominant antigen. When combined with complete surgical excision, neoadjuvant PD-1 immune checkpoint blockade led to formation of immunologic memory capable of preventing engraftment of tumors lacking the immunodominant but retaining subdominant antigen. CONCLUSIONS: Together, these results implicate PD-1 expression by T cells in the mechanism of functional immunodominance among independent T-cell clones within a progressing tumor and support the use of neoadjuvant PD-1 immune checkpoint blockade in patients with surgically resectable carcinomas.

Complete genome sequence of the extremophilic Bacillus cereus strain Q1 with industrial applications.

Bacillus cereus strain Q1 was isolated from a deep-subsurface oil reservoir in the Daqing oil field in northeastern China. This strain is able to produce biosurfactants and to survive in extreme environments. Here we report the finished and annotated genome sequence of this organism.