Molecular determinants and heterogeneity underlying host response to EV-A71 infection at single-cell resolution.

The pathogenic human enterovirus EV-A71 has raised serious public health concerns. A hallmark of EV-A71 infection is the distortion of host transcriptomes in favour of viral replication. While high-throughput approaches have been exploited to dissect these gene dysregulations, they do not fully capture molecular perturbations at the single-cell level and in a physiologically relevant context. In this study, we applied a single-cell RNA sequencing approach on infected differentiated enterocyte cells (C2BBe1), which model the gastrointestinal epithelium targeted initially by EV-A71. Our single-cell analysis of EV-A71-infected culture provided several lines of illuminating observations: 1) This systems approach demonstrated extensive cell-to-cell variation in a single culture upon viral infection and delineated transcriptomic differences between the EV-A71-infected and bystander cells. 2) By analysing expression profiles of known EV-A71 receptors and entry facilitation factors, we found that ANXA2 was closely correlated in expression with the viral RNA in the infected population, supporting its role in EV-A71 entry in the enteric cells. 3) We further catalogued dysregulated lncRNAs elicited by EV-A71 infection and demonstrated the functional implication of lncRNA CYTOR in promoting EV-A71 replication. Viewed together, our single-cell transcriptomic analysis illustrated at the single-cell resolution the heterogeneity of host susceptibility to EV-A71 and revealed the involvement of lncRNAs in host antiviral response.

A novel antisense RNA ASPACT confers multi-level suppression of PACT and associated signalling.

The innate immune system is the frontline host protection against pathogens. Effective antiviral immunity is elicited upon recognition of viral RNAs by the host pattern recognition receptors. One of the major viral RNA sensors is retinoic acid inducible gene-1, which triggers the production of interferons (IFNs). In turn, this protective response requires another viral sensor and immunity factor interferon-inducible protein kinase RNA activator (PACT/PRKRA). Here, we report the identification and characterization of a novel antisense PACT gene that expresses a non-coding RNA in a convergent and interferon-inducible manner. Publicly available gene structure and expression data revealed that this gene, that we termed ASPACT, overlaps with the 3' -end of the PACT locus and is highly expressed during viral infection. Our results confirm the IFN-ß-inducibility of ASPACT, which is dependent on STAT-1/2. We further discovered that downregulation of ASPACT impacts both the expression and localization of the PACT transcript. At the transcription level, ChIP and ChIRP assays demonstrated that the ASPACT non-coding RNA occupies distinct chromatin regions of PACT gene and is important for promoter recruitment of the epigenetic silencer HDAC1. In parallel, ASPACT was also found to mediate nuclear retention of the PACT mRNA via direct RNA-RNA interaction, as revealed by RNA antisense purification assay. In summary, our results support the model that the non-coding RNA ASPACT acts as a negative regulator of PACT at multiple levels, and reveal a novel regulator of the viral counteractive response.

Isoform Evolution in Primates through Independent Combination of Alternative RNA Processing Events.

Recent RNA-seq technology revealed thousands of splicing events that are under rapid evolution in primates, whereas the reliability of these events, as well as their combination on the isoform level, have not been adequately addressed due to its limited sequencing length. Here, we performed comparative transcriptome analyses in human and rhesus macaque cerebellum using single molecule long-read sequencing (Iso-seq) and matched RNA-seq. Besides 359 million RNA-seq reads, 4,165,527 Iso-seq reads were generated with a mean length of 14,875 bp, covering 11,466 human genes, and 10,159 macaque genes. With Iso-seq data, we substantially expanded the repertoire of alternative RNA processing events in primates, and found that intron retention and alternative polyadenylation are surprisingly more prevalent in primates than previously estimated. We then investigated the combinatorial mode of these alternative events at the whole-transcript level, and found that the combination of these events is largely independent along the transcript, leading to thousands of novel isoforms missed by current annotations. Notably, these novel isoforms are selectively constrained in general, and 1,119 isoforms have even higher expression than the previously annotated major isoforms in human, indicating that the complexity of the human transcriptome is still significantly underestimated. Comparative transcriptome analysis further revealed 502 genes encoding selectively constrained, lineage-specific isoforms in human but not in rhesus macaque, linking them to some lineage-specific functions. Overall, we propose that the independent combination of alternative RNA processing events has contributed to complex isoform evolution in primates, which provides a new foundation for the study of phenotypic difference among primates.

Epigenetic silencing of myogenic gene program by Myb-binding protein 1a suppresses myogenesis.

Skeletal myogenesis involves highly coordinated steps that integrate developmental cues at the chromatin of muscle progenitors. Here, we identify Myb-binding protein 1a (Mybbp1a) as a novel negative regulator of muscle-specific gene expression and myoblast differentiation. The mode of action of Mybbp1a was linked to promoter regulation as illustrated by its interaction with MyoD at the genomic regions of silent muscle-specific genes as well as its negative effect on MyoD-mediated transcriptional activity. We propose that Mybbp1a exerts its repressive role by inducing a less permissible chromatin structure following recruitment of negative epigenetic modifiers such as HDAC1/2 and Suv39h1. At the onset of differentiation, Mybbp1a undergoes a promoter disengagement that may be due to the differentiation-responsive, miR-546-mediated downregulation of Mybbp1a expression. Moreover, such alteration gave rise to promoter enrichment of activators and histone acetylation, an epigenetic status amenable to gene activation. Together, these findings unveil a hitherto unrecognized transcriptional co-repressor role of Mybbp1a in proliferating muscle progenitor cells, and highlight an epigenetic mechanism by which Mybbp1a and miR-546 interplay to control myoblast differentiation transition.

Comprehensive functional genomic analyses link APC somatic mutation and mRNA-miRNA networks to the clinical outcome of stage-III colorectal cancer patients.

BACKGROUND: Colorectal cancer (CRC) is a major health concern globally, but exhibits regional and/or environmental distinctions in terms of outcome especially for patients with stage III CRC. METHODS: From 2014 to 2016, matched pairs of tumor and adjacent normal tissue samples from 60 patients with stage I-IV CRC from Chang Gung Memorial Hospital in Taiwan were analyzed using next-generation sequencing. The DNA, mRNA, and miRNA sequences of paired tumor tissues were profiled. An observational study with survival analysis was done. Online datasets of The Cancer Genome Atlas (TCGA) and The International Cancer Genome Consortium (ICGC) were also integrated and compared. RESULTS: The gene that exhibited the highest mutation rate was adenomatous polyposis coli (APC) (75.0%), followed by TP53 (70.0%), KRAS (56.6%), and TTN (48.3%). APC was also the most frequently mutated gene in TCGA and ICGC datasets. Surprisingly, for non-metastatic cases (stages I-III), CRC patients with mutated APC had better outcome in terms of overall survival (p = 0.041) and recurrence free survival (p = 0.0048). Particularly for stage III CRC, the overall survival rate was 94.4% and 67.7%, respectively (p = 0.018), and the recurrence free survival rate was 94.4% and 16.7%, respectively (p = 0.00044). Further clinical and gene expression analyses revealed that the APC wt specimens to a greater extent exhibit poor differentiation state as well as EGFR upregulation, providing molecular basis for the poor prognosis of these patients. Finally, based on integrated transcriptome analysis, we constructed the mRNA-miRNA networks underlying disease recurrence of the stage III CRC and uncovered potential therapeutic targets for this clinical condition. CONCLUSION: For stage III CRC, patients with mutated APC had better overall and recurrence free survival.

Defects in CISD-1, a mitochondrial iron-sulfur protein, lower glucose level and ATP production in Caenorhabditis elegans.

BACKGROUND: CDGSH iron sulfur domain-containing protein 1 (CISD-1) belongs to the CISD protein family that is evolutionary conserved across different species. In mammals, CISD-1 protein has been implicated in diseases such as cancers and diabetes. As a tractable model organism to study disease-associated proteins, we employed Caenorhabditis elegans in this study with an aim to establish a model for interrogating the functional relevance of CISD-1 in human metabolic conditions. METHODS: We first bioinformatically identified the human Cisd-1 homologue in worms. We then employed N2 wild-type and cisd-1(tm4993) mutant to investigate the consequences of CISD-1 loss-of-function on: 1) the expression pattern of CISD-1, 2) mitochondrial morphology pattern, 3) mitochondrial function and bioenergetics, and 4) the effects of anti-diabetes drugs. RESULTS: We first identified C. elegans W02B12.15 gene as the human Cisd-1 homologous gene, and pinpointed the localization of CISD-1 to the outer membrane of mitochondria. As compared with the N2 wild-type worm, cisd-1(tm4993) mutant exhibited a higher proportion of hyperfused form of mitochondria. This structural abnormality was associated with the generation of higher levels of ROS and mitochondrial superoxide but lower ATP. These physiological changes in mutants did not result in discernable effects on animal motility and lifespan. Moreover, the amount of glucose in N2 wild-type worms treated with troglitazone and pioglitazone, derivatives of TZD, was reduced to a comparable level as in the mutant animals. CONCLUSIONS: By focusing on the Cisd-1 gene, our study established a C. elegans genetic system suitable for modeling human diabetes-related diseases.