Roles of canonical Wnt/TCF-1 signaling for hematopoiesis, lymphocyte development and function

Summary : The canonical Wnt signaling pathway plays key roles in the maintenance of self-renewing tissues, like the gut or the skin. In contrast, the role of this pathway in hematopoiesis remains poorly defined. Wnt ligands transmit signals through ß-catenin which activates gene transcription upon its association with Lymphoid Cell Enhancer/T Cell Factor (LEF/TCF). Currently, v-catenin is the only alternative factor known to transduce canonical Wnt signals. The ß-/γ-catenin bindiná domain in TCF-1 is required to partly rescue thymopoiesis and NK cell development in TCF-1-deficient mice. However, T cell development and hematopoiesis w-as normal in mice deficient of ß-catenin, or of γ-catenin. Surprisingly we found that hematopoiesis and thymopoiesis was also normal in the combined absence of ß- and γ-catenin. Reporter assays showed that double-deficient lymphocytes were still able to transduce canonical wnt signals. These data provided evidence that hematopoietic cells can transduce canonical Wnt signals in the combined absence of ß- and γ-catenin. There exist numerous TCF-1 isoforrns including those that harbor the N-terminal ß-/y-catenin binding domain or that contains a C-terminal CRARF domain whose role in vivo has not been previously tested. We found that the CRARF domain influences lymphocyte development in conjunction with the N-treminal ß-/γ-catenin binding. The presence of the two domains directs thymocytes to the CD8+ T cell lineage whereas NK cell development is abolished. Roles of the canonical Wnt/TCF-1 pathway for lymphocyte function have not been defined. We demonstrate that TCF-1 deficient CDBT T cells mount a normal primary response to viral infection but these T cells fail to expand upon restimulation. The failure of CD8+ T cells to respond to IL-2 during primary infection seems to account for this phenotype. Thus, TCF-1 is essential for programming functional CD8+ T cell memory. Collectively, these data provide significant new insights into the role of Wnt/TCF-1 pathway for lymphocyte development and function and suggest a novel mechanism of Wnt signal transuction in hematopoietic cells. Résumé : La voie de signalisation canonique Wnt joue un rôle prépondérant dans le renouvellement de tissus, comme l'intestin ou la peau. Son rôle dans l'hématopoïèse est quant à lui mal défini. Le ligand Wnt transmet le signal via la ß-catenin qui active la transcription de gènes cibles quand il est associé avec Lymphoid Cell Enhancer,~T Cell Factor (LEF/TCF). Actuellement, la γ-catenin est le seul autre facteur connu pouvant se substituer à la fonction de la ß-catenin. Un variant de TCF-1 contenant le domaine liant ß-/,~-catenin est capable de restaurer le développement des lymphocytes T et NK en l'absence de TCF-1. Cependant la thymopoïèse et l'hématopoïèse sont normales dans les souris déficientes pour la ß-catenin ou la γ-catenin. De façon surprenante, nous avons trouvé que l'hématopoïèse et le développement des lymphocytes sont normaux lors de l'absence combinée de ß-/γ-catenin. De plus, la transduction des signaux de la voie de signalisation Wnt est maintenue dans des lymphocytes déficients pour ß-/γ-catenin. Ces résultats démontrent que les cellules hématopoïétiques peuvent transmettre les signaux de la voie canonique Wnt lors de l'absence combinée de la ß et la γ -catenin. Il existe de nombreuses isofonnes de TCF-1, y compris certaines qui comprennent un domaine qui lie ß-/γ-catenin du côté N-terminus ou qui contiennent un domaine CRARF du côté C-terminus. Nous montrons ici que le domaine CRARF influence le développement des lymphocytes en conjonction avec le domaine liant ß-/γ-catenin. La présence des deux domaines dirige les thymocytes vers la lignée de cellules T CD8, alors que le développement des cellules NK est aboli. Au-delà de sa fonction sur le développement des lymphocytes, le rôle de la soie de signalisation canonique Wnt/TCF-1 lors d'une infection n'a pas été défini. Nous avons montré que les cellules T CD8, déficientes pour TCF-1, développent une réponse primaire normale à une infection virale, mais qu'elles ne s'accumulent pas après restimulation. L'incapacité des cellules TCD8 à répondre à l'IL-2 durant la réponse primaire peut expliquer ce phénotype. Ainsi; TCF-1 est essentiel pour la programmation de cellules T CD8 mémoires fonctionnelles. L'ensemble de ces résultats fournit de nouveaux aperçus du rôle de la voie de signalisation Wnt/TCF-1 pour le développement et la fonction des lymphocytes et suggèrent un nouveau mécanisme de transduction du signal Wnt dans les cellules hématopoïétiques.

Mechanisms of reproductive isolation between an ant species of hybrid origin and one of its parents.

The establishment of new species by hybridization is difficult because it requires the development of reproductive isolation (RI) in sympatry to escape the homogenizing effects of gene flow from the parental species. Here we investigated the role of two pre- and two postzygotic mechanisms of RI in a system comprising two interdependent Pogonomyrmex harvester ant lineages (the H1 and H2 lineages) of hybrid origin and one of their parental species (P. rugosus). Similar to most other ants, P. rugosus is characterized by an environmental system of caste determination with female brood developing either into queens or workers depending on nongenetic factors. By contrast, there is a strong genetic component to caste determination in the H1 and H2 lineages because the developmental fate of female brood depends on the genetic origin of the parents, with interlineage eggs developing into workers and intralineage eggs developing into queens. The study of a mixed mating aggregation revealed strong differences in mating flight timing between P. rugosus and the two lineages as a first mechanism of RI. A second important prezygotic mechanism was assortative mating. Laboratory experiments also provided support for one of the two investigated mechanisms of postzygotic isolation. The majority of offspring produced from the few matings between P. rugosus and the lineages aborted at the egg stage. This hybrid inviability was under maternal influence, with hybrids produced by P. rugosus queens being always inviable whereas a small proportion of H2 lineage queens produced large numbers of adult hybrid offspring. Finally, we found no evidence that genetic caste determination acted as a second postzygotic mechanism reducing gene flow between P. rugosus and the H lineages. The few viable P. rugosus-H hybrids were not preferentially shunted into functionally sterile workers but developed into both workers and queens. Overall, these results reveal that the nearly complete (99.5%) RI between P. rugosus and the two hybrid lineages stems from the combination of two typical prezygotic mechanisms (mating time divergence and assortative mating) and one postzygotic mechanism (hybrid inviability).

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures.

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions-regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots-are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

RNA-based gene duplication sheds new light on mammalian sex chromosome evolution

ABSTRACT : Gene duplication is a fundamental source of raw material for the origin of genetic novelty. It has been assumed for a long time that DNA-based gene duplication was the only source of new genes. Recently however, RNA-based gene duplication (retroposition) was shown in multiple organisms to contribute significantly to their genetic diversity. This mechanism produces intronless gene copies (retrocopies) that are inserted in random genomic position, independent of the position of the parental source genes. In human, mouse and fruit fly, it was demonstrated that the X-linked genes spawned an excess of functional retroposed gene copies (retrogenes). In human and mouse, the X chromosome also recruited an excess of retrogenes. Here we further characterized these interesting biases related to the X chromosome in mammals. Firstly, we have confirmed presence of the aforementioned biases in dog and opossum genome. Then based on the expression profile of retrogenes during various spermatogenetic stages, we have provided solid evidence that meiotic sex chromosome inactivation (MSCI) is responsible for an excess of retrogenes stemming from the X chromosome. Moreover, we showed that the X-linked genes started to export an excess of retrogenes just after the split of eutherian and marsupial mammalian lineages. This suggests that MSCI has originated around this time as well. More fundamentally, as MSCI reflects the spread of recombination barrier between the X and Y chromosomes during their evolution, our observation allowed us to re-estimate the age of mammalian sex chromosomes. Previous estimates suggested that they emerged in the common ancestor of all mammals (before the split of monotreme lineage); whereas, here we showed that they originated around the split of marsupial and eutherian lineages, after the divergence of monotremes. Thus, the therian (marsupial and eutherian) sex chromosomes are younger than previously thought. Thereafter, we have characterized the bias related to the recruitment of genes to the X chromosome. Sexually antagonistic forces are most likely driving this pattern. Using our limited retrogenes expression data, it is difficult to determine the exact nature of these forces but some conclusions have been made. Lastly, we looked at the history of this biased recruitment: it commenced around the split of marsupial and eutherian lineages (akin to the biased export of genes out of the X). In fact, the sexually antagonistic forces are predicted to appear just around that time as well. Thereby, the history of the recruitment of genes to the X, provides an indirect evidence that these forces are responsible for this bias.

Genetic diversity and genetic exchange in Trypanosoma cruzi: dual drug-resistant "progeny" from episomal transformants

Extensive characterisation of Trypanosoma cruzi by isoenzyme phenotypes has separated the species into three principal zymodeme groups, Z1, Z2 and Z3, and into many individual zymodemes. There is marked diversity within Z2. A strong correlation has been demonstrated between the strain clusters determined by isoenzymes and those obtained using random amplified polymorphic DNA (RAPD) profiles. Polymorphisms in ribosomal RNA genes, in mini-exon genes, and microsatellite fingerprinting indicate the presence of at least two principal T. cruzi genetic lineages. Lineage 1 appears to correspond with Z2 and lineage 2 with Z1. Z1 (lineage 2) is associated with Didelphis. Z2 (lineage 1) may be associated with a primate host. Departures from Hardy-Weinberg equilibrium and linkage disequilibrium indicate that propagation of T. cruzi is predominantly clonal. Nevertheless, two studies show putative homozygotes and heterozygotes circulating sympatrically: the allozyme frequencies for phosphoglucomutase, and hybrid RAPD profiles suggest that genetic exchange may be a current phenomenon in some T. cruzi transmission cycles. We were able to isolate dual drug-resistant T. cruzi biological clones following copassage of putative parents carrying single episomal drug-resistant markers. A multiplex PCR confirmed that dual drug-resistant clones carried both episomal plasmids. Preliminary karyotype analysis suggests that recombination may not be confined to the extranuclear genome.

Protein-binding site at the immunoglobulin mu membrane polyadenylylation signal: possible role in transcription termination.

mRNAs specifying immunoglobulin mu and delta heavy chains are encoded by a single large, complex transcription unit (mu + delta gene). The transcriptional activity of delta gene segments in terminally differentiated, IgM-secreting B lymphocytes is 10-20 times lower than in earlier B-lineage cells expressing delta mRNA. We find that transcription of the mu + delta gene in IgM-secreting murine myeloma cells terminates within a region of 500-1000 nucleotides immediately following the mu membrane (mu m) polyadenylylation site. Transcription decreases only minimally through this region in murine cell lines representative of earlier stages in B-cell development. A DNA fragment containing the mu m polyadenylylation signal gives protein-DNA complexes with different mobilities in gel retardation assays with nuclear extracts from myeloma cells than with nuclear extracts from earlier B-lineage cells. However, using a recently developed "footprinting" procedure in which protein-DNA complexes resolved in gel retardation assays are subjected to nucleolytic cleavage while still in the polyacrylamide gel, we find that the DNA sequences protected by factors from the two cell types are indistinguishable. The factor-binding site on the DNA is located 5' of the mu m polyadenylylation signal AATAAA and includes the 15-nucleotide-long A + T-rich palindrome CTGTAAACAAATGTC. This type of palindromic binding site exhibits orientation-dependent activity consistent with the reported properties of polymerase II termination signals. This binding site is followed by two sets of directly repeated DNA sequences with different helical conformation as revealed by their reactivity with the chemical nuclease 1,10-phenanthroline-copper. The close proximity of these features to the signals for mu m mRNA processing may reflect a linkage of the processes of developmentally regulated mu m polyadenylylation and transcription termination.

E2F1 mediates DNA damage and apoptosis through HCF-1 and the MLL family of histone methyltransferases.

E2F1 is a key positive regulator of human cell proliferation and its activity is altered in essentially all human cancers. Deregulation of E2F1 leads to oncogenic DNA damage and anti-oncogenic apoptosis. The molecular mechanisms by which E2F1 mediates these two processes are poorly understood but are important for understanding cancer progression. During the G1-to-S phase transition, E2F1 associates through a short DHQY sequence with the cell-cycle regulator HCF-1 together with the mixed-lineage leukaemia (MLL) family of histone H3 lysine 4 (H3K4) methyltransferases. We show here that the DHQY HCF-1-binding sequence permits E2F1 to stimulate both DNA damage and apoptosis, and that HCF-1 and the MLL family of H3K4 methyltransferases have important functions in these processes. Thus, HCF-1 has a broader role in E2F1 function than appreciated earlier. Indeed, sequence changes in the E2F1 HCF-1-binding site can modulate both up and down the ability of E2F1 to induce apoptosis indicating that HCF-1 association with E2F1 is a regulator of E2F1-induced apoptosis.

Non-nest mate discrimination and clonal colony structure in the parthenogenetic ant Cerapachys biroi

Understanding the interplay between cooperation and conflict in social groups is a major goal of biology. One important factor is genetic relatedness, and animal societies are usually composed of related but genetically different individuals, setting the stage for conflicts over reproductive allocation. Recently, however, it has been found that several ant species reproduce predominantly asexually. Although this can potentially give rise to clonal societies, in the few well-studied cases, colonies are often chimeric assemblies of different genotypes, due to worker drifting or colony fusion. In the ant Cerapachys biroi, queens are absent and all individuals reproduce via thelytokous parthenogenesis, making this species an ideal study system of asexual reproduction and its consequences for social dynamics. Here, we show that colonies in our study population on Okinawa, Japan, recognize and effectively discriminate against foreign workers, especially those from unrelated asexual lineages. In accord with this finding, colonies never contained more than a single asexual lineage and average pairwise genetic relatedness within colonies was extremely high (r = 0.99). This implies that the scope for social conflict in C. biroi is limited, with unusually high potential for cooperation and altruism.

Isoenzyme profiles of Trypanosoma cruzi stocks from different areas of Paraguay

Twenty one Trypanosoma cruzi stocks from humans, domiciliary triatomines and one sylvatic animal of different areas of Paraguay were subjected to isoenzyme analysis. Thirteen enzyme systems (15 loci in total) were studied. MN cl2 (clonets 39) and SO34 cl4 (clonets 20) were used as references. Relationships between stocks were depicted by an UPGMA dendrogram constructed using the Jaccard´s distances matrix. Among the Paraguayan stocks 14 zymodemes were identified (Par1 to Par14), Par 5 being the most frequent. Polymorphism rate and clonal diversity were 0.73 and 0.93, respectively. Average number of alleles per polymorphic locus was 2.5 (range 2-4). These measurements show a high diversity, which is confirmed by the dendrogram topology. All stocks belong to the same lineage, as MN cl2 reference strain (T. cruzi II). Moreover three distinct subgroups were identified and two of them correspond to Brazilian and Bolivian zymodemes, respectively. The third subgroup, the most common in Paraguay, is related to Tulahuen stock. The large geographical distribution of some zymodemes agrees with the hypothesis of clonality for T. cruzi populations. However sample size was not adequate to detect genetic recombination in any single locality.

Deciphering the unusual HLA-A2/Melan-A/MART-1-specific TCR repertoire in humans.

The Melan-A/MART-1(26-35) antigenic peptide is one of the best studied human tumor-associated antigens. It is expressed in healthy melanocytes and malignant melanoma and is recognized by CD8(+) T cells in the context of the MHC class I molecule HLA-A*0201. While an unusually large repertoire of CD8(+) T cells specific for this antigen has been documented, the reasons for its generation have remained elusive. In this issue of the European Journal of Immunology, Pinto et al. [Eur. J. Immunol. 2014. 44: 2811-2821] uncover one important mechanism by comparing the thymic expression of the Melan-A gene to that in the melanocyte lineage. This study shows that medullary thymic epithelial cells (mTECs) dominantly express a truncated Melan-A transcript, the product of misinitiation of transcription. Consequently, the protein product in mTECs lacks the immunodominant epitope spanning residues 26-35, thus precluding central tolerance to this antigen. In contrast, melanocytes and melanoma tumor cells express almost exclusively the full-length Melan-A transcript, thus providing the target antigen for efficient recognition by HLA-A2-restricted CD8(+) T cells. The frequency of these alternative gene transcription modes may be more common than previously appreciated and may represent an important factor modulating the efficiency of central tolerance induction in the thymus.

The early use of yellow fever virus strain 17D for vaccine production in Brazil - a review

The use of yellow fever (YF) virus 17D strain for vaccine production adapted in Brazil since its introduction in 1937 was reviewed. This was possible due to the availability of official records of vaccine production. The retrieved data highlight the simultaneous use of several serially passaged 17D substrain viruses for both inocula and vaccine preparation that allowed uninterrupted production. Substitution of these substrain viruses became possible with the experience gained during quality control and human vaccination. Post-vaccinal complications in humans and the failure of some viruses in quality control tests (neurovirulence for monkeys) indicated that variables needed to be reduced during vaccine production, leading to the development of the seed lot system. The 17DD substrain, still used today, was the most frequently used substrain and the most reliable in terms of safety and efficacy. For this reason, it is possible to derive an infectious cDNA clone of this substrain combined with production in cell culture that could be used to direct the expression of heterologous antigens and lead to the development of new live vaccines.

Deficient T cell fate specification in mice with an induced inactivation of Notch1.

Notch proteins are cell surface receptors that mediate developmental cell specification events. To explore the function of murine Notch1, an essential portion of the gene was flanked with loxP sites and inactivation induced via interferon-regulated Cre recombinase. Mice with a neonatally induced loss of Notch1 function were transiently growth retarded and had a severe deficiency in thymocyte development. Competitive repopulation of lethally irradiated wild-type hosts with wild-type- and Notch1-deficient bone marrow revealed a cell autonomous blockage in T cell development at an early stage, before expression of T cell lineage markers. Notch1-deficient bone marrow did, however, contribute normally to all other hematopoietic lineages. These findings suggest that Notch1 plays an obligatory and selective role in T cell lineage induction.

E2F activation of S phase promoters via association with HCF-1 and the MLL family of histone H3K4 methyltransferases.

E2F transcriptional regulators control human-cell proliferation by repressing and activating the transcription of genes required for cell-cycle progression, particularly the S phase. E2F proteins repress transcription in association with retinoblastoma pocket proteins, but less is known about how they activate transcription. Here, we show that the human G1 phase regulator HCF-1 associates with both activator (E2F1 and E2F3a) and repressor (E2F4) E2F proteins, properties that are conserved in insect cells. Human HCF-1-E2F interactions are versatile: their associations and binding to E2F-responsive promoters are cell-cycle selective, and HCF-1 displays coactivator properties when bound to the E2F1 activator and corepressor properties when bound to the E2F4 repressor. During the G1-to-S phase transition, HCF-1 recruits the mixed-lineage leukemia (MLL) and Set-1 histone H3 lysine 4 methyltransferases to E2F-responsive promoters and induces histone methylation and transcriptional activation. These results suggest that HCF-1 induces cell-cycle-specific transcriptional activation by E2F proteins to promote cell proliferation.

New insights into the long-term evolutionary history of the avian MHC

SummaryGene duplication and neofunctidnalization are important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family with a central role in vertebrates' adaptive immune system. Multigene families, which evolved in large part through duplication, represent promising systems to study the still strongly depbated relative roles of neutral and adaptive processes in the evolution of phenotypic complexity. Detailed knowledge on ecological function and a well-characterized evolutionary history place the mammals' MHC amongst ideal study systems. However mammalian MHCs usually encompass several million base pairs and hold a large number of functional and non-functional duplicate genes, which makes their study complex. Avian MHCs on the other hand are usually way more compact, but the reconstruction of. their evolutionary history has proven notoriously difficult. However, no focused attempt has been undertaken so far to study the avian MHC evolutionary history in a broad phylogenetic context and using adequate gene regions.In the present PhD, we were able to make important contributions to the understanding of the long-term evolution of the avian MHC class II Β (MHCI1B). First, we isolated and characterized MHCIIB genes in barn owl (Tyto alba?, Strigiformes, Tytonidae), a species from an avian lineage in which MHC has not been studied so far. Our results revealed that with only two functional MHCIIB genes the MHC organization of barn owl may be similar to the 'minimal essential' MHC of chicken (Gallus gallus), indicating that simple MHC organization may be ancestral to birds. Taking advantage of the sequence information from barn owl, we studied the evolution of MHCIIB genes in 13 additional species of 'typical' owls (Strigiformes, Strigidae). Phylogenetic analyses revealed that according to their function, in owls the peptide-binding region (PBR) encoding exon 2 and the non-PBR encoding exon 3 evolve by different patterns. Exon 2 exhibited an evolutionary history of positive selection and recombination, while exon 3 traced duplication history and revealed two paralogs evolving divergently from each other in owls, and in a shorebird, the great snipe {Gallinago media). The results from exon 3 were the first ever from birds to demonstrate gene orthology in species that diverged tens of millions of years ago, and strongly questioned whether the taxa studied before provided an adequate picture of avian MHC evolution. In a follow-up study, we aimed at explaining a striking pattern revealed by phylogenetic trees analyzing the owl sequences along with MHCIIB sequences from other birds: One owl paralog (termed DAB1) grouped with sequences of passerines and falcons, while the other (DAB2) grouped with wildfowl, penguins and birds of prey. This could be explained by either a duplication event preceding the evolution of these bird orders, or by convergent evolution of similar sequences in a number of orders. With extensive phylogenetic analyses we were able to show, that indeed a duplication event preceeded the major avian radiation -100 my ago, and that following this duplication, the paralogs evolved under positive selection. Furthermore, we showed that the divergently evolving amino acid residues in the MHCIIB-encoded β-chain potentially interact with the MHCI I α-chain, and that molecular coevolution of the interacting residues may have been involved in the divergent evolution of the MHCIIB paralogs.The findings of this PhD are of particular interest to the understanding of the evolutionary history of the avian MHC and, by providing essential information on long-term gene history in the avian MHC, open promising perspectives for advances in the understanding of the evolution of multigene families in general, and for avian MHC organization in particular. Amongst others I discuss the importance of including protein structure in the phylogenetic study of multigene families, and the roles of ecological versus molecular selection pressures. I conclude by providing a population genomic perspective on avian MHC, which may serve as a basis for future research to investigate the relative roles of neutral processes involving effective population size effects and of adaptation in the evolution of avian MHC diversity and organization.RésuméLa duplication de gènes et leur néo-fonctionnalisation sont des processus importants dans l'évolution de la complexité phénotypique. Ils sont impliqués dans l'apparition d'importantes nouveautés évolutives favorisant l'adaptation écologique, comme c'est le cas pour le complexe majeur d'histocompatibilité

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.