A phase 1 trial of human telomerase reverse transcriptase (hTERT) vaccination combined with therapeutic strategies to control immune-suppressor mechanisms.

The presence of inhibitory immune cells and difficulty in generating activated effector T cells remain obstacles to development of effective cancer vaccines. We designed a vaccine regimen combining human telomerase reverse transcriptase (hTERT) peptides with concomitant therapies targeting regulatory T cells (Tregs) and cyclooxygenase-2 (COX2)-mediated immunosuppression. This Phase 1 trial combined an hTERT-derived 7-peptide library, selected to ensure presentation by both HLA class-I and class-II in 90% of patients, with oral low-dose cyclophosphamide (to modulate Tregs) and the COX2 inhibitor celecoxib. Adjuvants were Montanide and topical TLR-7 agonist, to optimise antigen presentation. The primary objective was determination of the safety and tolerability of this combination therapy, with anti-cancer activity, immune response and detection of antigen-specific T cells as additional endpoints. Twenty-nine patients with advanced solid tumours were treated. All were multiply-pretreated, and the majority had either colorectal or prostate cancer. The most common adverse events were injection-site reactions, fatigue and nausea. Median progression-free survival was 9 weeks, with no complete or partial responses, but 24% remained progression-free for ≥6 months. Immunophenotyping showed post-vaccination expansion of CD4+ and CD8+ T cells with effector phenotypes. The in vitro re-challenge of T cells with hTERT peptides, TCR sequencing, and TCR similarity index analysis demonstrated the expansion following vaccination of oligoclonal T cells with specificity for hTERT. However, a population of exhausted PD-1+ cytotoxic T cells was also expanded in vaccinated patients. This vaccine combination regimen was safe and associated with antigen-specific immunological responses. Clinical activity could be improved in future by combination with anti-PD1 checkpoint inhibition to address the emergence of an exhausted T cell population.

Evolutionary remodelling of N-terminal domain loops fine-tunes SARS-CoV-2 spike.

The emergence of SARS-CoV-2 variants has exacerbated the COVID-19 global health crisis. Thus far, all variants carry mutations in the spike glycoprotein, which is a critical determinant of viral transmission being responsible for attachment, receptor engagement and membrane fusion, and an important target of immunity. Variants frequently bear truncations of flexible loops in the N-terminal domain (NTD) of spike; the functional importance of these modifications has remained poorly characterised. We demonstrate that NTD deletions are important for efficient entry by the Alpha and Omicron variants and that this correlates with spike stability. Phylogenetic analysis reveals extensive NTD loop length polymorphisms across the sarbecoviruses, setting an evolutionary precedent for loop remodelling. Guided by these analyses, we demonstrate that variations in NTD loop length, alone, are sufficient to modulate virus entry. We propose that variations in NTD loop length act to fine-tune spike; this may provide a mechanism for SARS-CoV-2 to navigate a complex selection landscape encompassing optimisation of essential functionality, immune-driven antigenic variation and ongoing adaptation to a new host.

An entropic safety catch controls hepatitis C virus entry and antibody resistance.

E1 and E2 (E1E2), the fusion proteins of Hepatitis C Virus (HCV), are unlike that of any other virus yet described, and the detailed molecular mechanisms of HCV entry/fusion remain unknown. Hypervariable region-1 (HVR-1) of E2 is a putative intrinsically disordered protein tail. Here, we demonstrate that HVR-1 has an autoinhibitory function that suppresses the activity of E1E2 on free virions; this is dependent on its conformational entropy. Thus, HVR-1 is akin to a safety catch that prevents premature triggering of E1E2 activity. Crucially, this mechanism is turned off by host receptor interactions at the cell surface to allow entry. Mutations that reduce conformational entropy in HVR-1, or genetic deletion of HVR-1, turn off the safety catch to generate hyper-reactive HCV that exhibits enhanced virus entry but is thermally unstable and acutely sensitive to neutralising antibodies. Therefore, the HVR-1 safety catch controls the efficiency of virus entry and maintains resistance to neutralising antibodies. This discovery provides an explanation for the ability of HCV to persist in the face of continual immune assault and represents a novel regulatory mechanism that is likely to be found in other viral fusion machinery.

Mechanistic diversity in MHC class I antigen recognition.

Throughout its evolution, the human immune system has developed a plethora of strategies to diversify the antigenic peptide sequences that can be targeted by the CD8+ T cell response against pathogens and aberrations of self. Here we provide a general overview of the mechanisms that lead to the diversity of antigens presented by MHC class I complexes and their recognition by CD8+ T cells, together with a more detailed analysis of recent progress in two important areas that are highly controversial: the prevalence and immunological relevance of unconventional antigen peptides; and cross-recognition of antigenic peptides by the T cell receptors of CD8+ T cells.

Characterization of the MurT/GatD complex in Mycobacterium tuberculosis towards validating a novel anti-tubercular drug target.

OBJECTIVES: Identification and validation of novel therapeutic targets is imperative to tackle the rise of drug resistance in tuberculosis. An essential Mur ligase-like gene (Rv3712), expected to be involved in cell-wall peptidoglycan (PG) biogenesis and conserved across mycobacteria, including the genetically depleted Mycobacterium leprae, was the primary focus of this study. METHODS: Biochemical analysis of Rv3712 was performed using inorganic phosphate release assays. The operon structure was identified using reverse-transcriptase PCR and a transcription/translation fusion vector. In vivo mycobacterial protein fragment complementation assays helped generate the interactome. RESULTS: Rv3712 was found to be an ATPase. Characterization of its operon revealed a mycobacteria-specific promoter driving the co-transcription of Rv3712 and Rv3713. The two gene products were found to interact with each other in vivo. Sequence-based functional assignments reveal that Rv3712 and Rv3713 are likely to be the mycobacterial PG precursor-modifying enzymes MurT and GatD, respectively. An in vivo network involving Mtb-MurT, regulatory proteins and cell division proteins was also identified. CONCLUSIONS: Understanding the role of the enzyme complex in the context of PG metabolism and cell division, and the implications for antimicrobial resistance and host immune responses will facilitate the design of therapeutics that are targeted specifically to M. tuberculosis.

Characterization of the UK anthrax vaccine and human immunogenicity.

The manufacture of the UK Anthrax vaccine (AVP) focuses on the production of Protective Antigen (PA) from the Bacillus anthracis Sterne strain. Although used for decades, several of AVP's fundamental properties are poorly understood, including its exact composition, the extent to which proteins other than PA may contribute to protection, and whether the degree of protection varies between individuals.This study involved three innovative investigations. Firstly, the composition of AVP was analyzed using liquid chromatography tandem mass-spectrometry (LC-MS/MS), requiring the development of a novel desorption method for releasing B. anthracis proteins from the vaccine's aluminum-containing adjuvant. Secondly, computational MHC-binding predictions using NetMHCIIpan were made for the eight most abundant proteins of AVP, for the commonest HLA alleles in multiple ethnic groups, and for multiple B. anthracis strains. Thirdly, antibody levels and toxin neutralizing antibody (TNA) levels were measured in sera from AVP human vaccinees for both PA and Lethal Factor (LF).It was demonstrated that AVP is composed of at least 138 B. anthracis proteins, including PA (65%), LF (8%) and Edema Factor (EF) (3%), using LC-MS/MS. NetMHCIIpan predicted that peptides from all eight abundant proteins are likely to be presented to T cells, a pre-requisite for protection; however, the number of such peptides varied considerably between different HLA alleles.These analyses highlight two important properties of the AVP vaccine that have not been established previously. Firstly, the effectiveness of AVP within humans may not depend on PA alone; there is compelling evidence to suggest that LF has a protective role, with computational predictions suggesting that additional proteins may be important for individuals with specific HLA allele combinations. Secondly, in spite of differences in the sequences of key antigenic proteins from different B. anthracis strains, these are unlikely to affect the cross-strain protection afforded by AVP.

Cholesterol sensing by CD81 is important for hepatitis C virus entry.

CD81 plays a central role in a variety of physiological and pathological processes. Recent structural analysis of CD81 indicates that it contains an intramembrane cholesterol-binding pocket and that interaction with cholesterol may regulate a conformational switch in the large extracellular domain of CD81. Therefore, CD81 possesses a potential cholesterol-sensing mechanism; however, its relevance for protein function is thus far unknown. In this study we investigate CD81 cholesterol sensing in the context of its activity as a receptor for hepatitis C virus (HCV). Structure-led mutagenesis of the cholesterol-binding pocket reduced CD81-cholesterol association but had disparate effects on HCV entry, both reducing and enhancing CD81 receptor activity. We reasoned that this could be explained by alterations in the consequences of cholesterol binding. To investigate this further we performed molecular dynamic simulations of CD81 with and without cholesterol; this identified a potential allosteric mechanism by which cholesterol binding regulates the conformation of CD81. To test this, we designed further mutations to force CD81 into either the open (cholesterol-unbound) or closed (cholesterol-bound) conformation. The open mutant of CD81 exhibited reduced HCV receptor activity, whereas the closed mutant enhanced activity. These data are consistent with cholesterol sensing switching CD81 between a receptor active and inactive state. CD81 interactome analysis also suggests that conformational switching may modulate the assembly of CD81-partner protein networks. This work furthers our understanding of the molecular mechanism of CD81 cholesterol sensing, how this relates to HCV entry, and CD81's function as a molecular scaffold; these insights are relevant to CD81's varied roles in both health and disease.

Flexibility and intrinsic disorder are conserved features of hepatitis C virus E2 glycoprotein.

The glycoproteins of hepatitis C virus, E1E2, are unlike any other viral fusion machinery yet described, and are the current focus of immunogen design in HCV vaccine development; thus, making E1E2 both scientifically and medically important. We used pre-existing, but fragmentary, structures to model a complete ectodomain of the major glycoprotein E2 from three strains of HCV. We then performed molecular dynamic simulations to explore the conformational landscape of E2, revealing a number of important features. Despite high sequence divergence, and subtle differences in the models, E2 from different strains behave similarly, possessing a stable core flanked by highly flexible regions, some of which perform essential functions such as receptor binding. Comparison with sequence data suggest that this consistent behaviour is conferred by a network of conserved residues that act as hinge and anchor points throughout E2. The variable regions (HVR-1, HVR-2 and VR-3) exhibit particularly high flexibility, and bioinformatic analysis suggests that HVR-1 is a putative intrinsically disordered protein region. Dynamic cross-correlation analyses demonstrate intramolecular communication and suggest that specific regions, such as HVR-1, can exert influence throughout E2. To support our computational approach we performed small-angle X-ray scattering with purified E2 ectodomain; this data was consistent with our MD experiments, suggesting a compact globular core with peripheral flexible regions. This work captures the dynamic behaviour of E2 and has direct relevance to the interaction of HCV with cell-surface receptors and neutralising antibodies.

Germline immunoglobulin genes: disease susceptibility genes hidden in plain sight?

Immunoglobulin genes are rarely considered as disease susceptibility genes despite their obvious and central contributions to immune function. This appears to be a consequence of historical views on antibody repertoire formation that no longer stand, and of difficulties that until recently surrounded the documentation of the suite of antibody genes in any individual. If these important genes are to be accessible to GWAS studies, allelic variation within the human population needs to be better documented, and a curated set of genomic variations associated with antibody genes needs to be formulated. Repertoire studies arising from the COVID-19 pandemic provide an opportunity to meet these needs, and may provide insights into the profound variability that is seen in outcomes to this infection.

OGRDB: a reference database of inferred immune receptor genes.

High-throughput sequencing of the adaptive immune receptor repertoire (AIRR-seq) is providing unprecedented insights into the immune response to disease and into the development of immune disorders. The accurate interpretation of AIRR-seq data depends on the existence of comprehensive germline gene reference sets. Current sets are known to be incomplete and unrepresentative of the degree of polymorphism and diversity in human and animal populations. A key issue is the complexity of the genomic regions in which they lie, which, because of the presence of multiple repeats, insertions and deletions, have not proved tractable with short-read whole genome sequencing. Recently, tools and methods for inferring such gene sequences from AIRR-seq datasets have become available, and a community approach has been developed for the expert review and publication of such inferences. Here, we present OGRDB, the Open Germline Receptor Database (https://ogrdb.airr-community.org), a public resource for the submission, review and publication of previously unknown receptor germline sequences together with supporting evidence.

VDJbase: an adaptive immune receptor genotype and haplotype database.

VDJbase is a publicly available database that offers easy searching of data describing the complete sets of gene sequences (genotypes and haplotypes) inferred from adaptive immune receptor repertoire sequencing datasets. VDJbase is designed to act as a resource that will allow the scientific community to explore the genetic variability of the immunoglobulin (Ig) and T cell receptor (TR) gene loci. It can also assist in the investigation of Ig- and TR-related genetic predispositions to diseases. Our database includes web-based query and online tools to assist in visualization and analysis of the genotype and haplotype data. It enables users to detect those alleles and genes that are significantly over-represented in a particular population, in terms of genotype, haplotype and gene expression. The database website can be freely accessed at https://www.vdjbase.org/, and no login is required. The data and code use creative common licenses and are freely downloadable from https://bitbucket.org/account/user/yaarilab/projects/GPHP.

sumrep: A Summary Statistic Framework for Immune Receptor Repertoire Comparison and Model Validation.

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

Factor VIII cross-matches to the human proteome reduce the predicted inhibitor risk in missense mutation hemophilia A.

Single missense mutations in the F8 gene encoding the coagulation protein factor VIII give rise predominantly to non-severe hemophilia A. Despite only a single amino acid sequence difference between the replacement, therapeutic factor VIII and the patient's endogenous factor VIII, therapeutic factor VIII may still be perceived as foreign by the recipient's immune system and trigger an immune response (inhibitor). Inhibitor formation is a life-long risk for patients with non-severe hemophilia A treated with therapeutic factor VIII, but remains difficult to predict. The aim of this study was to understand whether fortuitous, primary sequence cross-matches between therapeutic factor VIII and proteins in the human proteome are the reason why certain F8 mutations are not associated with inhibitor formation. We predicted which therapeutic factor VIII differences are potentially perceived as foreign by helper T cells - a necessary precursor to inhibitor development - and then scanned potentially immunogenic peptides against more than 100,000 proteins in the proteome. As there are hundreds of disease-causing F8 missense mutations and the human leukocyte antigen gene complex governing peptide presentation to helper T cells is highly polymorphic, these calculations pose a huge combinatorial challenge that we addressed computationally. We found that cross-matches between therapeutic factor VIII and the human proteome are commonplace and have a profound impact on the predicted risk of inhibitor development. Our results emphasize the importance of knowing both the F8 missense mutation and the human leukocyte antigen alleles of a patient with missense mutation hemophilia A if his underlying risk of inhibitor development is to be estimated.

Increased Epitope Complexity Correlated with Antibody Affinity Maturation and a Novel Binding Mode Revealed by Structures of Rabbit Antibodies against the Third Variable Loop (V3) of HIV-1 gp120.

The third variable (V3) loop of HIV-1 gp120 is an immunodominant region targeted by neutralizing antibodies (nAbs). Despite limited breadth, better characterization of the structural details of the interactions between these nAbs and their target epitopes would enhance our understanding of the mechanism of neutralization and facilitate designing better immunogens to induce nAbs with greater breadth. Recently, we isolated two anti-V3 neutralizing monoclonal antibodies (MAbs), 10A3 and 10A37, from a rabbit immunized with gp120 of the M group consensus sequence. In this study, crystal structures of these MAbs bound to target epitopes were determined. 10A3 binds to the V3 crown (303TRKSIHIGPGRAF317) using the cradle binding mode, similar to human V3 MAbs encoded by IGHV5-51 germ line genes, and its epitope structure resembles that bound to the human antibodies. In contrast, 10A37, which exhibits greater breadth and potency than 10A3, binds the V3 crown and the succeeding stem region (308HIGPGRAFYTTGEI323). Unexpectedly, the 315RAFYTT320 portion of the epitope existed as helical turns, a V3 structure that has not been observed previously. Its main chain-dominated antigen-antibody interactions not only explain the broad neutralization of 10A37 but also show that its epitope is a potential vaccine target to be further evaluated. In conclusion, our study provides novel insights about neutralization-susceptible epitope structures of the V3 loop of HIV-1 gp120 and demonstrates that, despite low amino acid sequence similarity to human antibody germ line genes, rabbits can serve as a useful animal model to evaluate human vaccine candidates.IMPORTANCE The apex crown of V3 of HIV-1 gp120 is the most immunogenic region of the surface glycoprotein, and many MAbs targeting this region have been developed. Structural understanding of V3 crown MAbs not only can help understand how antibody responses target this unique region but also contribute to immunogen design for vaccine development. We present here crystal structures of two neutralizing V3 MAbs, 10A3 and 10A37, developed from a rabbit immunized with gp120. Our analysis of 10A3 in complex with V3 provided a detailed example of how epitope complexity can evolve with affinity maturation, while that of 10A37 revealed a novel V3 binding mode targeting the C-terminal side of the V3 crown and showed that this region can form a helical structure. Our study provides novel insights about neutralization-susceptible V3 epitope structures and demonstrates that rabbits can serve as a useful animal model to evaluate human vaccine candidates.

Investigating Substitutions in Antibody-Antigen Complexes Using Molecular Dynamics: A Case Study with Broad-spectrum, Influenza A Antibodies.

In studying the binding of host antibodies to the surface antigens of pathogens, the structural and functional characterization of antibody-antigen complexes by X-ray crystallography and binding assay is important. However, the characterization requires experiments that are typically time consuming and expensive: thus, many antibody-antigen complexes are under-characterized. For vaccine development and disease surveillance, it is often vital to assess the impact of amino acid substitutions on antibody binding. For example, are there antibody substitutions capable of improving binding without a loss of breadth, or antigen substitutions that lead to antigenic escape? The questions cannot be answered reliably from sequence variation alone, exhaustive substitution assays are usually impractical, and alanine scans provide at best an incomplete identification of the critical residue-residue interactions. Here, we show that, given an initial structure of an antibody bound to an antigen, molecular dynamics simulations using the energy method molecular mechanics with Generalized Born surface area (MM/GBSA) can model the impact of single amino acid substitutions on antibody-antigen binding energy. We apply the technique to three broad-spectrum antibodies to influenza A hemagglutinin and examine both previously characterized and novel variant strains observed in the human population that may give rise to antigenic escape. We find that in some cases the impact of a substitution is local, while in others it causes a reorientation of the antibody with wide-ranging impact on residue-residue interactions: this explains, in part, why the change in chemical properties of a residue can be, on its own, a poor predictor of overall change in binding energy. Our estimates are in good agreement with experimental results-indeed, they approximate the degree of agreement between different experimental techniques. Simulations were performed on commodity computer hardware; hence, this approach has the potential to be widely adopted by those undertaking infectious disease research. Novel aspects of this research include the application of MM/GBSA to investigate binding between broadly binding antibodies and a viral glycoprotein; the development of an approach for visualizing substrate-ligand interactions; and the use of experimental assay data to rescale our predictions, allowing us to make inferences about absolute, as well as relative, changes in binding energy.

Evaluation of a novel multi-immunogen vaccine strategy for targeting 4E10/10E8 neutralizing epitopes on HIV-1 gp41 membrane proximal external region.

The membrane proximal external region (MPER) of HIV-1 gp41 is targeted by broadly neutralizing antibodies (bnAbs) 4E10 and 10E8. In this proof-of-concept study, we evaluated a novel multi-immunogen vaccine strategy referred to as Incremental, Phased Antigenic Stimulation for Rapid Antibody Maturation (IPAS-RAM) to induce 4E10/10E8-like bnAbs. Rabbits were immunized sequentially, but in a phased manner, with three immunogens that are progressively more native (gp41-28×3, gp41-54CT, and rVV-gp160DH12). Although nAbs were not induced, epitope-mapping analyses indicated that IPAS-RAM vaccination was better able to target antibodies towards the 4E10/10E8 epitopes than homologous prime-boost immunization using gp41-28×3 alone. MPER-specific rabbit monoclonal antibodies were generated, including 9F6. Although it lacked neutralizing activity, the target epitope profile of 9F6 closely resembled those of 4E10 and 10E8 (671NWFDITNWLWYIK683). B-cell repertoire analyses suggested the importance of co-immunizations for maturation of 9F6, which warrants further evaluation of our IPAS-RAM vaccine strategy using an improved priming immunogen.

Metabolic Pathway Mining.

Understanding metabolic pathways is one of the most important fields in bioscience in the post-genomic era, but curating metabolic pathways requires considerable man-power. As such there is a lack of reliable, experimentally verified metabolic pathways in databases and databases are forced to predict all but the most immediately useful pathways.Text-mining has the potential to solve this problem, but while sophisticated text-mining methods have been developed to assist the curation of many types of biomedical networks, such as protein-protein interaction networks, the mining of metabolic pathways from the literature has been largely neglected by the text-mining community. In this chapter we describe a pipeline for the extraction of metabolic pathways built on freely available open-source components and a heuristic metabolic reaction extraction algorithm.

Studying Antibody Repertoires with Next-Generation Sequencing.

Next-generation sequencing is making it possible to study the antibody repertoire of an organism in unprecedented detail, and, by so doing, to characterize its behavior in the response to infection and in pathological conditions such as autoimmunity and cancer. The polymorphic nature of the repertoire poses unique challenges that rule out the use of many commonly used NGS methods and require tradeoffs to be made when considering experimental design.We outline the main contexts in which antibody repertoire analysis has been used, and summarize the key tools that are available. The humoral immune response to vaccination has been a particular focus of repertoire analyses, and we review the key conclusions and methods used in these studies.

HVint: A Strategy for Identifying Novel Protein-Protein Interactions in Herpes Simplex Virus Type 1.

Human herpesviruses are widespread human pathogens with a remarkable impact on worldwide public health. Despite intense decades of research, the molecular details in many aspects of their function remain to be fully characterized. To unravel the details of how these viruses operate, a thorough understanding of the relationships between the involved components is key. Here, we present HVint, a novel protein-protein intraviral interaction resource for herpes simplex virus type 1 (HSV-1) integrating data from five external sources. To assess each interaction, we used a scoring scheme that takes into consideration aspects such as the type of detection method and the number of lines of evidence. The coverage of the initial interactome was further increased using evolutionary information, by importing interactions reported for other human herpesviruses. These latter interactions constitute, therefore, computational predictions for potential novel interactions in HSV-1. An independent experimental analysis was performed to confirm a subset of our predicted interactions. This subset covers proteins that contribute to nuclear egress and primary envelopment events, including VP26, pUL31, pUL40, and the recently characterized pUL32 and pUL21. Our findings support a coordinated crosstalk between VP26 and proteins such as pUL31, pUS9, and the CSVC complex, contributing to the development of a model describing the nuclear egress and primary envelopment pathways of newly synthesized HSV-1 capsids. The results are also consistent with recent findings on the involvement of pUL32 in capsid maturation and early tegumentation events. Further, they open the door to new hypotheses on virus-specific regulators of pUS9-dependent transport. To make this repository of interactions readily accessible for the scientific community, we also developed a user-friendly and interactive web interface. Our approach demonstrates the power of computational predictions to assist in the design of targeted experiments for the discovery of novel protein-protein interactions.

Missed, Not Missing: Phylogenomic Evidence for the Existence of Avian FoxP3.

The Forkhead box transcription factor FoxP3 is pivotal to the development and function of regulatory T cells (Tregs), which make a major contribution to peripheral tolerance. FoxP3 is believed to perform a regulatory role in all the vertebrate species in which it has been detected. The prevailing view is that FoxP3 is absent in birds and that avian Tregs rely on alternative developmental and suppressive pathways. Prompted by the automated annotation of foxp3 in the ground tit (Parus humilis) genome, we have questioned this assumption. Our analysis of all available avian genomes has revealed that the foxp3 locus is missing, incomplete or of poor quality in the relevant genomic assemblies for nearly all avian species. Nevertheless, in two species, the peregrine falcon (Falco peregrinus) and the saker falcon (F. cherrug), there is compelling evidence for the existence of exons showing synteny with foxp3 in the ground tit. A broader phylogenomic analysis has shown that FoxP3 sequences from these three species are similar to crocodilian sequences, the closest living relatives of birds. In both birds and crocodilians, we have also identified a highly proline-enriched region at the N terminus of FoxP3, a region previously identified only in mammals.