Evolutionary constraints in the β-globin cluster: the signature of purifying selection at the δ-globin (HBD) locus and its role in developmental gene regulation

CRISPR-directed gene-editing induces genetic rearrangement within the human globin gene locus

CRISPR-Cas is a revolutionary technology but has already demonstrated significant feasibility for clinical and non-clinical applications. While the efficiency and precision of this remarkable genetic tool is unprecedented, unfortunately, a series of collateral genetic rearrangement have been reported in response to double-stranded DNA breakage. Once these molecular scissions occur, the cascade of DNA repair reactions can lead to genomic rearrangements especially if breakage takes place within a family of sequence related genes. Here, we demonstrate that CRISPR- directed gene editing near the sickle cell mutation site generates a curious genetic outcome; a footprint of the δ globin gene proximal to the CRISPR/Cas cut site(s). This rearrangement is not dependent on the presence of an exogenously added DNA template but is apparently dependent on a double strand break. Our results the highlight recombinational capacity of double strand breaks in human chromosomes where the aim is to edit a human gene.

Secondary structure analysis of proteins within the same topology group

The native conformation of a protein plays a decisive role in ensuring its functionality. It is established that the spatial structure of proteins may exhibit a greater degree of conservation than the corresponding amino acid sequences. This study aims to clarify structural distinctions between homologous and non-homologous proteins with identical topology. The analysis focuses on secondary structures with special emphasis on their fraction, distribution along the polypeptide chain, and chirality. Three different groups of proteins with identical topology were considered according to the CATH database: a homologous group of Globins, a group of Phycocyanins, which is often considered as a potential relative of globins, and a diverse assembly of other globin-like proteins. Some structural patterns in the distribution of secondary structure have been identified within Globins. A similar profile was observed in Phycocyanins, in contrast to the third group. In addition, a distinguishable structural motif, including structures such as 310-helix and irregular structure, has been found in both Globins and Phycocyanins, which can be proposed as an evolutionary imprint.

Precise correction of a spectrum of β-thalassemia mutations in coding and non-coding regions by base editors

β-thalassemia/HbE results from mutations in the β-globin locus that impede the production of functional adult hemoglobin. Base editors (BEs) could facilitate the correction of the point mutations with minimal or no indel creation, but its efficiency and bystander editing for the correction of β-thalassemia mutations in coding and non-coding regions remains unexplored. Here, we screened BE variants in HUDEP-2 cells for their ability to correct a spectrum of β-thalassemia mutations that were integrated into the genome as fragments of HBB. The identified targets were introduced into their endogenous genomic location using BEs and Cas9/homology-directed repair (HDR) to create cellular models with β-thalassemia/HbE. These β-thalassemia/HbE models were then used to assess the efficiency of correction in the native locus and functional β-globin restoration. Most bystander edits produced near target sites did not interfere with adult hemoglobin expression and are not predicted to be pathogenic. Further, the effectiveness of BE was validated for the correction of the pathogenic HbE variant in severe β0/βE-thalassaemia patient cells. Overall, our study establishes a novel platform to screen and select optimal BE tools for therapeutic genome editing by demonstrating the precise, efficient, and scarless correction of pathogenic point mutations spanning multiple regions of HBB including the promoter, intron, and exons.

Engineering of the endogenous HBD promoter increases HbA2

The β-hemoglobinopathies, such as sickle cell disease and β-thalassemia, are one of the most common genetic diseases worldwide and are caused by mutations affecting the structure or production of β-globin subunits in adult hemoglobin. Many gene editing efforts to treat the β-hemoglobinopathies attempt to correct β-globin mutations or increase γ-globin for fetal hemoglobin production. δ-globin, the subunit of adult hemoglobin A2, has high homology to β-globin and is already pan-cellularly expressed at low levels in adult red blood cells. However, upregulation of δ-globin is a relatively unexplored avenue to increase the amount of functional hemoglobin. Here, we use CRISPR-Cas9 to repair non-functional transcriptional elements in the endogenous promoter region of δ-globin to increase overall expression of adult hemoglobin 2 (HbA2). We find that insertion of a KLF1 site alone is insufficient to upregulate δ-globin. Instead, multiple transcription factor elements are necessary for robust upregulation of δ-globin from the endogenous locus. Promoter edited HUDEP-2 immortalized erythroid progenitor cells exhibit striking increases of HBD transcript, from less than 5% to over 20% of total β-like globins in clonal populations. Edited CD34 +hematopoietic stem and progenitors (HSPCs) differentiated to primary human erythroblasts express up to 46% HBD in clonal populations. These findings add mechanistic insight to globin gene regulation and offer a new therapeutic avenue to treat β-hemoglobinopathies.

Ancestral reconstruction and correlation of the frequencies of the hemoglobin S allele and the Duffy blood group alleles in human populations

OBJECTIVES

Malaria is an important selective force for human genetic adaptation due to the sustained, lethal impact it has had on populations worldwide. High frequencies of both hemoglobin S and the null allele FYB^ES of the Duffy blood group have been found in areas where this disease is endemic, attributed to the protective action of the carriers of these variants against malaria infection. The objective of this work was to perform ancestral reconstruction and analyze the correlation of the frequencies of these alleles throughout the phylogeny of 24 human populations.

METHODS

A tree topology and the allelic frequencies reported in the literature for the 24 populations were used. The ancestral frequencies for the two alleles were reconstructed using the maximum likelihood method and the Brownian model of evolution (CI = 95%), and the correlation analysis was performed using phylogenetically independent contrasts (PICs). Statistical analyses were performed with the statistical software R version 3.4.1.

RESULTS

For both alleles, a correspondence was found in the reconstruction of the ancestral frequencies, and a significant statistical correlation (p = .001) was observed between the S and FYB^ES alleles.

CONCLUSIONS

These results provide evidence of an epistatic relationship between the two alleles, which may influence the fitness of the individuals who present with them when they are subjected to a selective force such as malaria.

Hb F Levels in β-Thalassemia Carriers and Normal Individuals: Known and Unknown Quantitative Trait Loci in the β-Globin Gene Cluster

In the already identified quantitative trait loci (QTL), modulating Hb F levels are cis-acting haplotypes of the β-globin gene cluster itself, although the single nucleotide polymorphisms (SNPs) accounting more for the association, remain uncertain. In this study, the role in Hb F production of previously reported candidate SNPs within the β-globin gene cluster was reexamined, along with a yet poorly studied variation in the BGLT3 gene. In a sample of β-thalassemia (β-thal) carriers, we succeeded in replicating the significant association between increased Hb F levels and rs7482144 (C>T) (HBG2 XmnI), which is the most well-established variation in the cluster influencing the trait. This SNP was found to be in strong linkage disequilibrium (LD) with a variation in the HBBP1 gene [rs10128556 (G>A)], which consistently revealed a similar association signal. Remarkably, much stronger than the latter associations were those involving both rs968857 (T allele) (3' HBBP1) and rs7924684 (G allele) (BGLT3), two SNPs that were also in strong LD. As the pattern of LD detected in the β-globin gene cluster does not correlate with a tight linkage between markers, complex interactions between SNPs at the cluster seem to modulate Hb F. Seeing that no such associations were detected in normal subjects, the question can be raised on whether, under erythropoiesis stress, epigenetic mechanisms contribute to change the regulation of the entire β-globin gene cluster. In conclusion, we provide statistical evidence for a new player within the β-globin gene cluster, BGLT3, that in cooperation with other regions influences Hb F levels in β-thal carriers.

Molecular genetics of β-thalassemia: A narrative review

ABSTRACT

β-thalassemia is a hereditary hematological disease caused by over 350 mutations in the β-globin gene (HBB). Identifying the genetic variants affecting fetal hemoglobin (HbF) production combined with the α-globin genotype provides some prediction of disease severity for β-thalassemia. However, the generation of an additive composite genetic risk score predicts prognosis, and guide management requires a larger panel of genetic modifiers yet to be discovered.Presently, using data from prior clinical trials guides the design of further research and academic studies based on gene augmentation, while fundamental insights into globin switching and new technology developments have inspired the investigation of novel gene therapy approaches.Genetic studies have successfully characterized the causal variants and pathways involved in HbF regulation, providing novel therapeutic targets for HbF reactivation. In addition to these HBB mutation-independent strategies involving HbF synthesis de-repression, the expanding genome editing toolkit provides increased accuracy to HBB mutation-specific strategies encompassing adult hemoglobin restoration for personalized treatment of hemoglobinopathies. Allogeneic hematopoietic stem cell transplantation was, until very recently, the curative option available for patients with transfusion-dependent β-thalassemia. Gene therapy currently represents a novel therapeutic promise after many years of extensive preclinical research to optimize gene transfer protocols.We summarize the current state of developments in the molecular genetics of β-thalassemia over the last decade, including the mechanisms associated with ineffective erythropoiesis, which have also provided valid therapeutic targets, some of which have been shown as a proof-of-concept.

Characterization of Gene Alterations following Editing of the β-Globin Gene Locus in Hematopoietic Stem/Progenitor Cells

The use of engineered nucleases combined with a homologous DNA donor template can result in targeted gene correction of the sickle cell disease mutation in hematopoietic stem and progenitor cells. However, because of the high homology between the adjacent human β- and δ-globin genes, off-target cleavage is observed at δ-globin when using some endonucleases targeted to the sickle mutation in β-globin. Introduction of multiple double-stranded breaks by endonucleases has the potential to induce intergenic alterations. Using a novel droplet digital PCR assay and high-throughput sequencing, we characterized the frequency of rearrangements between the β- and δ-globin paralogs when delivering these nucleases. Pooled CD34+ cells and colony-forming units from sickle bone marrow were treated with nuclease only or including a donor template and then analyzed for potential gene rearrangements. It was observed that, in pooled CD34+ cells and colony-forming units, the intergenic β-δ-globin deletion was the most frequent rearrangement, followed by inversion of the intergenic fragment, with the inter-chromosomal translocation as the least frequent. No rearrangements were observed when endonuclease activity was restricted to on-target β-globin cleavage. These findings demonstrate the need to develop site-specific endonucleases with high specificity to avoid unwanted gene alterations.

CKS1BP7, a Pseudogene of CKS1B, is Co-Amplified with IGF1R in Breast Cancers

Pseudogenes have been reported to exhibit functional roles. Amplification or overexpression of CDC28 protein kinase regulatory subunit 1B (CKS1B) was found in various human cancers. But it was known little about CKS1B pseudogene 7 (CKS1BP7), a pseudogene sharing considerable sequence identity with CKS1B. The aim of this study was to evaluate copy number alterations (CNAs) of CKS1BP7 and address its potential roles in breast cancer. We detected copy numbers of CKS1BP7 and insulin-like growth factor 1 receptor (IGF1R) using quantitative multi-gene fluorescence in situ hybridization (QM-FISH) technique, compared their status in both invasive carcinoma and ductal carcinoma in situ (DCIS) components within the same tumors, and investigated the associations of CNAs with tumor features and patients outcomes. Amplification of CKS1BP7 (dot-like pattern) was found in 28.8% of all cases, while amplified IGF1R (cluster pattern) was identified in 24.2% of all patients. The two events often co-existed (p = 0.01). Within the same tumors, identical CNAs of CKS1BP7 and IGF1R were found in DCIS and invasive carcinoma. Moreover, amplification of both genes was more frequent in aneuploidy tumors and the tumors with high ki67, but wasn't associated with patients' outcome. In summary, CKS1BP7 amplification is a frequent event in breast cancer and often co-occurs with amplified IGF1R, which provides evidence supporting the interactions between CKS1BP7 and IGF1R during mammary carcinogenesis. Our findings suggest that CKS1BP7 as well as IGF1R may serve as potential biomarkers for early detection and predict prognosis in breast cancer.

Comparative genomics of canine hemoglobin genes reveals primacy of beta subunit delta in adult carnivores

Background

The main function of hemoglobin (Hb) is to transport oxygen in the circulation. It is among the most highly studied proteins due to its roles in physiology and disease, and most of our understanding derives from comparative research. There is great diversity in Hb gene evolution in placental mammals, mostly in the repertoire and regulation of the β-globin subunits. Dogs are an ideal model in which to study Hb genes because: 1) they are members of Laurasiatheria, our closest relatives outside of Euarchontoglires (including primates, rodents and rabbits), 2) dog breeds are isolated populations with their own Hb-associated genetics and diseases, and 3) their high level of health care allows for development of biomedical investigation and translation.

Results

We established that dogs have a complement of five α and five β-globin genes, all of which can be detected as spliced mRNA in adults. Strikingly, HBD, the allegedly-unnecessary adult β-globin protein in humans, is the primary adult β-globin in dogs and other carnivores; moreover, dogs have two active copies of the HBD gene. In contrast, the dominant adult β-globin of humans, HBB, has high sequence divergence and is expressed at markedly lower levels in dogs. We also showed that canine HBD and HBB genes are complex chimeras that resulted from multiple gene conversion events between them. Lastly, we showed that the strongest signal of evolutionary selection in a high-altitude breed, the Bernese Mountain Dog, lies in a haplotype block that spans the β-globin locus.

Conclusions

We report the first molecular genetic characterization of Hb genes in dogs. We found important distinctions between adult β-globin expression in carnivores compared to other members of Laurasiatheria. Our findings are also likely to raise new questions about the significance of human HBD. The comparative genomics of dog hemoglobin genes sets the stage for diverse research and translation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3513-0) contains supplementary material, which is available to authorized users.

Extracellular hemoglobin: the case of a friend turned foe

Hemoglobin (Hb) is a highly conserved molecule present in all life forms and functionally tied to the complexity of aerobic organisms on earth in utilizing oxygen from the atmosphere and delivering to cells and tissues. This primary function sustains the energy requirements of cells and maintains cellular homeostasis. Decades of intensive research has presented a paradigm shift that shows how the molecule also functions to facilitate smooth oxygen delivery through the cardiovascular system for cellular bioenergetic homeostasis and signaling for cell function and defense. These roles are particularly highlighted in the binding of Hb to gaseous molecules carbon dioxide (CO₂), nitric oxide (NO) and carbon monoxide (CO), while also serving indirectly or directly as sources of these signaling molecules. The functional activities impacted by Hb outside of bioenergetics homeostasis, include fertilization, signaling functions, modulation of inflammatory responses for defense and cell viability. These activities are efficiently executed while Hb is sequestered safely within the confines of the red blood cell (rbc). Outside of rbc confines, Hb disaggregates and becomes a danger molecule to cell survival. In these perpectives, Hb function is broadly dichotomous, either a friend in its natural environment providing and facilitating the means for cell function or foe when dislocated from its habitat under stress or pathological condition disrupting cell function. The review presents insights into how this dichotomy in function manifests.

Distinctive patterns of evolution of the δ-globin gene (HBD) in primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.

DivStat: a user-friendly tool for single nucleotide polymorphism analysis of genomic diversity

Recent developments have led to an enormous increase of publicly available large genomic data, including complete genomes. The 1000 Genomes Project was a major contributor, releasing the results of sequencing a large number of individual genomes, and allowing for a myriad of large scale studies on human genetic variation. However, the tools currently available are insufficient when the goal concerns some analyses of data sets encompassing more than hundreds of base pairs and when considering haplotype sequences of single nucleotide polymorphisms (SNPs). Here, we present a new and potent tool to deal with large data sets allowing the computation of a variety of summary statistics of population genetic data, increasing the speed of data analysis.

Repeated evolution of chimeric fusion genes in the β-globin gene family of laurasiatherian mammals

The evolutionary fate of chimeric fusion genes may be strongly influenced by their recombinational mode of origin and the nature of functional divergence between the parental genes. In the β-globin gene family of placental mammals, the two postnatally expressed δ- and β-globin genes (HBD and HBB, respectively) have a propensity for recombinational exchange via gene conversion and unequal crossing-over. In the latter case, there are good reasons to expect differences in retention rates for the reciprocal HBB/HBD and HBD/HBB fusion genes due to thalassemia pathologies associated with the HBD/HBB "Lepore" deletion mutant in humans. Here, we report a comparative genomic analysis of the mammalian β-globin gene cluster, which revealed that chimeric HBB/HBD fusion genes originated independently in four separate lineages of laurasiatherian mammals: Eulipotyphlans (shrews, moles, and hedgehogs), carnivores, microchiropteran bats, and cetaceans. In cases where an independently derived "anti-Lepore" duplication mutant has become fixed, the parental HBD and/or HBB genes have typically been inactivated or deleted, so that the newly created HBB/HBD fusion gene is primarily responsible for synthesizing the β-type subunits of adult and fetal hemoglobin (Hb). Contrary to conventional wisdom that the HBD gene is a vestigial relict that is typically inactivated or expressed at negligible levels, we show that HBD-like genes often encode a substantial fraction (20-100%) of β-chain Hbs in laurasiatherian taxa. Our results indicate that the ascendancy or resuscitation of genes with HBD-like coding sequence requires the secondary acquisition of HBB-like promoter sequence via unequal crossing-over or interparalog gene conversion.