Effects of soil ingestion on nutrient digestibility and rumen bacterial diversity of Tibetan sheep.

Tibetan sheep (Ovis aries) are the most numerous livestock in Tibet Plateau pasture ecosystem and have strong ecological adaptability. In the natural grazing system, soil as a natural nutrient carrier and involuntarily or intentionally ingested by Tibetan sheep contribute as an important feed approach. However, quantifying the dosages of soil ingestion for the Tibetan sheep still needs to be clarified. This study aims to characterize nutrient digestibility and rumen bacterial communities by Tibetan sheep in response to different levels of soil ingestion. Thirty sheep were selected and divided into five treatments with soil ingestion (0%, 5%, 10%, 15%, and 20%). The conclusion demonstrated that soil ingestion improved the dry matter digestibility (59.3-62.97%), ether extract (59.79-67.87%) and crude protein (59.81-66.47%) digestibility, particularly 10% soil ingestion has highest nutrient digestibility. The rumen fermentation environment adjusted after soil ingestion by improvement of pH, ammonia nitrogen and volatile fatty acids. Appropriate soil ingestion reduced the bacterial diversity ranged from 946 to 1000 OUTs as compared control (1012), and the rumen bacterial community dominant by typical fiber digestion associated Firmicutes (47.48-53.56%), Bacteroidetes (34.93-40.02%) and Fibrobacteres (4.36-9.27%). Especially, the highest digestible feed capacity and stronger environment adaptability present in 10% soil ingestion Tibetan sheep. Overall, soil ingestion stimulates rumen metabolism by creating a favorable environment for microbial fermentation, improved bacterial community abundance associated with cellulose and saccharide degradation, contribute nutrient digestibility and growth performance of Tibetan sheep.

Self-Organization Formation of Multicellular Spheroids Mediated by Mechanically Tunable Hydrogel Platform: Toward Revealing the Synergy of Chemo- and Noninvasive Photothermal Therapy against Colon Microtumor.

Three-dimensional (3D) tumor cell culture offers a more tissue-recapitulating model in cancer treatment evaluation. However, conventional models based on cell-substrate adhesion deprivation are still of insufficient real tumor mimic. In this work, a novel method is proposed for inducing multicellular spheroids (MCSs) formation based on hydrogel with tunable microenvironmental properties. Colon tumor cells DLD1 cultured on hydrogel substrate with proper physical stimulation form MCSs via self-organization. Chemotherapy based on clinical drug and far-infrared photothermal therapy is evaluated with DLD1 MCSs obtained by this method. The synergism of chemotherapy and noninvasive photothermal therapy based on graphene device is further verified in MCSs model and it is believed this method holds potential in in vitro anti-tumor strategies evaluation for precision medicine.

A Combined Effect of Expression Levels of Obesity-Related Genes and Clinical Factors on Cancer Survival Rate.

Obesity is directly associated with the risk of cancer in different organs, including breast, colon, and kidney. However, adipocytes could be utilized to control progression for some types of cancer, such as leukemia and breast cancer. To explore the potential correlation between adipocytes and cancer, the combined effect of expression levels of obesity-related genes and clinical factors (i.e., gender, race, menopausal status, history of smoking, tumor grade, body mass index (BMI), and history of drinking) on cancer survival rate was systemically studied. The expression levels of obesity-related genes in cancer tissues and normal tissues were downloaded from The Cancer Genome Atlas (TCGA). Kaplan-Meier curves were plotted using R programming language. The log-rank test was applied to explore the correlation between different clinical subgroups. The overexpression of the nine obesity-related genes (MC4R, TMEM18, KCTD15, GNPDA2, SH2B1, MTCH2, FTO, PCSK1, and GPR120) may associate with tumor-promoting factors in some organs (head and neck, gastrointestinal tract, liver, and gallbladder). Underexpressed LEPR, NEGR1, TMEM18, and SH2B1 genes prevented the progression and metastasis of kidney cancer. The combined effect of clinical factors and the expression levels of obesity-related genes on patients' survival was found to be significant. Our outcomes suggested that the alternations of DNA methylation patterns could result in the changes of expression levels of obesity-related genes, playing a critical role in tumor progression. The results of the current study may be utilized to supplement precision and personalized medicine, as well as provide novel insights for the development of treatment approaches for cancer.

Cooperative binding of transcription factors in the human genome.

Transcription factors (TFs) cooperatively bind to specific DNA sequences to control chromatin and gene transcription in eukaryotes. Here, we searched canonical binding, co-binding and tethered binding regions of a TF within ChIP-seq peaks, and investigated the effect of TF-TF cooperation in human GM12878 and K562 cells. We found that TFs except for CTCF and SPI1 showed a large proportion of co-binding and tethered binding regions, and TFs frequently co-binding with other TFs would also frequently tether other TFs to their binding positions. We further observed lower in vivo nucleosome occupancy, higher in vitro nucleosome occupancy and higher levels of H2A.Z, H3K27ac, H3K9ac, H3K4me1, H3K4me2 and H3K4me3 within distal co-binding regions where other TFs were recruited. In addition, target genes for proximal co-binding regions where other TFs were recruited showed significantly higher expression levels. These results indicated that TF-TF cooperation directly associates with the chromatin structure and gene transcription.

miRNA Mediated Noise Making of 3'UTR Mutations in Cancer.

Somatic mutations in 3'-untranslated regions (3'UTR) do not alter amino acids and are considered to be silent in cancers. We found that such mutations can promote tumor progression by altering microRNA (miRNA) targeting efficiency and consequently affecting miRNA⁻mRNA interactions. We identified 67,159 somatic mutations located in the 3'UTRs of messenger RNAs (mRNAs) which can alter miRNA⁻mRNA interactions (functional somatic mutations, funcMutations), and 69.3% of these funcMutations (the degree of energy change > 12 kcal/mol) were identified to significantly promote loss of miRNA-mRNA binding. By integrating mRNA expression profiles of 21 cancer types, we found that the expression of target genes was positively correlated with the loss of absolute affinity level and negatively correlated with the gain of absolute affinity level. Functional enrichment analysis revealed that genes carrying funcMutations were significantly enriched in the MAPK and WNT signaling pathways, and analysis of regulatory modules identified eighteen miRNA modules involved with similar cellular functions. Our findings elucidate a complex relationship between miRNA, mRNA, and mutations, and suggest that 3'UTR mutations may play an important role in tumor development.

Nucleosome Positioning of Intronless Genes in the Human Genome.

Nucleosomes, the basic units of chromatin, are involved in transcription regulation and DNA replication. Intronless genes, which constitute 3 percent of the human genome, differ from intron-containing genes in evolution and function. Our analysis reveals that nucleosome positioning shows a distinct pattern in intronless and intron-containing genes. The nucleosome occupancy upstream of transcription start sites of intronless genes is lower than that of intron-containing genes. In contrast, high occupancy and well positioned nucleosomes are observed along the gene body of intronless genes, which is perfectly consistent with the barrier nucleosome model. Intronless genes have a significantly lower expression level than intron-containing genes and most of them are not expressed in CD4+ T cell lines and GM12878 cell lines, which results from their tissue specificity. However, the highly expressed genes are at the same expression level between the two types of genes. The highly expressed intronless genes require a higher density of RNA Pol II in an elongating state to compensate for the lack of introns. Additionally, 5' and 3' nucleosome depleted regions of highly expressed intronless genes are deeper than those of highly expressed intron-containing genes.

An Integrating Approach for Genome-Wide Screening of MicroRNA Polymorphisms Mediated Drug Response Alterations.

MicroRNAs (miRNAs) are a class of evolutionarily conserved small noncoding RNAs, ~22 nt in length, and found in diverse organisms and play important roles in the regulation of mRNA translation and degradation. It was shown that miRNAs were involved in many key biological processes through regulating the expression of targets. Genetic polymorphisms in miRNA target sites may alter miRNA regulation and therefore result in the alterations of the drug targets. Recent studies have demonstrated that SNPs in miRNA target sites can affect drug efficiency. However, there are still a large number of specific genetic variants related to drug efficiency that are yet to be discovered. We integrated large scale of genetic variations, drug targets, gene interaction networks, biological pathways, and seeds region of miRNA to identify miRNA polymorphisms affecting drug response. In addition, harnessing the abundant high quality biological network/pathways, we evaluated the cascade distribution of tarSNP impacts. We showed that the predictions can uncover most of the known experimentally supported cases as well as provide informative candidates complementary to existing methods/tools. Although there are several existing databases predicting the gain or loss of targeting function of miRNA mediated by SNPs, such as PolymiRTS, miRNASNP, MicroSNiPer, and MirSNP, none of them evaluated the influences of tarSNPs on drug response alterations. We developed a user-friendly online database of this approach named Mir2Drug.

Nucleosome organization in the vicinity of transcription factor binding sites in the human genome.

BACKGROUND: The binding of transcription factors (TFs) to specific DNA sequences is an initial and crucial step of transcription. In eukaryotes, this process is highly dependent on the local chromatin state, which can be modified by recruiting chromatin remodelers. However, previous studies have focused mainly on nucleosome occupancy around the TF binding sites (TFBSs) of a few specific TFs. Here, we investigated the nucleosome occupancy profiles around computationally inferred binding sites, based on 519 TF binding motifs, in human GM12878 and K562 cells. RESULTS: Although high nucleosome occupancy is intrinsically encoded at TFBSs in vitro, nucleosomes are generally depleted at TFBSs in vivo, and approximately a quarter of TFBSs showed well-positioned in vivo nucleosomes on both sides. RNA polymerase near the transcription start site (TSS) has a large effect on the nucleosome occupancy distribution around the binding sites located within one kilobase to the nearest TSS; fuzzier nucleosome positioning was thus observed around these sites. In addition, in contrast to yeast, repressors, rather than activators, were more likely to bind to nucleosomal DNA in the human cells, and nucleosomes around repressor sites were better positioned in vivo. Genes with repressor sites exhibiting well-positioned nucleosomes on both sides, and genes with activator sites occupied by nucleosomes had significantly lower expression, suggesting that actions of activators and repressors are associated with the nucleosome occupancy around their binding sites. It was also interesting to note that most of the binding sites, which were not in the DNase I-hypersensitive regions, were cell-type specific, and higher in vivo nucleosome occupancy were observed at these binding sites. CONCLUSIONS: This study demonstrated that RNA polymerase and the functions of bound TFs affected the local nucleosome occupancy around TFBSs, and nucleosome occupancy patterns around TFBSs were associated with the expression levels of target genes.

Identification of novel microRNAs in primates by using the synteny information and small RNA deep sequencing data.

Current technologies that are used for genome-wide microRNA (miRNA) prediction are mainly based on BLAST tool. They often produce a large number of false positives. Here, we describe an effective approach for identifying orthologous pre-miRNAs in several primates based on syntenic information. Some of them have been validated by small RNA high throughput sequencing data. This approach uses the synteny information and experimentally validated miRNAs of human, and incorporates currently available algorithms and tools to identify the pre-miRNAs in five other primates. First, we identified 929 potential pre-miRNAs in the marmoset in which miRNAs have not yet been reported. Then, we predicted the miRNAs in other primates, and we successfully re-identified most of the published miRNAs and found 721, 979, 650 and 639 new potential pre-miRNAs in chimpanzee, gorilla, orangutan and rhesus macaque, respectively. Furthermore, the miRNA transcriptome in the four primates have been re-analyzed and some novel predicted miRNAs have been supported by the small RNA sequencing data. Finally, we analyzed the potential functions of those validated miRNAs and explored the regulatory elements and transcription factors of some validated miRNA genes of interest. The results show that our approach can effectively identify novel miRNAs and some miRNAs that supported by small RNA sequencing data maybe play roles in the nervous system.

The patterns of histone modifications in the vicinity of transcription factor binding sites in human lymphoblastoid cell lines.

Transcription factor (TF) binding at specific DNA sequences is the fundamental step in transcriptional regulation and is highly dependent on the chromatin structure context, which may be affected by specific histone modifications and variants, known as histone marks. The lack of a global binding map for hundreds of TFs means that previous studies have focused mainly on histone marks at binding sites for several specific TFs. We therefore studied 11 histone marks around computationally-inferred and experimentally-determined TF binding sites (TFBSs), based on 164 and 34 TFs, respectively, in human lymphoblastoid cell lines. For H2A.Z, methylation of H3K4, and acetylation of H3K27 and H3K9, the mark patterns exhibited bimodal distributions and strong pairwise correlations in the 600-bp region around enriched TFBSs, suggesting that these marks mainly coexist within the two nucleosomes proximal to the TF sites. TFs competing with nucleosomes to access DNA at most binding sites, contributes to the bimodal distribution, which is a common feature of histone marks for TF binding. Mark H3K79me2 showed a unimodal distribution on one side of TFBSs and the signals extended up to 4000 bp, indicating a longer-distance pattern. Interestingly, H4K20me1, H3K27me3, H3K36me3 and H3K9me3, which were more diffuse and less enriched surrounding TFBSs, showed unimodal distributions around the enriched TFBSs, suggesting that some TFs may bind to nucleosomal DNA. Besides, asymmetrical distributions of H3K36me3 and H3K9me3 indicated that repressors might establish a repressive chromatin structure in one direction to repress gene expression. In conclusion, this study demonstrated the ranges of histone marks associated with TF binding, and the common features of these marks around the binding sites. These findings have epigenetic implications for future analysis of regulatory elements.