Search

Background Reducing abdominal fat (AF) accumulation and increasing the level of intramuscular fat (IMF) simultaneously is a major breeding goal in poultry industry. Results To explore the different molecular mechanisms between AF and IMF, gene expression profiles in the breast muscle (BM) and AF among three chicken breeds were analyzed. A total of 4737 shared DEGs were identified between BM and AF, of which 2602 DEGs were upregulated and 2135 DEGs were downregulated in the BM groups compared with that in AF groups. DEGs involved in glycerophospholipid metabolism and glycerolipid metabolism were potential regulators resulting in the difference of lipid metabolite accumulation in IMF and AF. PPAR signaling pathway was the most important pathway involved in tissue-specific lipid deposition. Correlation analysis showed that most representative DEGs enriched in PPAR signaling pathway, such as FABP5, PPARG, ACOX1, GK2, were negatively correlated with PUFA-enriched glycerophospholipid molecules. Most DEGs related to glycerophospholipid metabolism, such as GPD2, GPD1, PEMT, CRLS1 and GBGT1, were positively correlated with glycerophospholipid molecules, especially DHA- and arachidonic acid (ARA)- containing glycerophospholipid molecules. Conclusion This study elucidated the molecular mechanism of tissue-specific lipid deposition and poultry meat quality.

SUMMARYSystematic understanding of immune dynamics across the entire human lifespan at single-cell resolution is currently lacking. Here, we performed single-cell RNA sequencing (scRNA-seq) and single-cell T cell receptor (TCR)/B cell receptor (BCR) sequencing (scTCR/BCR-seq) on over 380,000 peripheral blood mononuclear cells collected from 45 healthy participants aged 0 to over 90 years. We revealed that the functions of T cell subsets were most susceptible to senescence among all PBMCs, featured by increased NF-κB signaling and IFN-γ responses pathways, while reduced telomere maintenance and energy metabolism. We subsequently explored the rewiring of cell-cell interactions among different immune cells across the lifespan and revealed the major alteration of immune checkpoints in T cells within the cellular interaction networks. By combined analysis of scRNA-seq and scTCR-seq, we revealed that 1) GNLY+CD8 Effector T cells exhibited a high clonal expansion with distinct functional signatures in children and the elderly; 2) Naïve CD4+T and naïve CD8+T cells displayed different aging patterns in both transcriptomes and immune repertoires; and 3) CD8+MAIT cells showed a higher cell abundance and clonal diversity in adolescents. Furthermore, we identified a unique cytotoxic B cell subset enriched in children by scRNA-seq and scBCR-seq analysis and experimental validations. In summary, our work provided valuable insights and rich resources for understanding the development and aging of the human immune system across the lifespan.

Small RNAs have important functions. However, small RNAs in primate oocytes remain unexplored. Herein, we develop CAS-seq, a single-cell small RNA sequencing method, and profile the small RNAs in human oocytes and embryos. We discover a class of ~20-nt small RNAs that are predominantly expressed in human and monkey oocytes, but not in mouse oocytes. They are specifically associated with HIWI3 (PIWIL3), whereas significantly shorter than the commonly known PIWI-interacting RNAs (piRNAs), designated as oocyte short piRNAs (os-piRNAs). Notably, the os-piRNAs in human oocytes lack 2’-O-methylation at the 3’ end, a hallmark of the classic piRNAs. In addition, the os-piRNAs have a strong 1U/10 A bias and are enriched on the antisense strands of recently evolved transposable elements (TEs), indicating the potential function of silencing TEs by cleavage. Therefore, our study has identified an oocyte-specific piRNA family with distinct features and provides valuable resources for studying small RNAs in primate oocytes., PIWI-interacting RNAs (piRNAs) are ~25–33 nt small RNAs expressed in animal germ cells. Here, the authors develop a single-cell small RNA sequencing method and report that a class of ~20-nt piRNAs lacking 3′ end 2′-O-methylation are associated with PIWIL3 protein and predominantly expressed in human and monkey oocytes.

We systematically investigated the transcriptomes of the peripheral immune cells from 6 inactivated vaccine, BBIBP-CorV recipients at 4 pivotal time points using single-cell RNA-seq technique. First, the significant variation of the canonical immune-responsive signals of both humoral and cellular immunity, as well as other possible symptom-driver signals were evaluated in the specific cell types. Second, we described and compared the common and distinct variation trends across COVID-19 vaccination, disease progression, and flu vaccination to achieve in-depth understandings of the manifestation of immune response in peripheral blood under different stimuli. Third, the expanded T cell and B cell clones were correlated to the specific phenotypes which allowed us to characterize the antigen-specific ones much easier in the future. At last, other than the coagulopathy, the immunogenicity of megakaryocytes in vaccination were highlighted in this study. In brief, our study provided a rich data resource and the related methodology to explore the details of the classical immunity scenarios.

With the discovery of magnetoreceptor mechanisms in animals, it materialized the novel applications of controlling cell and animal behaviors using magnetic fields. T cells have shown to be sensitive to magnetic fields. Here, we reported that exposure to moderate SMFs (static magnetic fields) led to increased granule and cytokine secretion as well as ATP production and mitochondrial respiration from CD8+ T cells. These effects were inhibited by knocking down the Uqcrb and Ndufs6 genes of mitochondrial respiratory chain, whose transcriptions were regulated by candidate magnetoreceptor genes Isca1 and Cry1/Cry2. SMF exposure also promoted CD8+ T cell granule and cytokine secretion and repressed tumor growth in vivo. SMFs enhanced CD8+ T cell cytotoxicity, and the adoptive transfer into tumor-bearing mice resulted in enhanced antitumor effects. Collectively, our study suggests that moderate SMFs enhance CD8+ T cell cytotoxicity by promoting mitochondrial respiration and promoted the antitumor function of CD8+ T cells.

The Tudor domain-containing (Tdrd) family proteins play a critical role in transposon silencing in animal gonads by recognizing the symmetrically dimethylated arginine (sDMA) on the (G/A)R motif of the N-terminal of PIWI family proteins via the eTud domains. Papi, also known as “Tdrd2,” is involved in Zucchini-mediated PIWI-interacting RNA (piRNA) 3′-end maturation. Intriguingly, a recent study showed that, in papi mutant flies, only Piwi-bound piRNAs increased in length, and not Ago3-bound or Aub-bound piRNAs. However, the molecular and structural basis of the Papi–Piwi complex is still not fully understood, which limits mechanistic understanding of the function of Papi in piRNA biogenesis. In the present study, we determined the crystal structures of Papi-eTud in the apo form and in complex with a peptide containing unmethylated or dimethylated R10 residues. Structural and biochemical analysis showed that the Papi interaction region on the Drosophila Piwi contains an RGRRR motif (R7–R11) distinct from the consensus (G/A)R motif recognized by canonical eTud. Mass spectrometry results indicated that Piwi is the major binding partner of Papi in vivo. The papi mutant flies suffered from both fertility and transposon-silencing defects, supporting the important role conferred to Papi in piRNA 3′ processing through direct interaction with Piwi proteins.

In this study, we identified the direct interaction region between Drosophila Piwi and Papi. We further determined the crystal structures of Papi-eTud in the apo form, in complex with unmethylated Piwi peptide, and in complex with symmetrically dimethylated Piwi peptide at arginine-10, which demonstrated that Papi interacts with an RGRRR motif on the N terminus of Piwi in a sequence-specific manner both in vitro and in vivo. This recognition sequence, which determines the specificity of Papi–Piwi interactions, is different from all previously reported (G/A)R repeats. Our studies provide mechanistic insights into the important role of Papi–Piwi interactions in the 3′ end-trimming process of PIWI-interacting RNA biogenesis and facilitate the identification of new PIWI-interacting partners of Tudor domain-containing proteins.

Invariant natural killer T cells (iNKT cells) are a specific subset of T cells that recognize glycolipid antigens and upon activation rapidly exert effector functions. This unique function is established during iNKT cell development; the detailed mechanisms of this process, however, remain to be elucidated. Here the authors show that deletion of the mediator subunit Med23 in CD4+CD8+ double positive (DP) thymocytes completely blocks iNKT cell development at stage 2. This dysregulation is accompanied by a bias in the expression of genes related to the regulation of transcription and metabolism, and functional impairment of the cells including the loss of NK cell characteristics, reduced ability to secrete cytokines and attenuated recruitment capacity upon activation. Moreover, Med23-deficient iNKT cells exhibit impaired anti-tumor activity. Our study identifies Med23 as an essential transcriptional regulator that controls iNKT cell differentiation and terminal maturation.

Recent studies have suggested that long non-coding RNAs (lncRNAs) can interact with microRNAs (miRNAs) and indirectly regulate miRNA targets though competing interactions. However, the molecular mechanisms underlying these interactions are still largely unknown. In this study, these lncRNA-miRNA-gene interactions were defined as lncRNA-associated competing triplets (LncACTs), and an integrated pipeline was developed to identify lncACTs that are active in cancer. Competing lncRNAs had sponge features distinct from non-competing lncRNAs. In the lncACT cross-talk network, disease-associated lncRNAs, miRNAs and coding-genes showed specific topological patterns indicative of their competence and control of communication within the network. The construction of global competing activity profiles revealed that lncACTs had high activity specific to cancers. Analyses of clustered lncACTs revealed that they were enriched in various cancer-related biological processes. Based on the global cross-talk network and cluster analyses, nine cancer-specific sub-networks were constructed. H19- and BRCA1/2-associated lncACTs were able to discriminate between two groups of patients with different clinical outcomes. Disease-associated lncACTs also showed variable competing patterns across normal and cancer patient samples. In summary, this study uncovered and systematically characterized global properties of human lncACTs that may have prognostic value for predicting clinical outcome in cancer patients.

Single cell RNA-seq is a powerful and sensitive way to capture the genome-wide gene expression. Here, single cell RNA-seq was utilized to study the transcriptomic profile of early zebrafish PGCs (primordial germ cells) at three different developmental stages. The three stages were 6, 11 and 24 hpf (hours post fertilization). For each developmental stage, three zebrafish PGCs from one embryo were collected, and 9 samples in total were used in this experiment. Single cell RNA-seq results showed that 5099-7376 genes were detected among the 9 samples, and the number of expressed genes decreased as development progressed. Based on the gene expression pattern, samples from 6 and 11 hpf clustered closely, while samples from 24 hpf were more dispersed. By WGCNA (weighted gene co-expression network analysis), the two biggest modules that had inverse gene expression patterns were found to be related to PGC formation or migration. Functional enrichment analysis for these two modules showed that PGCs mainly conducted migration and cell division in early development (6/11 hpf) and translation activity became active in late development (24 hpf). Differentially expressed gene analyses showed that more genes were downregulated than upregulated between two adjacent stages, and genes related to PGC formation or migration reported by previous studies decreased significantly from 11 to 24 hpf. Our results provide base knowledge about zebrafish PGC development at the single cell level and can be further studied by other researchers interested in biological development.

Identification of unknown metabolites is a major challenge in metabolomics. Without the identities of the metabolites, the metabolome data generated from a biological sample cannot be readily linked with the proteomic and genomic information for studies in systems biology and medicine. We have developed a web-based metabolite identification tool ( http://www.mycompoundid.org ) that allows searching and interpreting mass spectrometry (MS) data against a newly constructed metabolome library composed of 8,021 known human endogenous metabolites and their predicted metabolic products (375,809 compounds from one metabolic reaction and 10,583,901 from two reactions). As an example, in the analysis of a simple extract of human urine or plasma and the whole human urine by liquid chromatography-mass spectrometry and MS/MS, we are able to identify at least two times more metabolites in these samples than by using a standard human metabolome library. In addition, it is shown that the evidence-based metabolome library (EML) provides a much superior performance in identifying putative metabolites from a human urine sample, compared to the use of the ChemPub and KEGG libraries.

Expression, modification, and activity of microRNAs are dynamically regulated in early mouse embryos unveiled by sensitive sequencing., Small RNAs play important roles in early embryonic development. However, their expression dynamics and modifications are poorly understood because of the scarcity of RNA that is obtainable for sequencing analysis. Using an improved deep sequencing method that requires as little as 10 ng of total RNA or 50 oocytes, we profile small RNAs in mouse oocytes and early embryos. We find that microRNA (miRNA) expression starts soon after fertilization, and the mature miRNAs carried into the zygote by sperm during fertilization are relatively rare compared to the oocyte miRNAs. Intriguingly, the zygotic miRNAs display a marked increase in 3′ mono- and oligoadenylation in one- to two-cell embryos, which may protect the miRNAs from the massive degradation taking place during that time. Moreover, bioinformatics analyses show that the function of miRNA is suppressed from the oocyte to the two-cell stage and appears to be reactivated after the two-cell stage to regulate genes important in embryonic development. Our study thus provides a highly sensitive profiling method and valuable data sets for further examination of small RNAs in early embryos.

The significant roles of genetic variants in myasthenia gravis (MG) pathogenesis have been demonstrated in many studies, and recently it has been revealed that aberrant level/regulation of microRNAs (miRNAs) might contribute to the initiation and progression of MG. However, the dysfunction of miRNA associated with single nucleotide polymorphisms (miRSNPs) has not been well investigated in MG. In this study, we created a contemporary catalog of 89 MG risk genes via manual literature-mining. Based on this risk gene catalog, we obtained 18 MG risk pathways. Furthermore, we identified 93 miRNAs that target MG risk pathways and revealed the miRSNPs 'switches' in miRNA regulation in the MG risk pathways by integrating the database information of miRSNPs. We also constructed a miRNA-mediated SNP switching pathway network (MSSPN) to intuitively analyze miRNA regulation of MG risk pathways and the relationship of the polymorphism 'switch' with these changes in regulation. Moreover, we carried out in-depth dissection on the correlation between hsa05200 (pathway in cancer) and MG development, and elaborated the significance of 4 high-risk genes. By network analysis and literature mining, we proposed a potential mechanism of miRSNPs→gene→pathway effects on MG pathogenesis, especially for rs28457673 (miR-15/16/195/424/497 family)→IGF1R→hsa05200 (pathway in cancer). Therefore, our studies have revealed a functional role for genetic modulators in MG pathogenesis at a systemic level, which could be informative for further miRNA and miRSNPs studies in MG.

Small peptides, such as dipeptides and tripeptides, are naturally present in many biological samples (e.g., human biofluids and cell extracts). They have attracted great attention in many research fields because of their important biological functions as well as potential roles as disease biomarkers. Tandem mass spectrometry (MS/MS) can be used to profile these small peptides. However, the type and number of fragment ions generated in MS/MS are often limited for unambiguous identification. Herein we report a novel database-search strategy based on the use of MS/MS spectra of both unlabeled and dimethyl labeled peptides to identify and confirm amino acid sequences of di/tripeptides that are separated using hydrophilic interaction (HILIC) liquid chromatography (LC). To facilitate the di/tripeptide identification, a database consisting of all the predicted MS/MS spectra from 400 dipeptides and 8000 tripeptides was created, and a search tool, PEP Search, was developed and housed at the MyCompoundID website ( www.mycompoundid.org/PEP). To evaluate the identification specificity of this method, we used acid hydrolysis to degrade a standard protein, cytochrome c, to produce many di/tripeptides with known sequences for LC/MS/MS. The resultant MS/MS spectra were searched against the database to generate a list of matches which were compared to the known sequences. We correctly identified the di/tripeptides in the protein hydrolysate. We then applied this method to detect and identify di/tripeptides naturally present in human urine samples with high confidence. We envisage the use of this method as a complementary tool to various LC/MS techniques currently available for small molecule or metabolome profiling with an added benefit of covering all di/tripeptide chemical space.

Genome-wide association studies (GWAS) have successfully identified a large number of single nucleotide polymorphisms (SNPs) that are associated with a wide range of human diseases. However, many of these disease-associated SNPs are located in non-coding regions and have remained largely unexplained. Recent findings indicate that disease-associated SNPs in human large intergenic non-coding RNA (lincRNA) may lead to susceptibility to diseases through their effects on lincRNA expression. There is, therefore, a need to specifically record these SNPs and annotate them as potential candidates for disease. We have built LincSNP, an integrated database, to identify and annotate disease-associated SNPs in human lincRNAs. The current release of LincSNP contains approximately 140,000 disease-associated SNPs (or linkage disequilibrium SNPs), which can be mapped to around 5,000 human lincRNAs, together with their comprehensive functional annotations. The database also contains annotated, experimentally supported SNP-lincRNA-disease associations and disease-associated lincRNAs. It provides flexible search options for data extraction and searches can be performed by disease/phenotype name, SNP ID, lincRNA name and chromosome region. In addition, we provide users with a link to download all the data from LincSNP and have developed a web interface for the submission of novel identified SNP-lincRNA-disease associations. The LincSNP database aims to integrate disease-associated SNPs and human lincRNAs, which will be an important resource for the investigation of the functions and mechanisms of lincRNAs in human disease. The database is available at http://bioinfo.hrbmu.edu.cn/LincSNP .

MicroRNAs (miRNAs) are a class of small (19–25 nt) non-coding RNAs. This important class of gene regulator downregulates gene expression through sequence-specific binding to the 3′untranslated regions (3′UTRs) of target mRNAs. Several computational target prediction approaches have been developed for predicting miRNA targets. However, the predicted target lists often have high false positive rates. To construct a workable target list for subsequent experimental studies, we need novel approaches to properly rank the candidate targets from traditional methods. We performed a systematic analysis of experimentally validated miRNA targets using functional genomics data, and found significant functional associations between genes that were targeted by the same miRNA. Based on this finding, we developed a miRNA target prioritization method named mirTarPri to rank the predicted target lists from commonly used target prediction methods. Leave-one-out cross validation has proved to be successful in identifying known targets, achieving an AUC score up to 0. 84. Validation in high-throughput data proved that mirTarPri was an unbiased method. Applying mirTarPri to prioritize results of six commonly used target prediction methods allowed us to find more positive targets at the top of the prioritized candidate list. In comparison with other methods, mirTarPri had an outstanding performance in gold standard and CLIP data. mirTarPri was a valuable method to improve the efficacy of current miRNA target prediction methods. We have also developed a web-based server for implementing mirTarPri method, which is freely accessible at http://bioinfo.hrbmu.edu.cn/mirTarPri.

Drug design is crucial for the effective discovery of anti-cancer drugs. The success or failure of drug design often depends on the leading compounds screened in pre-clinical studies. Many efforts, such as in vivo animal experiments and in vitro drug screening, have improved this process, but these methods are usually expensive and laborious. In the post-genomics era, it is possible to seek leading compounds for large-scale candidate small-molecule screening with multiple OMICS datasets. In the present study, we developed a computational method of prioritizing small molecules as leading compounds by integrating transcriptomics and toxicogenomics data. This method provides priority lists for the selection of leading compounds, thereby reducing the time required for drug design. We found 11 known therapeutic small molecules for breast cancer in the top 100 candidates in our list, 2 of which were in the top 10. Furthermore, another 3 of the top 10 small molecules were recorded as closely related to cancer treatment in the DrugBank database. A comparison of the results of our approach with permutation tests and shared gene methods demonstrated that our OMICS data-based method is quite competitive. In addition, we applied our method to a prostate cancer dataset. The results of this analysis indicated that our method surpasses both the shared gene method and random selection. These analyses suggest that our method may be a valuable tool for directing experimental studies in cancer drug design, and we believe this time- and cost-effective computational strategy will be helpful in future studies in cancer therapy.

Large intergenic noncoding RNAs (lincRNAs) are emerging as key factors of multiple cellular processes. Cumulative evidence has linked lincRNA polymorphisms to diverse diseases. However, the global properties of lincRNA polymorphisms and their implications for human disease remain largely unknown. Here we performed a systematic analysis of naturally occurring variants in human lincRNAs, with a particular focus on lincRNA polymorphism as novel risk factor of disease etiology. We found that lincRNAs exhibited a relatively low level of polymorphisms, and low single-nucleotide polymorphism (SNP) density lincRNAs might have a broad range of functions. We also found that some polymorphisms in evolutionarily conserved regions of lincRNAs had significant effects on predicted RNA secondary structures, indicating their potential contribution to diseases. We mapped currently available phenotype-associated SNPs to lincRNAs and found that lincRNAs were associated with a wide range of human diseases. Some lincRNAs could be responsible for particular diseases. Our results provided not only a global perspective on genetic variants in human lincRNAs but also novel insights into the function and etiology of lincRNA. All the data in this study can be accessed and retrieved freely via a web server at http://bioinfo.hrbmu.edu.cn/lincPoly.

The application of microRNAs (miRNAs) in the therapeutics of glioma and other human diseases is an area of intense interest. However, it’s still a great challenge to interpret the functional consequences of using miRNAs in glioma therapy. Here, we examined paired deep sequencing expression profiles of miRNAs and mRNAs from human glioma cell lines after manipulating the levels of miRNAs miR-181d, -21, and -23b, as well as transcriptional regulators β-catenin, CBP, and STAT3. An integrated approach was used to identify functional miRNA-pathway regulatory networks (MPRNs) responding to each manipulation. MiRNAs were identified to regulate glioma related biological pathways collaboratively after manipulating the level of either post-transcriptional or transcriptional regulators, and functional synergy and crosstalk was observed between different MPRNs. MPRNs responsive to multiple interventions were found to occupy central positions in the comprehensive MPRN (cMPRN) generated by integrating all the six MPRNs. Finally, we identified a core module comprising 14 miRNAs and five pathways that could predict the survival of glioma patients and represent potential targets for glioma therapy. Our results provided novel insight into miRNA regulatory mechanisms implicated in therapeutic interventions and could offer more inspiration to miRNA-based glioma therapy.