Evolution of the Beta Defensin 2 Gene in Primates
Total Page:16
File Type:pdf, Size:1020Kb
Genes and Immunity (2003) 4, 251–257 & 2003 Nature Publishing Group All rights reserved 1466-4879/03 $25.00 www.nature.com/gene Evolution of the beta defensin 2 gene in primates M Boniotto1, A Tossi2, M DelPero3, S Sgubin1, N Antcheva2, D Santon1, J Masters4 and S Crovella1 1Dipartimento di Scienze della Riproduzione e dello Sviluppo, Universita` di Trieste, Trieste, Italy; 2Dipartimento di Biochimica, Biofisica e Chimica delle Macromolecole, Universita` di Trieste, Trieste, Italy; 3Dipartimento di Biologia Animale e dell’Uomo, Universita` di Torino, Torino, Italy; 4Natal Museum, Pietermaritzburg, South Africa With the aim of further investigating the molecular evolution of beta defensin genes, after having analysed beta defensin 1 (DEFB1) in humans and several nonhuman primate species, we have studied the evolution of the beta defensin 2 gene (DEFB2), which codifies for a peptide with antimicrobial and chemoattractant activity, in humans and 16 primate species. We have found evidence of positive selection during the evolution of orthologous DEFB2 genes at two points on a phylogenetic tree relating these primates: during the divergence of the platyrrhines from the catarrhines and during the divergence of the Cercopithecidae from the Hylobatidae, Great Apes and humans. Furthermore, amino acid variations in Old World Monkeys seem to centre either on residues that are involved in oligomerisation in the human molecule, or that are conserved (40–80%) in beta-defensins in general. It is thus likely that these variations affect the biological function of the molecules and suggest that their synthesis and functional analysis might reveal interesting new information as to their role in innate immunity. Genes and Immunity (2003) 4, 251–257. doi:10.1038/sj.gene.6363958 Keywords: antimicrobial peptides; beta defensins; innate immunity; primates; nucleotide sequences; amino-acid sequences Introduction Results The human beta defensin 2 (hBD2) was originally The alignment of nucleotide sequences of the primate isolated from the skin of psoriatic patients, and was DEFB2 coding region is shown in Table 1. We calculated found to be constitutively expressed in foreskin, lung the number of synonymous nucleotide substitutions per and trachea using the RT-PCR method.1 Moreover, synonymous site (dS) and the number of nonsynon- different authors have demonstrated that the hBD2 gene ymous nucleotide substitutions per nonsynonymous is upregulated upon inflammatory stimuli and its site (dN) by pairwise comparison among the 17 primate expression can subsequently be detected in many other DEFB2 sequences (Table 2). We restricted the analysis tissues.2–5 The peptide has a broad antimicrobial activity by studying only the sequence corresponding to the ranging from Gram negative bacteria to yeasts.1 It antimicrobially active, mature DEFB2. We hypothesised appears that hBD2 kills the pathogens via an electrostatic that only this region has been subject to a positive interaction between the oligomeric form with and a selection caused by direct interaction with pathogens negatively charged microbial surface.6 It has also been and cellular receptors. We found that the dN value was demonstrated that hBD2 shows a chemoattractant higher than dS in several of the nucleotide sequence activity for immature dendritic cells and memory T cells.7 comparisons, indicating that a positive selection pressure Different genes coding for proteins involved, directly acts at the level of this region (see Table 2). or indirectly, in the innate immune response have In order to identify at which node the active portion evolved rapidly by positive selection in response to of the molecule was altered under selection, we have changes in the pathogens to which the host is exposed.8,9 computed the number of synonymous (bS) and non- In this study, we investigate the molecular evolution of synonymous (bN) substitutions for each branch of the beta defensin 2 (DEFB2) gene in humans and 16 primate phylogenetic tree depicting the evolutionary relation species. Comparisons among orthologous DEFB2 genes among the primate species under investigation. For this from several primate species give new insights into the analysis, we used a phylogenetic tree that is widely evolution of the DEFB2 gene and the relation between accepted by primatologists (Figure 1).10 Several branches the structure of the protein and its function. showed greater values of bN than bS, but bN was significantly greater than bS in branch x (Z ¼ 2.56; Po0.05) and in branch y (Z ¼ 2.75; Po 0.05). These results indicate that positive selection has acted in the lineage leading to the clade including the homind and hylobatid families. The deduced amino-acid sequences for the primate Correspondence: Dr S Crovella, Dipartimento di Scienze della Riprodu- zione e dello Sviluppo, Universita` di Trieste, Via dell’Istria, 65/1- 34100 species are reported in Figure 2. Humans and Great Apes Trieste, Italy. E-mail: [email protected] (GA) showed identical amino-acid sequences, with the DEFB2 evolution in primate M Boniotto et al 252 Table 1 Alignment of nucleotide sequences from coding region of orthologous primate DEFB2 genes HSS ATGAGGGTCTTGTATCTCCTCTTCTCGTTCCTCTTCATATTCCTGATGCCTCTTCCAGGTGTTTTTGGTGGTATAGGCGA PTR ...........................................................................A.... PPY ................................................................................ GGO ................................................................................ HMO ........GC.............T..A....C......................................A....A.GA. HCO ......................................................................A....A.GA. HLA .........C................C....C......................................A....A.GA. PME .........C............................................................A....A.GA. PCR ....A....C............................................................A....A.GA. POB ....A....C............................................................A....A.GA. MMU ....A....C..........G.................................................A....A.GA. MFA ....A....C..........G.................................................A....A.GA. PPA ....A....C............................................................A....A.GA. CAE ....A....C...................................................C.C......A....A.GA. CER ....A....C...................................................C........A....A.GA. CPR ....A....C...................................................C........A....C.GA. CJA ....A....C..........G.................................................A..C.A.GA. HSS TCCTGTTACCTGCCTTAAGAGTGGAGCCATATGTCATCCAGTCTTTTGCCCTAGAAGGTATAAACAAATTGGCACCTGTG PTR .........................................................A...................... PPY ................................................................................ GGO ................................................................................ HMO ........................T............................................C.......... HCO ........................T............................................C.......... HLA ........................T............................................C.......... PME ....A...................T................G..........C.............T..C...GT..... PCR ....A...................T................G..........G.............T..C...GT..... POB ....A...................T................G..........C.............T..C...GT..... MMU .............G...G......T..........T.....G........................T......GT..... MFA .............G...G......T..........T.....G........................T......GT..... PPA .............G...G......T..........T.....G........................T......GT..... CAE .................G......T................G........................T..C...GT..... CER .................G......T................G........................T..C...GT..... CPR .................G......T................G........................T..C..TGT..... CJA ..T.T........A..TG......T................G..........G...C.........T..C....T..... HSS GTCTCCCTGGAACAAAATGCTGCAAAAAGCCATGA PTR ................................... PPY ................................... GGO ................................... HMO ................................... HCO .......................CG.......... HLA .......................CG.......... PME ..G..T...C..T...................... PCR ..G..T...C..T...................... POB ..G..T...C..T...................... MMU ..G..T...C..T...................... MFA ..G..T...C..T...................... PPA ..G..T...C..T...................... CAE .CG..T...C..T...................... CER ..G..T...C..T...................... CPR ..G..T...C..T...................... CJA ..G.....CT..T...................... exception of a Gly3-Ser substitution in the N-terminal based on the hBD2 sequence, relating these residues region seen only in P. troglodytes. The hylobatid sequences to the secondary structure elements is also provided show very little variation with respect to the hBD2. in Figure 2. Amino-acid changes found in the gibbon Three amino acid substitutions were observed in the (Hylobatidae) sequences, with respect to the human and N-terminal region of the active peptide, while only one GA sequences, are located exclusively at the N-terminus, conservative substitution was present in the C-terminal preceding the short a-helical stretch. Cercopithecid region. The Old World Monkeys (OWM, Cercopitheci- sequence variations are more extensive and mainly dae) are characterised by very similar amino acid involve the