Molecular evolution of GYPC: evidence for recent structural innovation and positive selection in humans

Genetic variation of glycophorins and infectious disease

Glycophorins are transmembrane proteins of red blood cells (RBCs), heavily glycosylated on their external-facing surface. In humans, there are four glycophorin proteins, glycophorins A, B, C and D. Glycophorins A and B are encoded by two similar genes GYPA and GYPB, and glycophorin C and glycophorin D are encoded by a single gene, GYPC. The exact function of glycophorins remains unclear. However, given their abundance on the surface of RBCs, it is likely that they serve as a substrate for glycosylation, giving the RBC a negatively charged, complex glycan “coat”. GYPB and GYPE (a closely related pseudogene) were generated from GYPA by two duplication events involving a 120-kb genomic segment between 10 and 15 million years ago. Non-allelic homologous recombination between these 120-kb repeats generates a variety of duplication alleles and deletion alleles, which have been systematically catalogued from genomic sequence data. One allele, called DUP4, encodes the Dantu NE blood type and is strongly protective against malaria as it alters the surface tension of the RBC membrane. Glycophorins interact with other infectious pathogens, including viruses, as well as the malarial parasite Plasmodium falciparum, but the role of glycophorin variation in mediating the effects of these pathogens remains underexplored.

Readthrough Errors Purge Deleterious Cryptic Sequences, Facilitating the Birth of Coding Sequences

De novo protein-coding innovations sometimes emerge from ancestrally noncoding DNA, despite the expectation that translating random sequences is overwhelmingly likely to be deleterious. The "preadapting selection" hypothesis claims that emergence is facilitated by prior, low-level translation of noncoding sequences via molecular errors. It predicts that selection on polypeptides translated only in error is strong enough to matter and is strongest when erroneous expression is high. To test this hypothesis, we examined noncoding sequences located downstream of stop codons (i.e., those potentially translated by readthrough errors) in Saccharomyces cerevisiae genes. We identified a class of "fragile" proteins under strong selection to reduce readthrough, which are unlikely substrates for co-option. Among the remainder, sequences showing evidence of readthrough translation, as assessed by ribosome profiling, encoded C-terminal extensions with higher intrinsic structural disorder, supporting the preadapting selection hypothesis. The cryptic sequences beyond the stop codon, rather than spillover effects from the regular C-termini, are primarily responsible for the higher disorder. Results are robust to controlling for the fact that stronger selection also reduces the length of C-terminal extensions. These findings indicate that selection acts on 3' UTRs in Saccharomyces cerevisiae to purge potentially deleterious variants of cryptic polypeptides, acting more strongly in genes that experience more readthrough errors.

X chromosome-dependent disruption of placental regulatory networks in hybrid dwarf hamsters

Embryonic development in mammals is highly sensitive to changes in gene expression within the placenta. The placenta is also highly enriched for genes showing parent-of-origin or imprinted expression, which is predicted to evolve rapidly in response to parental conflict. However, little is known about the evolution of placental gene expression, or if divergence of placental gene expression plays an important role in mammalian speciation. We used crosses between two species of dwarf hamsters (Phodopus sungorus and Phodopus campbelli) to examine the genetic and regulatory underpinnings of severe placental overgrowth in their hybrids. Using quantitative genetic mapping and mitochondrial substitution lines, we show that overgrowth of hybrid placentas was primarily caused by genetic differences on the maternally inherited P. sungorus X chromosome. Mitochondrial interactions did not contribute to abnormal hybrid placental development, and there was only weak correspondence between placental disruption and embryonic growth. Genome-wide analyses of placental transcriptomes from the parental species and first- and second-generation hybrids revealed a central group of co-expressed X-linked and autosomal genes that were highly enriched for maternally biased expression. Expression of this gene network was strongly correlated with placental size and showed widespread misexpression dependent on epistatic interactions with X-linked hybrid incompatibilities. Collectively, our results indicate that the X chromosome is likely to play a prominent role in the evolution of placental gene expression and the accumulation of hybrid developmental barriers between mammalian species.

Unannotated single nucleotide polymorphisms in the TATA box of erythropoiesis genes show in vitro positive involvements in cognitive and mental disorders

BACKGROUND

Hemoglobin is a tetramer consisting of two α-chains and two β-chains of globin. Hereditary aberrations in the synthesis of one of the globin chains are at the root of thalassemia, one of the most prevalent monogenic diseases worldwide. In humans, in addition to α- and β-globins, embryonic zeta-globin and fetal γ-globin are expressed. Immediately after birth, the expression of fetal Aγ- and Gγ-globin ceases, and then adult β-globin is mostly expressed. It has been shown that in addition to erythroid cells, hemoglobin is widely expressed in nonerythroid cells including neurons of the cortex, hippocampus, and cerebellum in rodents; embryonic and adult brain neurons in mice; and mesencephalic dopaminergic brain cells in humans, mice, and rats. Lately, there is growing evidence that different forms of anemia (changes in the number and quality of blood cells) may be involved in (or may accompany) the pathogenesis of various cognitive and mental disorders, such as Alzheimer's and Parkinson's diseases, depression of various severity levels, bipolar disorders, and schizophrenia. Higher hemoglobin concentrations in the blood may lead to hyperviscosity, hypovolemia, and lung diseases, which may cause brain hypoxia and anomalies of brain function, which may also result in cognitive deficits.

METHODS

In this study, a search for unannotated single-nucleotide polymorphisms (SNPs) of erythroid genes was initially performed using our previously created and published SNP-TATA_Z-tester, which is a Web service for computational analysis of a given SNP for in silico estimation of its influence on the affinity of TATA-binding protein (TBP) for TATA and TATA-like sequences. The obtained predictions were finally verified in vitro by an electrophoretic mobility shift assay (EMSA).

RESULTS

On the basis of these experimental in vitro results and literature data, we studied TATA box SNPs influencing both human erythropoiesis and cognitive abilities. For instance, TBP-TATA affinity in the HbZ promoter decreases 6.6-fold as a result of a substitution in the TATA box (rs113180943), thereby possibly disrupting stage-dependent events of "switching" of hemoglobin genes and thus causing erythroblastosis. Therefore, rs113180943 may be a candidate marker of severe hemoglobinopathies with comorbid cognitive and mental disorders associated with cerebral blood flow disturbances.

CONCLUSIONS

The literature data and experimental and computations results suggest that the uncovered candidate SNP markers of erythropoiesis anomalies may also be studied in cohorts of patients with cognitive and/or mental disorders with comorbid erythropoiesis diseases in comparison to conventionally healthy volunteers. Research into the regulatory mechanisms by which the identified SNP markers contribute to the development of hemoglobinopathies and of the associated cognitive deficits will allow physicians not only to take timely and adequate measures against hemoglobinopathies but also to implement strategies preventing cognitive and mental disorders.

Phospho-proteomic discovery of novel signal transducers including thioredoxin-interacting protein as mediators of erythropoietin-dependent human erythropoiesis

Erythroid cell formation critically depends on signals transduced via erythropoietin (EPO)/EPO receptor (EPOR)/JAK2 complexes. This includes not only core response modules (e.g., JAK2/STAT5, RAS/MEK/ERK), but also specialized effectors (e.g., erythroferrone, ASCT2 glutamine transport, Spi2A). By using phospho-proteomics and a human erythroblastic cell model, we identify 121 new EPO target proteins, together with their EPO-modulated domains and phosphosites. Gene ontology (GO) enrichment for "Molecular Function" identified adaptor proteins as one top EPO target category. This includes a novel EPOR/JAK2-coupled network of actin assemblage modifiers, with adaptors DLG-1, DLG-3, WAS, WASL, and CD2AP as prime components. "Cellular Component" GO analysis further identified 19 new EPO-modulated cytoskeletal targets including the erythroid cytoskeletal targets spectrin A, spectrin B, adducin 2, and glycophorin C. In each, EPO-induced phosphorylation occurred at pY sites and subdomains, which suggests coordinated regulation by EPO of the erythroid cytoskeleton. GO analysis of "Biological Processes" further revealed metabolic regulators as a likewise unexpected EPO target set. Targets included aldolase A, pyruvate dehydrogenase α1, and thioredoxin-interacting protein (TXNIP), with EPO-modulated p-Y sites in each occurring within functional subdomains. In TXNIP, EPO-induced phosphorylation occurred at novel p-T349 and p-S358 sites, and was paralleled by rapid increases in TXNIP levels. In UT7epo-E and primary human stem cell (HSC)-derived erythroid progenitor cells, lentivirus-mediated short hairpin RNA knockdown studies revealed novel pro-erythropoietic roles for TXNIP. Specifically, TXNIP's knockdown sharply inhibited c-KIT expression; compromised EPO dose-dependent erythroblast proliferation and survival; and delayed late-stage erythroblast formation. Overall, new insight is provided into EPO's diverse action mechanisms and TXNIP's contributions to EPO-dependent human erythropoiesis.

Few Plasmodium falciparum merozoite ligand and erythrocyte receptor pairs show evidence of balancing selection

Erythrocyte surface proteins have been identified as receptors of Plasmodium falciparum merozoite proteins. The ligand-receptor interactions enable the parasite to invade human erythrocytes, initiating the clinical symptoms of malaria. These interactions are likely to have had an evolutionary impact on the genes that encode the ligand and receptor proteins. We used sequence data from Kilifi, Kenya to detect departures from neutrality in a paired analysis of P. falciparum merozoite ligands and their erythrocyte receptor genes from the same population. We genotyped parasite and human DNA obtained from 93 individuals with severe malaria. We examined six merozoite ligands EBA175, EBL1, EBA140, MSP1, Rh4 and Rh5, and their corresponding erythrocyte receptors, glycophorin (Gyp) A, GypB, GypC, band 3, complement receptor (CR) 1 and basigin, focusing on the regions involved in the ligand-receptor interactions. Positive Tajima's D values (>1) were observed only in the MSP1 C-terminal region and EBA175 region II, while negative values (<-1) were observed in EBL-1 region II, Rh4, basigin exons 3 and 5, CR1 exon 5, Gyp B exons 2, 3 and 4 and Gyp C exon 2. Additionally, ebl-1 region II and basigin exon 3 showed extreme negative values in all three tests, Tajima's D, Fu & Li D* and F*, ≤ - 2. A large majority of the erythrocyte receptor and merozoite genes have a negative Tajima's D even when compared with previously published whole genome data. Thus, highlighting EBA175 region II and MSP1-33, as outlier genes with a positive Tajima's D (>1). Both these genes contain multiple polymorphisms, which in the case of EBA175 may counteract receptor polymorphisms and/or evade host immune responses and in MSP1 the polymorphisms may primarily evade host immune responses.

The Gerbich blood group system: old knowledge, new importance

Antigens of the Gerbich blood group system are expressed on glycophorin C (GPC) and glycophorin D (GPD), minor sialoglycoproteins of human erythrocytes. GPC and GPD help maintain erythrocyte shape of and contributes to the stability of its membrane. There are six high-prevalence Gerbich antigens: Ge2, Ge3, Ge4, GEPL (GE10), GEAT (GE11), GETI (GE12) and five low-prevalence Gerbich antigens: Wb (GE5), Ls^a (GE6), An^a (GE7), Dh^a (GE8), GEIS (GE9). Some Gerbich antigens (Ge4, Wb, Dh^a, GEAT) are expressed only on GPC, two (Ge2, An^a) are expressed only on GPD, while others (Ge3, Ls^a, GEIS, GEPL, GETI) are expressed on both GPC and GPD. Antibodies recognizing GPC/GPD may arise naturally (so-called "naturally-occurring RBC antibodies") or as the result of alloimmunization, and some of them may be clinically relevant. Gerbich antibodies usually do not cause serious hemolytic transfusion reactions (HTR); autoantibodies of anti-Ge2- or anti-Ge3 specificity can cause autoimmune hemolytic anemia (AIHA).

Plasmodium reichenowi EBA-140 merozoite ligand binds to glycophorin D on chimpanzee red blood cells, shedding new light on origins of Plasmodium falciparum

Background

All symptoms of malaria are caused by the intraerythrocytic proliferation of Plasmodium merozoites. Merozoites invade erythrocytes using multiple binding ligands that recognise specific surface receptors. It has been suggested that adaptation of Plasmodium parasites to infect specific hosts is driven by changes in genes encoding Plasmodium erythrocyte-binding ligands (EBL) and reticulocyte-binding ligands (RBL). Homologs of both EBL and RBL, including the EBA-140 merozoite ligand, have been identified in P. falciparum and P. reichenowi, which infect humans and chimpanzees, respectively. The P. falciparum EBA-140 was shown to bind human glycophorin C, a minor erythrocyte sialoglycoprotein. Until now, the erythrocyte receptor for the P. reichenowi EBA-140 remained unknown.

Methods

The baculovirus expression vector system was used to obtain the recombinant EBA-140 Region II, and flow cytometry and immunoblotting methods were applied to characterise its specificity.

Results

We showed that the chimpanzee glycophorin D is the receptor for the P. reichenowi EBA-140 ligand on chimpanzee red blood cells.

Conclusions

We propose that the development of glycophorin C specificity is spurred by the P. falciparum lineage. We speculate that the P. falciparum EBA-140 evolved to hijack GPC on human erythrocytes during divergence from its ape ancestor.

The Recent De Novo Origin of Protein C-Termini

Protein-coding sequences can arise either from duplication and divergence of existing sequences, or de novo from noncoding DNA. Unfortunately, recently evolved de novo genes can be hard to distinguish from false positives, making their study difficult. Here, we study a more tractable version of the process of conversion of noncoding sequence into coding: the co-option of short segments of noncoding sequence into the C-termini of existing proteins via the loss of a stop codon. Because we study recent additions to potentially old genes, we are able to apply a variety of stringent quality filters to our annotations of what is a true protein-coding gene, discarding the putative proteins of unknown function that are typical of recent fully de novo genes. We identify 54 examples of C-terminal extensions in Saccharomyces and 28 in Drosophila, all of them recent enough to still be polymorphic. We find one putative gene fusion that turns out, on close inspection, to be the product of replicated assembly errors, further highlighting the issue of false positives in the study of rare events. Four of the Saccharomyces C-terminal extensions (to ADH1, ARP8, TPM2, and PIS1) that survived our quality filters are predicted to lead to significant modification of a protein domain structure.

Pathogen-driven selection in the human genome

Infectious diseases and epidemics have always accompanied and characterized human history, representing one of the main causes of death. Even today, despite progress in sanitation and medical research, infections are estimated to account for about 15% of deaths. The hypothesis whereby infectious diseases have been acting as a powerful selective pressure was formulated long ago, but it was not until the availability of large-scale genetic data and the development of novel methods to study molecular evolution that we could assess how pervasively infectious agents have shaped human genetic diversity. Indeed, recent evidences indicated that among the diverse environmental factors that acted as selective pressures during the evolution of our species, pathogen load had the strongest influence. Beside the textbook example of the major histocompatibility complex, selection signatures left by pathogen-exerted pressure can be identified at several human loci, including genes not directly involved in immune response. In the future, high-throughput technologies and the availability of genetic data from different populations are likely to provide novel insights into the evolutionary relationships between the human host and its pathogens. Hopefully, this will help identify the genetic determinants modulating the susceptibility to infectious diseases and will translate into new treatment strategies.

Putatively noncoding transcripts show extensive association with ribosomes

There have been recent surprising reports that whole genes can evolve de novo from noncoding sequences. This would be extraordinary if the noncoding sequences were random with respect to amino acid identity. However, if the noncoding sequences were previously translated at low rates, with the most strongly deleterious cryptic polypeptides purged by selection, then de novo gene origination would be more plausible. Here we analyze Saccharomyces cerevisiae data on noncoding transcripts found in association with ribosomes. We find many such transcripts. Although their average ribosomal densities are lower than those of protein-coding genes, a significant proportion of noncoding transcripts nevertheless have ribosomal densities comparable to those of coding genes. Most show increased ribosomal association in response to starvation, as has been previously reported for other noncoding sequences such as untranslated regions and introns. In rich media, ribosomal association is correlated with start codons but is not usually consistent and contiguous beyond that, suggesting that translation occurs only at low rates. One transcript contains a 28-codon open reading frame, which we name RDT1, which shows evidence of translation, and may be a new protein-coding gene that originated de novo from noncoding sequence. But the bulk of the ribosomal association cannot be attributed to unannotated protein-coding genes. Our primary finding of extensive ribosome association shows that a necessary precondition for selective purging is met, making de novo gene evolution more plausible. Our analysis is also proof of principle of the utility of ribosomal profiling data for the purpose of gene annotation.

Genome-wide association study SNPs in the human genome diversity project populations: does selection affect unlinked SNPs with shared trait associations?

Genome-wide association studies (GWAS) have identified more than 2,000 trait-SNP associations, and the number continues to increase. GWAS have focused on traits with potential consequences for human fitness, including many immunological, metabolic, cardiovascular, and behavioral phenotypes. Given the polygenic nature of complex traits, selection may exert its influence on them by altering allele frequencies at many associated loci, a possibility which has yet to be explored empirically. Here we use 38 different measures of allele frequency variation and 8 iHS scores to characterize over 1,300 GWAS SNPs in 53 globally distributed human populations. We apply these same techniques to evaluate SNPs grouped by trait association. We find that groups of SNPs associated with pigmentation, blood pressure, infectious disease, and autoimmune disease traits exhibit unusual allele frequency patterns and elevated iHS scores in certain geographical locations. We also find that GWAS SNPs have generally elevated scores for measures of allele frequency variation and for iHS in Eurasia and East Asia. Overall, we believe that our results provide evidence for selection on several complex traits that has caused changes in allele frequencies and/or elevated iHS scores at a number of associated loci. Since GWAS SNPs collectively exhibit elevated allele frequency measures and iHS scores, selection on complex traits may be quite widespread. Our findings are most consistent with this selection being either positive or negative, although the relative contributions of the two are difficult to discern. Our results also suggest that trait-SNP associations identified in Eurasian samples may not be present in Africa, Oceania, and the Americas, possibly due to differences in linkage disequilibrium patterns. This observation suggests that non-Eurasian and non-East Asian sample populations should be included in future GWAS.

Evolution of prokaryotic genes by shift of stop codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene's reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene's end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3'UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3'UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.

Alternative translation start sites are conserved in eukaryotic genomes

Alternative start AUG codons within a single transcript can contribute to diversity of the proteome; however, their functional significance remains controversial. Here, we provide comparative genomics evidence that alternative start codons are under negative selection in vertebrates, insects and yeast. In genes where the annotated start codon (sAUG) resides within the suboptimal nucleotide context, the downstream in-frame AUG codons (dAUG) among the first ∼30 codon sites are significantly more conserved between species than in genes where the sAUG resides within the optimal context. Proteomics data show that this difference is not an annotation artifact and that dAUGs are in fact under selection as alternative start sites. The key optimal, and sometimes suboptimal, context-determining nucleotides of both the sAUG and dAUGs are conserved. Selection for secondary start sites is stronger in genes with the weak primary start site. Genes with multiple conserved start sites are enriched for transcription factors, and tend to have longer 5'UTRs and higher degree of alternative splicing. Together, these results imply that the use of alternative start sites by means of leaky mRNA scanning is a functional mechanism under selection for increased efficiency of translation and/or for translation of different N-terminal protein variants.

From evolutionary genetics to human immunology: how selection shapes host defence genes

Pathogens have always been a major cause of human mortality, so they impose strong selective pressure on the human genome. Data from population genetic studies, including genome-wide scans for selection, are providing important insights into how natural selection has shaped immunity and host defence genes in specific human populations and in the human species as a whole. These findings are helping to delineate genes that are important for host defence and to increase our understanding of how past selection has had an impact on disease susceptibility in modern populations. A tighter integration between population genetic studies and immunological phenotype studies is now necessary to reveal the mechanisms that have been crucial for our past and present survival against infection.