Dissecting embryonic and extraembryonic lineage crosstalk with stem cell co-culture.

Embryogenesis necessitates harmonious coordination between embryonic and extraembryonic tissues. Although stem cells of both embryonic and extraembryonic origins have been generated, they are grown in different culture conditions. In this study, utilizing a unified culture condition that activates the FGF, TGF-ß, and WNT pathways, we have successfully derived embryonic stem cells (FTW-ESCs), extraembryonic endoderm stem cells (FTW-XENs), and trophoblast stem cells (FTW-TSCs) from the three foundational tissues of mouse and cynomolgus monkey (Macaca fascicularis) blastocysts. This approach facilitates the co-culture of embryonic and extraembryonic stem cells, revealing a growth inhibition effect exerted by extraembryonic endoderm cells on pluripotent cells, partially through extracellular matrix signaling. Additionally, our cross-species analysis identified both shared and unique transcription factors and pathways regulating FTW-XENs. The embryonic and extraembryonic stem cell co-culture strategy offers promising avenues for developing more faithful embryo models and devising more developmentally pertinent differentiation protocols.

Single-Cell Sequencing of Peripheral Mononuclear Cells Reveals Distinct Immune Response Landscapes of COVID-19 and Influenza Patients.

The coronavirus disease 2019 (COVID-19) pandemic poses a current world-wide public health threat. However, little is known about its hallmarks compared to other infectious diseases. Here, we report the single-cell transcriptional landscape of longitudinally collected peripheral blood mononuclear cells (PBMCs) in both COVID-19- and influenza A virus (IAV)-infected patients. We observed increase of plasma cells in both COVID-19 and IAV patients and XIAP associated factor 1 (XAF1)-, tumor necrosis factor (TNF)-, and FAS-induced T cell apoptosis in COVID-19 patients. Further analyses revealed distinct signaling pathways activated in COVID-19 (STAT1 and IRF3) versus IAV (STAT3 and NFκB) patients and substantial differences in the expression of key factors. These factors include relatively increase of interleukin (IL)6R and IL6ST expression in COVID-19 patients but similarly increased IL-6 concentrations compared to IAV patients, supporting the clinical observations of increased proinflammatory cytokines in COVID-19 patients. Thus, we provide the landscape of PBMCs and unveil distinct immune response pathways in COVID-19 and IAV patients.

Cell competition constitutes a barrier for interspecies chimerism.

Cell competition involves a conserved fitness-sensing process during which fitter cells eliminate neighbouring less-fit but viable cells1. Cell competition has been proposed as a surveillance mechanism to ensure normal development and tissue homeostasis, and has also been suggested to act as a barrier to interspecies chimerism2. However, cell competition has not been studied in an interspecies context during early development owing to the lack of an in vitro model. Here we developed an interspecies pluripotent stem cell (PSC) co-culture strategy and uncovered a previously unknown mode of cell competition between species. Interspecies competition between PSCs occurred in primed but not naive pluripotent cells, and between evolutionarily distant species. By comparative transcriptome analysis, we found that genes related to the NF-κB signalling pathway, among others, were upregulated in less-fit 'loser' human cells. Genetic inactivation of a core component (P65, also known as RELA) and an upstream regulator (MYD88) of the NF-κB complex in human cells could overcome the competition between human and mouse PSCs, thereby improving the survival and chimerism of human cells in early mouse embryos. These insights into cell competition pave the way for the study of evolutionarily conserved mechanisms that underlie competitive cell interactions during early mammalian development. Suppression of interspecies PSC competition may facilitate the generation of human tissues in animals.

Plasma cell-free RNA characteristics in COVID-19 patients.

The pathogenesis of COVID-19 is still elusive, which impedes disease progression prediction, differential diagnosis, and targeted therapy. Plasma cell-free RNAs (cfRNAs) carry unique information from human tissue and thus could point to resourceful solutions for pathogenesis and host-pathogen interactions. Here, we performed a comparative analysis of cfRNA profiles between COVID-19 patients and healthy donors using serial plasma. Analyses of the cfRNA landscape, potential gene regulatory mechanisms, dynamic changes in tRNA pools upon infection, and microbial communities were performed. A total of 380 cfRNA molecules were up-regulated in all COVID-19 patients, of which seven could serve as potential biomarkers (AUC > 0.85) with great sensitivity and specificity. Antiviral (NFKB1A, IFITM3, and IFI27) and neutrophil activation (S100A8, CD68, and CD63)-related genes exhibited decreased expression levels during treatment in COVID-19 patients, which is in accordance with the dynamically enhanced inflammatory response in COVID-19 patients. Noncoding RNAs, including some microRNAs (let 7 family) and long noncoding RNAs (GJA9-MYCBP) targeting interleukin (IL6/IL6R), were differentially expressed between COVID-19 patients and healthy donors, which accounts for the potential core mechanism of cytokine storm syndromes; the tRNA pools change significantly between the COVID-19 and healthy group, leading to the accumulation of SARS-CoV-2 biased codons, which facilitate SARS-CoV-2 replication. Finally, several pneumonia-related microorganisms were detected in the plasma of COVID-19 patients, raising the possibility of simultaneously monitoring immune response regulation and microbial communities using cfRNA analysis. This study fills the knowledge gap in the plasma cfRNA landscape of COVID-19 patients and offers insight into the potential mechanisms of cfRNAs to explain COVID-19 pathogenesis.

Plasma cell-free DNA promise monitoring and tissue injury assessment of COVID-19.

Coronavirus 2019 (COVID-19) is a complex disease that affects billions of people worldwide. Currently, effective etiological treatment of COVID-19 is still lacking; COVID-19 also causes damages to various organs that affects therapeutics and mortality of the patients. Surveillance of the treatment responses and organ injury assessment of COVID-19 patients are of high clinical value. In this study, we investigated the characteristic fragmentation patterns and explored the potential in tissue injury assessment of plasma cell-free DNA in COVID-19 patients. Through recruitment of 37 COVID-19 patients, 32 controls and analysis of 208 blood samples upon diagnosis and during treatment, we report gross abnormalities in cfDNA of COVID-19 patients, including elevated GC content, altered molecule size and end motif patterns. More importantly, such cfDNA fragmentation characteristics reflect patient-specific physiological changes during treatment. Further analysis on cfDNA tissue-of-origin tracing reveals frequent tissue injuries in COVID-19 patients, which is supported by clinical diagnoses. Hence, our work demonstrates and extends the translational merit of cfDNA fragmentation pattern as valuable analyte for effective treatment monitoring, as well as tissue injury assessment in COVID-19.

Noninvasive preeclampsia prediction using plasma cell-free RNA signatures.

BACKGROUND: Preeclampsia, especially preterm preeclampsia and early-onset preeclampsia, is a life-threating pregnancy disorder, and the heterogeneity and complexity of preeclampsia make it difficult to predict risk and to develop treatments. Plasma cell-free RNA carries unique information from human tissue and may be useful for noninvasive monitoring of maternal, placental, and fetal dynamics during pregnancy. OBJECTIVE: This study aimed to investigate various RNA biotypes associated with preeclampsia in plasma and to develop classifiers to predict preterm preeclampsia and early-onset preeclampsia before diagnosis. STUDY DESIGN: We performed a novel, cell-free RNA sequencing method termed polyadenylation ligation-mediated sequencing to investigate the cell-free RNA characteristics of 715 healthy pregnancies and 202 pregnancies affected by preeclampsia before symptom onset. We explored differences in the abundance of different RNA biotypes in plasma between healthy and preeclampsia samples and built preterm preeclampsia and early-onset preeclampsia prediction classifiers using machine learning methods. Furthermore, we validated the performance of the classifiers using the external and internal validation cohorts and assessed the area under the curve and positive predictive value. RESULTS: We detected 77 genes, including messenger RNA (44%) and microRNA (26%), that were differentially expressed in healthy mothers and mothers with preterm preeclampsia before symptom onset, which could separate participants with preterm preeclampsia from healthy samples and that played critical functional roles in preeclampsia physiology. We developed 2 classifiers for predicting preterm preeclampsia and early-onset preeclampsia before diagnosis based on 13 cell-free RNA signatures and 2 clinical features (in vitro fertilization and mean arterial pressure), respectively. Notably, both classifiers showed enhanced performance when compared with the existing methods. The preterm preeclampsia prediction model achieved 81% area under the curve and 68% positive predictive value in an independent validation cohort (preterm, n=46; control, n=151); the early-onset preeclampsia prediction model had an area under the curve of 88% and a positive predictive value of 73% in an external validation cohort (early-onset preeclampsia, n=28; control, n=234). Furthermore, we demonstrated that downregulation of microRNAs may play vital roles in preeclampsia through the upregulation of preeclampsia-relevant target genes. CONCLUSION: In this cohort study, a comprehensive transcriptomic landscape of different RNA biotypes in preeclampsia was presented and 2 advanced classifiers with substantial clinical importance for preterm preeclampsia and early-onset preeclampsia prediction before symptom onset were developed. We demonstrated that messenger RNA, microRNA, and long noncoding RNA can simultaneously serve as potential biomarkers of preeclampsia, holding the promise of prevention of preeclampsia in the future. Abnormal cell-free messenger RNA, microRNA, and long noncoding RNA molecular changes may help to elucidate the pathogenic determinants of preeclampsia and open new therapeutic windows to effectively reduce pregnancy complications and fetal morbidity.

Quadruple-editing of the MAPK and PI3K pathways effectively blocks the progression of KRAS-mutated colorectal cancer cells.

Mutated KRAS promotes the activation of the MAPK pathway and the progression of colorectal cancer (CRC) cells. Aberrant activation of the PI3K pathway strongly attenuates the efficacy of MAPK suppression in KRAS-mutated CRC. The development of a novel strategy targeting a dual pathway is therefore highly essential for the therapy of KRAS-mutated CRC. In this study, a quadruple-depleting system for the KRAS, MEK1, PIK3CA, and MTOR genes based on CRISPR/SaCas9 was developed. Adenovirus serotype 5 (ADV5) was integrated with two engineered proteins, an adaptor and a protector, to form ADV-protein complex (APC) for systemic delivery of the CRISPR system. Quadruple-editing could significantly inhibit the MAPK and PI3K pathways in CRC cells with oncogenic mutations of KRAS and PIK3CA or with KRAS mutation and compensated PI3K activation. Compared with MEK and PI3K/MTOR inhibitors, quadruple-editing induced more significant survival inhibition on primary CRC cells with oncogenic mutations of KRAS and PIK3CA. The adaptor specifically targeting EpCAM and the hexon-shielding protector could dramatically enhance ADV5 infection efficiency to CRC cells and significantly reduce off-targeting tropisms to many organs except the colon. Moreover, quadruple-editing intravenously delivered by APC significantly blocked the dual pathway and tumor growth of KRAS-mutated CRC cells, without influencing normal tissues in cell- and patient-derived xenograft models. Therefore, APC-delivered quadruple-editing of the MAPK and PI3K pathways shows a promising therapeutic potential for KRAS-mutated CRC.

Prophage Hunter: an integrative hunting tool for active prophages.

Identifying active prophages is critical for studying coevolution of phage and bacteria, investigating phage physiology and biochemistry, and engineering designer phages for diverse applications. We present Prophage Hunter, a tool aimed at hunting for active prophages from whole genome assembly of bacteria. Combining sequence similarity-based matching and genetic features-based machine learning classification, we developed a novel scoring system that exhibits higher accuracy than current tools in predicting active prophages on the validation datasets. The option of skipping similarity matching is also available so that there's higher chance for novel phages to be discovered. Prophage Hunter provides a one-stop web service to extract prophage genomes from bacterial genomes, evaluate the activity of the prophages, identify phylogenetically related phages, and annotate the function of phage proteins. Prophage Hunter is freely available at https://pro-hunter.bgi.com/.

Dehydration stress extends mRNA 3' untranslated regions with noncoding RNA functions in Arabidopsis.

The 3' untranslated regions (3' UTRs) of mRNAs play important roles in the regulation of mRNA localization, translation, and stability. Alternative cleavage and polyadenylation (APA) generates mRNAs with different 3' UTRs, but the involvement of this process in stress response has not yet been clarified. Here, we report that a subset of stress-related genes exhibits 3' UTR extensions of their mRNAs during dehydration stress. These extended 3' UTRs have characteristics of long noncoding RNAs and likely do not interact with miRNAs. Functional studies using T-DNA insertion mutants reveal that they can act as antisense transcripts to repress expression levels of sense genes from the opposite strand or can activate the transcription or lead to read-through transcription of their downstream genes. Further analysis suggests that transcripts with 3' UTR extensions have weaker poly(A) signals than those without 3' UTR extensions. Finally, we show that their biogenesis is partially dependent on a trans-acting factor FPA. Taken together, we report that dehydration stress could induce transcript 3' UTR extensions and elucidate a novel function for these stress-induced 3' UTR extensions as long noncoding RNAs in the regulation of their neighboring genes.

ELF18-INDUCED LONG-NONCODING RNA Associates with Mediator to Enhance Expression of Innate Immune Response Genes in Arabidopsis.

The plant immune response is a complex process involving transcriptional and posttranscriptional regulation of gene expression. Responses to plant immunity are initiated upon the perception of pathogen-associated molecular patterns, including peptide fragment of bacterial flagellin (flg22) or translation elongation factor Tu (elf18). Here, we identify an Arabidopsis thaliana long-noncoding RNA, designated ELF18-INDUCED LONG-NONCODING RNA1 (ELENA1), as a factor enhancing resistance against Pseudomonas syringe pv tomato DC3000. ELENA1 knockdown plants show decreased expression of PATHOGENESIS-RELATED GENE1 (PR1) and the plants are susceptible to pathogens. By contrast, plants overexpressing ELENA1 show elevated PR1 expression after elf18 treatment and display a pathogen resistance phenotype. RNA-sequencing analysis of ELENA1-overexpressing plants after elf18 treatment confirms increased expression of defense-related genes compared with the wild type. ELENA1 directly interacts with Mediator subunit 19a (MED19a) and affects enrichment of MED19a on the PR1 promoter. These results show that MED19a regulates PR1 expression through ELENA1. Our findings uncover an additional layer of complexity, implicating long-noncoding RNAs in the transcriptional regulation of plant innate immunity.

CURLY LEAF Regulates Gene Sets Coordinating Seed Size and Lipid Biosynthesis.

CURLY LEAF (CLF), a histone methyltransferase of Polycomb Repressive Complex 2 (PRC2) for trimethylation of histone H3 Lys 27 (H3K27me3), has been thought as a negative regulator controlling mainly postgermination growth in Arabidopsis (Arabidopsis thaliana). Approximately 14% to 29% of genic regions are decorated by H3K27me3 in the Arabidopsis genome; however, transcriptional repression activities of PRC2 on a majority of these regions remain unclear. Here, by analysis of transcriptome profiles, we found that approximately 11.6% genes in the Arabidopsis genome were repressed by CLF in various organs. Unexpectedly, approximately 54% of these genes were preferentially repressed in siliques. Further analyses of 118 transcriptome datasets uncovered a group of genes that was preferentially expressed and repressed by CLF in embryos at the mature-green stage. This observation suggests that CLF mediates a large-scale H3K27me3 programming/reprogramming event during embryonic development. Plants of clf-28 produced bigger and heavier seeds with higher oil content, larger oil bodies, and altered long-chain fatty acid composition compared with wild type. Around 46% of CLF-repressed genes were associated with H3K27me3 marks; moreover, we verified histone modification and transcriptional repression by CLF on regulatory genes. Our results suggest that CLF silences specific gene expression modules. Genes operating within a module have various molecular functions, but they cooperate to regulate a similar physiological function during embryo development.

Noncoding and coding transcriptome responses of a marine diatom to phosphate fluctuations.

Phosphorus (P) is an essential element to all living cells, yet fluctuations in P concentrations are recurrent in the marine environment. Diatoms are amongst the most successful phytoplankton groups, adapting to and surviving periods of suboptimal conditions and resuming growth as soon as nutrient concentrations permit. A knowledge of the molecular underpinnings of diatom ecological success is, however, still very incomplete. By strand-specific RNA sequencing, we analyzed the global transcriptome changes of the diatom Phaeodactylum tricornutum in response to P fluctuations over a course of 8 d, defining five distinct physiological states. This study reports previously unidentified genes highly responsive to P stress in P. tricornutum. Our data also uncover the complexity of the P. tricornutum P-responsive sensory and signaling system that combines bacterial two-component systems with more complex pathways reminiscent of metazoans. Finally, we identify a multitude of novel long intergenic nonprotein coding RNAs (lincRNAs) specifically responsive to P depletion, suggesting putative regulatory roles in the regulation of P homeostasis. Our work provides additional molecular insights into the resilience of diatoms and their ecological success, and opens up novel routes to address and explore the function and regulatory roles of P. tricornutum lincRNAs in the context of nutrient stress.

Insights into salt tolerance from the genome of Thellungiella salsuginea.

Thellungiella salsuginea, a close relative of Arabidopsis, represents an extremophile model for abiotic stress tolerance studies. We present the draft sequence of the T. salsuginea genome, assembled based on ~134-fold coverage to seven chromosomes with a coding capacity of at least 28,457 genes. This genome provides resources and evidence about the nature of defense mechanisms constituting the genetic basis underlying plant abiotic stress tolerance. Comparative genomics and experimental analyses identified genes related to cation transport, abscisic acid signaling, and wax production prominent in T. salsuginea as possible contributors to its success in stressful environments.

Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis in Salvia miltiorrhiza.

BACKGROUND: Plant natural products have been co-opted for millennia by humans for various uses such as flavor, fragrances, and medicines. These compounds often are only produced in relatively low amounts and are difficult to chemically synthesize, limiting access. While elucidation of the underlying biosynthetic processes might help alleviate these issues (e.g., via metabolic engineering), investigation of this is hindered by the low levels of relevant gene expression and expansion of the corresponding enzymatic gene families. However, the often-inducible nature of such metabolic processes enables selection of those genes whose expression pattern indicates a role in production of the targeted natural product. RESULTS: Here, we combine metabolomics and transcriptomics to investigate the inducible biosynthesis of the bioactive diterpenoid tanshinones from the Chinese medicinal herb, Salvia miltiorrhiza (Danshen). Untargeted metabolomics investigation of elicited hairy root cultures indicated that tanshinone production was a dominant component of the metabolic response, increasing at later time points. A transcriptomic approach was applied to not only define a comprehensive transcriptome (comprised of 20,972 non-redundant genes), but also its response to induction, revealing 6,358 genes that exhibited differential expression, with significant enrichment for up-regulation of genes involved in stress, stimulus and immune response processes. Consistent with our metabolomics analysis, there appears to be a slower but more sustained increased in transcript levels of known genes from diterpenoid and, more specifically, tanshinone biosynthesis. Among the co-regulated genes were 70 transcription factors and 8 cytochromes P450, providing targets for future investigation. CONCLUSIONS: Our results indicate a biphasic response of Danshen terpenoid metabolism to elicitation, with early induction of sesqui- and tri- terpenoid biosynthesis, followed by later and more sustained production of the diterpenoid tanshinones. Our data provides a firm foundation for further elucidation of tanshinone and other inducible natural product metabolism in Danshen.

Mutation characteristics of cancer susceptibility genes in Chinese ovarian cancer patients.

Introduction: The association between mutations in susceptibility genes and the occurrence of ovarian cancer has been extensively studied. Previous research has primarily concentrated on genes involved in the homologous recombination repair pathway, particularly BRCA1 and BRCA2. However, a wider range of genes related to the DNA damage response pathways has not been fully explored. Methods: To investigate the mutation characteristics of cancer susceptibility genes in the Chinese ovarian cancer population and the associations between gene mutations and clinical data, this study initially gathered a total of 1171 Chinese ovarian cancer samples and compiled a dataset of germline mutations in 171 genes. Results: In this study, it was determined that MC1R and PRKDC were high-frequency ovarian cancer susceptibility genes in the Chinese population, exhibiting notable distinctions from those in European and American populations; moreover high-frequency mutation genes, such as MC1R: c.359T>C and PRKDC: c.10681T>A, typically had high-frequency mutation sites. Furthermore, we identified c.8187G>T as a characteristic mutation of BRCA2 in the Chinese population, and the CHEK2 mutation was significantly associated with the early onset of ovarian cancer, while the CDH1 and FAM175A mutations were more prevalent in Northeast China. Additionally, Fanconi anemia pathway-related genes were significantly associated with ovarian carcinogenesis. Conclusion: In summary, this research provided fundamental data support for the optimization of ovarian cancer gene screening policies and the determination of treatment, and contributed to the precise intervention and management of patients.

Single-cell transcriptomic analysis reveals transcriptional and cell subpopulation differences between human and pig immune cells.

BACKGROUND: The pig is a promising donor candidate for xenotransplantation. Understanding the differences between human and swine immune systems is critical for addressing xenotransplant rejection and hematopoietic reconstitution. The gene transcriptional profile differences between human and pig immune cell subpopulations have not been studied. To assess the similarities and differences between pigs and humans at the levels of gene transcriptional profiles or cell subpopulations are important for better understanding the cross-species similarity of humans and pigs, and it would help establish the fundamental principles necessary to genetically engineer donor pigs and improve xenotransplantation. OBJECTIVE: To assess the gene transcriptional similarities and differences between pigs and humans. METHODS: Two pigs and two healthy humans' PBMCs were sorted for 10 × genomics single-cell sequence. We generated integrated human-pig scRNA-seq data from human and pig PBMCs and defined the overall gene expression landscape of pig peripheral blood immune cell subpopulations by updating the set of human-porcine homologous genes. The subsets of immune cells were detected by flow cytometry. RESULTS: There were significantly less T cells, NK cells and monocytes but more B cells in pig peripheral blood than those in human peripheral blood. High oxidative phosphorylation, HIF-1, glycolysis, and lysosome-related gene expressions in pig CD14+ monocytes were observed, whereas pig CD14+ monocytes exhibited lower levels of cytokine receptors and JAK-STAT-related genes. Pig activated CD4+T cells decreased cell adhesion and inflammation, while enriched for migration and activation processes. Porcine GNLY+CD8+T cells reduced cytotoxicity and increased proliferation compared with human GNLY+CD8+T cells. Pig CD2+CD8+γδT cells were functionally homologous to human CD2+CD4+ γδT cells. Pig CD2-CD8-γδT cells expressed genes with quiescent and precursor characteristics, while CD2-CD8+γδT cells expressed migration and memory-related molecules. Pig CD24+ and CD5+B cells are associated with inflammatory responses. CONCLUSION: Our research with integrated scRNA-seq assays identified the different distribution of pig immune cell subpopulations and the different transcriptional profiles of human and pig immune cells. This study enables a deeper understanding of the development and function of porcine immune cells.

Multiome in the Same Cell Reveals the Impact of Osmotic Stress on Arabidopsis Root Tip Development at Single-Cell Level.

Cell-specific transcriptional regulatory networks (TRNs) play vital roles in plant development and response to environmental stresses. However, traditional single-cell mono-omics techniques are unable to directly capture the relationships and dynamics between different layers of molecular information within the same cells. While advanced algorithm facilitates merging scRNA-seq and scATAC-seq datasets, accurate data integration remains a challenge, particularly when investigating cell-type-specific TRNs. By examining gene expression and chromatin accessibility simultaneously in 16,670 Arabidopsis root tip nuclei, the TRNs are reconstructed that govern root tip development under osmotic stress. In contrast to commonly used computational integration at cell-type level, 12,968 peak-to-gene linkage is captured at the bona fide single-cell level and construct TRNs at an unprecedented resolution. Furthermore, the unprecedented datasets allow to more accurately reconstruct the coordinated changes of gene expression and chromatin states during cellular state transition. During root tip development, chromatin accessibility of initial cells precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for subsequent differentiation steps. Pseudo-time trajectory analysis reveal that osmotic stress can shift the functional differentiation of trichoblast. Candidate stress-related gene-linked cis-regulatory elements (gl-cCREs) as well as potential target genes are also identified, and uncovered large cellular heterogeneity under osmotic stress.

Computational Approaches and Challenges in Spatial Transcriptomics.

The development of spatial transcriptomics (ST) technologies has transformed genetic research from a single-cell data level to a two-dimensional spatial coordinate system and facilitated the study of the composition and function of various cell subsets in different environments and organs. The large-scale data generated by these ST technologies, which contain spatial gene expression information, have elicited the need for spatially resolved approaches to meet the requirements of computational and biological data interpretation. These requirements include dealing with the explosive growth of data to determine the cell-level and gene-level expression, correcting the inner batch effect and loss of expression to improve the data quality, conducting efficient interpretation and in-depth knowledge mining both at the single-cell and tissue-wide levels, and conducting multi-omics integration analysis to provide an extensible framework toward the in-depth understanding of biological processes. However, algorithms designed specifically for ST technologies to meet these requirements are still in their infancy. Here, we review computational approaches to these problems in light of corresponding issues and challenges, and present forward-looking insights into algorithm development.

A novel computational method enables RNA editome profiling during human hematopoiesis from scRNA-seq data.

RNA editing is a post-transcriptional modification with a cell-specific manner and important biological implications. Although single-cell RNA-seq (scRNA-seq) is an effective method for studying cellular heterogeneity, it is difficult to detect and study RNA editing events from scRNA-seq data because of the low sequencing coverage. To overcome this, we develop a computational method to systematically identify RNA editing sites of cell types from scRNA-seq data. To demonstrate its effectiveness, we apply it to scRNA-seq data of human hematopoietic stem/progenitor cells (HSPCs) with an annotated lineage differentiation relationship according to previous research and study the impacts of RNA editing on hematopoiesis. The dynamic editing patterns reveal the relevance of RNA editing on different HSPCs. For example, four microRNA (miRNA) target sites on 3' UTR of EIF2AK2 are edited across all HSPC populations, which may abolish the miRNA-mediated inhibition of EIF2AK2. Elevated EIF2AK2 may thus activate the integrated stress response (ISR) pathway to initiate global translational attenuation as a protective mechanism to maintain cellular homeostasis during HSPCs' differentiation. Besides, our findings also indicate that RNA editing plays an essential role in the coordination of lineage commitment and self-renewal of hematopoietic stem cells (HSCs). Taken together, we demonstrate the capacity of scRNA-seq data to exploit RNA editing events of cell types, and find that RNA editing may exert multiple modules of regulation in hematopoietic processes.

Establishment of bovine trophoblast stem cells.

Here, we report that a chemical cocktail (LCDM: leukemia inhibitory factor [LIF], CHIR99021, dimethinedene maleate [DiM], minocycline hydrochloride), previously developed for extended pluripotent stem cells (EPSCs) in mice and humans, enables de novo derivation and long-term culture of bovine trophoblast stem cells (TSCs). Bovine TSCs retain developmental potency to differentiate into mature trophoblast cells and exhibit transcriptomic and epigenetic (chromatin accessibility and DNA methylome) features characteristic of trophectoderm cells from early bovine embryos. The bovine TSCs established in this study will provide a model to study bovine placentation and early pregnancy failure.