Objective: Programmed death 1 and its ligand 1 (PD-1/PD-L1) immunotherapy is promising for late-stage lung cancer treatment, however, the response rate needs to be improved. Gut microbiota plays a crucial role in immunotherapy sensitisation and Panax ginseng has been shown to possess immunomodulatory potential. In this study, we aimed to investigate whether the combination treatment of ginseng polysaccharides (GPs) and αPD-1 monoclonal antibody (mAb) could sensitise the response by modulating gut microbiota. Design: Syngeneic mouse models were administered GPs and αPD-1 mAb, the sensitising antitumour effects of the combination therapy on gut microbiota were assessed by faecal microbiota transplantation (FMT) and 16S PacBio single-molecule real-time (SMRT) sequencing. To assess the immune-related metabolites, metabolomics analysis of the plasma samples was performed. Results: We found GPs increased the antitumour response to αPD-1 mAb by increasing the microbial metabolites valeric acid and decreasing L-kynurenine, as well as the ratio of Kyn/Trp, which contributed to the suppression of regulatory T cells and induction of T _eff cells after combination treatment. Besides, the microbial analysis indicated that the abundance of Parabacteroides distasonis and Bacteroides vulgatus was higher in responders to anti-PD-1 blockade than non-responders in the clinic. Furthermore, the combination therapy sensitised the response to PD-1 inhibitor in the mice receiving microbes by FMT from six non-responders by reshaping the gut microbiota from non-responders towards that of responders. Conclusion: Our results demonstrate that GPs combined with αPD-1 mAb may be a new strategy to sensitise non-small cell lung cancer patients to anti-PD-1 immunotherapy. The gut microbiota can be used as a novel biomarker to predict the response to anti-PD-1 immunotherapy.

Background: MicroRNAs are a family of ~22 nt small RNAs that can regulate gene expression at the post-transcriptional level. Identification of these molecules and their targets can aid understanding of regulatory processes. Recently, HTS has become a common identification method but there are two major limitations associated with the technique. Firstly, the method has low efficiency, with typically less than 1 in 10,000 sequences representing miRNA reads and secondly the method preferentially targets highly expressed miRNAs. If sequences are available, computational methods can provide a screening step to investigate the value of an HTS study and aid interpretation of results. However, current methods can only predict miRNAs for short fragments and have usually been trained against small datasets which don't always reflect the diversity of these molecules. Results: We have developed a software tool, miRPara, that predicts most probable mature miRNA coding regions from genome scale sequences in a species specific manner. We classified sequences from miRBase into animal, plant and overall categories and used a support vector machine to train three models based on an initial set of 77 parameters related to the physical properties of the pre-miRNA and its miRNAs. By applying parameter filtering we found a subset of ~25 parameters produced higher prediction ability compared to the full set. Our software achieves an accuracy of up to 80% against experimentally verified mature miRNAs, making it one of the most accurate methods available. Conclusions: miRPara is an effective tool for locating miRNAs coding regions in genome sequences and can be used as a screening step prior to HTS experiments. It is available at http://www.whiov.ac.cn/bioinformatics/mirpara

The spread of SARS-CoV-2 in Beijing before May, 2020 resulted from transmission following both domestic and global importation of cases. Here we present genomic surveillance data on 102 imported cases, which account for 17.2% of the total cases in Beijing. Our data suggest that all of the cases in Beijing can be broadly classified into one of three groups: Wuhan exposure, local transmission and overseas imports. We classify all sequenced genomes into seven clusters based on representative high-frequency single nucleotide polymorphisms (SNPs). Genomic comparisons reveal higher genomic diversity in the imported group compared to both the Wuhan exposure and local transmission groups, indicating continuous genomic evolution during global transmission. The imported group show region-specific SNPs, while the intra-host single nucleotide variations present as random features, and show no significant differences among groups. Epidemiological data suggest that detection of cases at immigration with mandatory quarantine may be an effective way to prevent recurring outbreaks triggered by imported cases. Notably, we also identify a set of novel indels. Our data imply that SARS-CoV-2 genomes may have high mutational tolerance.

Bats have been identified as a natural reservoir of a variety of coronaviruses (CoVs). Several of them have caused diseases in humans and domestic animals by interspecies transmission. Considering the diversity of bat coronaviruses, bat species and populations, we expect to discover more bat CoVs through virus surveillance. In this study, we described a new member of alphaCoV (BtCoV/Rh/YN2012) in bats with unique genome features. Unique accessory genes, ORF4a and ORF4b were found between the spike gene and the envelope gene, while ORF8 gene was found downstream of the nucleocapsid gene. All the putative genes were further confirmed by reverse-transcription analyses. One unique gene at the 3’ end of the BtCoV/Rh/YN2012 genome, ORF9, exhibits ~30% amino acid identity to ORF7a of the SARS-related coronavirus. Functional analysis showed ORF4a protein can activate IFN-β production, whereas ORF3a can regulate NF-κB production. We also screened the spike-mediated virus entry using the spike-pseudotyped retroviruses system, although failed to find any fully permissive cells. Our results expand the knowledge on the genetic diversity of bat coronaviruses. Continuous screening of bat viruses will help us further understand the important role played by bats in coronavirus evolution and transmission.

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

Since 2016, the highly pathogenic avian influenza H5N8 virus has emerged in the Central Asian flyway and Europe, causing massive deaths in poultry and wild birds. In this study, we isolated and ide...