Identification of the transcriptome signatures and immune-inflammatory responses in postmenopausal osteoporosis.

Postmenopausal osteoporosis is the most common type of osteoporosis in women. To date, little is known about their transcriptome signatures, although biomarkers from peripheral blood mononuclear cells are attractive for postmenopausal osteoporosis diagnoses. Here, we performed bulk RNA sequencing of 206 samples (124 postmenopausal osteoporosis and 82 normal samples) and described the clinical phenotypic characteristics of postmenopausal women. We then highlighted the gene set enrichment analyses between the extreme T-score group and the heathy control group, revealing that some immune-inflammatory responses were enhanced in postmenopausal osteoporosis, with representative pathways including the mitogen-activated protein kinase (NES = 1.6, FDR <0.11) pathway and B_CELL_RECEPTOR (NES = 1.69, FDR <0.15) pathway. Finally, we developed a combined risk prediction model based on lasso-logistic regression to predict postmenopausal osteoporosis, which combined eleven genes (PTGS2, CXCL16, NECAP1, RPS23, SSR3, CD74, IL4R, BTBD2, PIGS, LILRA2, MAP3K11) and three pieces of clinical information (age, procollagen I N-terminal propeptide, ß isomer of C-terminal telopeptide of type I) and provided the best prediction ability (AUC = 0.97). Taken together, this study filled a gap in the large-scale transcriptome signature profiles and revealed the close relationship between immune-inflammatory responses and postmenopausal osteoporosis, providing a unique perspective for understanding the occurrence and development of postmenopausal osteoporosis.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development.

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

Genome-Wide Diversity Analysis of African Swine Fever Virus Based on a Curated Dataset.

African swine fever (ASF) is a lethal contagious viral disease of domestic pigs and wild boars caused by the African swine fever virus (ASFV). The pandemic spread of ASF has had serious effects on the global pig industry. Virus genome sequencing and comparison play an important role in tracking the outbreaks of the disease and tracing the transmission of the virus. Although more than 140 ASFV genome sequences have been deposited in the public databases, the genome-wide diversity of ASFV remains unclear. Here we prepared a curated dataset of ASFV genome sequences by filtering genomes with sequencing errors as well as duplicated genomes. A total of 123 ASFV genome sequences were included in the dataset, representing 10 genotypes collected between 1949 and 2020. Phylogenetic analysis based on whole-genome sequences provided high-resolution topology in differentiating closely related ASFV isolates, and drew new clues in the classification of some ASFV isolates. Genome-wide diversity of ASFV genomes was explored by pairwise sequence similarity comparison and ORF distribution comparison. Tandem repeat sequences were found widely distributed and highly varied in ASFV genomes. Structural variation and highly variable poly G or poly C tracts also contributed to the genome diversity. This study expanded our knowledge on the patterns of genetic diversity and evolution of ASFV, and provided valuable information for diagnosis improvement and vaccine development.

Massively targeted evaluation of therapeutic CRISPR off-targets in cells.

Methods for sensitive and high-throughput evaluation of CRISPR RNA-guided nucleases (RGNs) off-targets (OTs) are essential for advancing RGN-based gene therapies. Here we report SURRO-seq for simultaneously evaluating thousands of therapeutic RGN OTs in cells. SURRO-seq captures RGN-induced indels in cells by pooled lentiviral OTs libraries and deep sequencing, an approach comparable and complementary to OTs detection by T7 endonuclease 1, GUIDE-seq, and CIRCLE-seq. Application of SURRO-seq to 8150 OTs from 110 therapeutic RGNs identifies significantly detectable indels in 783 OTs, of which 37 OTs are found in cancer genes and 23 OTs are further validated in five human cell lines by targeted amplicon sequencing. Finally, SURRO-seq reveals that thermodynamically stable wobble base pair (rG•dT) and free binding energy strongly affect RGN specificity. Our study emphasizes the necessity of thoroughly evaluating therapeutic RGN OTs to minimize inevitable off-target effects.

Effect of high-concentrate diets on mRNA expression of genes related to muscle fiber type and metabolism of psoas major muscle in goats.

In the process of modern breeding, high-concentrate diets are widely used to meet the high energy nutritional requirements of animals but change the form of access to energy and nutrients and the way the organism metabolizes them. Goat psoas major (PM) muscle is a hybrid skeletal muscle whose characteristics are important for the motility and meat quality of goats. However, there are few studies on the effects of high-concentrate diets on the muscle type and metabolic characteristics of PM in goats. In this study, two treatment groups were set up: high concentrate group (HC) and control group (C). The expression of genes related to muscle type and metabolism of the PM was examined by quantitative PCR. The results showed that high concentrate promoted the conversion of PM fibers from intermediate to slow type at the mRNA level, improved the absorption, transport, and oxidation of fat by PM, and upregulated the expression of calpain system. These changes may be regulated by the involvement of differential expression of MSTN, Myf-5, and IGF-2. These results suggest that high concentrate may exert a positive effect on skeletal muscle function, metabolism, and meat quality in goats by affecting the expression of muscle type and metabolism-related genes.

Single-cell RNA sequencing technologies and applications: A brief overview.

Single-cell RNA sequencing (scRNA-seq) technology has become the state-of-the-art approach for unravelling the heterogeneity and complexity of RNA transcripts within individual cells, as well as revealing the composition of different cell types and functions within highly organized tissues/organs/organisms. Since its first discovery in 2009, studies based on scRNA-seq provide massive information across different fields making exciting new discoveries in better understanding the composition and interaction of cells within humans, model animals and plants. In this review, we provide a concise overview about the scRNA-seq technology, experimental and computational procedures for transforming the biological and molecular processes into computational and statistical data. We also provide an explanation of the key technological steps in implementing the technology. We highlight a few examples on how scRNA-seq can provide unique information for better understanding health and diseases. One important application of the scRNA-seq technology is to build a better and high-resolution catalogue of cells in all living organism, commonly known as atlas, which is key resource to better understand and provide a solution in treating diseases. While great promises have been demonstrated with the technology in all areas, we further highlight a few remaining challenges to be overcome and its great potentials in transforming current protocols in disease diagnosis and treatment.

A Brief Overview of Global Trends in MSC-Based Cell Therapy.

Human mesenchymal stem cells (MSCs), also known as mesenchymal stromal cells or medicinal signaling cells, are important adult stem cells for regenerative medicine, largely due to their regenerative characteristics such as self-renewal, secretion of trophic factors, and the capability of inducing mesenchymal cell lineages. MSCs also possess homing and trophic properties modulating immune system, influencing microenvironment around damaged tissues and enhancing tissue repair, thus offering a broad perspective in cell-based therapies. Therefore, it is not surprising that MSCs have been the broadly used adult stem cells in clinical trials. To gain better insights into the current applications of MSCs in clinical applications, we perform a comprehensive review of reported data of MSCs clinical trials conducted globally. We summarize the biological effects and mechanisms of action of MSCs, elucidating recent clinical trials phases and findings, highlighting therapeutic effects of MSCs in several representative diseases, including neurological, musculoskeletal diseases and most recent Coronavirus infectious disease. Finally, we also highlight the challenges faced by many clinical trials and propose potential solutions to streamline the use of MSCs in routine clinical applications and regenerative medicine.

Genome-wide annotation of protein-coding genes in pig.

BACKGROUND: There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs. RESULTS: An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented. CONCLUSIONS: Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Systematical analysis reveals a strong cancer relevance of CREB1-regulated genes.

The transcription factor cyclic-AMP response element-binding protein 1 (CREB1) responds to cAMP level and controls the expression of target genes, which regulates nutrition partitioning. The promoters of CREB1-targeted genes responsive to cAMP have been extensively investigated and characterized with the presence of both cAMP response element and TATA box. Compelling evidence demonstrates that CREB1 also plays an essential role in promoting tumor development. However, only very few genes required for cell survival, proliferation and migration are known to be constitutively regulated by CREB1 in tumors. Their promoters mostly do not harbor any cAMP response element. Thus, it is very likely that CREB1 regulates the expressions of distinct sets of target genes in normal tissues and tumors. The whole gene network constitutively regulated by CREB1 in tumors has remained unrevealed. Here, we employ a systematical and integrative approach to decipher this gene network in the context of both tissue cultured cancer cells and patient samples. We combine transcriptomic, Rank-Rank Hypergeometric Overlap, and Chipseq analysis, to define and characterize CREB1-regulated genes in a multidimensional fashion. A strong cancer relevance of those top-ranked targets, which meet the most stringent criteria, is eventually verified by overall survival analysis of cancer patients. These findings strongly suggest the importance of genes constitutively regulated by CREB1 for their implicative involvement in promoting tumorigenesis.

The Dynamic Changes of Transcription Factors During the Development Processes of Human Biparental and Uniparental Embryos.

Previous studies have revealed that transcription factors (TFs) play important roles in biparental (BI) early human embryogenesis. However, the contribution of TFs during early uniparental embryo development is still largely unknown. Here we systematically studied the expression profiles of transcription factors in early embryonic development and revealed the dynamic changes of TFs in human biparental and uniparental embryogenesis by single-cell RNA sequencing (scRNA-seq). In general, the TF expression model of uniparental embryos showed a high degree of conformity with biparental embryos. The detailed network analysis of three different types of embryos identified that 10 out of 17 hub TFs were shared or specifically owned, such as ZNF480, ZNF581, PHB, and POU5F1, were four shared TFs, ZFN534, GTF3A, ZNF771, TEAD4, and LIN28A, were androgenic (AG) specific TFs, and ZFP42 was the only one parthenogenetic (PG) specific TF. All the four shared TFs were validated using human embryonic stem cell (hESC) differentiation experiments; most of their target genes are responsible for stem cell maintenance and differentiation. We also found that Zf-C2H2, HMG, and MYB were three dominant transcription factor families that appeared in early embryogenesis. Altogether, our work provides a comprehensive regulatory framework and better understanding of TF function in human biparental and uniparental embryogenesis.

Enhancing CRISPR-Cas9 gRNA efficiency prediction by data integration and deep learning.

The design of CRISPR gRNAs requires accurate on-target efficiency predictions, which demand high-quality gRNA activity data and efficient modeling. To advance, we here report on the generation of on-target gRNA activity data for 10,592 SpCas9 gRNAs. Integrating these with complementary published data, we train a deep learning model, CRISPRon, on 23,902 gRNAs. Compared to existing tools, CRISPRon exhibits significantly higher prediction performances on four test datasets not overlapping with training data used for the development of these tools. Furthermore, we present an interactive gRNA design webserver based on the CRISPRon standalone software, both available via https://rth.dk/resources/crispr/ . CRISPRon advances CRISPR applications by providing more accurate gRNA efficiency predictions than the existing tools.

Feeder-free generation and transcriptome characterization of functional mesenchymal stromal cells from human pluripotent stem cells.

Induced mesenchymal stromal cells (iMSCs) derived from human pluripotent stem cells (PSCs) are attractive cells for regenerative medicine. However, the transcriptome of iMSCs and signature genes that can distinguish MSCs from fibroblasts and other cell types are rarely explored. In this study, we reported an optimized feeder-free method for the generation of iMSCs from human pluripotent stem cells. These iMSCs display a typical MSC morphology, express classic MSC markers (CD29, CD44, CD73, CD90, CD105, CD166), are negative for lymphocyte markers (CD11b, CD14, CD31, CD34, CD45, HLA-DR), and are potent for osteogenic and chondrogenic differentiation. Using genome-wide transcriptome profiling, we created an easily accessible transcriptome reference for the process of differentiating PSCs into iMSCs. The iMSC transcriptome reference revealed clear patterns in the silencing of pluripotency genes, activation of lineage commitment genes, and activation of mesenchymal genes during iMSC generation. All previously known positive and negative markers for MSCs were confirmed by our iMSC transcriptomic reference, and most importantly, gene classification and time course analysis identified 52 genes including FN1, TGFB1, TAGLN and SERPINE1, which showed significantly higher expression in MSCs (over 3 folds) than fibroblasts and other cell types. Taken together, these results provide a useful method and important resources for developing and understanding iMSCs in regenerative medicine.

Haplotyping by CRISPR-mediated DNA circularization (CRISPR-hapC) broadens allele-specific gene editing.

Allele-specific protospacer adjacent motif (asPAM)-positioning SNPs and CRISPRs are valuable resources for gene therapy of dominant disorders. However, one technical hurdle is to identify the haplotype comprising the disease-causing allele and the distal asPAM SNPs. Here, we describe a novel CRISPR-based method (CRISPR-hapC) for haplotyping. Based on the generation (with a pair of CRISPRs) of extrachromosomal circular DNA in cells, the CRISPR-hapC can map haplotypes from a few hundred bases to over 200 Mb. To streamline and demonstrate the applicability of the CRISPR-hapC and asPAM CRISPR for allele-specific gene editing, we reanalyzed the 1000 human pan-genome and generated a high frequency asPAM SNP and CRISPR database (www.crispratlas.com/knockout) for four CRISPR systems (SaCas9, SpCas9, xCas9 and Cas12a). Using the huntingtin (HTT) CAG expansion and transthyretin (TTR) exon 2 mutation as examples, we showed that the asPAM CRISPRs can specifically discriminate active and dead PAMs for all 23 loci tested. Combination of the CRISPR-hapC and asPAM CRISPRs further demonstrated the capability for achieving highly accurate and haplotype-specific deletion of the HTT CAG expansion allele and TTR exon 2 mutation in human cells. Taken together, our study provides a new approach and an important resource for genome research and allele-specific (haplotype-specific) gene therapy.

Genome Sequencing Explores Complexity of Chromosomal Abnormalities in Recurrent Miscarriage.

Recurrent miscarriage (RM) affects millions of couples globally, and half of them have no demonstrated etiology. Genome sequencing (GS) is an enhanced and novel cytogenetic tool to define the contribution of chromosomal abnormalities in human diseases. In this study we evaluated its utility in RM-affected couples. We performed low-pass GS retrospectively for 1,090 RM-affected couples, all of whom had routine chromosome analysis. A customized sequencing and interpretation pipeline was developed to identify chromosomal rearrangements and deletions/duplications with confirmation by fluorescence in situ hybridization, chromosomal microarray analysis, and PCR studies. Low-pass GS yielded results in 1,077 of 1,090 couples (98.8%) and detected 127 chromosomal abnormalities in 11.7% (126/1,077) of couples; both members of one couple were identified with inversions. Of the 126 couples, 39.7% (50/126) had received former diagnostic results by karyotyping characteristic of normal human male or female karyotypes. Low-pass GS revealed additional chromosomal abnormalities in 50 (4.0%) couples, including eight with balanced translocations and 42 inversions. Follow-up studies of these couples showed a higher miscarriage/fetal-anomaly rate of 5/10 (50%) compared to 21/93 (22.6%) in couples with normal GS, resulting in a relative risk of 2.2 (95% confidence interval, 1.1 to 4.6). In these couples, this protocol significantly increased the diagnostic yield of chromosomal abnormalities per couple (11.7%) in comparison to chromosome analysis (8.0%, chi-square test p = 0.000751). In summary, low-pass GS identified underlying chromosomal aberrations in 1 in 9 RM-affected couples, enabling identification of a subgroup of couples with increased risk of subsequent miscarriage who would benefit from a personalized intervention.

Single-Cell Transcriptome Analysis of Uniparental Embryos Reveals Parent-of-Origin Effects on Human Preimplantation Development.

To investigate the contribution of parental genomes to early embryogenesis, we profiled the single-cell transcriptomes of human biparental and uniparental embryos systematically from the 1-cell to the morula stage. We observed that uniparental embryos exhibited variable and decreased embryonic genome activation (EGA). Comparative transcriptome analysis identified 807 maternally biased expressed genes (MBGs) and 581 paternally biased expressed genes (PBGs) in the preimplantation stages. MBGs became apparent at the 4-cell stage and contributed to the initiation of EGA, whereas PBGs preferentially appeared at the 8-cell stage and might affect embryo compaction and trophectoderm specification. Regulatory network analysis revealed that DUX4, EGR2, and DUXA are key transcription factors in MBGs' expression; ZNF263 and KLF3 are important for PBGs' expression. We demonstrated that parent-specific DNA methylation might account for the expression of most PBGs. Our results provide a valuable resource to understand parental genome activation and might help to elucidate parent-of-origin effects in early human development.

STAT3 associates with vacuolar H+-ATPase and regulates cytosolic and lysosomal pH.

Dysregulated intracellular pH is emerging as a hallmark of cancer. In spite of their acidic environment and increased acid production, cancer cells maintain alkaline intracellular pH that promotes cancer progression by inhibiting apoptosis and increasing glycolysis, cell growth, migration, and invasion. Here we identify signal transducer and activator of transcription-3 (STAT3) as a key factor in the preservation of alkaline cytosol. STAT3 associates with the vacuolar H+-ATPase in a coiled-coil domain-dependent manner and increases its activity in living cells and in vitro. Accordingly, STAT3 depletion disrupts intracellular proton equilibrium by decreasing cytosolic pH and increasing lysosomal pH, respectively. This dysregulation can be reverted by reconstitution with wild-type STAT3 or STAT3 mutants unable to activate target genes (Tyr705Phe and DNA-binding mutant) or to regulate mitochondrial respiration (Ser727Ala). Upon cytosolic acidification, STAT3 is transcriptionally inactivated and further recruited to lysosomal membranes to reestablish intracellular proton equilibrium. These data reveal STAT3 as a regulator of intracellular pH and, vice versa, intracellular pH as a regulator of STAT3 localization and activity.

Genome-wide determination of on-target and off-target characteristics for RNA-guided DNA methylation by dCas9 methyltransferases.

Background: Fusion of DNA methyltransferase domains to the nuclease-deficient clustered regularly interspaced short palindromic repeat (CRISPR) associated protein 9 (dCas9) has been used for epigenome editing, but the specificities of these dCas9 methyltransferases have not been fully investigated. Findings: We generated CRISPR-guided DNA methyltransferases by fusing the catalytic domain of DNMT3A or DNMT3B to the C terminus of the dCas9 protein from Streptococcus pyogenes and validated its on-target and global off-target characteristics. Using targeted quantitative bisulfite pyrosequencing, we prove that dCas9-BFP-DNMT3A and dCas9-BFP-DNMT3B can efficiently methylate the CpG dinucleotides flanking its target sites at different genomic loci (uPA and TGFBR3) in human embryonic kidney cells (HEK293T). Furthermore, we conducted whole genome bisulfite sequencing (WGBS) to address the specificity of our dCas9 methyltransferases. WGBS revealed that although dCas9-BFP-DNMT3A and dCas9-BFP-DNMT3B did not cause global methylation changes, a substantial number (more than 1000) of the off-target differentially methylated regions (DMRs) were identified. The off-target DMRs, which were hypermethylated in cells expressing dCas9 methyltransferase and guide RNAs, were predominantly found in promoter regions, 5΄ untranslated regions, CpG islands, and DNase I hypersensitivity sites, whereas unexpected hypomethylated off-target DMRs were significantly enriched in repeated sequences. Through chromatin immunoprecipitation with massive parallel DNA sequencing analysis, we further revealed that these off-target DMRs were weakly correlated with dCas9 off-target binding sites. Using quantitative polymerase chain reaction, RNA sequencing, and fluorescence reporter cells, we also found that dCas9-BFP-DNMT3A and dCas9-BFP-DNMT3B can mediate transient inhibition of gene expression, which might be caused by dCas9-mediated de novo DNA methylation as well as interference with transcription. Conclusion: Our results prove that dCas9 methyltransferases cause efficient RNA-guided methylation of specific endogenous CpGs. However, there is significant off-target methylation indicating that further improvements of the specificity of CRISPR-dCas9 based DNA methylation modifiers are required.

Central and Peripheral Nervous System Progenitors Derived from Human Pluripotent Stem Cells Reveal a Unique Temporal and Cell-Type Specific Expression of PMCAs.

The P-type ATPases family consists of ion and lipid transporters. Their unique diversity in function and expression is critical for normal development. In this study we investigated human pluripotent stem cells (hPSC) and different neural progenitor states to characterize the expression of the plasma membrane calcium ATPases (PMCAs) during human neural development and in mature mesencephalic dopaminergic (mesDA) neurons. Our RNA sequencing data identified a dynamic change in ATPase expression correlating with the differentiation time of the neural progenitors, which was independent of the neuronal progenitor type. Expression of ATP2B1 and ATP2B4 were the most abundantly expressed, in accordance with their main role in Ca2+ regulation and we observed all of the PMCAs to have a subcellular punctate localization. Interestingly in hPSCs ATP2B1 and ATP2B3 were highly expressed in a cell cycle specific manner and ATP2B2 and ATP2B4 were highly expressed in a hPSC sub-population. In neural rosettes a strong apical PMCA expression was identified in the luminal region. Lastly, we confirmed all PMCAs to be expressed in mesDA neurons, however at varying levels. Our results reveal that PMCA expression dynamically changes during stem cell differentiation and highlights the diverging needs of cell populations to regulate and properly integrate Ca2+ changes, which can ultimately correspond to changes in specific stem cell transcription states.

BS-virus-finder: virus integration calling using bisulfite sequencing data.

Background: DNA methylation plays a key role in the regulation of gene expression and carcinogenesis. Bisulfite sequencing studies mainly focus on calling single nucleotide polymorphism, different methylation region, and find allele-specific DNA methylation. Until now, only a few software tools have focused on virus integration using bisulfite sequencing data. Findings: We have developed a new and easy-to-use software tool, named BS-virus-finder (BSVF, RRID:SCR_015727), to detect viral integration breakpoints in whole human genomes. The tool is hosted at https://github.com/BGI-SZ/BSVF. Conclusions: BS-virus-finder demonstrates high sensitivity and specificity. It is useful in epigenetic studies and to reveal the relationship between viral integration and DNA methylation. BS-virus-finder is the first software tool to detect virus integration loci by using bisulfite sequencing data.