Testes-specific hemoglobins in Drosophila evolved by a combination of sub- and neofunctionalization after gene duplication

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Androglobin gene expression patterns and FOXJ1-dependent regulation indicate its functional association with ciliogenesis

Androglobin (ADGB) represents the latest addition to the globin superfamily in metazoans. The chimeric protein comprises a calpain domain and a unique circularly permutated globin domain. ADGB expression levels are most abundant in mammalian testis, but its cell-type-specific expression, regulation, and function have remained unexplored. Analyzing bulk and single-cell mRNA-Seq data from mammalian tissues, we found that—in addition to the testes—ADGB is prominently expressed in the female reproductive tract, lungs, and brain, specifically being associated with cell types forming motile cilia. Correlation analysis suggested coregulation of ADGB with FOXJ1, a crucial transcription factor of ciliogenesis. Investigating the transcriptional regulation of the ADGB gene, we characterized its promoter using epigenomic datasets, exogenous promoter-dependent luciferase assays, and CRISPR/dCas9-VPR-mediated activation approaches. Reporter gene assays revealed that FOXJ1 indeed substantially enhanced luciferase activity driven by the ADGB promoter. ChIP assays confirmed binding of FOXJ1 to the endogenous ADGB promoter region. We dissected the minimal sequence required for FOXJ1-dependent regulation and fine mapped the FOXJ1 binding site to two evolutionarily conserved regions within the ADGB promoter. FOXJ1 overexpression significantly increased endogenous ADGB mRNA levels in HEK293 and MCF-7 cells. Similar results were observed upon RFX2 overexpression, another key transcription factor in ciliogenesis. The complex transcriptional regulation of the ADGB locus was illustrated by identifying a distal enhancer, responsible for synergistic regulation by RFX2 and FOXJ1. Finally, cell culture studies indicated an ADGB-dependent increase in the number of ciliated cells upon overexpression of the full-length protein, confirming a ciliogenesis-associated role of ADGB in mammals.

The Globin Gene Family in Arthropods: Evolution and Functional Diversity

Globins are small heme-proteins that reversibly bind oxygen. Their most prominent roles in vertebrates are the transport and storage of O₂ for oxidative energy metabolism, but recent research has suggested alternative, non-respiratory globin functions. In the species-rich and ecologically highly diverse taxon of arthropods, the copper-containing hemocyanin is considered the main respiratory protein. However, recent studies have suggested the presence of globin genes and their proteins in arthropod taxa, including model species like Drosophila. To systematically assess the taxonomic distribution, evolution and diversity of globins in arthropods, we systematically searched transcriptome and genome sequence data and found a conserved, widespread occurrence of three globin classes in arthropods: hemoglobin-like (HbL), globin X (GbX), and globin X-like (GbXL) protein lineages. These globin types were previously identified in protostome and deuterostome animals including vertebrates, suggesting their early ancestry in Metazoa. The HbL genes show multiple, lineage-specific gene duplications in all major arthropod clades. Some HbL genes (e.g., Glob2 and 3 of Drosophila) display particularly fast substitution rates, possibly indicating the evolution of novel functions, e.g., in spermatogenesis. In contrast, arthropod GbX and GbXL globin genes show high evolutionary stability: GbXL is represented by a single-copy gene in all arthropod groups except Brachycera, and representatives of the GbX clade are present in all examined taxa except holometabolan insects. GbX and GbXL both show a brain-specific expression. Most arthropod GbX and GbXL proteins, but also some HbL variants, include sequence motifs indicative of potential N-terminal acylation (i.e., N-myristoylation, 3C-palmitoylation). All arthropods except for the brachyceran Diptera harbor at least one such potentially acylated globin copy, confirming the hypothesis of an essential, conserved globin function associated with the cell membrane. In contrast to other animals, the fourth ancient globin lineage, represented by neuroglobin, appears to be absent in arthropods, and the putative arthropod orthologs of the fifth metazoan globin lineage, androglobin, lack a recognizable globin domain. Thus, the remarkable evolutionary stability of some globin variants is contrasted by occasional dynamic gene multiplication or even loss of otherwise strongly conserved globin lineages in arthropod phylogeny.

Transcriptomes reveal expression of hemoglobins throughout insects and other Hexapoda

Insects have long been thought to largely not require hemoglobins, with some notable exceptions like the red hemolymph of chironomid larvae. The tubular, branching network of tracheae in hexapods is traditionally considered sufficient for their respiration. Where hemoglobins do occur sporadically in plants and animals, they are believed to be either convergent, or because they are ancient in origin and their expression is lost in many clades. Our comprehensive analysis of 845 Hexapod transcriptomes, totaling over 38 Gbases, revealed the expression of hemoglobins in all 32 orders of hexapods, including the 29 recognized orders of insects. Discovery and identification of 1333 putative hemoglobins were achieved with target-gene BLAST searches of the NCBI TSA database, verifying functional residues, secondary- and tertiary-structure predictions, and localization predictions based on machine learning. While the majority of these hemoglobins are intracellular, extracellular ones were recovered in 38 species. Gene trees were constructed via multiple-sequence alignments and phylogenetic analyses. These indicate duplication events within insects and a monophyletic grouping of hemoglobins outside other globin clades, for which we propose the term insectahemoglobins. These hemoglobins are phylogenetically adjacent and appear structurally convergent with the clade of chordate myoglobins, cytoglobins, and hemoglobins. Their derivation and co-option from early neuroglobins may explain the widespread nature of hemoglobins in various kingdoms and phyla. These results will guide future work involving genome comparisons to transcriptome results, experimental investigations of gene expression, cell and tissue localization, and gas binding properties, all of which are needed to further illuminate the complex respiratory adaptations in insects.

Divergent roles of the Drosophila melanogaster globins

In contrast to long-held assumptions, the gene repertoire of most insects includes hemoglobins. Analyses of the genome of the fruitfly Drosophila melanogaster identified three distinct hemoglobin genes (glob1, glob2, and glob3). While glob1 is predominantly associated with the tracheal system and fat body, glob2 and glob3 are almost exclusively expressed in the testis. The physiological role of globins in Drosophila is uncertain. Here, we studied the functions of the three globins in a cell culture system. Drosophila Schneider 2 (S2) cells were stably transfected with each of the three globins and the empty vector as control. Under hypoxia (1% atmospheric O₂), only glob1 overexpression enhanced the activity of mitochondrial oxidases and the ATP content. However, the positive effect of glob1 expression disappeared after 24h hypoxia, suggesting metabolic adaptations of the S2 cells. glob2 and glob3 had no positive effect on hypoxia-survival. After application of oxidative stress by H₂O₂, glob2 dramatically enhanced the viability of S2 cells. Evaluation of the intracellular localization of the globins using specific antibodies and green fluorescent protein-fusion constructs suggested that glob1 and glob2 most likely reside in the cytoplasm, while glob3 is associated with structures that may represent parts of the intracellular transport machinery. In silico analyses of public RNA-Seq data from different developmental stages provided that glob1 is co-expressed with genes of the aerobic energy metabolism, while glob2 and glob3 expression can be related to spermatogenesis and reproduction. Together, the results indicate divergent functions of the Drosophila globins: glob1 may play a role in the O₂-dependent metabolism while glob2 may protect spermatogenesis from reactive oxygen species.

Analysis of Mitochondrial haemoglobin in Parkinson's disease brain

Mitochondrial dysfunction is an early feature of neurodegeneration. We have shown there are mitochondrial haemoglobin changes with age and neurodegeneration. We hypothesised that altered physiological processes are associated with recruitment and localisation of haemoglobin to these organelles. To confirm a dynamic localisation of haemoglobin we exposed Drosophila melanogaster to cyclical hypoxia with recovery. With a single cycle of hypoxia and recovery we found a relative accumulation of haemoglobin in the mitochondria compared with the cytosol. An additional cycle of hypoxia and recovery led to a significant increase of mitochondrial haemoglobin (p<0.05). We quantified ratios of human mitochondrial haemoglobin in 30 Parkinson's and matched control human post-mortem brains. Relative mitochondrial/cytosolic quantities of haemoglobin were obtained for the cortical region, substantia nigra and cerebellum. In age matched post-mortem brain mitochondrial haemoglobin ratios change, decreasing with disease duration in female cerebellum samples (n=7). The change is less discernible in male cerebellum (n=18). In cerebellar mitochondria, haemoglobin localisation in males with long disease duration shifts from the intermembrane space to the outer membrane of the organelle. These new data illustrate dynamic localisation of mitochondrial haemoglobin within the cell. Mitochondrial haemoglobin should be considered in the context of gender differences characterised in Parkinson's disease. It has been postulated that cerebellar circuitry may be activated to play a protective role in individuals with Parkinson's. The changing localisation of intracellular haemoglobin in response to hypoxia presents a novel pathway to delineate the role of the cerebellum in Parkinson's disease.

Knockdown of Drosophila hemoglobin suggests a role in O2 homeostasis

Almost all insects are equipped with a tracheal system, which appears to be sufficient for O2 supply even in phases of high metabolic activity. Therefore, with the exception of a few species dwelling in hypoxic habitats, specialized respiratory proteins had been considered unnecessary in insects. The recent discovery and apparently universal presence of intracellular hemoglobins in insects has remained functionally unexplained. The fruitfly Drosophila melanogaster harbors three different globin genes (referred to as glob1-3). Glob1 is the most highly expressed globin and essentially occurs in the tracheal system and the fat body. To better understand the functions of insect globins, the levels of glob1 were modulated in Drosophila larvae and adults by RNAi-mediated knockdown and transgenic over-expression. No effects on the development were observed in flies with manipulated glob1 levels. However, the knockdown of glob1 led to a significantly reduced survival rate of adult flies under hypoxia (5% and 1.5% O2). Surprisingly, the glob1 knockdown flies also displayed increased resistance towards the reactive oxygen species-forming agent paraquat, which may be explained by a restricted availability of O2 resulting in decreased formation of harmful O2(-). In summary, our results suggest an important functional role of glob1 in O2 homeostasis, possibly by enhancing O2 supply.

Faster evolving Drosophila paralogs lose expression rate and ubiquity and accumulate more non-synonymous SNPs

Background

Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates.

Results

We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs’ substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart.

Conclusions

Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles.

Reviewers

This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf.

Characterization of the hemoglobin of the backswimmer Anisops deanei (Hemiptera)

While O(2)-binding hemoglobin-like proteins are present in many insects, prominent amounts of hemoglobin have only been found in a few species. Backswimmers of the genera Anisops and Buenoa (Notonectidae) have high concentrations of hemoglobin in the large tracheal cells of the abdomen. Oxygen from the hemoglobin is delivered to a gas bubble and controls the buoyant density, which enables the bugs to maintain their position without swimming and to remain stationary in the mid-water zone where they hunt for prey. We have obtained the cDNA sequences of three Anisops deanei hemoglobin chains by RT-PCR and RACE techniques. The deduced amino acid sequences show an unusual insertion of a single amino acid in the conserved helix E, but this does not affect protein stability or ligand binding kinetics. Recombinant A. deanei hemoglobin has an oxygen affinity of P(50) = 2.4 kPa (18 torr) and reveals the presence of a dimeric fraction or two different conformations. The absorption spectra demonstrate that the Anisops hemoglobin is a typical pentacoordinate globin. Phylogenetic analyses show that the backswimmer hemoglobins evolved within Heteroptera and most likely originated from an intracellular hemoglobin with divergent function.