Genome-wide DNA methylation profiles regulate distinct heat stress response in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) cattle.

Epigenetic variations result from long-term adaptation to environmental factors. The Bos indicus (zebu) adapted to tropical conditions, whereas Bos taurus adapted to temperate conditions; hence native zebu cattle and its crossbred (B indicus × B taurus) show differences in responses to heat stress. The present study evaluated genome-wide DNA methylation profiles of these two breeds of cattle that may explain distinct heat stress responses. Physiological responses to heat stress and estimated values of Iberia heat tolerance coefficient and Benezra's coefficient of adaptability revealed better relative thermotolerance of Hariana compared to the Vrindavani cattle. Genome-wide DNA methylation patterns were different for Hariana and Vrindavani cattle. The comparison between breeds indicated the presence of 4599 significant differentially methylated CpGs with 756 hypermethylated and 3845 hypomethylated in Hariana compared to the Vrindavani cattle. Further, we found 79 genes that showed both differential methylation and differential expression that are involved in cellular stress response functions. Differential methylations in the microRNA coding sequences also revealed their functions in heat stress responses. Taken together, epigenetic differences represent the potential regulation of long-term adaptation of Hariana (B indicus) cattle to the tropical environment and relative thermotolerance.

Genome-wide expression analysis reveals different heat shock responses in indigenous (Bos indicus) and crossbred (Bos indicus X Bos taurus) cattle.

Environmental heat stress in dairy cattle leads to poor health, reduced milk production and decreased reproductive efficiency. Multiple genes interact and coordinate the response to overcome the impact of heat stress. The present study identified heat shock regulated genes in the peripheral blood mononuclear cells (PBMC). Genome-wide expression patterns for cellular stress response were compared between two genetically distinct groups of cattle viz., Hariana (B. indicus) and Vrindavani (B. indicus X B. taurus). In addition to major heat shock response genes, oxidative stress and immune response genes were also found to be affected by heat stress. Heat shock proteins such as HSPH1, HSPB8, FKB4, DNAJ4 and SERPINH1 were up-regulated at higher fold change in Vrindavani compared to Hariana cattle. The oxidative stress response genes (HMOX1, BNIP3, RHOB and VEGFA) and immune response genes (FSOB, GADD45B and JUN) were up-regulated in Vrindavani whereas the same were down-regulated in Hariana cattle. The enrichment analysis of dysregulated genes revealed the biological functions and signaling pathways that were affected by heat stress. Overall, these results show distinct cellular responses to heat stress in two different genetic groups of cattle. This also highlight the long-term adaptation of B. indicus (Hariana) to tropical climate as compared to the crossbred (Vrindavani) with mixed genetic makeup (B. indicus X B. taurus).

The role of continental evapotranspiration on water vapour isotopic variability in the troposphere.

Rainforests play an important role in hydrological and carbon cycles, both at regional and global scales. They pump large quantities of moisture from the soil to the atmosphere and are major rainfall hotspots of the world. Satellite-observed stable water isotope ratios have played an essential role in determining sources of moisture in the atmosphere. Satellites provide information about the processes involving vapour transport in different zones of the world, identifying sources of rainfall and distinguishing moisture transport in monsoonal systems. This paper focuses on major rainforests of the world (Southern Amazon, Congo and Northeast India) to understand the role of continental evapotranspiration in influencing tropospheric water vapour. We have used satellite measurements of 1H2H16O/1H216O from Atmospheric InfraRed Sounder (AIRS), evapotranspiration (ET), solar-induced fluorescence (SIF), precipitation (P), atmospheric reanalysis-derived moisture flux convergence (MFC) and wind to discern the role of ET in influencing water vapour isotopes. A global map of the correlation between δ2Hv and ET-P flux indicates that densely vegetated regions in the tropics show the highest positive correlation (r > 0.5). Using mixing models and observations of specific humidity and isotopic ratio over these forested regions, we discern the source of moisture in pre-wet and wet seasons.

N6-Methyladenosine RNA Modification in Host Cells Regulates Peste des Petits Ruminants Virus Replication.

N6-methyladenosine (m6A) modification is a major RNA epigenetic regulatory mechanism. The dynamics of m6A levels in viral genomic RNA and their mRNAs have been shown to have either pro- or antiviral functions, and therefore, m6A modifications influence virus-host interactions. Currently, no reports are available on the effect of m6A modifications in the genome of Peste des petits ruminants virus (PPRV). In the present study, we took PPRV as a model for nonsegmented negative-sense single-stranded RNA viruses and elucidate the role of m6A modification on viral replication. We detected m6A-modified sites in the mRNA of the virus and host cells, as well as the PPRV RNA genome. Further, it was found that the level of m6A modification in host cells alters the viral gene expression. Knockdown of the METTL3 and FTO genes (encoding the m6A RNA modification writer and eraser proteins, respectively) results in alterations of the levels of m6A RNA modifications in the host cells. Experiments using these genetically modified clones of host cells infected with PPRV revealed that both higher and lower m6A RNA modification in the host cells negatively affect PPRV replication. We found that m6A-modified viral transcripts had better stability and translation efficiency compared to the unmodified mRNA. Altogether, from these data, we conclude that the m6A modification of RNA regulates PPRV replication. These findings contribute toward a way forward for developing novel antiviral strategies against PPRV by modulating the dynamics of host m6A RNA modification. IMPORTANCE Peste des petits ruminants virus (PPRV) causes a severe disease in sheep and goats. PPRV infection is a major problem, causing significant economic losses to small ruminant farmers in regions of endemicity. N6-methyladenosine (m6A) is an important RNA modification involved in various functions, including virus-host interactions. In the present study, we used stable clones of Vero cells, having knocked down the genes encoding proteins involved in dynamic changes of the levels of m6A modification. We also used small-molecule compounds that interfere with m6A methylation. This resulted in a platform of host cells with various degrees of m6A RNA modification. The host cells with these different microenvironments were useful for studying the effect of m6A RNA modification on the expression of viral genes and viral replication. The results pinpoint the level of m6A modifications that facilitate the maximum replication of PPRV. These findings will be useful in increasing the virus titers in cultured cells needed for the economical development of the vaccine. Furthermore, the findings have guiding significance for the development of novel antiviral strategies for limiting PPRV replication in infected animals.

Design, Synthesis, and Biological Evaluation of Novel Dihydropyrimidinone Derivatives as Potential Anticancer Agents and Tubulin Polymerization Inhibitors.

The severity and prevalence of cancer in modern time are a huge global health burden. Continuous efforts are being made toward the development of newer therapeutic candidates to treat and manage this ailment. The dihydropyrimidinone scaffold is one of the key nuclei that have been highly explored and investigated against cancer. It has the potential to combat the consequences of cancer by interacting with several biological targets. Tubulin polymerization inhibition is one such strategy to prevent the progression of cancer. In the presented work, we have synthesized a series of sixteen dihydropyrimidinone derivatives by following a rational drug design. The synthesized compounds have been characterized by 1H NMR and 13C NMR and were further evaluated for cytotoxic activity against breast cancer cell lines (MCF-7 and MDA-MB-231), lung cancer cell lines (A549), and colon cancer cell lines (HCT-116). Compounds 5D and 5P were found most potent and revealed a better cytotoxic activity compared with the standard drug colchicine. Furthermore, the tubulin polymerization inhibition assay revealed that compound 5D showed better inhibition than colchicines, whereas compound 5P revealed an almost equal inhibition to that of colchicine. Furthermore, to investigate the possible mode of action and binding patterns, compounds 5P and 5D were subjected to molecular docking against tubulin (Protein Data Bank ID: ISA0). The results showed that compounds revealed significant interactions and were well occupied inside the cavity of tubulin. The compounds 5D and 5P may serve as new leads in drug development against cancer.

Deriving wetland-cover types (WCTs) from integration of multispectral indices based on Earth observation data.

The wetland cover is defined as the spatially homogenous region of a wetland attributed to the underlying biophysical conditions such as vegetation, turbidity, hydric soil, and the amount of water. Here, we present a novel method to derive the wetland-cover types (WCTs) combining three commonly used multispectral indices, NDVI, MNDWI, and NDTI, in three large Ramsar wetlands located in different geomorphic and climatic settings across India. These wetlands include the Kaabar Tal, a floodplain wetland in east Ganga Plains, Chilika Lagoon, a coastal wetland in eastern India, and Nal Sarovar in semi-arid western India. The novelty of our approach is that the derived WCTs are stable in space and time, and therefore, a given WCT across different wetlands or within different zones of a large wetland will imply similar underlying biophysical attributes. The WCTs can therefore provide a novel tool for monitoring and change detection of wetland cover types. We have automated the proposed WCT algorithm using the Google Earth Engine (GEE) environment and by developing ArcGIS tools. The method can be implemented on any wetland and using any multispectral imagery dataset with visible and NIR bands. The proposed methodology is simple yet robust and easy to implement and, therefore, holds significant importance in wetland monitoring and management.

Integrated analysis of long-noncoding RNA and circular RNA expression in Peste-des-Petits-Ruminants Virus (PPRV) infected marmoset B lymphocyte (B95a) cells.

Peste-des-Petits-Ruminants (PPR) or goat plague is an important viral disease of sheep and goats caused by the small ruminant morbilli virus or PPR virus (PPRV). Long non coding RNAs (lncRNA) and circular RNAs (circRNA) play a pivotal role in several biological processes including regulation of virus-host interactions. The present study explored the expression of lncRNA, circRNA and their functions in PPRV infected B-lymphocyte (B95a) cells. The results revealed a total of 4531 lncRNA and 2348 circRNA expression in both mock and PPRV infected samples. Analysis of differentially expressed (DE) RNA identified 123 DE-lncRNA and 39 DE-circRNA as significantly dysregulated. Functional analysis of cis-target genes of DE-lncRNA indicated activation of TCF dependent WNT signaling and PKN1 stimulated transcription process. Interactions (sponging) of microRNA (miRNA) revealed 344 DE-lncRNA-miRNA and 93 DE-circRNA-miRNA pairs. The competing endogenous RNA (ceRNA) network of lncRNA/circRNA-miRNA-mRNA in PPRV infected B95a cells was represented by 69 ceRNA pairs. We validated the DE-circRNA by targeted amplification and sequencing of back spliced junctions (BSJs). The present study revealed a profile of lncRNA, circRNA and their potential ceRNA network in PPRV infection. The results provide insight for better understanding of PPRV-host interactions.

Effect of acute heat shock on stress gene expression and DNA methylation in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) dairy cattle.

Environmental temperature is one of the major factors to affect health and productivity of dairy cattle. Gene expression networks within the cells and tissues coordinate stress response, metabolism, and milk production in dairy cattle. Epigenetic DNA methylations were found to mediate the effect of environment by regulating gene expression patterns. In the present study, we compared three Indian native zebu cattle, Bos indicus (Sahiwal, Tharparkar, and Hariana) and one crossbred Bos indicus × Bos taurus (Vrindavani) for stress gene expression and differences in the DNA methylation patterns. The results indicated acute heat shock to cultured PBMC affected their proliferation, stress gene expression, and DNA methylation. Interestingly, expressions of HSP70, HSP90, and STIP1 were found more pronounced in zebu cattle than the crossbred cattle. However, no significant changes were observed in global DNA methylation due to acute heat shock, even though variations were observed in the expression patterns of DNA methyltransferases (DNMT1, DNMT3a) and demethylases (TET1, TET2, and TET3) genes. The treatment 5-AzaC (5-azacitidine) that inhibit DNA methylation in proliferating PBMC caused significant increase in heat shock-induced HSP70 and STIP1 expression indicating that hypomethylation facilitated stress gene expression. Further targeted analysis DNA methylation in the promoter regions revealed no significant differences for HSP70, HSP90, and STIP1. However, there was a significant hypomethylation for BDNF in both zebu and crossbred cattle. Similarly, NR3C1 promoter region showed hypomethylation alone in crossbred cattle. Overall, the results indicated that tropically adapted zebu cattle had comparatively higher expression of stress genes than the crossbred cattle. Furthermore, DNA methylation may play a role in regulating expression of certain genes involved in stress response pathways.

Predictor of in-stent restenosis in patients with drug-eluting stent (PRIDE)- a retrospective cohort study.

BACKGROUND: It is a fact that coronary artery disease (CAD) is more prevalent in India as compared to western countries. The major risk factors associated with the early CAD are a high prevalence of diabetes mellitus, atherogenic lipid profile, smoking habits, sedentary lifestyle, low socioeconomic condition and high prevalence of obesity. Is this true for restenosis after drug-eluting stent (DES) implantation and factors associated with it? The main objective of the study was to determine the rate of in-stent restenosis (ISR) in patients with DES and risk factors associated with it from our region. METHODS: It was a single-center, retrospective cohort study in which 550 patients who underwent DES implantation were included. Patient's demographic data, coronary angiography findings, procedural characteristics and development of ISR were noted. RESULTS: Out of 550 patients, 31 developed ISR with a rate of restenosis of 5.63% and target lesion revascularization (TLR) of 5.63%. On multiple Cox-regression analysis, only diabetes mellitus (DM) (p=0.008, adjusted hazard ratio (HR): 2.757, 95% confidence interval (CI): 1.296-5.863), deployment of stent in the left anterior descending (LAD) artery (p=0.031, adjusted HR: 3.342, 95% CI: 1.115-10.017) and periprocedural complication during percutaneous coronary intervention (p=0.040, adjusted HR: 2.824, 95% CI: 1.049-7.603) were found to be significantly associated with increased risk of ISR. Kaplan-Meier survival analysis of event-free survival for restenosis showed patients with DM had significantly lower event-free survival compared to patients without DM (p=0.005 by log-rank test). CONCLUSIONS: In our study, the rate of restenosis after DES implantation was 5.63%. The presence of DM, the stent in the LAD territory and the periprocedural complication is strongly associated with the development of ISR.

Functional characterization of partial recombinant goat conglutinin: Its role as innate immunity marker and use as antigen in sandwich ELISA.

Conglutinin, a liver synthesized versatile innate immune marker consisting C-type lectin domain belongs to collectin superfamily of proteins. The protein, first detected in bovine serum as soluble pattern recognition receptor (PRR) has wide range of antimicrobial activities. In the present study, open reading frame (ORF) encoding neck and carbohydrate recognition domain (NCRD) of goat conglutinin gene ligated to the vector pRSET-A was expressed in E. coli BL-21(pLys) cells. The 27 kDa recombinant protein (rGCGN) purified by single step Ni+2 -NTA affinity chromatography was found to cross-react with recombinant anti-buffalo conglutinin antibody raised in poultry. Further, it displayed calcium-dependant sugar binding activity towards yeast mannan and calcium-independent binding activity towards LPS. The mannan binding activity of rGCGN was inhibited in the presence of N-acetyl-glucosamine because of higher affinity towards this sugar. The recombinant protein was found to stimulate production of superoxide ions and hydrogen peroxide in goat neutrophils, which are instrumental in stimulating phagocytic activity of cells. When used as antigen in Sandwich ELISA, straight line (Y = 0.299x + 0.067, R2 = 0.997) was observed within the concentration range of 200-1000 ng/100 µl of rGCGN. Using this equation, the native conglutinin concentration in goat sera was estimated to be 0.5-7.5 µg/ml. The results indicated that prokaryotically expressed functionally active rGCGN can be used as antigen to assess native serum conglutinin levels in Sandwich ELISA and as immunomodulator in therapeutic applications to sequester unwanted immune complexes from the circulation.

Lentiviral-mediated delivery of classical swine fever virus Erns gene into porcine kidney-15 cells for production of recombinant ELISA diagnostic antigen.

Classical swine fever virus (CSFV), a member of the Pestivirus genus within the Flaviviridae family causes contagious fatal disease in swine. Antibodies against E2, Erns and NS3 proteins of virus can be detected in infected animals. Development of an ELISA coating antigen to improve the sensitivity of detecting Erns-specific antibodies in pig sera is always desirable for diagnosis as well as for differentiation of infected from vaccinated animals. In present study, a lentivirus-based gene delivery system was used to develop a stable PK-15 cell line expressing Erns (PK-Erns) for production of diagnostic antigen. The Lenti-Erns virus was purified from the supernatant of co-transfected 293LTV cells and used to transduce PK-15 cells. The homogenous PK-Erns cell line was produced by single cell cloning by monitoring eGFP expression. The Erns gene in the genomic DNA and RNA transcripts in total RNA isolated from PK-Erns cells were detected by PCR and RT-PCR, respectively. Expression of 45 kDa Erns glycoprotein was detected in western blot using CSFV-specific hyperimmune sera. The use of PK-Erns cell lysate as antigen in serial dilution and single dilution ELISAs with known positive and negative pig sera was investigated. The PK-Erns ELISA revealed sensitivity equivalent to commercial HerdChek ELISA kit. The sensitivity, specificity and accuracy of the PK-Erns ELISA was 95%, 100% and 96.66%, respectively compared to ELISA using purified CSFV as coating antigen. When field pig sera (n = 69) were tested in PK-Erns ELISA, a significant correlation between the titers from serial dilution and single dilution ELISA was observed. This indicated that PK-Erns cell line can serve as continuous source of ELISA diagnostic antigen for detection of CSFV-specific antibodies in pig sera.

Multidrug resistance protein 4 (MRP4) is expressed as transcript variants in both Gallus domesticus and Gyps himalyanesis.

The non-steroidal anti-inflammatory drug (NSAID) diclofenac, known to cause hyperuricemia and concomitant visceral gout in Gyps vultures is suggested to be a result of interference with renal uric acid excretion. Three species of Gyps vultures are on the verge of extinction due to nephrotoxic veterinary diclofenac having entered the food chain, notwithstanding the fact that the toxicity of different avian species to the NSAIDs like diclofenac varies. The multidrug resistance protein 4 (MRP4), an organic anion transporter in birds has unique role in unidirectional efflux of urate into proximal renal tubular lumen for excretion and maintenance of homeostasis. We characterized MRP4 channel at molecular level to predict its structural based ligand binding activity in Gallus domesticus (Indian domestic chicken) and Gyps himalayensis (Himalayan griffon vulture). MRP4 gene was amplified using reverse transcribed cDNA from renal tissue sample in overlapping fragments. The obtained amplicons were cloned, sequenced, assembled and analyzed. Multiple alignment and blast analysis revealed point variations and presence of additional stretch of 57 bp towards the 3' end which was confirmed in Real time PCR. Predicted MRP4 polypeptides revealed presence of characteristic 12 transmembrane helices (TMH) with two nucleotide binding domains (NBD). Additional 19 amino acids in transcript variant was found to be localized in NBD2 that might influence the transporter function. The homology modeling and pocket identification throws ample light on varying transport efficacy and paves the way for depicting its role of these amino acids in effect of diclofenac on urate transport in further studies.

Rapid Veterinary Diagnosis of Bovine Reproductive Infectious Diseases from Semen Using Paper-Origami DNA Microfluidics.

The health and well-being of cattle is an important issue in maintaining and increasing global agricultural output. In dairy production within low and middle income countries (LMICs), there is a significant biosensing challenge in detecting sexually transmitted infection (STI) pathogens during animal husbandry, due in part to difficulties associated with the limited infrastructure for veterinary medicine. Here we demonstrate low-cost, multiplexed, and sample-to-answer paper-origami tests for the detection of three bovine infectious reproductive diseases in semen samples, collected at a test site in rural India. Pathogen DNA from one viral pathogen, bovine herpes virus-1 (BoHV-1), and two bacteria (Brucella and Leptospira) was extracted, amplified (using loop-mediated isothermal amplification, LAMP), and detected fluorescently, enabling <1 pg (∼ from 115 to 274 copies per reaction) of target genomic DNA to be measured. Data was collected as a fluorescence signal either visually, using a low-cost hand-held torch, or digitally with a mobile-phone camera. Limits of detection and sensitivities of the paper-origami device for the three pathogens were also evaluated using pathogen-inoculated semen samples and were as few as 50 Leptospira organisms, 50 CFU Brucella, and 1 TCID50 BoHV-1. Semen samples from elite bulls at a germplasm center were also tested in double-blind tests, as a demonstrator for a low-cost, user-friendly point-of-care sensing platform, for in-the-field resource-limited regions. The sensors showed excellent levels of sensitivity and specificity, and for the first time a demonstrated ability of the application of paper microfluidics devices for the diagnosis multiple infectious diseases from semen samples.

Visual detection of Brucella in bovine biological samples using DNA-activated gold nanoparticles.

Brucellosis is a bacterial disease, which, although affecting cattle primarily, has been associated with human infections, making its detection an important challenge. The existing gold standard diagnosis relies on the culture of bacteria which is a lengthy and costly process, taking up to 45 days. New technologies based on molecular diagnosis have been proposed, either through dip-stick, immunological assays, which have limited specificity, or using nucleic acid tests, which enable to identify the pathogen, but are impractical for use in the field, where most of the reservoir cases are located. Here we demonstrate a new test based on hybridization assays with metal nanoparticles, which, upon detection of a specific pathogen-derived DNA sequence, yield a visual colour change. We characterise the components used in the assay with a range of analytical techniques and show sensitivities down to 1000 cfu/ml for the detection of Brucella. Finally, we demonstrate that the assay works in a range of bovine samples including semen, milk and urine, opening up the potential for its use in the field, in low-resource settings.

EvaGreen-based Multiplex Real-time PCR Assay for Rapid Differentiation of Wild-Type and Glycoprotein E-Deleted Bovine Herpesvirus-1 Strains.

Bovine herpesvirus-1 (BoHV-1) is an important viral pathogen causing significant economic losses to the cattle industry. Glycoprotein E-deleted marker vaccines form the basis for BoHV-1 control programs widely, wherein detection and differentiation of wild-type and gE-deleted vaccine strains is of crucial importance for proper disease management. In the present study, we report an EvaGreen-based multiplex real-time polymerase chain reaction (EGRT-PCR) assay for rapid differentiation of wild-type and glycoprotein E-deleted strains of BoHV-1. The EGRT-PCR assay could simultaneously detect two viral genes (glycoprotein B and E) and an internal positive control gene (bovine growth hormone- bGH), in a single-tube reaction. The analytical sensitivity of the EGRT-PCR assay was as little as 10 copies of the BoHV-1 DNA per reaction. The modified real-time PCR assay could successfully differentiate wild-type and gE-deleted BoHV-1 strains based on gene specific melting temperatures (Tm) peaks. Our results have shown that the EGRT-PCR developed in this study might prove to be a promising tool in disease management by enabling rapid differentiation of wild-type and gE-deleted strains of BoHV-1.

Visual Detection of Brucella spp. in Spiked Bovine Semen Using Loop-Mediated Isothermal Amplification (LAMP) Assay.

Several pathogens including Brucella spp. are shed in semen of infected bulls and can be transmitted to cows through contaminated semen during artificial insemination. The present study reports omp2a and bcsp31 gene based loop-mediated isothermal amplification (LAMP) assays for detection of Brucella genomic DNA in semen from infected bulls. The positive results could be interpreted visually by change in colour of reaction mixture containing hydroxyl naphthol blue (HNB) dye from violet to sky blue. LAMP assays based on omp2a and bcsp31 could detect as little as 10 and 100 fg of B. abortus S19 genomic DNA, respectively. Sensitivity of omp2a and bcsp31 LAMP assays for direct detection of organisms in bovine semen was 2.28 × 10(1) CFU and 2.28 × 10(2) CFU of B. abortus S19 in spiked bovine semen, respectively. The omp2a LAMP assay was found equally sensitive to TaqMan probe based real-time PCR and 100 times more sensitive than conventional PCR in identifying Brucella in spiked semen. The diagnostic applicability of the omp2a LAMP assay was evaluated with seventy-nine bovine semen samples and results were re-evaluated through TaqMan probe based real-time PCR and conventional PCR. Taken together, the omp2a LAMP assay is easy to perform, rapid and sensitive in diagnosis of Brucella spp. in bovine semen.

Characterization of codon usage pattern and influencing factors in Japanese encephalitis virus.

Recently, several outbreaks of Japanese encephalitis (JE), caused by Japanese encephalitis virus (JEV), have been reported and it has become cause of concern across the world. In this study, detailed analysis of JEV codon usage pattern was performed. The relative synonymous codon usage (RSCU) values along with mean effective number of codons (ENC) value of 55.30 indicated the presence of low codon usages bias in JEV. The effect of mutational pressure on codon usage bias was confirmed by significant correlations of A3s, U3s, G3s, C3s, GC3s, ENC values, with overall nucleotide contents (A%, U%, G%, C%, and GC%). The correlation analysis of A3s, U3s, G3s, C3s, GC3s, with axis values of correspondence analysis (CoA) further confirmed the role of mutational pressure. However, the correlation analysis of Gravy values and Aroma values with A3s, U3s, G3s, C3s, and GC3s, indicated the presence of natural selection on codon usage bias in addition to mutational pressure. The natural selection was further confirmed by codon adaptation index (CAI) analysis. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern to some extent.