BackgroundIn combination with gene expression profiles, the protein interaction network (PIN) constructs a dynamic network that includes multiple functional modules. Previous studies have demonstrated that rifampin can influence drug metabolism by regulating drug-metabolizing enzymes, transporters, and microRNAs (miRNAs). Rifampin induces gene expression, at least in part, by activating the pregnane X receptor (PXR), which induces gene expression; however, the impact of rifampin on global gene regulation has not been examined under the molecular network frameworks.MethodsIn this study, we extracted rifampin-induced significant differentially expressed genes (SDG) based on the gene expression profile. By integrating the SDG and human protein interaction network (HPIN), we constructed the rifampin-regulated protein interaction network (RrPIN). Based on gene expression measurements, we extracted a subnetwork that showed enriched changes in molecular activity. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG), we identified the crucial rifampin-regulated biological pathways and associated genes. In addition, genes targeted by miRNAs that were significantly differentially expressed in the miRNA expression profile were extracted based on the miRNA-gene prediction tools. The miRNA-regulated PIN was further constructed using associated genes and miRNAs. For each miRNA, we further evaluated the potential impact by the gene interaction network using pathway analysis.Results and DisccussionWe extracted the functional modules, which included 84 genes and 89 interactions, from the RrPIN, and identified 19 key rifampin-response genes that are associated with seven function pathways that include drug response and metabolism, and cancer pathways; many of the pathways were supported by previous studies. In addition, we identified that a set of 6 genes (CAV1, CREBBP, SMAD3, TRAF2, KBKG, and THBS1) functioning as gene hubs in the subnetworks that are regulated by rifampin. It is also suggested that 12 differentially expressed miRNAs were associated with 6 biological pathways.ConclusionsOur results suggest that rifampin contributes to changes in the expression of genes by regulating key molecules in the protein interaction networks. This study offers valuable insights into rifampin-induced biological mechanisms at the level of miRNAs, genes and proteins.

BackgroundBidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.ResultsHere, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.ConclusionsOur analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

BackgroundIn combination with gene expression profiles, the protein interaction network (PIN) constructs a dynamic network that includes multiple functional modules. Previous studies have demonstrated that rifampin can influence drug metabolism by regulating drug-metabolizing enzymes, transporters, and microRNAs (miRNAs). Rifampin induces gene expression, at least in part, by activating the pregnane X receptor (PXR), which induces gene expression; however, the impact of rifampin on global gene regulation has not been examined under the molecular network frameworks.MethodsIn this study, we extracted rifampin-induced significant differentially expressed genes (SDG) based on the gene expression profile. By integrating the SDG and human protein interaction network (HPIN), we constructed the rifampin-regulated protein interaction network (RrPIN). Based on gene expression measurements, we extracted a subnetwork that showed enriched changes in molecular activity. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG), we identified the crucial rifampin-regulated biological pathways and associated genes. In addition, genes targeted by miRNAs that were significantly differentially expressed in the miRNA expression profile were extracted based on the miRNA-gene prediction tools. The miRNA-regulated PIN was further constructed using associated genes and miRNAs. For each miRNA, we further evaluated the potential impact by the gene interaction network using pathway analysis.Results and DisccussionWe extracted the functional modules, which included 84 genes and 89 interactions, from the RrPIN, and identified 19 key rifampin-response genes that are associated with seven function pathways that include drug response and metabolism, and cancer pathways; many of the pathways were supported by previous studies. In addition, we identified that a set of 6 genes (CAV1, CREBBP, SMAD3, TRAF2, KBKG, and THBS1) functioning as gene hubs in the subnetworks that are regulated by rifampin. It is also suggested that 12 differentially expressed miRNAs were associated with 6 biological pathways.ConclusionsOur results suggest that rifampin contributes to changes in the expression of genes by regulating key molecules in the protein interaction networks. This study offers valuable insights into rifampin-induced biological mechanisms at the level of miRNAs, genes and proteins.

BackgroundBidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.ResultsHere, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.ConclusionsOur analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

BackgroundThe P1B-ATPase subfamily is an important group involved in transporting heavy metals and has been extensively studied in model plants, such as Arabidopsis thaliana and Oryza sativa. Emerging evidence indicates that one homolog in Glycine max is also involved in cadmium (Cd) stress, but the gene family has not been fully investigated in soybean.ResultsHere, we identified 20 heavy metal ATPase (HMA) family members in the soybean genome, presented as 10 paralogous pairs, which is significantly greater than the number in Arabidopsis or rice, and was likely caused by the latest whole genome duplication event in soybean. A phylogenetic analysis divided the 20 members into six groups, each having conserved or divergent gene structures and protein motif patterns. The integration of RNA-sequencing and qRT-PCR data from multiple tissues provided an overall expression pattern for the HMA family in soybean. Further comparisons of expression patterns and the single nucleotide polymorphism distribution between paralogous pairs suggested functional conservation and the divergence of HMA genes during soybean evolution. Finally, analyses of the HMAs expressed in response to Cd stress provided evidence on how plants manage Cd tolerance, at least in the two contrasting soybean genotypes examined.ConclusionsThe genome-wide identification, chromosomal distribution, gene structures, and evolutionary and expression analyses of the 20 HMA genes in soybean provide an overall insight into their potential involvement in Cd responses. These results will facilitate further research on the HMA gene family, and their conserved and divergent biological functions in soybean.

BackgroundThe P1B-ATPase subfamily is an important group involved in transporting heavy metals and has been extensively studied in model plants, such as Arabidopsis thaliana and Oryza sativa. Emerging evidence indicates that one homolog in Glycine max is also involved in cadmium (Cd) stress, but the gene family has not been fully investigated in soybean.ResultsHere, we identified 20 heavy metal ATPase (HMA) family members in the soybean genome, presented as 10 paralogous pairs, which is significantly greater than the number in Arabidopsis or rice, and was likely caused by the latest whole genome duplication event in soybean. A phylogenetic analysis divided the 20 members into six groups, each having conserved or divergent gene structures and protein motif patterns. The integration of RNA-sequencing and qRT-PCR data from multiple tissues provided an overall expression pattern for the HMA family in soybean. Further comparisons of expression patterns and the single nucleotide polymorphism distribution between paralogous pairs suggested functional conservation and the divergence of HMA genes during soybean evolution. Finally, analyses of the HMAs expressed in response to Cd stress provided evidence on how plants manage Cd tolerance, at least in the two contrasting soybean genotypes examined.ConclusionsThe genome-wide identification, chromosomal distribution, gene structures, and evolutionary and expression analyses of the 20 HMA genes in soybean provide an overall insight into their potential involvement in Cd responses. These results will facilitate further research on the HMA gene family, and their conserved and divergent biological functions in soybean.