Digestive system is a potential route of COVID-19: an analysis of single-cell coexpression pattern of key proteins in viral entry process

ObjectiveSince December 2019, a newly identified coronavirus (severe acute respiratory syndrome coronavirus (SARS-CoV-2)) has caused outbreaks of pneumonia in Wuhan, China. SARS-CoV-2 enters host cells via cell receptor ACE II (ACE2) and the transmembrane serine protease 2 (TMPRSS2). In order to identify possible prime target cells of SARS-CoV-2 by comprehensive dissection of ACE2 and TMPRSS2 coexpression pattern in different cell types, five datasets with single-cell transcriptomes of lung, oesophagus, gastric mucosa, ileum and colon were analysed.DesignFive datasets were searched, separately integrated and analysed. Violin plot was used to show the distribution of differentially expressed genes for different clusters. The ACE2-expressing and TMPRRSS2-expressing cells were highlighted and dissected to characterise the composition and proportion.ResultsCell types in each dataset were identified by known markers. ACE2 and TMPRSS2 were not only coexpressed in lung AT2 cells and oesophageal upper epithelial and gland cells but also highly expressed in absorptive enterocytes from the ileum and colon. Additionally, among all the coexpressing cells in the normal digestive system and lung, the expression of ACE2 was relatively highly expressed in the ileum and colon.ConclusionThis study provides the evidence of the potential route of SARS-CoV-2 in the digestive system along with the respiratory tract based on single-cell transcriptomic analysis. This finding may have a significant impact on health policy setting regarding the prevention of SARS-CoV-2 infection. Our study also demonstrates a novel method to identify the prime cell types of a virus by the coexpression pattern analysis of single-cell sequencing data.

Genomewide Expression and Functional Interactions of Genes under Drought Stress in Maize

A genomewide transcriptome assay of two subtropical genotypes of maize was used to observe the expression of genes at seedling stage of drought stress. The number of genes expressed differentially was greater in HKI1532 (a drought tolerant genotype) than in PC3 (a drought sensitive genotype), indicating primary differences at the transcriptional level in stress tolerance. The global coexpression networks of the two genotypes differed significantly with respect to the number of modules and the coexpression pattern within the modules. A total of 174 drought-responsive genes were selected from HKI1532, and their coexpression network revealed key correlations between different adaptive pathways, each cluster of the network representing a specific biological function. Transcription factors related to ABA-dependent stomatal closure, signalling, and phosphoprotein cascades work in concert to compensate for reduced photosynthesis. Under stress, water balance was maintained by coexpression of the genes involved in osmotic adjustments and transporter proteins. Metabolism was maintained by the coexpression of genes involved in cell wall modification and protein and lipid metabolism. The interaction of genes involved in crucial biological functions during stress was identified and the results will be useful in targeting important gene interactions to understand drought tolerance in greater detail.

Coexpression Pattern Analysis of NPM1-Associated Genes in Chronic Myelogenous Leukemia

Background. Nucleophosmin 1 (NPM1) plays an important role in ribosomal synthesis and malignancies, but NPM1 mutations occur rarely in the blast-crisis and chronic-phase chronic myelogenous leukemia (CML) patients. The NPM1-associated gene set (GCM_NPM1), in total 116 genes including NPM1, was chosen as the candidate gene set for the coexpression analysis. We wonder if NPM1-associated genes can affect the ribosomal synthesis and translation process in CML.Results. We presented a distribution-based approach for gene pair classification by identifying a disease-specific cutoff point that classified the coexpressed gene pairs into strong and weak coexpression structures. The differences in the coexpression patterns between the normal and the CML groups were reflected from the overall structure by performing two-sample Kolmogorov-Smirnov test. Our developed method effectively identified the coexpression pattern differences from the overall structure:P value=1.71×10-22<0.05for the maximum deviationD=0.109. Moreover, we found that genes involved in the ribosomal synthesis and translation process tended to be coexpressed in the CML group.Conclusion. Our developed method can identify the coexpression difference between two different groups. Dysregulation of ribosomal synthesis and translation process may be related to the CML disease. Our significant findings may provide useful information for the novel CML mechanism exploration and cancer treatment.

Stratification of Gene Coexpression Patterns and GO Function Mining for a RNA-Seq Data Series

RNA-Seq is emerging as an increasingly important tool in biological research, and it provides the most direct evidence of the relationship between the physiological state and molecular changes in cells. A large amount of RNA-Seq data across diverse experimental conditions have been generated and deposited in public databases. However, most developed approaches for coexpression analyses focus on the coexpression pattern mining of the transcriptome, thereby ignoring the magnitude of gene differences in one pattern. Furthermore, the functional relationships of genes in one pattern, and notably among patterns, were not always recognized. In this study, we developed an integrated strategy to identify differential coexpression patterns of genes and probed the functional mechanisms of the modules. Two real datasets were used to validate the method and allow comparisons with other methods. One of the datasets was selected to illustrate the flow of a typical analysis. In summary, we present an approach to robustly detect coexpression patterns in transcriptomes and to stratify patterns according to their relative differences. Furthermore, a global relationship between patterns and biological functions was constructed. In addition, a freely accessible web toolkit “coexpression pattern mining and GO functional analysis” (COGO) was developed.

A Major Lineage of Enteroendocrine Cells Coexpress CCK, Secretin, GIP, GLP-1, PYY, and Neurotensin but Not Somatostatin

Enteroendocrine cells such as duodenal cholecystokinin (CCK cells) are generally thought to be confined to certain segments of the gastrointestinal (GI) tract and to store and release peptides derived from only a single peptide precursor. In the current study, however, transgenic mice expressing enhanced green fluorescent protein (eGFP) under the control of the CCK promoter demonstrated a distribution pattern of CCK-eGFP positive cells that extended throughout the intestine. Quantitative PCR and liquid chromatography-mass spectrometry proteomic analyses of isolated, FACS-purified CCK-eGFP-positive cells demonstrated expression of not only CCK but also glucagon-like peptide 1 (GLP-1), gastric inhibitory peptide (GIP), peptide YY (PYY), neurotensin, and secretin, but not somatostatin. Immunohistochemistry confirmed this expression pattern. The broad coexpression phenomenon was observed both in crypts and villi as demonstrated by immunohistochemistry and FACS analysis of separated cell populations. Single-cell quantitative PCR indicated that approximately half of the duodenal CCK-eGFP cells express one peptide precursor in addition to CCK, whereas an additional smaller fraction expresses two peptide precursors in addition to CCK. The coexpression pattern was further confirmed through a cell ablation study based on expression of the human diphtheria toxin receptor under the control of the proglucagon promoter, in which activation of the receptor resulted in a marked reduction not only in GLP-1 cells, but also PYY, neurotensin, GIP, CCK, and secretin cells, whereas somatostatin cells were spared. Key elements of the coexpression pattern were confirmed by immunohistochemical double staining in human small intestine. It is concluded that a lineage of mature enteroendocrine cells have the ability to coexpress members of a group of functionally related peptides: CCK, secretin, GIP, GLP-1, PYY, and neurotensin, suggesting a potential therapeutic target for the treatment and prevention of diabetes and obesity.