Evolution of hemoglobin and its genes

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells

Background

Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis.

Methods

Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells.

Results

Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding.

Conclusion

Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

C2H2 Zinc Finger Transcription Factors Associated with Hemoglobinopathies

In humans, specific aberrations in β-globin results in sickle cell disease and β-thalassemia, symptoms of which can be ameliorated by increased expression of fetal globin (HbF). Two recent CRISPR-Cas9 screens, centered on ∼1500 annotated sequence-specific DNA binding proteins and performed in a human erythroid cell line that expresses adult hemoglobin, uncovered four groups of candidate regulators of HbF gene expression. They are (1) members of the nucleosome remodeling and deacetylase (NuRD) complex proteins that are already known for HbF control; (2) seven C2H2 zinc finger (ZF) proteins, including some (ZBTB7A and BCL11A) already known for directly silencing the fetal γ-globin genes in adult human erythroid cells; (3) a few other transcription factors of different structural classes that might indirectly influence HbF gene expression; and (4) DNA methyltransferase 1 (DNMT1) that maintains the DNA methylation marks that attract the MBD2-associated NuRD complex to DNA as well as associated histone H3 lysine 9 methylation. Here we briefly discuss the effects of these regulators, particularly C2H2 ZFs, in inducing HbF expression for treating β-hemoglobin disorders, together with recent advances in developing safe and effective small-molecule therapeutics for the regulation of this well-conserved hemoglobin switch.

Local energetic frustration conservation in protein families and superfamilies

Energetic local frustration offers a biophysical perspective to interpret the effects of sequence variability on protein families. Here we present a methodology to analyze local frustration patterns within protein families and superfamilies that allows us to uncover constraints related to stability and function, and identify differential frustration patterns in families with a common ancestry. We analyze these signals in very well studied protein families such as PDZ, SH3, ɑ and β globins and RAS families. Recent advances in protein structure prediction make it possible to analyze a vast majority of the protein space. An automatic and unsupervised proteome-wide analysis on the SARS-CoV-2 virus demonstrates the potential of our approach to enhance our understanding of the natural phenotypic diversity of protein families beyond single protein instances. We apply our method to modify biophysical properties of natural proteins based on their family properties, as well as perform unsupervised analysis of large datasets to shed light on the physicochemical signatures of poorly characterized proteins such as the ones belonging to emergent pathogens.

Comparison of Evolutionary Relationships between Branchiostoma floridae, Ciona intestinalis, and Homo sapiens Globins Provide Evidence of Gene Co-Option and Convergent Evolution

Globins have been studied as model proteins to elucidate the principles of protein evolution. This was achieved by understanding the relationship between amino acid sequence, three-dimensional structure, physicochemical properties, and physiological function. Previous molecular phylogenies of chordate globin genes revealed the monophyletic evolution of urochordate globins and suggested convergent evolution. However, to provide evidence of convergent evolution, it is necessary to determine the physicochemical and functional similarities between vertebrates and urochordate globins. In this study, we determined the expression patterns of Ciona globin genes using real-time RT-PCR. Two genes (Gb-1 and Gb-2) were predominantly expressed in the branchial sac, heart, and hemocytes and were induced under hypoxia. Combined with the sequence analysis, our findings suggest that Gb-1/-2 correspond to vertebrate hemoglobin-α/-β. However, we did not find a robust similarity between Gb-3, Gb-4, and vertebrate globins. These results suggested that, even though Ciona globins obtained their unique functions differently from vertebrate globins, the two of them shared some physicochemical features and physiological functions. Our findings offer a good example for understanding the molecular mechanisms underlying gene co-option and convergence, which could lead to evolutionary innovations.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

Characterization of transcriptional enhancers in the chicken genome using CRISPR-mediated activation

DNA regulatory elements intricately control when, where, and how genes are activated. Therefore, understanding the function of these elements could unveil the complexity of the genetic regulation network. Genome-wide significant variants are predominantly found in non-coding regions of DNA, so comprehending the predicted functional regulatory elements is crucial for understanding the biological context of these genomic markers, which can be incorporated into breeding programs. The emergence of CRISPR technology has provided a powerful tool for studying non-coding regulatory elements in genomes. In this study, we leveraged epigenetic data from the Functional Annotation of Animal Genomes project to identify promoter and putative enhancer regions associated with three genes (HBBA, IRF7, and PPARG) in the chicken genome. To identify the enhancer regions, we designed guide RNAs targeting the promoter and candidate enhancer regions and utilized CRISPR activation (CRISPRa) with dCas9-p300 and dCas9-VPR as transcriptional activators in chicken DF-1 cells. By comparing the expression levels of target genes between the promoter activation and the co-activation of the promoter and putative enhancers, we were able to identify functional enhancers that exhibited augmented upregulation. In conclusion, our findings demonstrate the remarkable efficiency of CRISPRa in precisely manipulating the expression of endogenous genes by targeting regulatory elements in the chicken genome, highlighting its potential for functional validation of non-coding regions.

IFF: Identifying key residues in intrinsically disordered regions of proteins using machine learning

Conserved residues in protein homolog sequence alignments are structurally or functionally important. For intrinsically disordered proteins or proteins with intrinsically disordered regions (IDRs), however, alignment often fails because they lack a steric structure to constrain evolution. Although sequences vary, the physicochemical features of IDRs may be preserved in maintaining function. Therefore, a method to retrieve common IDR features may help identify functionally important residues. We applied unsupervised contrastive learning to train a model with self-attention neuronal networks on human IDR orthologs. Parameters in the model were trained to match sequences in ortholog pairs but not in other IDRs. The trained model successfully identifies previously reported critical residues from experimental studies, especially those with an overall pattern (e.g., multiple aromatic residues or charged blocks) rather than short motifs. This predictive model can be used to identify potentially important residues in other proteins, improving our understanding of their functions. The trained model can be run directly from the Jupyter Notebook in the GitHub repository using Binder (mybinder.org). The only required input is the primary sequence. The training scripts are available on GitHub (https://github.com/allmwh/IFF). The training datasets have been deposited in an Open Science Framework repository (https://osf.io/jk29b).

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Point-of-care paper-based analytical device for potentiometric detection of myoglobin as a cardiovascular disease biomarker

One of the cardiac biomarkers, myoglobin (Mb), is important in the rapid identification of cardio-vascular disorders. Therefore, point-of-care monitoring is essential. Pursuing this goal, a robust, reliable, and affordable paper-based analytical apparatus for potentiometric sensing has been developed and characterized. The molecular imprint technique was used to create a customized biomimetic antibody for myoglobin (Mb) on the surface of carboxylated multiwalled carbon nanotubes (MWCNT-COOH). This was accomplished by attaching Mb to carboxylated MWCNTs' surfaces and then filling the empty spaces through the mild polymerization of acrylamide in N,N-methylenebisacrylamide and ammonium persulphate. The modification of the MWCNTs' surface was verified by SEM and FTIR analysis. A hydrophobic paper substrate coated with fluorinated alkyl silane (CF₃(CF₂)₇CH₂CH₂SiCl₃, CF¹⁰) has been coupled with a printed all-solid-state Ag/AgCl reference electrode. The presented sensors showed a linear range of 5.0 × 10^-8 to 1.0 × 10^-4 M with a potentiometric slope of -57.1 ± 0.3 mV decade^-1 (R² = 0.9998) and a detection limit of 28 nM at pH 4. Compared to creatinine, sucrose, fructose, galactose, sodium glutamate, thiamine, alanine, ammonium, uric acid, albumin, glutamine, guanine, troponine T, and glucose, the sensor showed good selectivity for Mb. It demonstrated a good recovery for the detection of Mb in several fake serum samples (93.0-103.3%), with an average relative standard deviation of 4.5%. The current approach might be viewed as a potentially fruitful analytical tool for obtaining disposable, cost-effective paper-based potentiometric sensing devices. These types of analytical devices can be potentially manufacturable at large scales in clinical analysis.

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Thalassemia and Nanotheragnostics: Advanced Approaches for Diagnosis and Treatment

Thalassemia is a monogenic autosomal recessive disorder caused by mutations, which lead to abnormal or reduced production of hemoglobin. Ineffective erythropoiesis, hemolysis, hepcidin suppression, and iron overload are common manifestations that vary according to genotypes and dictate, which diagnosis and therapeutic modalities, including transfusion therapy, iron chelation therapy, HbF induction, gene therapy, and editing, are performed. These conventional therapeutic methods have proven to be effective, yet have several disadvantages, specifically iron toxicity, associated with them; therefore, there are demands for advanced therapeutic methods. Nanotechnology-based applications, such as the use of nanoparticles and nanomedicines for theragnostic purposes have emerged that are simple, convenient, and cost-effective methods. The therapeutic potential of various nanoparticles has been explored by developing artificial hemoglobin, nano-based iron chelating agents, and nanocarriers for globin gene editing by CRISPR/Cas9. Au, Ag, carbon, graphene, silicon, porous nanoparticles, dendrimers, hydrogels, quantum dots, etc., have been used in electrochemical biosensors development for diagnosis of thalassemia, quantification of hemoglobin in these patients, and analysis of conventional iron chelating agents. This review summarizes the potential of nanotechnology in the development of various theragnostic approaches to determine thalassemia-causing gene mutations using various nano-based biosensors along with the employment of efficacious nano-based therapeutic procedures, in contrast to conventional therapies.

A specific G9a inhibitor unveils BGLT3 lncRNA as a universal mediator of chemically induced fetal globin gene expression

Sickle cell disease (SCD) is a heritable disorder caused by β-globin gene mutations. Induction of fetal γ-globin is an established therapeutic strategy. Recently, epigenetic modulators, including G9a inhibitors, have been proposed as therapeutic agents. However, the molecular mechanisms whereby these small molecules reactivate γ-globin remain unclear. Here we report the development of a highly selective and non-genotoxic G9a inhibitor, RK-701. RK-701 treatment induces fetal globin expression both in human erythroid cells and in mice. Using RK-701, we find that BGLT3 long non-coding RNA plays an essential role in γ-globin induction. RK-701 selectively upregulates BGLT3 by inhibiting the recruitment of two major γ-globin repressors in complex with G9a onto the BGLT3 gene locus through CHD4, a component of the NuRD complex. Remarkably, BGLT3 is indispensable for γ-globin induction by not only RK-701 but also hydroxyurea and other inducers. The universal role of BGLT3 in γ-globin induction suggests its importance in SCD treatment.

CRISPR Activator Approaches to Study Endogenous Androglobin Gene Regulation

Androglobin (ADGB), the most recently identified member of the mammalian globin family, is a chimeric protein with an unusual, embedded globin domain that is circularly permutated and exhibits hallmarks of a hexacoordinated heme-binding scheme. Whereas abundant expression of ADGB was initially found to be mainly restricted to cells in the postmeiotic stages of spermatogenesis, more recent RNA-Seq-based expression analysis data revealed that ADGB is detectable in cells carrying motile cilia or flagella. This very tight regulation of ADGB gene expression urges the need for alternative techniques to study endogenous expression in classical mammalian cell models, which do not express ADGB. We describe here the use of CRISPR activation (CRISPRa) technology to induce endogenous ADGB gene expression in HEK293T, MCF-7, and HeLa cells from its promoter and illustrate how this method can be employed to validate putative regulatory DNA elements of ADGB in promoter and enhancer regions.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as “MNˮ and “LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

Segregation of α- and β-Globin Gene Cluster in Vertebrate Evolution: Chance or Necessity?

The review is devoted to the patterns of evolution of α- and β-globin gene domains. A hypothesis is presented according to which segregation of the ancestral cluster of α/β-globin genes in Amniota occurred due to the performance by α-globins and β-globins of non-canonical functions not related to oxygen transport.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Limitations of mouse models for sickle cell disease conferred by their human globin transgene configurations

We characterized the human β-like globin transgenes in two mouse models of sickle cell disease (SCD) and tested a genome-editing strategy to induce red blood cell fetal hemoglobin (HbF; α₂γ₂). Berkeley SCD mice contain four to 22 randomly arranged, fragmented copies of three human transgenes (HBA1, HBG2-HBG1-HBD-HBB^S and a mini-locus control region) integrated into a single site of mouse chromosome 1. Cas9 disruption of the BCL11A repressor binding motif in the γ-globin gene (HBG1 and HBG2; HBG) promoters of Berkeley mouse hematopoietic stem cells (HSCs) caused extensive death from multiple double-strand DNA breaks. Long-range sequencing of Townes SCD mice verified that the endogenous Hbb genes were replaced by single-copy segments of human HBG1 and HBB^S including proximal but not some distal gene-regulatory elements. Townes mouse HSCs were viable after Cas9 disruption of the HBG1 BCL11A binding motif but failed to induce HbF to therapeutic levels, contrasting with human HSCs. Our findings provide practical information on the genomic structures of two common mouse SCD models, illustrate their limitations for analyzing therapies to induce HbF and confirm the importance of distal DNA elements in human globin regulation.

This article has an associated First Person interview with the first author of the paper.

Editor's choice: This study describes the genomic structures of two common sickle cell disease mouse models, illustrates their limitations for analyzing some genetic therapies and confirms the importance of distal DNA elements in human globin gene regulation.

Physiological roles of hydrogen sulfide in mammalian cells, tissues and organs

H2S belongs to the class of molecules known as gasotransmitters, which also includes nitric oxide (NO) and carbon monoxide (CO). Three enzymes are recognized as endogenous sources of H2S in various cells and tissues: cystathionine g-lyase (CSE), cystathionine β-synthase (CBS) and 3-mercaptopyruvate sulfurtransferase (3-MST). The current article reviews the regulation of these enzymes as well as the pathways of their enzymatic and non-enzymatic degradation and elimination. The multiple interactions of H2S with other labile endogenous molecules (e.g. NO) and reactive oxygen species are also outlined. The various biological targets and signaling pathways are discussed, with special reference to H2S and oxidative posttranscriptional modification of proteins, the effect of H2S on channels and intracellular second messenger pathways, the regulation of gene transcription and translation and the regulation of cellular bioenergetics and metabolism. The pharmacological and molecular tools currently available to study H2S physiology are also reviewed, including their utility and limitations. In subsequent sections, the role of H2S in the regulation of various physiological and cellular functions is reviewed. The physiological role of H2S in various cell types and organ systems are overviewed. Finally, the role of H2S in the regulation of various organ functions is discussed as well as the characteristic bell-shaped biphasic effects of H2S. In addition, key pathophysiological aspects, debated areas, and future research and translational areas are identified A wide array of significant roles of H2S in the physiological regulation of all organ functions emerges from this review.

From Water to Land: The Structural Construction and Molecular Switches in Lungs during Metamorphosis of Microhyla fissipes

Simple Summary

The functionalization of lungs is a necessity for most anurans to breathe on land. Previous studies have focused on the morphological and physiological functions of amphibian lungs, while the microstructural changes and molecular mechanisms that underpin the functional maturation of lungs remain under-researched. We used integrated histology and transcriptomics to study the critical cytological and molecular events associated with lung maturation in Microhyla fissipes. The results illuminated the molecular processes and their coordination in lung development, providing an insight into the transition of amphibians from aquatic to terrestrial life stages.

Abstract

Most anurans must undergo metamorphosis to adapt to terrestrial life. This process enhances the air-breathing ability of the lungs to cope with the change in oxygen medium from water to air. Revealing the structural construction and molecular switches of lung organogenesis is essential to understanding the realization of the air-breathing function. In this study, histology and transcriptomics were conducted in combination to explore these issues in Microhyla fissipes’ lungs during metamorphosis. During the pro-metamorphic phase, histological structural improvement of the alveolar wall is accompanied by robust substrate metabolism and protein turnover. The lungs, at the metamorphic climax phase, are characterized by an increased number of cilia in the alveolar epithelial cells and collagenous fibers in the connective tissues, corresponding to the transcriptional upregulation of cilia and extracellular matrix-related genes. Post-metamorphic lungs strengthen their contracting function, as suggested by the thickened muscle layer and the upregulated expression of genes involved in muscle contraction. The blood–gas barrier is fully developed in adult lungs, the transcriptional features of which are tissue growth and regulation of differentiation and immunity. Importantly, significant transcriptional switches of pulmonary surfactant protein and hemoglobin facilitate air breathing. Our results illuminated four key steps of lung development for amphibians to transition from water to land.

Alterations to cavefish red blood cells provide evidence of adaptation to reduced subterranean oxygen

Animals inhabiting extreme environments allow the powerful opportunity to examine adaptive evolution in response to diverse pressures. One such pressure is reduced oxygen, commonly present at high-altitude and subterranean environments. Cave-dwelling animals must also deal with darkness and starvation, both of which have been rigorously studied as key forces driving the evolution of cave-associated traits. Interestingly, hypoxia as an environmental pressure has received less attention. Here we examined putatively adaptive phenotypes evolving in a freshwater teleost fish, Astyanax mexicanus, which includes both surface- and cave-dwelling forms. This model system also provides the opportunity to identify convergent responses to hypoxia, owing to the presence of numerous natural and independently-colonised cave populations, alongside closely-related surface conspecifics. The focus of this study is hemoglobin, an essential molecule for oxygen transport and delivery. We found that multiple cave populations harbor a higher concentration of hemoglobin in their blood, which is coincident with an increase in cave morph erythrocyte size compared to surface fish. Interestingly, both cave and surface morphs have comparable numbers of erythrocytes per unit of blood, suggesting elevated hemoglobin is not due to overproduction of red blood cells. Alternatively, owing to an increased cell area of erythrocytes in cavefish, we reason that they contain more hemoglobin per erythrocyte. These findings support the notion that cavefish have adapted to hypoxia in caves through modulation of both hemoglobin production and erythrocyte size. This work reveals an additional adaptive feature of Astyanax cavefish, and demonstrates that coordinated changes between cellular architecture and molecular changes are necessary for organisms evolving under intense environmental pressure.

Photoglobin, a distinct family of non-heme binding globins, defines a potential photosensor in prokaryotic signal transduction systems

Globins constitute an ancient superfamily of proteins, exhibiting enormous structural and functional diversity, as demonstrated by many heme-binding families and two non-heme binding families that were discovered in bacterial stressosome component RsbR and in light-harvesting phycobiliproteins (phycocyanin) in cyanobacteria and red algae. By comprehensively exploring the globin repertoire using sensitive computational analyses of sequences, structures, and genomes, we present the identification of the third family of non-heme binding globins—the photoglobin. By conducting profile-based comparisons, clustering analyses, and structural modeling, we demonstrate that photoglobin is related to, but distinct from, the phycocyanin family. Photoglobin preserves a potential ligand-binding pocket, whose residue configuration closely resembles that of phycocyanin, indicating that photoglobin potentially binds to a comparable linear tetrapyrrole. By exploring the contextual information provided by the photoglobin’s domain architectures and gene-neighborhoods, we found that photoglobin is frequently associated with the B12-binding light sensor domain and many domains typical of prokaryotic signal transduction systems. Structural modeling using AlphaFold2 demonstrated that photoglobin and B12-binding domains form a structurally conserved hub among different domain architecture contexts. Based on these strong associations, we predict that the coupled photoglobin and B12-binding domains act as a light-sensing regulatory bundle, with each domain sensing different wavelengths of light resulting in switch-like regulation of downstream signaling effectors. Thus, based on the above lines of evidence, we present a distinct non-heme binding globin family and propose that it may define a new type of light sensor, by means of a linear tetrapyrrole, in complex prokaryotic signal transduction systems.

GAPDH is involved in the heme-maturation of myoglobin and hemoglobin

GAPDH, a heme chaperone, has been previously implicated in the incorporation of heme into iNOS and soluble guanylyl cyclase (sGC). Since sGC is critical for myoglobin (Mb) heme-maturation, we investigated the role of GAPDH in the maturation of this globin, as well as hemoglobins α, β, and γ. Utilizing cell culture systems, we found that overexpression of wild-type GAPDH increased, whereas GAPDH mutants H53A and K227A decreased, the heme content of Mb and Hbα and Hbβ. Overexpression of wild-type GAPDH fully recovered the heme-maturation inhibition observed with the GAPDH mutants. Partial rescue was observed by overexpression of sGCβ1 but not by overexpression of a sGCΔβ1 deletion mutant, which is unable to bind the sGCα1 subunit required to form the active sGCα1β1 complex. Wild type and mutant GAPDH was found to be associated in a complex with each of the globins and Hsp90. GAPDH at endogenous levels was found to be associated with Mb in differentiating C2C12 myoblasts, and with Hbγ or Hbα in differentiating HiDEP-1 erythroid progenitor cells. Knockdown of GAPDH in C2C12 cells suppressed Mb heme-maturation. GAPDH knockdown in K562 erythroleukemia cells suppressed Hbα and Hbγ heme-maturation as well as Hb dimerization. Globin heme incorporation was not only dependent on elevated sGCα1β1 heterodimer formation, but also influenced by iron provision and magnitude of expression of GAPDH, d-aminolevulinic acid, and FLVCR1b. Together, our data support an important role for GAPDH in the maturation of myoglobin and γ, β, and α hemoglobins.

Modulation of Functional Characteristics of Mesenchymal Stromal Cells by Acellular Preparation of Porcine Hemoglobin

Exploring the potential usage of the acellular preparation of porcine hemoglobin (PHb) isolated from slaughterhouse blood as a cell culture media component, we have tested its effects on the functional characteristics of stromal cells of mesodermal origin. Human peripheral blood mesenchymal stromal cells (PB-MSCs) were used in this study as a primary cell model system, along with three mouse cell lines (ATDC5, MC3T3-E1, and 3T3-L1), which represent more uniform model systems. We investigated the effect of PHb at concentrations of 0.1, 1, and 10 μM on these cells’ proliferation, cycle, and clonogenic and migratory potential, and found that PHb’s effect depended on both the cell type and its concentration. At the lowest concentration used (0.1 μM), PHb showed the least evident impact on the cell growth and migration; hence, we analyzed its effect on mesenchymal cell multilineage differentiation capacity at this concentration. Even under conditions that induce a specific type of MSC differentiation (cultivation in particular differentiation media), PHb modulated chondrogenic, osteogenic, and adipogenic differentiation, making it a potential candidate for a supplement of MSC culture. Through a model of porcine hemoglobin, these findings also contribute to improving the knowledge of extracellular hemoglobin’s influence on MSCs >in vivo.

Molecular mechanisms of hematological and biochemical alterations in malaria: A review

Malaria is a dangerous disease that contributes to millions of hospital visits and hundreds of thousands of deaths, especially in children residing in sub-Saharan Africa. Although several interventions such as vector control, case detection, and treatment are already in place, there is no substantive reduction in the disease burden. Several studies in the past have reported the emergence of resistant strains of malaria parasites (MPs) and mosquitoes, and poor adherence and inaccessibility to effective antimalarial drugs as the major factors for this persistent menace of malaria infections. Moreover, victory against MP infections for many years has been hampered by an incomplete understanding of the complex nature of malaria pathogenesis. Very recent studies have identified different complex interactions and hematological alterations induced by malaria parasites. However, no studies have hybridized these alterations for a better understanding of Malaria pathogenesis. Hence, this review thoroughly discusses the molecular mechanisms of all reported hematological and biochemical alterations induced by MPs infections. Specifically, the mechanisms in which MP-infection induces anemia, thrombocytopenia, leukopenia, dyslipidemia, hypoglycemia, oxidative stress, and liver and kidney malfunctions were presented. The study also discussed how MPs evade the host's immune response and suggested strategies to limit evasion of the host's immune response to combat malaria and its complications.

VTR: A Web Tool for Identifying Analogous Contacts on Protein Structures and Their Complexes

Evolutionarily related proteins can present similar structures but very dissimilar sequences. Hence, understanding the role of the inter-residues contacts for the protein structure has been the target of many studies. Contacts comprise non-covalent interactions, which are essential to stabilize macromolecular structures such as proteins. Here we show VTR, a new method for the detection of analogous contacts in protein pairs. The VTR web tool performs structural alignment between proteins and detects interactions that occur in similar regions. To evaluate our tool, we proposed three case studies: we 1) compared vertebrate myoglobin and truncated invertebrate hemoglobin; 2) analyzed interactions between the spike protein RBD of SARS-CoV-2 and the cell receptor ACE2; and 3) compared a glucose-tolerant and a non-tolerant β-glucosidase enzyme used for biofuel production. The case studies demonstrate the potential of VTR for the understanding of functional similarities between distantly sequence-related proteins, as well as the exploration of important drug targets and rational design of enzymes for industrial applications. We envision VTR as a promising tool for understanding differences and similarities between homologous proteins with similar 3D structures but different sequences. VTR is available at http://bioinfo.dcc.ufmg.br/vtr.

The road (not) taken - Placental transfer and interspecies differences

Species differences are among the main reasons for the high failure rate of preclinical studies. A better awareness and understanding of these differences might help to improve the outcome of preclinical research. In reproduction, the placenta is the central organ regulating fetal exposure to a substance circulating in the maternal organism. Exact information about placental transfer can help to better estimate the toxic potential of a substance. From an evolutionary point of view, the chorioallantoic placenta is the organ with the highest anatomical diversity among species. Moreover, frequently used animal models in reproduction belong to rodents and lagomorphs, two groups that are characterized by the generation of an additional type of placenta, which is crucial for fetal development, but absent from humans: the inverted yolk sac placenta. Taken together, the translatability of placental transfer studies from laboratory animals to humans is challenging, which is supported by the fact that numerous species-dependent toxic effects are described in literature. Thus, reliable human-relevant data are frequently lacking and the toxic potential of chemicals and pharmaceuticals for humans can hardly be estimated, often resulting in recommendations that medical treatments or exposure to chemicals should be avoided for safety reasons. Although species differences of placental anatomy have been described frequently and the need for human-relevant research models has been emphasized, analyses of substances with species-dependent placental transfer have been performed only sporadically. Here, we present examples for species-specific placental transfer, including that of nanoparticles and pharmaceuticals, and discuss potential underlying mechanisms. With respect to the COVID 19-pandemic it might be of interest that some antiviral drugs are reported to feature species-specific placental transfer. Further, differences in placental structure and antibody transfer may affect placental transfer of ZIKA virus.

The role of globins in cardiovascular physiology

Globin proteins exist in every cell type of the vasculature, from erythrocytes to endothelial cells, vascular smooth muscle cells, and peripheral nerve cells. Many globin subtypes are also expressed in muscle tissues (including cardiac and skeletal muscle), in other organ-specific cell types, and in cells of the central nervous system. The ability of each of these globins to interact with molecular oxygen (O₂) and nitric oxide (NO) is preserved across these contexts. Endothelial α-globin is an example of extra-erythrocytic globin expression. Other globins, including myoglobin, cytoglobin, and neuroglobin are observed in other vascular tissues. Myoglobin is observed primarily in skeletal muscle and smooth muscle cells surrounding the aorta or other large arteries. Cytoglobin is found in vascular smooth muscle but can also be expressed in non-vascular cell types, especially in oxidative stress conditions after ischemic insult. Neuroglobin was first observed in neuronal cells, and its expression appears to be restricted mainly to the central and peripheral nervous systems. Brain and central nervous system neurons expressing neuroglobin are positioned close to many arteries within the brain parenchyma and can control smooth muscle contraction and, thus, tissue perfusion and vascular reactivity. Overall, reactions between NO and globin heme-iron contribute to vascular homeostasis by regulating vasodilatory NO signals and scaveging reactive species in cells of the mammalian vascular system. Here, we discuss how globin proteins affect vascular physiology with a focus on NO biology, and offer perspectives for future study of these functions.

Gene duplication and subsequent diversification strongly affect phenotypic evolvability and robustness

We study the effects of non-determinism and gene duplication on the structure of genotype-phenotype (GP) maps by introducing a non-deterministic version of the Polyomino self-assembly model. This model has previously been used in a variety of contexts to model the assembly and evolution of protein quaternary structure. Firstly, we show the limit of the current deterministic paradigm which leads to built-in anti-correlation between evolvability and robustness at the genotypic level. We develop a set of metrics to measure structural properties of GP maps in a non-deterministic setting and use them to evaluate the effects of gene duplication and subsequent diversification. Our generalized versions of evolvability and robustness exhibit positive correlation for a subset of genotypes. This positive correlation is only possible because non-deterministic phenotypes can contribute to both robustness and evolvability. Secondly, we show that duplication increases robustness and reduces evolvability initially, but that the subsequent diversification that duplication enables has a stronger, inverse effect, greatly increasing evolvability and reducing robustness relative to their original values.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

Integrating genomics with biomarkers and therapeutic targets to invigorate cardiovascular drug development

Drug development in cardiovascular disease is stagnating, with lack of efficacy and adverse effects being barriers to innovation. Human genetics can provide compelling evidence of causation through approaches such as Mendelian randomization, with genetic support for causation increasing the probability of a clinical trial succeeding. Mendelian randomization applied to quantitative traits can identify risk factors for disease that are both causal and amenable to therapeutic modification. However, important differences exist between genetic investigations of a biomarker (such as HDL cholesterol) and a drug target aimed at modifying the same biomarker of interest (such as cholesteryl ester transfer protein), with implications for the methodology, interpretation and application of Mendelian randomization to drug development. Differences include the comparative nature of the genetic architecture - that is, biomarkers are typically polygenic, whereas protein drug targets are influenced by either cis-acting or trans-acting genetic variants - and the potential for drug targets to show disease associations that might differ from those of the biomarker that they are intended to modify (target-mediated pleiotropy). In this Review, we compare and contrast the use of Mendelian randomization to evaluate potential drug targets versus quantitative traits. We explain how genetic epidemiological studies can be used to assess the aetiological roles of biomarkers in disease and to prioritize drug targets, including designing their evaluation in clinical trials.

Expression of Myoglobin in Normal and Cancer Brain Tissues: Correlation With Hypoxia Markers

Background

Myoglobin (MB) is increasingly recognized as a key player in cancer growth and metastasis. Low oxygen tensions, commonly associated with highly aggressive and recurrent cancers, have been shown to regulate its expression in several cancers such as lung, neck, prostate and breast cancer. However, it is not yet known whether it contributes to the growth and spread of brain cancers especially Glioblastoma multiforme (GBM).

Methods

Here we investigate the expression of MB, and its correlation with the hypoxia markers carbonic anhydrase IX (CAIX) and lactate dehydrogenase A (LDHA), in human tissue microarrays of multiple organ tumors, brain tumors, and GBM tumors, and their respective cancer-adjacent normal tissues. Correlation between MB protein expression and tumor grade was also assessed.

Results

We show that MB protein is expressed in a wide variety of cancers, benign tumors, cancer-adjacent normal tissues, hyperplastic tissue samples and normal brain tissue, and low oxygen tensions modulate MB protein expression in different brain cancers, including GBM. Enhanced nuclear LDHA immune-reactivity in GBM was also observed. Finally, we report for the first time a positive correlation between MB expression and brain tumor grade.

Conclusion

Our data suggest that hypoxia regulate MB expression in different brain cancers (including GBM) and that its expression is associated with a more aggressive phenotype as indicated by the positive correlation with the brain tumor grade. Additionally, a role for nuclear LDHA in promoting aggressive tumor phenotype is also suggested based on enhanced nuclear expression which was observed only in GBM.

Multimerization variants as potential drivers of neofunctionalization

Whole-genome duplications are common during evolution, creating genetic redundancy that can enable cellular innovations. Novel protein-protein interactions provide a route to diversified gene functions, but, at present, there is limited proteome-scale knowledge on the extent to which variability in protein complex formation drives neofunctionalization. Here, we used protein correlation profiling to test for variability in apparent mass among thousands of orthologous proteins isolated from diverse species and cell types. Variants in protein complex size were unexpectedly common, in some cases appearing after relatively recent whole-genome duplications or an allopolyploidy event. In other instances, variants such as those in the carbonic anhydrase orthologous group reflected the neofunctionalization of ancient paralogs that have been preserved in extant species. Our results demonstrate that homo- and heteromer formation have the potential to drive neofunctionalization in diverse classes of enzymes, signaling, and structural proteins.

Whole blood and blood components from vertebrates differentially affect egg formation in three species of anautogenous mosquitoes

Background

Most female mosquitoes are anautogenous and must blood feed on a vertebrate host to produce eggs. Prior studies show that the number of eggs females lay per clutch correlates with the volume of blood ingested and that protein is the most important macronutrient for egg formation. In contrast, how whole blood, blood fractions and specific blood proteins from different vertebrates affect egg formation is less clear. Since egg formation is best understood in Aedes aegypti, we examined how blood and blood components from different vertebrates affect this species and two others: the malaria vector Anopheles gambiae and arbovirus vector Culex quinquefasciatus.

Methods

Adult female mosquitoes were fed blood, blood fractions and purified major blood proteins from different vertebrate hosts. Markers of reproductive response including ovary ecdysteroidogenesis, yolk deposition into oocytes and number of mature eggs produced were measured.

Results

Ae. aegypti, An. gambiae and C. quinquefasciatus responded differently to meals of whole blood, plasma or blood cells from human, rat, chicken and turkey hosts. We observed more similarities between the anthropophiles Ae. aegypti and An. gambiae than the ornithophile C. quinquefasciatus. Focusing on Ae. aegypti, the major plasma-derived proteins (serum albumin, fibrinogen and globulins) differentially stimulated egg formation as a function of vertebrate host source. The major blood cell protein, hemoglobin, stimulated yolk deposition when from pigs but not humans, cows or sheep. Serum albumins from different vertebrates also variably affected egg formation. Bovine serum albumin (BSA) stimulated ovary ecdysteroidogenesis, but more weakly induced digestive enzyme activities than whole blood. In contrast, BSA-derived peptides and free amino acids had no stimulatory effects on ecdysteroidogenesis or yolk deposition into oocytes.

Conclusions

Whole blood, blood fractions and specific blood proteins supported egg formation in three species of anautogenous mosquitoes but specific responses varied with the vertebrate source of the blood components tested.

Genome-wide analysis of pseudogenes reveals HBBP1's human-specific essentiality in erythropoiesis and implication in β-thalassemia

The human genome harbors 14,000 duplicated or retroposed pseudogenes. Given their functionality as regulatory RNAs and low conservation, we hypothesized that pseudogenes could shape human-specific phenotypes. To test this, we performed co-expression analyses and found that pseudogene exhibited tissue-specific expression, especially in the bone marrow. By incorporating genetic data, we identified a bone-marrow-specific duplicated pseudogene, HBBP1 (η-globin), which has been implicated in β-thalassemia. Extensive functional assays demonstrated that HBBP1 is essential for erythropoiesis by binding the RNA-binding protein (RBP), HNRNPA1, to upregulate TAL1, a key regulator of erythropoiesis. The HBBP1/TAL1 interaction contributes to a milder symptom in β-thalassemia patients. Comparative studies further indicated that the HBBP1/TAL1 interaction is human-specific. Genome-wide analyses showed that duplicated pseudogenes are often bound by RBPs and less commonly bound by microRNAs compared with retropseudogenes. Taken together, we not only demonstrate that pseudogenes can drive human evolution but also provide insights on their functional landscapes.

Detection of Common Deletional of α-Thalassemia 3.7 Kb from Metropolitan Region of Manaus, Amazonas, Brazil

Background

Alpha Thalassemia (α-thal) is a heterogeneous group of hereditary alterations caused by deletions that affect alpha regulatory genes, and the 3.7Kb deletion is the most frequent worldwide. The prevalence ranges from 20% and 35% in Brazil, depending mainly on race, predominant in Afro-descendants.

Purpose

The aim was to determine α-thal -^α3.7Kb and -^α4.2Kb deletions, estimating their frequency in individuals from six regions of Amazonas State.

Methods

Volunteers age between 18-59 years old of both genders participated in the study. Blood was collected from March 2014 to September 2017 at the health centers of each participant city. α-thal^3.7Kb was performed by GAP-PCR, while α-thal^4.2Kb by Multiplex-PCR. The total samples collected from each city were: Manaus (capital), 356 (19.7%); Iranduba 232 (12.8%); Manacapuru, 287 (15.9%); Presidente Figueiredo, 370 (20.5%); Itacoatiara, 301 (16.6%); and Coari, 263 (14.5%).

Results

The average age among males was 35.3±14.8, while for females, it was 36.7±14.9 years old. Microcytosis (MCV <80fL) was found in 158 individuals (8,46%) and α-thal diagnosed in 143 individuals (7.9%), and all of these individuals carried the ^3.7Kb deletion 5.95% in heterozygous and 1.95% in homozygous. α-thal^4.2kb was not found in any volunteer. The association analyses to the α-thal^3.7kb genotypes were statistically significant for all hematological parameters (p<.001), except serum iron and serum ferritin analyses.

Conclusion

This study highlights α-thal 3.7kb deletion as an important public health problem, especially in a population not yet characterized about this disease. Thus, epidemiological studies using molecular tools become relevant in regions where the disease is underestimated, contributing to a better understanding of thalassemia incidence and iron deficiency anemias incidence of the participating cities. We reinforce that future molecular studies in North Region from Brazil can be utilized to describe other genetic anemias as structural hemoglobinopathies that have already proven to be highly prevalent in Brazil.

Globins in the marine annelid Platynereis dumerilii shed new light on hemoglobin evolution in bilaterians

BACKGROUND

How vascular systems and their respiratory pigments evolved is still debated. While many animals present a vascular system, hemoglobin exists as a blood pigment only in a few groups (vertebrates, annelids, a few arthropod and mollusk species). Hemoglobins are formed of globin sub-units, belonging to multigene families, in various multimeric assemblages. It was so far unclear whether hemoglobin families from different bilaterian groups had a common origin.

RESULTS

To unravel globin evolution in bilaterians, we studied the marine annelid Platynereis dumerilii, a species with a slow evolving genome. Platynereis exhibits a closed vascular system filled with extracellular hemoglobin. Platynereis genome and transcriptomes reveal a family of 19 globins, nine of which are predicted to be extracellular. Extracellular globins are produced by specialized cells lining the vessels of the segmental appendages of the worm, serving as gills, and thus likely participate in the assembly of a previously characterized annelid-specific giant hemoglobin. Extracellular globin mRNAs are absent in smaller juveniles, accumulate considerably in growing and more active worms and peak in swarming adults, as the need for O2 culminates. Next, we conducted a metazoan-wide phylogenetic analysis of globins using data from complete genomes. We establish that five globin genes (stem globins) were present in the last common ancestor of bilaterians. Based on these results, we propose a new nomenclature of globins, with five clades. All five ancestral stem-globin clades are retained in some spiralians, while some clades disappeared early in deuterostome and ecdysozoan evolution. All known bilaterian blood globin families are grouped in a single clade (clade I) together with intracellular globins of bilaterians devoid of red blood.

CONCLUSIONS

We uncover a complex "pre-blood" evolution of globins, with an early gene radiation in ancestral bilaterians. Circulating hemoglobins in various bilaterian groups evolved convergently, presumably in correlation with animal size and activity. However, all hemoglobins derive from a clade I globin, or cytoglobin, probably involved in intracellular O2 transit and regulation. The annelid Platynereis is remarkable in having a large family of extracellular blood globins, while retaining all clades of ancestral bilaterian globins.

MZF1 regulates α-globin gene transcription via long-range interactions in erythroid differentiation

Precise spatiotemporal gene expression regulation is crucial for human erythropoiesis. However, dramatic changes in the chromatin structure and transcriptome involved in α-globin gene expression during erythropoiesis still not fully understand. To identify candidate regulators for α-globin gene regulation, we carried out an integrated approach by integrating publicly available transcriptomic and epigenomic data. We computed active enhancers by overlapping enriched regions marked with H3K4me1 and H3K27ac and correlated their activity with mRNA expression. Next, we cataloged potential transcription factors via de novo motif analysis. We highlighted the discovery of potential novel transcription factor MZF1 of the α-globin gene in erythroid differentiation. To validate the role of MZF1, we quantified the expression level of MZF1 and α-globin gene in HSPCs, early erythroid progenitors and late erythroid precursors cells. Both the mRNA and protein expression patterns of MZF1 were consistent with the α-globin gene. Also, the qPCR result showed that the expression of the α-globin gene was significantly increased by the MZF1 overexpression. To further investigate the role of MZF1 regulating α-globin gene transcriptional activity during erythroid differentiation, we performed ChIP-qPCR at the α-globin locus. Our results showed that MZF1 recruitment both at 4 upstream HS sites and α-globin gene promoter in erythroid precursor cells. To determine the importance of the MZF1 to enhancer-promoter interaction at the α-globin locus, we compared interaction frequency before and after knockdown of MZF1 by chromosome conformation capture (3C) assay. Upon MZF1 depletion, both the expression of the α-globin gene and all 3C signals were significantly decreased. Taken together, MZF1 plays an important role in regulating α-globin gene expression by binding to long-region enhancers and α-globin gene promoter and facilitates the organization of specific 3D chromatin architecture in erythroid differentiation.

Mechanisms mediating suppression of globin gene transcription in Danio rerio nonerythroid cells

We studied the repression of adult and embryo-larval genes of the major globin gene locus in D. rerio fibroblasts. The results obtained suggest that at least some of the globin genes are repressed by Polycomb, similarly to human α-globin genes. Furthermore, within two α/β globin gene pairs, repression of α-type and β-type genes appears to be mediated by different mechanisms, as increasing the level of histone acetylation can activate transcription of only β-type genes.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

LRF/ZBTB7A conservation accentuates its potential as a therapeutic target for the hematopoietic disorders

Conserved sequences across species have always provided valuable insights to improve our understanding on the human genome's entity and the interplay among different loci. Lymphoma/leukemia related factor (LRF) is encoded by ZBTB7A gene and belongs to an evolutionarily conserved family of transcription factors, implicated in vital cellular functions. The present data, demonstrating the wide-spread and the high overlap of the LRF/ZBTB7A recognition sites with genomic segments identified as CpG islands in the human genome, suggest that its binding capacity strongly depends on a specific sequence-encoded feature within CpGs. We have previously shown that de-methylation of the CpG island 326 lying in the ZBTB7A gene promoter is associated with impaired pharmacological induction of fetal hemoglobin in β-type hemoglobinopathies patients. Within this context we aimed to investigate the extent of the LRF/ZBTB7A conservation among primates and mouse genome, focusing our interest also on the CpG island flanking the gene's promoter region, in an effort to further establish its epigenetic regulatory role in human hematopoiesis and pharmacological involvement in hematopoietic disorders. Comparative analysis of the human ZBTB7A nucleotide and amino acid sequences and orthologous sequences among non-human primates and mouse, exhibited high conservation scores. Pathway analysis, clearly indicated that LRF/ZBTB7A influences conserved cellular processes. These data in conjunction with the high levels of expression foremost in hematopoietic tissues, highlighted LRF/ZBTB7A as an essential factor operating indisputably during hematopoiesis.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Enhancing droplet-based single-nucleus RNA-seq resolution using the semi-supervised machine learning classifier DIEM

Single-nucleus RNA sequencing (snRNA-seq) measures gene expression in individual nuclei instead of cells, allowing for unbiased cell type characterization in solid tissues. We observe that snRNA-seq is commonly subject to contamination by high amounts of ambient RNA, which can lead to biased downstream analyses, such as identification of spurious cell types if overlooked. We present a novel approach to quantify contamination and filter droplets in snRNA-seq experiments, called Debris Identification using Expectation Maximization (DIEM). Our likelihood-based approach models the gene expression distribution of debris and cell types, which are estimated using EM. We evaluated DIEM using three snRNA-seq data sets: (1) human differentiating preadipocytes in vitro, (2) fresh mouse brain tissue, and (3) human frozen adipose tissue (AT) from six individuals. All three data sets showed evidence of extranuclear RNA contamination, and we observed that existing methods fail to account for contaminated droplets and led to spurious cell types. When compared to filtering using these state of the art methods, DIEM better removed droplets containing high levels of extranuclear RNA and led to higher quality clusters. Although DIEM was designed for snRNA-seq, our clustering strategy also successfully filtered single-cell RNA-seq data. To conclude, our novel method DIEM removes debris-contaminated droplets from single-cell-based data fast and effectively, leading to cleaner downstream analysis. Our code is freely available for use at https://github.com/marcalva/diem.

A Modular RNA Domain That Confers Differential Ligand Specificity

The modularity of protein domains is well-known, but the existence of independent domains that confer function in RNA is less established. Recently, a family of RNA aptamers termed ykkC was discovered; they bind at least four ligands of very different chemical composition, including guanidine, phosphoribosyl pyrophosphate (PRPP), and guanosine tetraphosphate (ppGpp) (graphical abstract). Structures of these aptamers revealed an architecture characterized by two coaxial helical stacks. The first helix appears to be a generic scaffold, while the second helix forms the most contacts to the ligands. To determine if these two regions within the aptamer are modular units for ligand recognition, we swapped the ligand-binding coaxial stacks of a guanidine aptamer and a PRPP aptamer. This operation, in combination with a single mutation in the scaffold domain, achieved full switching of ligand specificity. This finding suggests that the ligand-binding helix largely dictates the ligand specificity of ykkC RNAs and that the scaffold coaxial stack is generally compatible with various ykkC ligand-binding modules. This work presents an example of RNA domain modularity comparable to that of a ligand-binding protein, showcasing the versatility of RNA as an entity capable of molecular evolution through adaptation of existing motifs.