Single-cell RNA sequencing reveals peripheral blood leukocyte responses to spinal cord injury in mice with humanised immune systems.

Next-generation humanised mouse models and single-cell RNA sequencing (scRNAseq) approaches enable in-depth studies into human immune cell biology. Here we used NSG-SGM3 mice engrafted with human umbilical cord haematopoietic stem cells to investigate how human immune cells respond to and/or are changed by traumatic spinal cord injury (SCI). We hypothesised that the use of such mice could help advance our understanding of spinal cord injury-induced immune depression syndrome (SCI-IDS), and also how human leukocytes change as they migrate from the circulation into the lesion site. Our scRNAseq experiments, supplemented by flow cytometry, demonstrate the existence of up to 11 human immune cell (sub-) types and/or states across the blood and injured spinal cord (7 days post-SCI) of humanised NSG-SGM3 mice. Further comparisons of human immune cell transcriptomes between naïve, sham-operated and SCI mice identified a total of 579 differentially expressed genes, 190 of which were 'SCI-specific' (that is, genes regulated only in response to SCI but not sham surgery). Gene ontology analysis showed a prominent downregulation of immune cell function under SCI conditions, including for T cell receptor signalling and antigen presentation, confirming the presence of SCI-IDS and the transcriptional signature of human leukocytes in association with this phenomenon. We also highlight the activating influence of the local spinal cord lesion microenvironment by comparing the transcriptomes of circulating versus infiltrated human immune cells; those isolated from the lesion site were enriched for genes relating to both immune cell activity and function (e.g., oxidative phosphorylation, T cell proliferation and antigen presentation). We lastly applied an integrated bioinformatics approach to determine where immune responses in humanised NSG-SGM3 mice appear congruent to the native responses of human SCI patients, and where they diverge. Collectively, our study provides a valuable resource and methodological framework for the use of these mice in translational research.

Clinicopathologic Characteristics of Thyroid Microcarcinoma: Findings from a Hospital-Based Study in Vietnam.

BACKGROUND: Thyroid microcarcinoma (TMC) incidence has significantly increased in recent decades. The rates of lymph node metastasis extrathyroidal extension have been significantly different in patients with TMC ≤5 mm versus those with size >5 mm. The current analysis aimed to examine the clinicopathologic features of TMC measuring <5 mm and to compare them with those of TMC ≥5 mm. METHODS: A total of 273 patients with TMC confirmed by histological examination from December 2020 to May 2021 were enrolled in Bach Mai Hospital, Hanoi, Vietnam. Unconditional logistic regression models were used to determine the association between clinicopathological factors and tumor size, central lymph node metastasis and extrathyroidal extension. RESULTS: We found 212/273 patients (77.7%) were diagnosed incidentally. The majority of patients were female (87.5%) and had a mean age of 44.2 years. The mean tumor size (±standard deviation (SD)) was 5.72 ± 2.33 mm. Most of the patients were also diagnosed with papillary TMC. Multifocal and bilateral lesions accounted for 13.2% and 12.1%, respectively. The extrathyroidal invasion was observed in 14.7% (40 patients), while 24.5% (67 patients) were those with central lymph node metastases. The rate of extrathyroidal extension in patients with tumor size ≥5 mm was significantly higher than in patients with tumor size <5 mm (odds ratio (OR) = 4.98; 95% confidence interval (CI): 1.48-16.70; p = 0.004). Patients with body mass index (BMI) <23 kg/m2 were found to be protected against the odds of extrathyroidal extension (OR = 0.38, 95% CI: 0.19-0.75; p = 0.004) compared to those with BMI ≥23 kg/m2. In univariable mode, central lymph node metastasis was positively associated with the odds of the presence of extrathyroidal extension (OR = 2.70, 95% CI: 1.34-5.45; p = 0.004). In the multivariable model, central lymph node metastasis was also associated with the presence of extrathyroidal extension (OR = 2.507, 95% CI: 1.194-5.264; p = 0.017). Univariate analysis demonstrated that tumor size ≥5 mm (OR = 2.04; 95% CI: 1.01-4.17; p = 0.047) and extrathyroidal extension (OR = 2.71; 95% CI: 1.34-5.45; p = 0.004) were risk factors of central cervical lymph node metastasis. In multivariable models, the extrathyroidal extension was associated with central lymph metastasis. CONCLUSIONS: TMC <5 mm tumor size is less likely to have aggressive characteristics, including extrathyroidal extension, than a TMC ≥5 mm. Long-term follow-up studies are thus warranted to investigate the factors in the prognosis of TMC.

Robust mapping of spatiotemporal trajectories and cell-cell interactions in healthy and diseased tissues.

Spatial transcriptomics (ST) technologies generate multiple data types from biological samples, namely gene expression, physical distance between data points, and/or tissue morphology. Here we developed three computational-statistical algorithms that integrate all three data types to advance understanding of cellular processes. First, we present a spatial graph-based method, pseudo-time-space (PSTS), to model and uncover relationships between transcriptional states of cells across tissues undergoing dynamic change (e.g. neurodevelopment, brain injury and/or microglia activation, and cancer progression). We further developed a spatially-constrained two-level permutation (SCTP) test to study cell-cell interaction, finding highly interactive tissue regions across thousands of ligand-receptor pairs with markedly reduced false discovery rates. Finally, we present a spatial graph-based imputation method with neural network (stSME), to correct for technical noise/dropout and increase ST data coverage. Together, the algorithms that we developed, implemented in the comprehensive and fast stLearn software, allow for robust interrogation of biological processes within healthy and diseased tissues.

The Association between Peptic Ulcer Disease and Gastric Cancer: Results from the Stomach Cancer Pooling (StoP) Project Consortium.

BACKGROUND: Gastric cancer (GC) is the fifth most common type of cancer and the fourth most common cause of cancer-related mortality. Although the risk of GC and peptic ulcer disease (PUD) is known to be increased by H. pylori infection, evidence regarding the direct relationship between PUD and GC across ethnicities is inconclusive. Therefore, we investigated the association between PUD and GC in the Stomach cancer Pooling (StoP) consortium. METHODS: History of peptic ulcer disease was collected using a structured questionnaire in 11 studies in the StoP consortium, including 4106 GC cases and 6922 controls. The two-stage individual-participant data meta-analysis approach was adopted to generate a priori. Unconditional logistic regression and Firth's penalized maximum likelihood estimator were used to calculate study-specific odds ratios (ORs) and 95% confidence intervals (CIs) for the association between gastric ulcer (GU)/duodenal ulcer (DU) and risk of GC. RESULTS: History of GU and DU was thoroughly reported and used in association analysis, respectively, by 487 cases (12.5%) and 276 controls (4.1%), and 253 cases (7.8%) and 318 controls (6.0%). We found that GU was associated with an increased risk of GC (OR = 3.04, 95% CI: 2.07-4.49). No association between DU and GC risk was observed (OR = 1.03, 95% CI: 0.77-1.39). CONCLUSIONS: In the pooled analysis of 11 case-control studies in a large consortium (i.e., the Stomach cancer Pooling (StoP) consortium), we found a positive association between GU and risk of GC and no association between DU and GC risk.

Safety and immunogenicity of Nanocovax, a SARS-CoV-2 recombinant spike protein vaccine: Interim results of a double-blind, randomised controlled phase 1 and 2 trial.

Background: Nanocovax is a recombinant severe acute respiratory syndrome coronavirus 2 subunit vaccine composed of full-length prefusion stabilized recombinant SARS-CoV-2 spike glycoproteins (S-2P) and aluminium hydroxide adjuvant. Methods: We conducted a dose-escalation, open label trial (phase 1) and a randomized, double-blind, placebo-controlled trial (phase 2) to evaluate the safety and immunogenicity of the Nanocovax vaccine (in 25 mcg, 50 mcg, and 75 mcg doses, aluminium hydroxide adjuvanted (0·5 mg/dose) in 2-dose regime, 28 days apart (ClinicalTrials.gov number, NCT04683484). In phase 1, 60 participants received two intramuscular injection of the vaccine following dose-escalation procedure. The primary outcomes were reactogenicity and laboratory tests to evaluate the vaccine safety. In phase 2, 560 healthy adults received either vaccine doses similar in phase 1 (25 or 50 or 75 mcg S antigen in 0·5 mg aluminium per dose) or adjuvant (0·5 mg aluminium) in a ratio of 2:2:2:1. One primary outcome was the vaccine safety, including solicited adverse events for 7 day and unsolicited adverse events for 28 days after each injection as well as serious adverse event or adverse events of special interest throughout the study period. Another primary outcome was anti-S IgG antibody response (Index unit/ml). Secondary outcomes were surrogate virus neutralisation (inhibition percentage), wild-type SARS-CoV-2 neutralisation (dilution fold), and T-cell responses by intracellular staining for interferon gamma (IFNg). Anti-S IgG and neutralising antibody levels were compared with convalescent serum samples from symptomatic Covid-19 patients. Findings: For phase 1 study, no serious adverse events were observed for all 60 participants. Most adverse events were grade 1 and disappeared shortly after injection. For phase 2 study, after randomisation, 480 participants were assigned to receive the vaccine with adjuvant, and 80 participants were assigned to receive the placebo (adjuvant only). Reactogenicity was absent or mild in the majority of participants and of short duration (mean ≤3 days). Unsolicited adverse events were mild in most participants. There were no serious adverse events related to Nanocovax. Regarding the immunogenicity, Nanocovax induced robust anti-S antibody responses. In general, there humoral responses were similar among vaccine groups which reached their peaks at day 42 and declined afterward. At day 42, IgG levels of vaccine groups were 60·48 [CI95%: 51·12-71·55], 49·11 [41·26-58·46], 57·18 [48·4-67·5] compared to 7·10 [6·32-13·92] of convalescent samples. IgG levels reported here can be converted to WHO international standard binding antibody unit (BAU/ml) by multiplying them to a conversion factor of 21·8. Neutralising antibody titre of vaccine groups at day 42 were 89·2 [52·2-152·3], 80·0 [50·8-125.9] and 95·1 [63·1-143·6], compared to 55·1 [33·4-91·0] of the convalescent group. Interpretation: Up to day 90, Nanocovax was found to be safe, well tolerated, and induced robust immune responses. Funding: This work was funded by the Coalition for Epidemic Preparedness Innovations (CEPI), the Ministry of Science and Technology of Vietnam, and Nanogen Pharmaceutical Biotechnology JSC.

A model of impaired Langerhans cell maturation associated with HPV induced epithelial hyperplasia.

Langerhans cells (LC) are skin-resident antigen-presenting cells that regulate immune responses to epithelial microorganisms. Human papillomavirus (HPV) infection can promote malignant epithelial transformation. As LCs are considered important for controlling HPV infection, we compared the transcriptome of murine LCs from skin transformed by K14E7 oncoprotein and from healthy skin. We identified transcriptome heterogeneity at the single cell level amongst LCs in normal skin, associated with ontogeny, cell cycle, and maturation. We identified a balanced co-existence of immune-stimulatory and immune-inhibitory LC cell states in normal skin that was significantly disturbed in HPV16 E7-transformed skin. Hyperplastic skin was depleted of immune-stimulatory LCs and enriched for LCs with an immune-inhibitory gene signature, and LC-keratinocyte crosstalk was dysregulated. We identified reduced expression of interleukin (IL)-34, a critical molecule for LC homeostasis. Enrichment of an immune-inhibitory LC gene signature and reduced levels of epithelial IL-34 were also found in human HPV-associated cervical epithelial cancers.

Genome-wide analysis of SARS-CoV-2 strains circulating in Vietnam: Understanding the nature of the epidemic and role of the D614G mutation.

Genome-wide analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains is essential to better understand infectivity and virulence and to track coronavirus disease 2019 (COVID-19) cases and outbreaks. We performed whole-genome sequencing of 27 SARS-CoV-2 strains isolated between January 2020 and April 2020. A total of 54 mutations in different genomic regions was found. The D614G mutation, first detected in March 2020, was identified in 18 strains and was more likely associated with a lower cycle threshold (<25) in real-time reverse-transcription polymerase chain reaction diagnostic tests than the original D614 (prevalence ratio = 2.75; 95% confidence interval, 1.19-6.38). The integration of sequencing and epidemiological data suggests that SARS-CoV-2 transmission in both quarantine areas and in the community in Vietnam occur at the beginning of the epidemic although the country implemented strict quarantine quite early, with strict contact tracing, and testing. These findings provide insights into the nature of the epidemic, as well as shape strategies for COVID-19 prevention and control in Vietnam.

Survival after out-of-hospital cardiac arrest, Viet Nam: multicentre prospective cohort study.

OBJECTIVE: To investigate factors associated with survival after out-of-hospital cardiac arrest in Viet Nam. METHODS: We did a multicentre prospective observational study of people (> 18 years) presenting with out-of-hospital cardiac arrest (not caused by trauma) to three tertiary hospitals in Viet Nam from February 2014 to December 2018. We collected data on characteristics, management and outcomes of patients with out-of-hospital cardiac arrest and compared these data by type of transportation to hospital and survival to hospital admission. We assessed factors associated with survival to admission to and discharge from hospital using logistic regression analysis. FINDINGS: Of 590 eligible people with out-of-hospital cardiac arrest, 440 (74.6%) were male and the mean age was 56.1 years (standard deviation: 17.2). Only 24.2% (143/590) of these people survived to hospital admission and 14.1% (83/590) survived to hospital discharge. Most cardiac arrests (67.8%; 400/590) occurred at home, 79.4% (444/559) were witnessed by bystanders and 22.3% (124/555) were given cardiopulmonary resuscitation by a bystander. Only 8.6% (51/590) of the people were taken to hospital by the emergency medical services and 32.2% (49/152) received pre-hospital defibrillation. Pre-hospital defibrillation (odds ratio, OR: 3.90; 95% confidence interval, CI: 1.54-9.90) and return of spontaneous circulation in the emergency department (OR: 2.89; 95% CI: 1.03-8.12) were associated with survival to hospital admission. Hypothermia therapy during post-resuscitation care was associated with survival to discharge (OR: 5.44; 95% CI: 2.33-12.74). CONCLUSION: Improvements are needed in the emergency medical services in Viet Nam such as increasing bystander cardiopulmonary resuscitation and public access defibrillation, and improving ambulance and post-resuscitation care.

Prognostic value of early leukocyte fluctuations for recovery from traumatic spinal cord injury.

BACKGROUND: Acute traumatic spinal cord injury (SCI) induces a systemic immune response involving circulating white blood cells (WBCs). How this response is influenced by overall trauma severity, the neurological level of injury and/or correlates with patient outcomes is poorly understood. The objective of this study was to identify relationships between early changes in circulating WBCs, injury characteristics and long-term patient outcomes in individuals with traumatic SCI. METHODS: We retrospectively analysed data from 161 SCI patients admitted to Brisbane's Princess Alexandra Hospital (exploration cohort). Logistic regression models in conjunction with receiver operating characteristic (ROC) analyses were used to assess the strength of specific links between the WBC response, respiratory infection incidence and neurological outcomes (American Spinal Injury Association Impairment Scale (AIS) grade conversion). An independent validation cohort from the Trauma Hospital Berlin, Germany (n = 49) was then probed to assess the robustness of effects and disentangle centre effects. RESULTS: We find that the extent of acute neutrophilia in human SCI patients is positively correlated with New Injury Severity Scores but inversely with the neurological outcome (AIS grade). Multivariate analysis demonstrated that acute SCI-induced neutrophilia is an independent predictor of AIS grade conversion failure, with an odds ratio (OR) of 4.16 and ROC area under curve (AUC) of 0.82 (P < 0.0001). SCI-induced lymphopenia was separately identified as an independent predictor of better recovery (OR = 24.15; ROC AUC = 0.85, P < 0.0001). Acute neutrophilia and increased neutrophil-lymphocyte ratios were otherwise significantly associated with respiratory infection presentation in both patient cohorts. CONCLUSIONS: Our findings demonstrate the prognostic value of modelling early circulating neutrophil and lymphocyte counts with patient characteristics for predicting the longer term recovery after SCI.

Clinical features, isolation, and complete genome sequence of severe acute respiratory syndrome coronavirus 2 from the first two patients in Vietnam.

In January 2020, we identified two severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected patients in a familial cluster with one person coming from Wuhan, China. The complete genome sequences of two SARS-CoV-2 strains isolated from these patients were identical and 99.98% similar to strains isolated in Wuhan. This is genetically suggestive of human-to-human transmission of SARS-CoV-2 and indicates Wuhan as the most plausible origin of the early outbreak in Vietnam. The younger patient had a mild upper respiratory illness and a brief viral shedding, whereas the elderly with multi-morbidity had pneumonia, prolonged viral shedding, and residual lung damage. The evidence of nonsynonymous substitutions in the ORF1ab region of the viral sequence warrants further studies.

Repopulating Microglia Promote Brain Repair in an IL-6-Dependent Manner.

Cognitive dysfunction and reactive microglia are hallmarks of traumatic brain injury (TBI), yet whether these cells contribute to cognitive deficits and secondary inflammatory pathology remains poorly understood. Here, we show that removal of microglia from the mouse brain has little effect on the outcome of TBI, but inducing the turnover of these cells through either pharmacologic or genetic approaches can yield a neuroprotective microglial phenotype that profoundly aids recovery. The beneficial effects of these repopulating microglia are critically dependent on interleukin-6 (IL-6) trans-signaling via the soluble IL-6 receptor (IL-6R) and robustly support adult neurogenesis, specifically by augmenting the survival of newborn neurons that directly support cognitive function. We conclude that microglia in the mammalian brain can be manipulated to adopt a neuroprotective and pro-regenerative phenotype that can aid repair and alleviate the cognitive deficits arising from brain injury.

ascend: R package for analysis of single-cell RNA-seq data.

BACKGROUND: Recent developments in single-cell RNA sequencing (scRNA-seq) platforms have vastly increased the number of cells typically assayed in an experiment. Analysis of scRNA-seq data is multidisciplinary in nature, requiring careful consideration of the application of statistical methods with respect to the underlying biology. Few analysis packages exist that are at once robust, are computationally fast, and allow flexible integration with other bioinformatics tools and methods. FINDINGS: ascend is an R package comprising tools designed to simplify and streamline the preliminary analysis of scRNA-seq data, while addressing the statistical challenges of scRNA-seq analysis and enabling flexible integration with genomics packages and native R functions, including fast parallel computation and efficient memory management. The package incorporates both novel and established methods to provide a framework to perform cell and gene filtering, quality control, normalization, dimension reduction, clustering, differential expression, and a wide range of visualization functions. CONCLUSIONS: ascend is designed to work with scRNA-seq data generated by any high-throughput platform and includes functions to convert data objects between software packages. The ascend workflow is simple and interactive, as well as suitable for implementation by a broad range of users, including those with little programming experience.

Detection of HPV E7 Transcription at Single-Cell Resolution in Epidermis.

Persistent human papillomavirus (HPV) infection is responsible for at least 5% of human malignancies. Most HPV-associated cancers are initiated by the HPV16 genotype, as confirmed by detection of integrated HPV DNA in cells of oral and anogenital epithelial cancers. However, single-cell RNA sequencing may enable prediction of HPV involvement in carcinogenesis at other sites. We conducted single-cell RNA sequencing on keratinocytes from a mouse transgenic for the E7 gene of HPV16 and showed sensitive and specific detection of HPV16-E7 mRNA, predominantly in basal keratinocytes. We showed that increased E7 mRNA copy number per cell was associated with increased expression of E7 induced genes. This technique enhances detection of active viral transcription in solid tissue and may clarify possible linkage of HPV infection to development of squamous cell carcinoma.

Single-cell RNA-seq of human induced pluripotent stem cells reveals cellular heterogeneity and cell state transitions between subpopulations.

Heterogeneity of cell states represented in pluripotent cultures has not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC-CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method and, through this, identified four subpopulations distinguishable on the basis of their pluripotent state, including a core pluripotent population (48.3%), proliferative (47.8%), early primed for differentiation (2.8%), and late primed for differentiation (1.1%). For each subpopulation, we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets composed of 165 unique genes that denote the specific pluripotency states; using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to threefold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations and support our conclusions with results from two orthogonal pseudotime trajectory methods.

Mammalian genomic regulatory regions predicted by utilizing human genomics, transcriptomics, and epigenetics data.

Genome sequences for hundreds of mammalian species are available, but an understanding of their genomic regulatory regions, which control gene expression, is only beginning. A comprehensive prediction of potential active regulatory regions is necessary to functionally study the roles of the majority of genomic variants in evolution, domestication, and animal production. We developed a computational method to predict regulatory DNA sequences (promoters, enhancers, and transcription factor binding sites) in production animals (cows and pigs) and extended its broad applicability to other mammals. The method utilizes human regulatory features identified from thousands of tissues, cell lines, and experimental assays to find homologous regions that are conserved in sequences and genome organization and are enriched for regulatory elements in the genome sequences of other mammalian species. Importantly, we developed a filtering strategy, including a machine learning classification method, to utilize a very small number of species-specific experimental datasets available to select for the likely active regulatory regions. The method finds the optimal combination of sensitivity and accuracy to unbiasedly predict regulatory regions in mammalian species. Furthermore, we demonstrated the utility of the predicted regulatory datasets in cattle for prioritizing variants associated with multiple production and climate change adaptation traits and identifying potential genome editing targets.

Single cell RNA sequencing of stem cell-derived retinal ganglion cells.

We used single cell sequencing technology to characterize the transcriptomes of 1,174 human embryonic stem cell-derived retinal ganglion cells (RGCs) at the single cell level. The human embryonic stem cell line BRN3B-mCherry (A81-H7), was differentiated to RGCs using a guided differentiation approach. Cells were harvested at day 36 and prepared for single cell RNA sequencing. Our data indicates the presence of three distinct subpopulations of cells, with various degrees of maturity. One cluster of 288 cells showed increased expression of genes involved in axon guidance together with semaphorin interactions, cell-extracellular matrix interactions and ECM proteoglycans, suggestive of a more mature RGC phenotype.

The nematode Caenorhabditis elegans displays a chemotaxis behavior to tuberculosis-specific odorants.

A simple, affordable diagnostic test for pulmonary tuberculosis (TB) is urgently needed to improve detection of active Mycobacterium tuberculosis. Recently, it has been suggested that animal behavior can be used as a biosensor to signal the presence of human disease. For example, the giant African pouched rats can detect tuberculosis by sniffing sputum specimens while trained honeybees respond to three of the volatile organic compounds (VOCs) detected in the breath of TB positive patients by proboscis extension. However, both rats and honeybees require animal housing facilities and professional trainers, which are outside the scope of most disease testing facilities. Here, we report that the innate olfactory behavioral response of the roundworm nematode Caenorhabditis elegans can be used to detect the TB-specific VOCs methyl p-anisate, methyl nicotinate, methyl phenylacetate and o-phenylanisole, in chemotaxis assays. Dauer larvae, a long-lived stress resistant alternative development state of C. elegans in which the animals can survive for extended periods of time in dry conditions with no food, were also demonstrated to detect the VOCs. We propose that exposing naive dauer larvae to TB-related VOCs and recording their response in this behavioral assay could lead to the development of a new method for TB diagnostics using breath as the sample type.