RBPLight: a computational tool for discovery of plant-specific RNA-binding proteins using light gradient boosting machine and ensemble of evolutionary features

RNA-binding proteins (RBPs) are essential for post-transcriptional gene regulation in eukaryotes, including splicing control, mRNA transport and decay. Thus, accurate identification of RBPs is important to understand gene expression and regulation of cell state. In order to detect RBPs, a number of computational models have been developed. These methods made use of datasets from several eukaryotic species, specifically from mice and humans. Although some models have been tested on Arabidopsis, these techniques fall short of correctly identifying RBPs for other plant species. Therefore, the development of a powerful computational model for identifying plant-specific RBPs is needed. In this study, we presented a novel computational model for locating RBPs in plants. Five deep learning models and ten shallow learning algorithms were utilized for prediction with 20 sequence-derived and 20 evolutionary feature sets. The highest repeated five-fold cross-validation accuracy, 91.24% AU-ROC and 91.91% AU-PRC, was achieved by light gradient boosting machine. While evaluated using an independent dataset, the developed approach achieved 94.00% AU-ROC and 94.50% AU-PRC. The proposed model achieved significantly higher accuracy for predicting plant-specific RBPs as compared to the currently available state-of-art RBP prediction models. Despite the fact that certain models have already been trained and assessed on the model organism Arabidopsis, this is the first comprehensive computer model for the discovery of plant-specific RBPs. The web server RBPLight was also developed, which is publicly accessible at https://iasri-sg.icar.gov.in/rbplight/, for the convenience of researchers to identify RBPs in plants.

DBPMod: a supervised learning model for computational recognition of DNA-binding proteins in model organisms

DNA-binding proteins (DBPs) play critical roles in many biological processes, including gene expression, DNA replication, recombination and repair. Understanding the molecular mechanisms underlying these processes depends on the precise identification of DBPs. In recent times, several computational methods have been developed to identify DBPs. However, because of the generic nature of the models, these models are unable to identify species-specific DBPs with higher accuracy. Therefore, a species-specific computational model is needed to predict species-specific DBPs. In this paper, we introduce the computational DBPMod method, which makes use of a machine learning approach to identify species-specific DBPs. For prediction, both shallow learning algorithms and deep learning models were used, with shallow learning models achieving higher accuracy. Additionally, the evolutionary features outperformed sequence-derived features in terms of accuracy. Five model organisms, including Caenorhabditis elegans, Drosophila melanogaster, Escherichia coli, Homo sapiens and Mus musculus, were used to assess the performance of DBPMod. Five-fold cross-validation and independent test set analyses were used to evaluate the prediction accuracy in terms of area under receiver operating characteristic curve (auROC) and area under precision-recall curve (auPRC), which was found to be ~89-92% and ~89-95%, respectively. The comparative results demonstrate that the DBPMod outperforms 12 current state-of-the-art computational approaches in identifying the DBPs for all five model organisms. We further developed the web server of DBPMod to make it easier for researchers to detect DBPs and is publicly available at https://iasri-sg.icar.gov.in/dbpmod/. DBPMod is expected to be an invaluable tool for discovering DBPs, supplementing the current experimental and computational methods.

Population structure and genetic diversity of hatchery stocks as revealed by combined mtDNA fragment sequences in Indian major carp, Catla catla

Catla catla is the second most important Indian major carp due to high growth rate and acceptance to consumers for food value. It is widely cultured in the Indian subcontinent as monoculture or polyculture. In the present study, genetic diversity among hatchery stocks (total 218 samples of catla) collected from different geographical regions of India was examined using mtDNA fragment sequence of Cyt b (306 bp) and D loop (710 bp). High numbers (57) of population specific haplotypes were observed in the present study. The results revealed significant genetic heterogeneity for the sequence data (F_ST = 0.27546, p < .05). Analysis of molecular variance revealed significant genetic differentiation among different catla populations. The information generated in present study could be useful to develop broad genetic base populations of catla.

nifPred: Proteome-Wide Identification and Categorization of Nitrogen-Fixation Proteins of Diaztrophs Based on Composition-Transition-Distribution Features Using Support Vector Machine

As inorganic nitrogen compounds are essential for basic building blocks of life (e.g., nucleotides and amino acids), the role of biological nitrogen-fixation (BNF) is indispensible. All nitrogen fixing microbes rely on the same nitrogenase enzyme for nitrogen reduction, which is in fact an enzyme complex consists of as many as 20 genes. However, the occurrence of six genes viz., nifB, nifD, nifE, nifH, nifK, and nifN has been proposed to be essential for a functional nitrogenase enzyme. Therefore, identification of these genes is important to understand the mechanism of BNF as well as to explore the possibilities for improving BNF from agricultural sustainability point of view. Further, though the computational tools are available for the annotation and phylogenetic analysis of nifH gene sequences alone, to the best of our knowledge no tool is available for the computational prediction of the above mentioned six categories of nitrogen-fixation (nif) genes or proteins. Thus, we proposed an approach, which is first of its kind for the computational identification of nif proteins encoded by the six categories of nif genes. Sequence-derived features were employed to map the input sequences into vectors of numeric observations that were subsequently fed to the support vector machine as input. Two types of classifier were constructed: (i) a binary classifier for classification of nif and non-nitrogen-fixation (non-nif) proteins, and (ii) a multi-class classifier for classification of six categories of nif proteins. Higher accuracies were observed for the combination of composition-transition-distribution (CTD) feature set and radial kernel, as compared to the other feature-kernel combinations. The overall accuracies were observed >90% in both binary and multi-class classifications. The developed approach further achieved >92% accuracy, while evaluated with blind (independent) test datasets. The developed approach also produced higher accuracy in identifying nif proteins, while evaluated using proteome-wide datasets of several species. Furthermore, we established a prediction server nifPred (http://webapp.cabgrid.res.in/nifPred) to assist the scientific community for proteome-wide identification of six categories of nif proteins. Besides, the source code of nifPred is also available at https://github.com/PrabinaMeher/nifPred. The developed web server is expected to supplement the transcriptional profiling and comparative genomics studies for the identification and functional annotation of genes related to BNF.

ir-HSP: Improved Recognition of Heat Shock Proteins, Their Families and Sub-types Based On g-Spaced Di-peptide Features and Support Vector Machine

Heat shock proteins (HSPs) play a pivotal role in cell growth and variability. Since conventional approaches are expensive and voluminous protein sequence information is available in the post-genomic era, development of an automated and accurate computational tool is highly desirable for prediction of HSPs, their families and sub-types. Thus, we propose a computational approach for reliable prediction of all these components in a single framework and with higher accuracy as well. The proposed approach achieved an overall accuracy of ~84% in predicting HSPs, ~97% in predicting six different families of HSPs, and ~94% in predicting four types of DnaJ proteins, with bench mark datasets. The developed approach also achieved higher accuracy as compared to most of the existing approaches. For easy prediction of HSPs by experimental scientists, a user friendly web server ir-HSP is made freely accessible at http://cabgrid.res.in:8080/ir-hsp. The ir-HSP was further evaluated for proteome-wide identification of HSPs by using proteome datasets of eight different species, and ~50% of the predicted HSPs in each species were found to be annotated with InterPro HSP families/domains. Thus, the developed computational method is expected to supplement the currently available approaches for prediction of HSPs, to the extent of their families and sub-types.

Rapid recovery of complete mitogenome of Indian major carp, Catla catla from low depth paired end Illumina sequencing

Here we report the reconstruction of the catla (Catla catla) complete mitochondrial genome sequence from low depth paired end Illumina sequencing. The genome is of 16,597 bp in size. Similar to other vertebrate mtgenomes, it consists of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and a putative control region. The present mtgenome is 3 bp longer than the earlier reported catla mtgenome from our laboratory. Majority of the mitochondrial genes are encoded by the H-strand. Phylogenetics analysis revealed that Catla catla is closer to Labeo rohita than other labeo species. Present study demonstrated the power of next generation sequencing towards hassle free and rapid sequencing of mitochondrial genomes of non-model organisms.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

Low-depth shotgun sequencing resolves complete mitochondrial genome sequence of Labeo rohita

Labeo rohita, popularly known as rohu, is a widely cultured species in whole Indian subcontinent. In the present study, we used in-silico approach to resolve complete mitochondrial genome of rohu. Low-depth shotgun sequencing using Roche 454 GS FLX (Branford, Connecticut, USA) followed by de novo assembly in CLC Genomics Workbench version 7.0.4 (Aarhus, Denmark) revealed the complete mitogenome of L. rohita to be 16 606 bp long (accession No. KR185963). It comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 putative control region. The gene order and organization are similar to most vertebrates. The mitogenome in the present investigation has 99% similarity with that of previously reported mitogenomes of rohu and this is also evident from the phylogenetic study using maximum-likelihood (ML) tree method. This study was done to determine the feasibility, accuracy and reliability of low-depth sequence data obtained from NGS platform as compared to the Sanger sequencing. Thus, NGS technology has proven to be competent and a rapid in-silico alternative to resolve the complete mitochondrial genome sequence, thereby reducing labors and time.

New encoded single-indicator sequences based on physico-chemical parameters for efficient exon identification

The first step in gene identification problem based on genomic signal processing is to convert character strings into numerical sequences. These numerical sequences are then analysed spectrally or using digital filtering techniques for the period-3 peaks, which are present in exons (coding areas) and absent in introns (non-coding areas). In this paper, we have shown that single-indicator sequences can be generated by encoding schemes based on physico-chemical properties. Two new methods are proposed for generating single-indicator sequences based on hydration energy and dipole moments. The proposed methods produce high peak at exon locations and effectively suppress false exons (intron regions having greater peak than exon regions) resulting in high discriminating factor, sensitivity and specificity.

Derivation and characterization of embryonic stem-like cells of Indian major carp Catla catla

Embryonic stem (ES)-like cells were derived from mid-blastula stage embryos of a freshwater fish, catla Catla catla, under feeder-free condition and designated as CCES cells. The conditioned media was optimized with 10% foetal bovine serum (FBS), fish embryo extract (FEE) having 100 µg ml(-1) protein concentration, 15 ng ml(-1) basic fibroblast growth factor (bFGF) and basic media containing Leibovitz-15, DMEM with 4·5 g l(-1) glucose and Ham's F12 (LDF) in 2:1:1 ratio using a primary culture of CCES cells. Cells attached to gelatin-coated plates after 24 h of seeding and ES-like colonies were obtained at day 5 onwards. A stable cell culture was obtained after passage 10 and further maintained up to passage 44. These cells were characterized by their typical morphology, high alkaline phosphatase activity, positive expression of cell-surface antigen SSEA-1, transcription factor Oct4, germ cell marker vasa and consistent karyotype up to extended periods. The undifferentiated state was confirmed by their ability to form embryoid bodies and their differentiation potential.

Effect of endotoxin on the immunity of Indian major carp, Labeo rohita

Endotoxin, a lipopolysaccharide component of outer cell wall membrane of the Gram-negative bacteria is a factor responsible for a number of biological effects including immunostimulatory activities in different animal species including fish. In this study, L. rohita yearlings of weight ranging from 80 to 100g were injected intraperitoneally with 0.5, 1, 2, 5, 10 and 20 EU/fish dose of endotoxin to find out its effect on the immunity. The L. rohita yearlings were found to resist the endotoxin dose up to 20 EU/fish and at the lower doses, i.e., at 1 and 2 EU/fish; it acted as an immune potentiator. Different serum and immune parameters like protein, globulin, lysozyme, respiratory burst activity, myeloperoxidase activity, natural agglutination titre were found to be significantly high (p<0.01) at a dose of 1 EU/fish. While at 10 and 20 EU/fish, most of these parameters were lower thereby indicating the immuno-suppressive nature of the endotoxin at these higher doses.

Non-specific immune parameters of brood Indian major carp Labeo rohita and their seasonal variations

Different non-specific immune parameters and their seasonal changes in brood Indian major carp Labeo rohita reared in two major freshwater aquaculture regions of India viz. West Bengal and Orissa were investigated. It was undertaken for 2 consecutive years and included three main seasons of a year such as summer (March-May), rainy (July-September) and winter (November-January). Total serum protein, albumin and globulin levels were not significantly different throughout the year (p>0.01). Serum lysozyme and myeloperoxidase activities were lower (7.26+/-0.87mg/ml and, 0.54+/-0.11 OD, respectively) in winter as compared to any other season of the year. The bacterial agglutination titer was higher (p<0.01) in the rainy season (8.70+/-1.70) compared to summer and winter seasons (3.40+/-0.60 and 4.00+/-0.89, respectively). Haemagglutination and haemolytic activities did not vary (p>0.01) throughout the year. In blood smears, lymphocyte percentage was higher (75-80%) as compared to those of neutrophil (10-15%) and monocytes (5-10%) but eosinophilic granulocytes were present only in few cases. The differential leucocyte count did not vary significantly (p>0.05) in any season. This study indicated that certain non-specific immune parameters of this species can be modulated at certain times of the year.

The immunomodulatory effects of tuftsin on the non-specific immune system of Indian Major carp, Labeo rohita

The purpose of this study was to determine if injections of different dosages of tuftsin would enhance the immune response and disease resistance against the infections due to the opportunistic pathogens Aeromonas hydrophila and Edwardsiella tarda in Labeo rohita fingerlings. Hence, four different dosages of tuftsin in PBS suspension at the rate of 0, 5, 10, 15 mg kg(-1) body weight of fish were injected intraperitoneally to the fingerlings of L. rohita at 2-week intervals for four times. After every 2-week interval, different serum biochemical, haematological and immunological parameters of fish were evaluated. Biochemical and haematological parameters including serum total protein content, albumin content, globulin content, albulin:globulin ratio, glucose content, leucocyte counts etc.; cellular immune parameters including superoxide anion production, phagocytic activities, lymphokine production index etc.; humoral immune parameters including lysozyme activity, complement activity, serum bactericidal activity etc., in the fish were evaluated after every 2-week interval. After 56 days, fish were divided into two subgroups under each major treatment group for challenge with two pathogens A. hydrophila and E. tarda. The mortality (%) and agglutinating antibody titre was recorded on 28th day post challenge. Most of the immune parameters including leucocyte count, phagocytic ratio, phagocytic index, lysozyme activity, complement activity, and serum bactericidal activity were significantly (p<or=0.05) maximum on 42 days after three i.p. injections of 10 mg kg(-1) body weight of tuftsin. Challenge study indicated least mortality in the group of fish injected with 10 mg kg(-1) body weight of tuftsin for four times. Multiple injections of tuftsin might have maintained the activation of phagocytic cells for a long period, which in turn led to long-term protection in the fish. Thus, multiple injections of 10 mg kg(-1) body weight of tuftsin for three times can be advocated for enhancing the immune response of fish species under aquaculture.

Passive transfer of maternal antibodies and their existence in eggs, larvae and fry of Indian major carp, Labeo rohita (Ham.)

Lack of immune competence in the early stages of life leads to severe mortality in larval stages of different fish species including Indian major carp (IMC). Investigation through indirect enzyme linked immunosorbent assay (ELISA) and agglutination test revealed a significant increase in specific serum antibody response in the brood fish of Indian major carp, Labeo rohita (Ham.) following immunisation with a virulent Aeromonas hydrophila bacterin 1 month prior to breeding, which was transferred to larvae through the egg. No significant differences (P > 0.05) in mean antibody levels in larvae at the 1st and 2nd weeks post-hatch was recorded while a slight rise in antibody level was observed in 3-week-old fry, perhaps due to exposure to A. hydrophila present in the aquatic environment. Immunised brood fish serum, egg and larval extracts in non-reducing sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) and subsequent western blot analysis revealed an antibody molecule of approximate molecular weight 210 kDa. On challenge with virulent A. hydrophila, a significant reduction in mortality was recorded in immunised larvae and fry (58.0, 43.75 and 37.14% in the 1st, 2nd and 3rd week, respectively) relative to control fish (87.0, 79.0 and 76.4% in 1st, 2nd and 3rd week, respectively). The present study indicated the role of maternally derived antibody in protection of hatchlings of Indian major carp against specific pathogens.

High antigenic cross-reaction among the bacterial species responsible for diseases of cultured freshwater fishes and strategies to overcome it for specific serodiagnosis

Antigenic sharing among the most commonly bacterial pathogens such as Aeromonas hydrophila, Edwardsiella tarda and Pseudomonas fluorescens of Indian major carps has been studied using immunological reactions such as cross-agglutination, disc diffusion and indirect enzyme linked immunosorbent assay (ELISA). The data were analysed using statistical analysis (SAS), version 6.12. The results showed high antigenic similarities among the bacterial whole cells, whole cell lysates, somatic 'O' antigens, lipopolysaccharides (LPS) and extracellular products (ECP). However, few or no similarities were observed in ECP components of <20kD. The present study indicates a need to develop differential diagnostic methods based on serology by choosing the highly specific less cross-reactive ECP antigen.

Bath immunisation of spawn, fry and fingerlings of Indian major carps using a particulate bacterial antigen

Larval mortality in Indian major carps is one of the major problems encountered in the pond culture system. The present investigation was carried out to investigate the proper age, duration of exposure, and optimum bacterin concentration for vaccinating rohu (Labeo rohita) and catla (Catla catla) at their early stages with a formalin killed Edwardsiella tarda bacterin suspension. The development of immunological competence was recorded with spawn of rohu and catla of 3 weeks of age exposed to a bacterin at a concentration 10(9) cfu ml(-1) for 15 min, where it persisted up to 4 weeks post vaccination. They showed significant resistance against challenge with virulent E. tarda bacteria. Significant antibody titre could be recorded in advanced fries and fingerlings exposed to 10(9) cfu/ml(-1) bacterin concentration for 45 and 60 min, respectively.