Sustained release of naringin from silk-fibroin-nanohydroxyapatite scaffold for the enhancement of bone regeneration

Bone defects are a common challenge in the clinical setting. Bone tissue engineering (BTE) is an effective treatment for the clinical problem of large bone defects. In this study, we fabricated silk fibroin (SF)/hydroxyapatite (HAp) scaffolds inlaid with naringin poly lactic-co-glycolic acid (PLGA) microspheres, investigating the feasibility of their application in BTE. Naringin PLGA microspheres were manufactured and adhered to the SF/HAp scaffold. Bone mesenchymal stem cells (BMSCs) were inoculated onto the SF/HAp scaffold containing naringin PLGA microsphere to examine the biocompatibility of the SF/HAp scaffolds. A rabbit femoral distal bone defect model was used to evaluate the in vivo function of the SF/HAp scaffolds containing naringin-loaded PLGA microspheres. The current study demonstrated that SF/HAp scaffolds containing naringin-loaded PLGA microspheres show promise as osteo-modulatory biomaterials for bone regeneration.

Mendelian randomization analysis identified genes potentially pleiotropically associated with periodontitis

OBJECTIVE

To prioritize genes that were pleiotropically or potentially causally associated with periodontitis.

METHODS

We applied the summary data-based Mendelian randomization (SMR) method integrating genome-wide association study (GWAS) for periodontitis and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with periodontitis. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analysis were done for participants of European and East Asian ancestries, separately.

RESULTS

We identified multiple genes showing pleiotropic association with periodontitis in participants of European ancestry and participants of East Asian ancestry. PDCD2 (corresponding probe: ILMN_1758915) was the top hit showing pleotropic association with periodontitis in the participants of European ancestry using CAGE eQTL data, and BX093763 (corresponding probe: ILMN_1899903) and AC104135.3 (corresponding probe: ENSG00000204792.2) were the top hits in the participants of East Asian ancestry using CAGE eQTL data and GTEx eQTL data, respectively.

CONCLUSION

We identified multiple genes that may be involved in the pathogenesis of periodontitis in participants of European ancestry and participants of East Asian ancestry. Our findings provided important leads to a better understanding of the mechanisms underlying periodontitis and revealed potential therapeutic targets for the effective treatment of periodontitis.

Interpopulation differences of retroduplication variations (RDVs) in rice retrogenes and their phenotypic correlations

Retroduplication variation (RDV), a type of retrocopy polymorphism, is considered to have essential biological significance, but its effect on gene function and species phenotype is still poorly understood. To this end, we analyzed the retrocopies and RDVs in 3,010 rice genomes. We calculated the RDV frequencies in the genome of each rice population; detected the mutated, ancestral and expressed retrogenes in rice genomes; and analyzed their RDV influence on rice phenotypic traits. Collectively, 73 RDVs were identified, and 14 RDVs in ancestral retrogenes can significantly affect rice phenotypes. Our research reveals that RDV plays an important role in rice migration, domestication and evolution. We think that RDV is a good molecular breeding marker candidate. To our knowledge, this is the first study on the relationship between retrogene function, expression, RDV and species phenotype.

Identification of potential therapeutic target of naringenin in breast cancer stem cells inhibition by bioinformatics and in vitro studies

Cancer therapy is a strategic measure in inhibiting breast cancer stem cell (BCSC) pathways. Naringenin, a citrus flavonoid, was found to increase breast cancer cells' sensitivity to chemotherapeutic agents. Bioinformatics study and 3D tumorsphere in vitro modeling in breast cancer (mammosphere) were used in this study, which aims to explore the potential therapeutic targets of naringenin (PTTNs) in inhibiting BCSCs. Bioinformatic analyses identified direct target proteins (DTPs), indirect target proteins (ITPs), naringenin-mediated proteins (NMPs), BCSC regulatory genes, and PTTNs. The PTTNs were further analyzed for gene ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, protein-protein interaction (PPI) networks, and hub protein selection. Mammospheres were cultured in serum-free media. The effects of naringenin were measured by MTT-based cytotoxicity, mammosphere forming potential (MFP), colony formation, scratch wound-healing assay, and flow cytometry-based cell cycle analyses and apoptosis assays. Gene expression analysis was performed using real-time quantitative polymerase chain reaction (q-RT PCR). Bioinformatics analysis revealed p53 and estrogen receptor alpha (ERα) as PTTNs, and KEGG pathway enrichment analysis revealed that TGF-ß and Wnt/ß-catenin pathways are regulated by PTTNs. Naringenin demonstrated cytotoxicity and inhibited mammosphere and colony formation, migration, and epithelial to mesenchymal transition in the mammosphere. The mRNA of tumor suppressors P53 and ERα were downregulated in the mammosphere, but were significantly upregulated upon naringenin treatment. By modulating the P53 and ERα mRNA, naringenin has the potential of inhibiting BCSCs. Further studies on the molecular mechanism and formulation of naringenin in BCSCs would be beneficial for its development as a BCSC-targeting drug.

Bioinformatics analyses of significant genes, related pathways, and candidate diagnostic biomarkers and molecular targets in SARS-CoV-2/COVID-19

Severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) infection is a leading cause of pneumonia and death. The aim of this investigation is to identify the key genes in SARS-CoV-2 infection and uncover their potential functions. We downloaded the expression profiling by high throughput sequencing of GSE152075 from the Gene Expression Omnibus database. Normalization of the data from primary SARS-CoV-2 infected samples and negative control samples in the database was conducted using R software. Then, joint analysis of the data was performed. Pathway and Gene ontology (GO) enrichment analyses were performed, and the protein-protein interaction (PPI) network, target gene - miRNA regulatory network, target gene - TF regulatory network of the differentially expressed genes (DEGs) were constructed using Cytoscape software. Identification of diagnostic biomarkers was conducted using receiver operating characteristic (ROC) curve analysis. 994 DEGs (496 up regulated and 498 down regulated genes) were identified. Pathway and GO enrichment analysis showed up and down regulated genes mainly enriched in the NOD-like receptor signaling pathway, Ribosome, response to external biotic stimulus and viral transcription in SARS-CoV-2 infection. Down and up regulated genes were selected to establish the PPI network, modules, target gene - miRNA regulatory network, target gene - TF regulatory network revealed that these genes were involved in adaptive immune system, fluid shear stress and atherosclerosis, influenza A and protein processing in endoplasmic reticulum. In total, ten genes (CBL, ISG15, NEDD4, PML, REL, CTNNB1, ERBB2, JUN, RPS8 and STUB1) were identified as good diagnostic biomarkers. In conclusion, the identified DEGs, hub genes and target genes contribute to the understanding of the molecular mechanisms underlying the advancement of SARS-CoV-2 infection and they may be used as diagnostic and molecular targets for the treatment of patients with SARS-CoV-2 infection in the future.

Lingzhi and San-Miao-San with hyaluronic acid gel mitigate cartilage degeneration in anterior cruciate ligament transection induced osteoarthritis

Objective

To investigate the mitigate efficacy of Chinese medicine Lingzhi (LZ) and San-Miao-San (SMS) combined with hyaluronic acid (HA)-gel in attenuating cartilage degeneration in traumatic osteoarthritis (OA).

Methods

The standardized surgery of anterior cruciate ligament transection (ACLT) was made from the medial compartment of right hind limbs of 8-week-old female SD rats and resulted in a traumatic OA. Rats (n ​= ​5/group) were treated once intra-articular injection of 50 ​μl HA-gel, 50 ​μl HA-gel+50 ​μg LZ-SMS, 50 ​μl of saline+50 ​μg LZ-SMS and null (ACLT group) respectively, except sham group. Limbs were harvested for μCT scan and histopathological staining 3-month post-treatment. Inflammatory cytokines from plasma and synovial fluid were detected using Immunology Multiplex Assay kit. The putative targets of active compounds in LZ-SMS and known therapeutic targets for OA were combined to construct protein–protein interaction network. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was adopted to predict the potential targets and signaling pathway of LZ-SMS in OA through the tool of DAVID Bioinformatics.

Results

In vivo, HA-gel ​+ ​LZ-SMS treatment resulted in a higher volume ratio of hyaline cartilage (HC)/calcified cartilage (CC) and HC/Sum (total volume of cartilage), compared to ACLT and HA-gel groups. In addition, histological results showed the elevated cartilage matrix, chondrogenic and osteoblastic signals in HA-gel ​+ ​LZ-SMS treatment. Treatment also significantly altered subchondral bone (SCB) structure including an increase in BV/TV, Tb.Th, BMD, Conn.Dn, Tb.N, and DA, as well as a significant decrease in Tb.Sp and Po(tot), which implied a protective effect on maintaining the stabilization of tibial SCB microstructure. Furthermore, there was also a down-regulated inflammatory cytokines and upregulated anti-inflammatory cytokine IL-10 in HA+LZ-SMS group. Finally, 64 shared targets from 37 active compounds in LZ-SMS related to the core genes for the development of OA. LZ-SMS has a putative role in regulating inflammatory circumstance through influencing the MAPK signaling pathway.

Conclusion

Our study elucidated a protective effect of HA-gel ​+ ​LZ-SMS in mitigating cartilage degradation and putative interaction with targets and signaling pathway for the development of traumatic OA.

The translational potential of this article

Our results provide a biological rationale for the use of LZ-SMS as a potential candidate for OA treatment.

ceRNA network development and tumour-infiltrating immune cell analysis of metastatic breast cancer to bone

PURPOSE

Advanced breast cancer commonly metastasises to bone; however, the molecular mechanisms underlying the affinity for breast cancer cells to bone remains unclear. Thus, we developed nomograms based on a competing endogenous RNA (ceRNA) network and analysed tumour-infiltrating immune cells to elucidate the molecular pathways that may predict prognosis in patients with breast cancer.

METHODS

We obtained the RNA expression profile of 1091 primary breast cancer samples included in The Cancer Genome Atlas database, 58 of which were from patients with bone metastasis. We analysed the differential RNA expression patterns between breast cancer with and without bone metastasis and developed a ceRNA network. Cibersort was employed to differentiate between immune cell types based on tumour transcripts. Nomograms were then established based on the ceRNA network and immune cell analysis. The value of prognostic factors was evaluated by Kaplan-Meier survival analysis and a Cox proportional risk model.

RESULTS

We found significant differences in long non-coding RNAs (lncRNAs), 18 microRNAs (miRNAs), and 20 messenger RNAs (mRNAs) between breast cancer with and without bone metastasis, which were used to construct a ceRNA network. We found that the protein-coding genes GJB3, CAMMV, PTPRZ1, and FBN3 were significantly differentially expressed by Kaplan-Meier analysis. We also observed significant differences in the abundance of plasma cell and follicular helper T cell populations between the two groups. In addition, the proportion of mast cells, gamma delta T cells, and plasma cells differed depending on disease location and stage. Our analysis showed that a high proportion of follicular helper T cells and a low proportion of eosinophils promoted survival and that DLX6-AS1, Wnt6, and GABBR2 expression may be associated with bone metastasis in breast cancer.

CONCLUSIONS

We developed a bioinformatic tool for exploring the molecular mechanisms of bone metastasis in patients with breast cancer and identified factors that may predict the occurrence of bone metastasis.

Iron deficiency triggered transcriptome changes in bread wheat

A series of complex transport, storage and regulation mechanisms control iron metabolism and thereby maintain iron homeostasis in plants. Despite several studies on iron deficiency responses in different plant species, these mechanisms remain unclear in the allohexaploid wheat, which is the most widely cultivated commercial crop. We used RNA sequencing to reveal transcriptomic changes in the wheat flag leaves and roots, when subjected to iron limited conditions. We identified 5969 and 2591 differentially expressed genes (DEGs) in the flag leaves and roots, respectively. Genes involved in the synthesis of iron ligands i.e., nicotianamine (NA) and deoxymugineic acid (DMA) were significantly up-regulated during iron deficiency. In total, 337 and 635 genes encoding transporters exhibited altered expression in roots and flag leaves, respectively. Several genes related to MAJOR FACILITATOR SUPERFAMILY (MFS), ATP-BINDING CASSETTE (ABC) transporter superfamily, NATURAL RESISTANCE ASSOCIATED MACROPHAGE PROTEIN (NRAMP) family and OLIGOPEPTIDE TRANSPORTER (OPT) family were regulated, indicating their important roles in combating iron deficiency stress. Among the regulatory factors, the genes encoding for transcription factors of BASIC HELIX-LOOP-HELIX (bHLH) family were highly up-regulated in both roots and the flag leaves. The jasmonate biosynthesis pathway was significantly altered but with notable expression differences between roots and flag leaves. Homoeologs expression and induction bias analysis revealed subgenome specific differential expression. Our findings provide an integrated overview on regulated molecular processes in response to iron deficiency stress in wheat. This information could potentially serve as a guideline for breeding iron deficiency stress tolerant crops as well as for designing appropriate wheat iron biofortification strategies.

An Evaluation of Machine Learning Approaches for the Prediction of Essential Genes in Eukaryotes Using Protein Sequence-Derived Features

The availability of whole-genome sequences and associated multi-omics data sets, combined with advances in gene knockout and knockdown methods, has enabled large-scale annotation and exploration of gene and protein functions in eukaryotes. Knowing which genes are essential for the survival of eukaryotic organisms is paramount for an understanding of the basic mechanisms of life, and could assist in identifying intervention targets in eukaryotic pathogens and cancer. Here, we studied essential gene orthologs among selected species of eukaryotes, and then employed a systematic machine-learning approach, using protein sequence-derived features and selection procedures, to investigate essential gene predictions within and among species. We showed that the numbers of essential gene orthologs comprise small fractions when compared with the total number of orthologs among the eukaryotic species studied. In addition, we demonstrated that machine-learning models trained with subsets of essentiality-related data performed better than random guessing of gene essentiality for a particular species. Consistent with our gene ortholog analysis, the predictions of essential genes among multiple (including distantly-related) species is possible, yet challenging, suggesting that most essential genes are unique to a species. The present work provides a foundation for the expansion of genome-wide essentiality investigations in eukaryotes using machine learning approaches.

De novo sequencing of Eucommia ulmoides flower bud transcriptomes for identification of genes related to floral development

Eucommia ulmoides Oliver is a woody perennial dioecious species native to China and has great economic value. However, little is known about flower bud development in this species. In this study, the transcriptomes of female and male flower buds were sequenced using the Illumina platform, a next-generation sequencing technology that provides cost-effective, highly efficient transcriptome profiling. In total, 11,558,188,080 clean reads were assembled into 75,065 unigenes with an average length of 1011 bp by de novo assembly using Trinity software. Through similarity comparisons with known protein databases, 47,071 unigenes were annotated, 146 of which were putatively related to the floral development of E. ulmoides. Fifteen of the 146 unigenes had significantly different expression levels between the two samples. Additionally, 24,346 simple sequence repeats were identified in 18,565 unigenes with 12,793 sequences suitable for the designed primers. In total, 67,447 and 58,236 single nucleotide polymorphisms were identified in male and female buds, respectively. This study provides a valuable resource for further conservation genetics and functional genomics research on E. ulmoides.

Adaptation of the targeted capture Methyl-Seq platform for the mouse genome identifies novel tissue-specific DNA methylation patterns of genes involved in neurodevelopment

Methyl-Seq was recently developed as a targeted approach to assess DNA methylation (DNAm) at a genome-wide level in human. We adapted it for mouse and sought to examine DNAm differences across liver and 2 brain regions: cortex and hippocampus. A custom hybridization array was designed to isolate 99 Mb of CpG islands, shores, shelves, and regulatory elements in the mouse genome. This was followed by bisulfite conversion and sequencing on the Illumina HiSeq2000. The majority of differentially methylated cytosines (DMCs) were present at greater than expected frequency in introns, intergenic regions, near CpG islands, and transcriptional enhancers. Liver-specific enhancers were observed to be methylated in cortex, while cortex specific enhancers were methylated in the liver. Interestingly, commonly shared enhancers were differentially methylated between the liver and cortex. Gene ontology and pathway analysis showed that genes that were hypomethylated in the cortex and hippocampus were enriched for neuronal components and neuronal function. In contrast, genes that were hypomethylated in the liver were enriched for cellular components important for liver function. Bisulfite-pyrosequencing validation of 75 DMCs from 19 different loci showed a correlation of r = 0.87 with Methyl-Seq data. We also identified genes involved in neurodevelopment that were not previously reported to be differentially methylated across brain regions. This platform constitutes a valuable tool for future genome-wide studies involving mouse models of disease.

Gene expression profiling data of Schizosaccharomyces pombe under nitrosative stress using differential display

Excess production of nitric oxide (NO) and reactive nitrogen intermediates (RNIs) causes nitrosative stress on cells. Schizosaccharomyces pombe was used as a model to study nitrosative stress response. In the present data article, we have used differential display to identify the differentially expressed genes in the fission yeast under nitrosative stress conditions. We have used pure NO donor compound detaNONOate at final concentrations of 0.1 mM and 1 mM to treat the cells for 15 min alongside control before studying their gene expression profiles. At both the treated conditions, we identified genes which were commonly repressed while several genes were induced upon both 0.1 mM and 1 mM treatments. The differentially expressed genes were further analyzed in DAVID and categorized into several different pathways.