Evolutionary analysis of globin domains from kinetoplastids

Globin (Gb) domains function in sensing gaseous ligands like oxygen and nitric oxide. In recent years, Gb domain containing heme binding adenylate cyclases (OsAC or GbAC) emerged as significant modulator of Leishmania response to hypoxia and oxidative stress. During progression of life cycle stages, kinetoplastids experience altered condition in insect vectors or other hosts. Moreover, marked diversity in life style has been accounted among kinetoplastids. Distribution and abundance of Gb-domains vary between different groups of kinetoplastids. While in bodonoids, Gbs are not combined with any other functional domains, in trypanosomatids it is either fused with adenylate cyclase (AC) or oxidoreductase (OxR) domains. In salivarian trypanosomatids and Leishmania (Viannia) subtypes, no gene product featuring Gbs can be identified. In this context, evolution of Gb-domains in kinetoplastids was explored. GbOxR derived Gbs clustered with bacterial flavohemoglobins (fHb) including one fHb from Advenella , an endosymbiont of monoxeneous trypanosomatids. Codon adaptation and other evolutionary analysis suggested that OsAC (LmjF.28.0090), the solitary Gb-domain featuring gene product in Leishmania , was acquired via possible horizontal gene transfer. Substantial functional divergence was estimated between orthologues of genes encoding GbAC or GbOxR; an observation also reflected in structural alignment and heme-binding residue predictions. Orthologue–paralogue and synteny analysis indicated genomic reduction in GbOxR and GbAC loci for dixeneous trypanosomatids.