Context-dependent gene regulatory network reveals regulation dynamics and cell trajectories using unspliced transcripts.

Gene regulatory networks govern complex gene expression programs in various biological phenomena, including embryonic development, cell fate decisions and oncogenesis. Single-cell techniques are increasingly being used to study gene expression, providing higher resolution than traditional approaches. However, inferring a comprehensive gene regulatory network across different cell types remains a challenge. Here, we propose to construct context-dependent gene regulatory networks (CDGRNs) from single-cell RNA sequencing data utilizing both spliced and unspliced transcript expression levels. A gene regulatory network is decomposed into subnetworks corresponding to different transcriptomic contexts. Each subnetwork comprises the consensus active regulation pairs of transcription factors and their target genes shared by a group of cells, inferred by a Gaussian mixture model. We find that the union of gene regulation pairs in all contexts is sufficient to reconstruct differentiation trajectories. Functions specific to the cell cycle, cell differentiation or tissue-specific functions are enriched throughout the developmental process in each context. Surprisingly, we also observe that the network entropy of CDGRNs decreases along differentiation trajectories, indicating directionality in differentiation. Overall, CDGRN allows us to establish the connection between gene regulation at the molecular level and cell differentiation at the macroscopic level.

Thymic macrophages consist of two populations with distinct localization and origin.

Tissue-resident macrophages are essential to protect from pathogen invasion and maintain organ homeostasis. The ability of thymic macrophages to engulf apoptotic thymocytes is well appreciated, but little is known about their ontogeny, maintenance, and diversity. Here, we characterized the surface phenotype and transcriptional profile of these cells and defined their expression signature. Thymic macrophages were most closely related to spleen red pulp macrophages and Kupffer cells and shared the expression of the transcription factor (TF) SpiC with these cells. Single-cell RNA sequencing (scRNA-Seq) showed that the macrophages in the adult thymus are composed of two populations distinguished by the expression of Timd4 and Cx3cr1. Remarkably, Timd4+ cells were located in the cortex, while Cx3cr1+ macrophages were restricted to the medulla and the cortico-medullary junction. Using shield chimeras, transplantation of embryonic thymuses, and genetic fate mapping, we found that the two populations have distinct origins. Timd4+ thymic macrophages are of embryonic origin, while Cx3cr1+ macrophages are derived from adult hematopoietic stem cells. Aging has a profound effect on the macrophages in the thymus. Timd4+ cells underwent gradual attrition, while Cx3cr1+ cells slowly accumulated with age and, in older mice, were the dominant macrophage population in the thymus. Altogether, our work defines the phenotype, origin, and diversity of thymic macrophages.

Mechanistic analysis of carbon-carbon bond formation by deoxypodophyllotoxin synthase.

Deoxypodophyllotoxin contains a core of four fused rings (A to D) with three consecutive chiral centers, the last being created by the attachment of a peripheral trimethoxyphenyl ring (E) to ring C. Previous studies have suggested that the iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, deoxypodophyllotoxin synthase (DPS), catalyzes the oxidative coupling of ring B and ring E to form ring C and complete the tetracyclic core. Despite recent efforts to deploy DPS in the preparation of deoxypodophyllotoxin analogs, the mechanism underlying the regio- and stereoselectivity of this cyclization event has not been elucidated. Herein, we report 1) two structures of DPS in complex with 2OG and (±)-yatein, 2) in vitro analysis of enzymatic reactivity with substrate analogs, and 3) model reactions addressing DPS's catalytic mechanism. The results disfavor a prior proposal of on-pathway benzylic hydroxylation. Rather, the DPS-catalyzed cyclization likely proceeds by hydrogen atom abstraction from C7', oxidation of the benzylic radical to a carbocation, Friedel-Crafts-like ring closure, and rearomatization of ring B by C6 deprotonation. This mechanism adds to the known pathways for transformation of the carbon-centered radical in Fe/2OG enzymes and suggests what types of substrate modification are likely tolerable in DPS-catalyzed production of deoxypodophyllotoxin analogs.

cgMLST@Taiwan: A web service platform for Vibrio cholerae cgMLST profiling and global strain tracking.

BACKGROUND: Cholera, a rapidly dehydrating diarrheal disease caused by toxigenic Vibrio cholerae, is a leading cause of morbidity and mortality in some regions of the world. Core genome multilocus sequence typing (cgMLST) is a promising approach in generating genetic fingerprints from whole-genome sequencing (WGS) data for strain comparison among laboratories. METHODS: We constructed a V. cholerae core gene allele database using an in-house developed computational pipeline, a database with cgMLST profiles converted from genomic sequences from the National Center for Biotechnology Information, and built a REST-based web accessible via the Internet. RESULTS: We built a web service platform-cgMLST@Taiwan and installed a V. cholerae allele database, a cgMLST profile database, and computational tools for generating V. cholerae cgMLST profiles (based on 3,017 core genes), performing rapid global strain tracking, and clustering analysis of cgMLST profiles. This web-based platform provides services to researchers, public health microbiologists, and physicians who use WGS data for the investigation of cholera outbreaks and tracking of V. cholerae strain transmission across countries and geographic regions. The cgMLST@Taiwan is accessible at http://rdvd.cdc.gov.tw/cgMLST.

Identification of cell states using super-enhancer RNA.

BACKGROUND: A new class of regulatory elements called super-enhancers, comprised of multiple neighboring enhancers, have recently been reported to be the key transcriptional drivers of cellular, developmental, and disease states. RESULTS: Here, we defined super-enhancer RNAs as highly expressed enhancer RNAs that are transcribed from a cluster of localized genomic regions. Using the cap analysis of gene expression sequencing data from FANTOM5, we systematically explored the enhancer and messenger RNA landscapes in hundreds of different cell types in response to various environments. Applying non-negative matrix factorization (NMF) to super-enhancer RNA profiles, we found that different cell types were well classified. In addition, through the NMF of individual time-course profiles from a single cell-type, super-enhancer RNAs were clustered into several states with progressive patterns. We further investigated the enriched biological functions of the proximal genes involved in each pattern, and found that they were associated with the corresponding developmental process. CONCLUSIONS: The proposed super-enhancer RNAs can act as a good alternative, without the complicated measurement of histone modifications, for identifying important regulatory elements of cell type specification and identifying dynamic cell states.

Emergence of Vibrio cholerae O1 Sequence Type 75 in Taiwan.

We investigated the epidemiology of cholera in Taiwan during 2002-2018. Vibrio cholerae sequence type (ST) 75 clone emerged in 2009 and has since become more prevalent than the ST69 clone from a previous pandemic. Closely related ST75 strains have emerged in 4 countries and may now be widespread in Asia.

Emergence of Multidrug-Resistant Salmonella enterica Serovar Goldcoast Strains in Taiwan and International Spread of the ST358 Clone.

Salmonella enterica serovar Goldcoast infection was rare in Taiwan; it was not detected in routine surveillance from 2004 to 2013. This serovar was first identified in 2014, but the frequency of infection remained low until 2017. From 2014 to 2016, all but one isolate was pan-susceptible. S Goldcoast infections abruptly increased in 2018, and all isolates were multidrug-resistant (MDR). All MDR isolates harbored an IncHI2 plasmid, and the majority carried 14 antimicrobial resistance genes, aac(3)-IId, aadA22, aph(3')-Ia, aph(6)-Id, blaTEM-1B, blaCTX-M-55, lnu(F), floR, qnrS13, arr-2, sul2, sul3, tet(A), and dfrA14. S Goldcoast strains recovered in Taiwan and 96 of 99 strains from Germany, the Netherlands, the United Kingdom, and the United States belonged to sequence type 358 (ST358). Whole-genome single-nucleotide polymorphism and core genome multilocus sequence type analyses revealed that all strains of the ST358 clone shared a high degree of genetic relatedness. The present study highlighted that a dramatic increase in S Goldcoast infections followed the emergence of MDR strains and indicated that a genetically closely related S Goldcoast ST358 clone may have widespread significance internationally.

Reaction Mechanism of a Nonheme Iron Enzyme Catalyzed Oxidative Cyclization via C-C Bond Formation.

A complementary study including design of mechanistic probes, biochemical assays, model analysis, and liquid chromatography coupled mass spectrometry was conducted to establish the reaction mechanism for a nonheme iron enzyme catalyzed (-)-podophyllotoxin formation. Our results indicate that the originally proposed hydroxylated intermediate is unlikely to be involved in this reaction. Instead, the formation of benzylic radical/carbocation intermediate can be utilized to trigger the C-C bond formation to construct the C-ring of (-)-podophyllotoxin.

New Multidrug-Resistant Salmonella enterica Serovar Anatum Clone, Taiwan, 2015-2017.

In 2011, a Salmonella enterica serovar Anatum clone emerged in Taiwan. During 2016-2017, infections increased dramatically, strongly associated with emergence and spread of multidrug-resistant strains with a plasmid carrying 11 resistance genes, including blaDHA-1. Because these resistant strains infect humans and food animals, control measures are urgently needed.

Genetic Relationships among Multidrug-Resistant Salmonella enterica Serovar Typhimurium Strains from Humans and Animals.

We identified 20 to 22 resistance genes, carried in four incompatibility groups of plasmids, in each of five genetically closely related Salmonella enterica serovar Typhimurium strains recovered from humans, pigs, and chickens. The genes conferred resistance to aminoglycosides, chloramphenicol, sulfonamides, trimethoprim, tetracycline, fluoroquinolones, extended-spectrum cephalosporins and cefoxitin, and azithromycin. This study demonstrates the transmission of multidrug-resistant Salmonella strains among humans and food animals and may be the first identification of mphA in azithromycin-resistant Salmonella strains in Taiwan.