The Diverging Routes of BORIS and CTCF: An Interactomic and Phylogenomic Analysis.

The CCCTC-binding factor (CTCF) is multi-functional, ubiquitously expressed, and highly conserved from Drosophila to human. It has important roles in transcriptional insulation and the formation of a high-dimensional chromatin structure. CTCF has a paralog called "Brother of Regulator of Imprinted Sites" (BORIS) or "CTCF-like" (CTCFL). It binds DNA at sites similar to those of CTCF. However, the expression profiles of the two proteins are quite different. We investigated the evolutionary trajectories of the two proteins after the duplication event using a phylogenomic and interactomic approach. We find that CTCF has 52 direct interaction partners while CTCFL only has 19. Almost all interactors already existed before the emergence of CTCF and CTCFL. The unique secondary loss of CTCF from several nematodes is paralleled by a loss of two of its interactors, the polycomb repressive complex subunit SuZ12 and the multifunctional transcription factor TYY1. In contrast to earlier studies reporting the absence of BORIS from birds, we present evidence for a multigene synteny block containing CTCFL that is conserved in mammals, reptiles, and several species of birds, indicating that not the entire lineage of birds experienced a loss of CTCFL. Within this synteny block, BORIS and its genomic neighbors seem to be partitioned into two nested chromatin loops. The high expression of SPO11, RAE1, RBM38, and PMEPA1 in male tissues suggests a possible link between CTCFL, meiotic recombination, and fertility-associated phenotypes. Using the 65,700 exomes and the 1000 genomes data, we observed a higher number of intergenic, non-synonymous, and loss-of-function mutations in CTCFL than in CTCF, suggesting a reduced strength of purifying selection, perhaps due to less functional constraint.

Genome-wide identification and tissue-specific expression analysis of nucleotide binding site-leucine rich repeat gene family in Cicer arietinum (kabuli chickpea).

The nucleotide binding site-leucine rich repeat (NBS-LRR) proteins play an important role in the defense mechanisms against pathogens. Using bioinformatics approach, we identified and annotated 104 NBS-LRR genes in chickpea. Phylogenetic analysis points to their diversification into two families namely TIR-NBS-LRR and non-TIR-NBS-LRR. Gene architecture revealed intron gain/loss events in this resistance gene family during their independent evolution into two families. Comparative genomics analysis elucidated its evolutionary relationship with other fabaceae species. Around 50% NBS-LRRs reside in macro-syntenic blocks underlining positional conservation along with sequence conservation of NBS-LRR genes in chickpea. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. Four cis-regulatory elements namely WBOX, DRE, CBF, and GCC boxes, that commonly occur in resistance genes, were present in the promoter regions of these genes. Further, the findings will provide a strong background to use candidate disease resistance NBS-encoding genes and identify their specific roles in chickpea.

Genome-wide identification and structure-function studies of proteases and protease inhibitors in Cicer arietinum (chickpea).

BACKGROUND: Proteases are a family of enzymes present in almost all living organisms. In plants they are involved in many biological processes requiring stress response in situations such as water deficiency, pathogen attack, maintaining protein content of the cell, programmed cell death, senescence, reproduction and many more. Similarly, protease inhibitors (PIs) are involved in various important functions like suppression of invasion by pathogenic nematodes, inhibition of spores-germination and mycelium growth of Alternaria alternata and response to wounding and fungal attack. As much as we know, no genome-wide study of proteases together with proteinaceous PIs is reported in any of the sequenced genomes till now. METHODS: Phylogenetic studies and domain analysis of proteases were carried out to understand the molecular evolution as well as gene and protein features. Structural analysis was carried out to explore the binding mode and affinity of PIs for cognate proteases and prolyl oligopeptidase protease with inhibitor ligand. RESULTS: In the study reported here, a significant number of proteases and PIs were identified in chickpea genome. The gene expression profiles of proteases and PIs in five different plant tissues revealed a differential expression pattern in more than one plant tissue. Molecular dynamics studies revealed the formation of stable complex owing to increased number of protein-ligand and inter and intramolecular protein-protein hydrogen bonds. DISCUSSION: The genome-wide identification, characterization, evolutionary understanding, gene expression, and structural analysis of proteases and PIs provide a framework for future analysis when defining their roles in stress response and developing a more stress tolerant variety of chickpea.

Genome-wide identification and tissue-specific expression analysis of UDP-glycosyltransferases genes confirm their abundance in Cicer arietinum (Chickpea) genome.

UDP-glycosyltransferases (EC 2.4.1.x; UGTs) are enzymes coded by an important gene family of higher plants. They are involved in the modification of secondary metabolites, phytohormones, and xenobiotics by transfer of sugar moieties from an activated nucleotide molecule to a wide range of acceptors. This modification regulates various functions like detoxification of xenobiotics, hormone homeostasis, and biosynthesis of secondary metabolites. Here, we describe the identification of 96 UGT genes in Cicer arietinum (CaUGT) and report their tissue-specific differential expression based on publically available RNA-seq and expressed sequence tag data. This analysis has established medium to high expression of 84 CaUGTs and low expression of 12 CaUGTs. We identified several closely related orthologs of CaUGTs in other genomes and compared their exon-intron arrangement. An attempt was made to assign functional specificity to chickpea UGTs by comparing substrate binding sites with experimentally determined specificity. These findings will assist in precise selection of candidate genes for various applications and understanding functional genomics of chickpea.

An improved method for specificity annotation shows a distinct evolutionary divergence among the microbial enzymes of the cholylglycine hydrolase family.

Bile salt hydrolases (BSHs) are gut microbial enzymes that play a significant role in the bile acid modification pathway. Penicillin V acylases (PVAs) are enzymes produced by environmental microbes, having a possible role in pathogenesis or scavenging of phenolic compounds in their microbial habitats. The correct annotation of such physiologically and industrially important enzymes is thus vital. The current methods relying solely on sequence homology do not always provide accurate annotations for these two members of the cholylglycine hydrolase (CGH) family as BSH/PVA enzymes. Here, we present an improved method [binding site similarity (BSS)-based scoring system] for the correct annotation of the CGH family members as BSH/PVA enzymes, which along with the phylogenetic information incorporates the substrate specificity as well as the binding site information. The BSS scoring system was developed through the analysis of the binding sites and binding modes of the available BSH/PVA structures with substrates glycocholic acid and penicillin V. The 198 sequences in the dataset were then annotated accurately using BSS scores as BSH/PVA enzymes. The dataset presented contained sequences from Gram-positive bacteria, Gram-negative bacteria and archaea. The clustering obtained for the dataset using the method described above showed a clear distinction in annotation of Gram-positive bacteria and Gram-negative bacteria. Based on this clustering and a detailed analysis of the sequences of the CGH family in the dataset, we could infer that the CGH genes might have evolved in accordance with the hypothesis stating the evolution of diderms and archaea from the monoderms.

In-silico analysis of binding site features and substrate selectivity in plant flavonoid-3-O glycosyltransferases (F3GT) through molecular modeling, docking and dynamics simulation studies.

Flavonoids are a class of plant secondary metabolites that act as storage molecules, chemical messengers, as well as participate in homeostasis and defense processes. They possess pharmaceutical properties important for cancer treatment such as antioxidant and anti-tumor activities. The drug-related properties of flavonoids can be improved by glycosylation. The enzymes glycosyltransferases (GTs) glycosylate acceptor molecules in a regiospecific manner with the help of nucleotide sugar donor molecules. Several plant GTs have been characterized and their amino acid sequences determined. However, three-dimensional structures of only a few are reported. Here, phylogenetic analysis using amino acid sequences have identified a group of GTs with the same regiospecific activity. The structures of these closely related GTs were modeled using homologous GT structures. Their substrate binding sites were elaborated by docking flavonoid acceptor and UDP-sugar donor molecules in the modeled structures. Eight regions near the acceptor binding site in the N- and C- terminal domain of GTs have been identified that bind and specifically glycosylate the 3-OH group of acceptor flavonoids. Similarly, a conserved motif in the C-terminal domain is known to bind a sugar donor substrate. In certain GTs, the substitution of a specific glutamine by histidine in this domain changes the preference of sugar from glucose to galactose as a result of changed pattern of interactions. The molecular modeling, docking, and molecular dynamics simulation studies have revealed the chemical and topological features of the binding site and thus provided insights into the basis of acceptor and donor recognition by GTs.