# Molecular evolution of the human SRPX2 gene that causes brain disorders of the Rolandic and Sylvian speech areas

Barbara Royer, Dinesh C Soares, Paul N Barlow, Ronald E Bontrop, Patrice Roll, Andrée Robaglia-Schlupp, Antoine Blancher, Anthony Levasseur, Pierre Cau, Pierre Pontarotti, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Barbara Royer, Dinesh C Soares, Paul N Barlow, Ronald E Bontrop, Patrice Roll, et al. Molecular evolution of the human SRPX2 gene that causes brain disorders of the Rolandic and Sylvian speech areas. BMC Genetics, 2007, 8, pp.72. 10.1186/1471-2156-8-72. hal-03509486

HAL Id: hal-03509486

https://hal.science/hal-03509486

Submitted on 7 Jan 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
Research article

Open Access

# Molecular evolution of the human SRPX2 gene that causes brain disorders of the Rolandic and Sylvian speech areas

Barbara Royer¹, Dinesh C Soares², Paul N Barlow², Ronald E Bontrop³, Patrice Roll¹, Andrée Robaglia-Schlupp¹, Antoine Blancher⁴, Anthony Levasseur⁵, Pierre Cau¹, Pierre Pontarotti⁵ and Pierre Szepetowski*¹

Address: ¹INSERM UMR 491, Université de la Méditerranée, 13385 Marseille, Cedex 5, France, ²Edinburgh Biomolecular NMR Unit, The University of Edinburgh, Edinburgh EH9 3JJ, UK, ³Biomedical Primate Research Centre, 2280 GH Rijswijk, The Netherlands, ⁴LIGM, Université Paul Sabatier, 31403 Toulouse, Cedex 4, France and ⁵EA 3781, Université d'Aix-Marseille I, 13331 Marseille, Cedex 3, France

Email: Barbara Royer - barbara.royer@medecine.univ-mrs.fr; Dinesh C Soares - D.C.Soares@sms.ed.ac.uk; Paul N Barlow - paul.barlow@ed.ac.uk; Ronald E Bontrop - bontrop@bprc.nl; Patrice Roll - patrice-roll@medecine.univ-mrs.fr; Andrée Robaglia-Schlupp - robaglia@medecine.univ-mrs.fr; Antoine Blancher - blancher@easynet.fr; Anthony Levasseur - levassr@up.univ-mrs.fr; Pierre Cau - pierre.cau@medecine.univ-mrs.fr; Pierre Pontarotti - pierre.pontarotti@up.univ-mrs.fr; Pierre Szepetowski* - szepetowski@medecine.univ-mrs.fr

Published: 18 October 2007

Received: 29 March 2007

BMC Genetics 2007, 8:72 doi:10.1186/1471-2156-8-72

Accepted: 18 October 2007

This article is available from: http://www.biomedcentral.com/1471-2156/8/72

© 2007 Royer et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Abstract

**Background:** The X-linked SRPX2 gene encodes a Sushi Repeat-containing Protein of unknown function and is mutated in two disorders of the Rolandic/Sylvian speech areas. Since it is linked to defects in the functioning and the development of brain areas for speech production, SRPX2 may thus have participated in the adaptive organization of such brain regions. To address this issue, we have examined the recent molecular evolution of the SRPX2 gene.

**Results:** The complete coding region was sequenced in 24 human X chromosomes from worldwide populations and in six representative nonhuman primate species. One single, fixed amino acid change (R75K) has been specifically incorporated in human SRPX2 since the human-chimpanzee split. The R75K substitution occurred in the first sushi domain of SRPX2, only three amino acid residues away from a previously reported disease-causing mutation (Y72S). Three-dimensional structural modeling of the first sushi domain revealed that Y72 and K75 are both situated in the hypervariable loop that is usually implicated in protein-protein interactions. The side-chain of residue 75 is exposed, and is located within an unusual and SRPX-specific protruding extension to the hypervariable loop. The analysis of non-synonymous/synonymous substitution rate (Ka/Ks) ratio in primates was performed in order to test for positive selection during recent evolution. Using the branch models, the Ka/Ks ratio for the human branch was significantly different (p = 0.027) from that of the other branches. In contrast, the branch-site tests did not reach significance. Genetic analysis was also performed by sequencing 9,908 kilobases (kb) of intronic SRPX2 sequences. Despite low nucleotide diversity, neither the HKA (Hudson-Kreitman-Aguadé) test nor the Tajima's D test reached significance.

**Conclusion:** The R75K human-specific variation occurred in an important functional loop of the first sushi domain of SRPX2, indicating that this evolutionary mutation may have functional importance; however, positive selection for R75K could not be demonstrated. Nevertheless, our data contribute to the first understanding of molecular evolution of the human SPRX2 gene. Further experiments are now required in order to evaluate the possible consequences of R75K on SRPX2 interactions and functioning.

<footnote>* Corresponding author</footnote>
## Background

Evolution studies have been undertaken to identify those genetic changes that underlie human-specific features such as susceptibility to acquired immunodeficiency syndrome, bipedalism, a large brain, and higher-order cognitive functions. Several phenotypic differences distinguishing human from other great apes species obviously rely on cerebral activity. Large-scale studies in human and chimpanzee using either genome comparisons [1,2] or brain transcriptome analyses [3-5] have led to the identification of a subset of genes that may have contributed to the evolution of human brain anatomy and activity from a common primate ancestor. An important complementary approach has relied on the study of candidate genes selected on the basis of their importance in specific human phenotypes. Consequently, several genes involved in the structure and/or functioning of the human brain have been associated with recent positive selection: ASPM [6,7], MCPH1 [8-10], GLUD2 [11], MAOA [12,13], SHH [14], and the "speech gene" FOXP2 [15-17]. More recently, accelerated evolution of noncoding sequences has also been shown [18,19].

The Rolandic and Sylvian fissures divide the cortex hemispheres of primates into their main anatomical structures. In human, these areas participate in speech production under the control of the Broca's area. We recently identified the *SRPX2* gene as being responsible for two related disorders of the Rolandic and Sylvian speech areas [20,21]. Since it is linked to defects in the functioning and the development of such brain regions, such as epileptic seizures, oral and speech dyspraxia, or bilateral perisylvian polymicrogyria, *SRPX2* may be one of the specific genes whose evolution at the DNA-level may have participated in the recent emergence of higher-order cognitive functions, including the adaptive organization of brain areas for speech production.

In this study, we have examined the molecular evolution of the *SRPX2* gene. One single, fixed amino acid change occurred in the first sushi domain (also known as CCP - complement control protein - module, or short consensus repeat) of *SRPX2* after the human-chimpanzee split. Three-dimensional modeling showed that both this evolutionary mutation and a previously identified disease-associated mutation [20] lie within a hypervariable loop shared by all sushi modules and that has been implicated in some cases in protein-protein interactions [22]. Using the branch models, the synonymous/non-synonymous analysis was consistent with accelerated evolution in the human lineage but this could not be confirmed when the branch-site models were used. Population genetics tests did not reach statistical significance, indicating either that a selective sweep may have occurred more than 100 000-

200 000 years ago, or that there has been no episode of positive selection on *SRPX2*.

## Results and discussion

**One single amino acid substitution (R75K) has been fixed in human SRPX2 since the human-chimpanzee split**

*In silico* screening of publicly available databases followed by phylogenetic analysis showed that *SRPX2* belongs to a family of five genes: *SRPX2* itself (sushi-repeat-containing protein, X-linked 2), *SRPX* (sushi-repeat-containing protein, X-linked), *SELP* (selectin P precursor), *SELE* (selectin E precursor), and *SVEP1* (selectin-like protein) gene (Fig. 1; Additional file 1). This family emerged during vertebrate evolution. In order to trace the recent evolutionary history of *SRPX2*, all coding exons were sequenced from a subset of primate species that represent key evolutionary steps: human [Genbank: NM_014467], chimpanzee [Genbank: EF369515], gorilla [Genbank: EF369516], orangutan [Genbank: EF369517], gibbon [Genbank: EF369518], macaque [Genbank: EF369519], and baboon [Genbank: EF369520]. Sequence comparison (Fig. 2)

<figure><img src="image_3.png" /><figcaption>**Figure 1**</figcaption></figure>

Phylogenetic tree of the *SRPX2* family members. ag: Anopheles gambiae; bt: Bos taurus; cn: Canis familiaris; dr: Danio rerio; gg: Gallus gallus; hs: Homo sapiens; mm: Mus musculus; ma: Macaca mulatta; rn: Rattus norvegicus; pt: Pan troglodytes; tn: Takifugu rubripes; tr: Tetraodon nigroviridis. The list of the ENSEMBL gene accession numbers used to construct the phylogenetic tree is available in Additional file 1.
<table><caption>**Figure 2**</caption><tbody><tr><td>Human</td><td>MA$S$LTQRGALFLLFFLPTAVPTWYAGSGYYPDE$S$NEVYAEVEPQAPALDYRVPRWCY</td><td>60</td></tr><tr><td>Chimp</td><td></td><td>60</td></tr><tr><td>Gorilla</td><td></td><td>60</td></tr><tr><td>Orangutan</td><td></td><td>60</td></tr><tr><td>Gibbon</td><td></td><td>60</td></tr><tr><td>Macaque</td><td></td><td>60</td></tr><tr><td>Baboon</td><td></td><td>60</td></tr><tr><td>Human</td><td>T$N$IQDGEATCYSF$G$NYHSSLGTRELC$D$RGF$LIG$R$SVQCL$PSR$R$W$GTAYVC$RQ$M</td><td>120</td></tr><tr><td>Chimp</td><td></td><td>120</td></tr><tr><td>Gorilla</td><td></td><td>120</td></tr><tr><td>Orangutan</td><td></td><td>120</td></tr><tr><td>Gibbon</td><td></td><td>120</td></tr><tr><td>Macaque</td><td></td><td>120</td></tr><tr><td>Baboon</td><td></td><td>120</td></tr><tr><td>Human</td><td>R$C$HALPFITS$G$TYTCTNG$V$LLO$D$R$CDY$C$SSGYH$LE$D$G$R$D$R$C$M$E$D$G$R$W$S$G$O$G$E$P$VC$D$D</td><td>180</td></tr><tr><td>Chimp</td><td></td><td>180</td></tr><tr><td>Gorilla</td><td></td><td>180</td></tr><tr><td>Orangutan</td><td></td><td>180</td></tr><tr><td>Gibbon</td><td></td><td>180</td></tr><tr><td>Macaque</td><td></td><td>180</td></tr><tr><td>Baboon</td><td></td><td>180</td></tr><tr><td>Human</td><td>P$P$KIRCPHSRE$R$MAE$P$EK$LT$AR$V$WD$P$LVKDSAG$T$IT$R$VT$LRG$E$P$G$SH$F$E$G$EH$V$IR</td><td>240</td></tr><tr><td>Chimp</td><td></td><td>240</td></tr><tr><td>Gorilla</td><td></td><td>240</td></tr><tr><td>Orangutan</td><td></td><td>240</td></tr><tr><td>Gibbon</td><td></td><td>240</td></tr><tr><td>Macaque</td><td></td><td>240</td></tr><tr><td>Baboon</td><td></td><td>240</td></tr><tr><td>Human</td><td>Y$T$AY$D$RAY$N$R$AS$C$F$IV$K$V$Q$V$RC$P$T$K$P$P$Q$H$GY$L$T$C$T$AG$D$N$Y$G$AT$C$E$Y$H$C$D$G$O$G$Y$D$R$Q$G</td><td>300</td></tr><tr><td>Chimp</td><td></td><td>300</td></tr><tr><td>Gorilla</td><td></td><td>300</td></tr><tr><td>Orangutan</td><td></td><td>300</td></tr><tr><td>Gibbon</td><td></td><td>300</td></tr><tr><td>Macaque</td><td></td><td>300</td></tr><tr><td>Baboon</td><td></td><td>300</td></tr><tr><td>Human</td><td>T$P$SRVCQ$SSR$Q$W$SG$P$PIC$A$P$K$IM$N$V$N$NSA$G$ALL$D$Q$F$Y$E$K$Q$R$LL$I$S$A$P$D$PS$N$R$Y$K$M$Q$</td><td>360</td></tr><tr><td>Chimp</td><td></td><td>360</td></tr><tr><td>Gorilla</td><td></td><td>360</td></tr><tr><td>Orangutan</td><td></td><td>360</td></tr><tr><td>Gibbon</td><td></td><td>360</td></tr><tr><td>Macaque</td><td></td><td>360</td></tr><tr><td>Baboon</td><td></td><td>360</td></tr><tr><td>Human</td><td>$I$S$M$L$Q$Q$STC$G$LD$R$LV$N$VT$II$E$LV$G$OP$P$Q$E$V$G$RI$EQ$R$Q$LS$AN$I$E$L$R$P$Q$R$L$TR$SY$F$N$M$VL</td><td>420</td></tr><tr><td>Chimp</td><td></td><td>420</td></tr><tr><td>Gorilla</td><td></td><td>420</td></tr><tr><td>Orangutan</td><td></td><td>420</td></tr><tr><td>Gibbon</td><td></td><td>420</td></tr><tr><td>Macaque</td><td></td><td>420</td></tr><tr><td>Baboon</td><td></td><td>420</td></tr><tr><td>Human</td><td>$I$D$Q$G$ID$R$D$RY$ME$P$VT$P$E$I$FT$F$ID$Y$LL$SN$EL$T$Q$R$EQ$R$D$ICE</td><td>465</td></tr><tr><td>Chimp</td><td></td><td>465</td></tr><tr><td>Gorilla</td><td></td><td>465</td></tr><tr><td>Orangutan</td><td></td><td>465</td></tr><tr><td>Gibbon</td><td></td><td>465</td></tr><tr><td>Macaque</td><td></td><td>465</td></tr><tr><td>Baboon</td><td></td><td>465</td></tr></tbody></table>

Amino acid sequence alignment of SRPX2 from human and nonhuman primates. Dots represent residues identical to the human amino acid sequence. Black lines represent the three Sushi domains and the dotted line represents the HYR domain. The site of the human-specific variation (R75K) is boxed.

detected only five amino acid variant sites (1.1%; 5/465). One discrepancy between the protein [Genbank: ABN46998] predicted from the five chimpanzee SRPX2 sequences obtained here (see Methods section), and the chimpanzee SRPX2 protein [Genbank: XP_521170] predicted from the genomic sequence at the UCSC database [23], was detected (D429 in ABN46998; N429 in XP_521170). This likely corresponds to an error in the chimpanzee genome sequence previously available, although it may be due to chimpanzee polymorphism. The only amino acid substitution specific to the human lineage is an Arg to Lys change at position 75, corresponding to a 224A>G mutation within exon 4 of SRPX2. The possible importance of the R75K substitution in the evolution of the human species was indicated by the conservation of R75 in all SRPX2 orthologues, from nonhuman primates to ray-finned fishes, as well as by the fixation of K75 in human, as demonstrated by the lack of any variation at position 75 in silico, in the 24 worldwide X-chromosomes tested here, and in the 624 Caucasian control X-chromosomes previously screened [20].

**K75 is situated in the hypervariable loop that is usually implicated in protein-protein interactions**

The R75K human-specific modification occurred in the first sushi domain of the protein. Sushi domains have been identified in several proteins of the complement system and in the selectin family of proteins [24]. They may serve in protein-protein interactions [22], as demonstrated in the case of the neurocan-L1 interaction [25]. R75K occurs only three amino acid residues away from the tyrosine residue (Y72) that is mutated in a patient with rolandic seizures and bilateral perisylvian polymicrogyria and in his female relatives with mild mental retardation [20]. The amino acid residue at position 75 may thus participate in the proper function of the first sushi domain of SRPX2.

To further address this issue, three-dimensional modeling of the first sushi module of the human SRPX2 protein (i.e. sushi 1) was undertaken by homology with a known sushi domain (CCP module) structure (Fig. 3; Additional file 2). Sushi domains are characterised by a compact hydrophobic core, containing an almost invariant Trp residue, which is enclosed in a framework of five extended seg-

<figure><img src="image_4.png" /><figcaption>**Figure 3**</figcaption></figure>

Three-dimensional modeling of the first sushi domain of human SRPX2. A cartoon representation [65] is shown, highlighting the conserved Trp and Cys residues. The model reveals an additional putative disulfide bridge (Cys71-Cys87), atypical of the classical sushi (CCP) module fold. The hyper-variable loop, with the sites of disease-causing mutation (Tyr72) and evolutionary change (Lys75) and the protrusion specific to SRPX, are indicated.
ments that form β-strands for all or part of their lengths. The extended segments are aligned with the long axis of the module and are held together by two strictly conserved disulfide bridges [24]. In a previously performed cluster assignment exercise undertaken for more than 240 sushi domain sequences [26], the first sushi domain of SRPX (which is highly similar to sushi 1) fell into cluster-D, which also includes many modules found within the selectin family of proteins.

In almost all sushi domains, a region that is highly variable in length, sequence and conformation, and commonly referred to as the hypervariable loop [27], is inserted within the second extended region. Depending on its length this projects laterally from the module and forms an obvious candidate surface for protein-protein recognition. Indeed, the hypervariable loop has been implicated as a "hot-spot" for several protein-protein interactions and disease-causing mutations in CCP-containing complement proteins [22,28-30]. In the sushi 1 model (Fig. 3) an additional disulfide bridge (cysteine residues C71-C87), atypical of the classical sushi domain fold, is evident. Cys71 lies at the beginning of the hypervariable loop that is exceptionally long in sushi 1. The hypervariable loop is thus forced to form a prominent protrusion extending towards the N-terminus of the module. This feature has neither been seen in any experimentally determined sushi module structure to date, nor is predicted to occur in other members of the D-cluster except SRPX.

The Y72S mutation is largely solvent exposed and located within the hypervariable loop, adjacent to the cysteine residue (C71) that participates in the non-typical, third, disulfide bridge. This change from a large aromatic side-chain (tyrosine; Y) to a small, polar one (serine; S) at position 72 will have a profound effect on the surface properties of this region that is close to the aforementioned prominent protrusion. K75 is located nearby, within the protrusion, and its side-chain is exposed (Fig. 3). It is reasonable to suggest that the unique structure formed by the hypervariable loop of sushi 1 performs a role that is specific to the SRPX2 protein. Presumably, it is not a coincidence, that it is also the site of the human-specific change. An R75K substitution is a conservative one; the substitution of one exposed, positively charged residue for another can easily be accommodated by small atom shifts in surrounding side-chains, and is not likely to affect the structure of the sushi domain. Such a change might have a small but not a dramatic functional effect. Conservative substitutions are thought to play a role in adaptive change [31]. Moreover, conservative R-to-K and K-to-R substitutions can result in the altered properties of either secreted proteins [32] or the extracellular domains of some plasma

membrane proteins [33,34], including a member of the selectin family [35].

## Analysis of non-synonymous/synonymous substitution rate ratio in primates

Human-specific modifications with putative functional consequences may sustain positive selection since the human-chimpanzee split. The ratio of non-synonymous (Ka) to synonymous (Ks) substitution rates was then calculated for each branch on the primate phylogenetic tree (Fig. 4). Using the branch models, the Ka/Ks value in the human lineage differed from all other branches with statistical significance (p = 0.027) (Additional file 3). In this model, the calculated Ka/Ks ratio has an infinite value, as no synonymous change has accumulated between the human and the chimpanzee. However, this elevated Ka/Ks ratio was not statistically >1 (p = 0.371) (Additional file 3). Using the branch-site tests, model A indicated that the K75 site might be under positive selection in the human lineage. However, model A did not differ from the null hypothesis with statistical significance (p = 0.089) (Additional figure 3). Overall, positive selection in the human lineage could not be clearly demonstrated, as the branch-site analysis did not confirm the data obtained using the branch models. This may be due either to the actual absence of any positive selection, or to the lack of power of the tests due to the low number of genetic variations within SRPX2.

## Population genetics analysis

Evidence for positive selection can also be found by analyzing the intra-specific variation in DNA sequences subjected to a selective sweep. DNA sequences adjacent to an advantageous allele should display lower-than-expected levels of diversity. The levels of polymorphism and diver-

<figure><img src="image_5.png" /><figcaption>**Figure 4**
Numbers of non-synonymous (left) and synonymous (right) substitutions of different evolutionary lineages in primates.</figcaption></figure>
gence were estimated by sequencing a total of 9,908 kb of intronic sequences surrounding exon 4 of SRPX2 in the 24 X-chromosomes mentioned above. Polymorphic sites are shown in Table 1. Nucleotide diversity at the SRPX2 locus (π = 0.00036) was 76% of the estimated nucleotide diversity on chromosome X (π = 0.00047) [36]. The divergence between humans and chimpanzee for SRPX2 (D = 0.0075) was nearly identical to the average divergence calculated for X-linked intronic regions (D = 0.0072) [37].

The McDonald-Kreitman test [38] that measures the fraction of site under positive selection pressure by comparing the ratio of nonsynonymous to synonymous divergence and the ratio of nonsynonymous to synonymous polymorphism, was useless in the case of SRPX2 because of the very few number of inter- and intra-species DNA variations: no synonymous mutation and only one non-synonymous mutation has occurred since the human-chimpanzee split, and neither non-synonymous nor synonymous polymorphisms were detected when 24 human SRPX2 entire coding sequences from two sub-Saharan African, three Asian, three Maghrebian, and four Caucasian women were tested.

Moreover, no polymorphism had been detected in the 198 unrelated patients previously screened for disease-causing mutations in the coding sequence of SRPX2 [20]. It is also noteworthy that while no synonymous polymorphism was found in silico, one non-synonymous variation was detected in silico [dbSNP:rs17851822]. Although the existence of a rare DNA variation cannot be ruled out, this change was detected in one single IMAGE clone only [Genbank: BC020733; IMAGE:4769946]. This is more likely to correspond to a clone artifact, as we indeed detected a large proportion of clones carrying artifactual mutations within SRPX2 inserts, using various *Escherichia coli* strains (unpublished data).

The HKA (Hudson-Kreitman-Aguadé) test [39] was then used in order to compare the intra- and inter-specific variations between the SRPX2 locus and control loci assumed to be under neutral selection, but did not yield significant results (Table 2). The Tajima's D test was also applied to the SRPX2 intron data. Tajima's D was clearly negative (D = -0.646) but did not reach significance (p > 0.05). As in

<table><caption>Table 1: Polymorphic sites of SRPX2 introns 3 and 4 in 12 women worldwide.</caption><thead><tr><th rowspan="2">Site*</th><th colspan="12">Intron 3</th><th colspan="3">Intron 4</th></tr><tr><th>3478</th><th>3749</th><th>4292</th><th>4756</th><th>5115</th><th>6076</th><th>6783</th><th>7862</th><th>8433</th><th>9691</th><th>9700</th><th>10308</th><th>10446</th><th>10732</th><th>11286</th><th>2112</th><th>2202</th><th>2369</th></tr><tr><th>SNPs (ID)</th><th>a</th><th>b</th><th>c</th><th>d</th><th>e</th><th>f</th><th>g</th><th>h</th><th>i</th><th>j</th><th>k</th><th>l</th><th>m</th><th>n</th><th>o</th><th>p</th><th>q</th><th>r</th></tr></thead><tbody><tr><td>Caucasian 1</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Caucasian 2</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/T</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/A</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Caucasian 3</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/C</td><td>G/G</td><td>A/G</td><td>T/C</td><td>T/T</td><td>T/T</td><td>C/T</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Caucasian 4</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Maghrebia 1</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Maghrebia 2</td><td>T/G</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/C</td><td>G/G</td><td>A/G</td><td>T/C</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Maghrebia 3</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Asian 1</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/C</td><td>G/A</td><td>A/G</td><td>T/C</td><td>T/T</td><td>T/T</td><td>C/T</td><td>G/G</td><td>C/C</td><td>C/C</td><td>A/A</td></tr><tr><td>Asian 2</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/A</td><td>G/G</td><td>A/A</td><td>T/T</td><td>T/T</td><td>T/T</td><td>C/C</td><td>G/G</td><td>T/T</td><td>C/C</td><td>A/A</td></tr><tr><td>Asian 3</td><td>T/T</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>G/G</td><td>A/C</td><td>G/A</td><td>A/G</td><td>T/C</td><td>T/T</td><td>T/T</td><td>C/T</td><td>G/G</td><td>T/C</td><td>C/C</td><td>A/A</td></tr><tr><td>African 1</td><td>G/G</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/C</td><td>C/T</td><td>G/G</td><td>C/C</td><td>G/G</td><td>G/G</td><td>C/C</td><td>T/T</td><td>T/T</td><td>C/T</td><td>T/G</td><td>T/T</td><td>C/T</td><td>A/G</td></tr><tr><td>African 2</td><td>G/G</td><td>T/T</td><td>G/G</td><td>C/C</td><td>T/T</td><td>C/C</td><td>A/A</td><td>C/C</td><td>G/G</td><td>G/G</td><td>C/C</td><td>T/T</td><td>C/C</td><td>T/T</td><td>G/G</td><td>T/T</td><td>C/C</td><td>G/G</td></tr></tbody></table>

<footnote>* Numbers refer to the positions in the sequences of human introns 3 and 4, respectively. ** For those single nucleotide polymorphisms (SNPs) that were not reported in current databases [23,66] at the time of analysis and that appeared in one single individual only, PCRs were redone and both strands resequenced in order to exclude PCR artifacts. SNPs IDs at NCBI [66] are: a: rs6616171. b: ss52048551. c: ss52048547. d: ss52085994. e: ss52048549. f: ss52084549. g: ss52084548. h: rs5966711. i: ss52048552. j: ss52048553. k: ss52048554. l: ss52048546. m: rs5921619. n: rs5920840. o: ss52084547. p: rs5921620. q: ss52048548. r: rs2022475.</footnote>
<table><caption>Table 2: Nucleotide polymorphism and divergence at SRPX2</caption><thead><tr><th></th><th>Sequence length (nt)</th><th>π (%)</th><th>θ (%)</th><th>Divergence D (%)</th><th>θ/D</th><th>HKA probability</th></tr></thead><tbody><tr><td>β-globin initiation at 11p15 [67]</td><td>6076</td><td>0.129</td><td>0.107</td><td>1.284</td><td>0.083</td><td>0.696</td></tr><tr><td>Noncoding region at 22q11 [68]</td><td>9901</td><td>0.088</td><td>0.139</td><td>1.353</td><td>0.103</td><td>0.605</td></tr><tr><td>Dystrophin intron-dys44 at Xp21 [69]</td><td>7475</td><td>0.135</td><td>0.102</td><td>0.604</td><td>0.169</td><td>0.499</td></tr><tr><td>PDHA1 introns at Xp22 [70]</td><td>3530</td><td>0.225</td><td>0.211</td><td>0.992</td><td>0.213</td><td>0.168</td></tr><tr><td>Noncoding region at Xq13.3 [71]</td><td>10200</td><td>0.045</td><td>0.083</td><td>0.922</td><td>0.090</td><td>0.655</td></tr><tr><td>SRPX2 introns at Xq22 (this study)</td><td>9908</td><td>0.048 (0.036)</td><td>0.062 (0.047)</td><td>0.750</td><td>0.084</td><td></td></tr></tbody></table>

nt: nucleotide; π: nucleotide diversity per site; for X chromosome data, π is corrected by multiplication by 4/3 and the uncorrected value is in brackets. θ: Watterson's estimate of polymorphism per site; for X chromosome data, θ is corrected by multiplication by 4/3 and the uncorrected value is in brackets. D: number of nucleotide differences per site between human and chimpanzee sequences. HKA probability: probability from the HKA test, with comparison to the SRPX2 intron data.

the case of the ASPM gene that is involved in human brain expansion [6], the present intra-specific analyses did not show significant evidence against neutral expectations. The lack of statistical significance may indicate that no selective sweep has ever occurred in the human lineage. However, it should be mentioned that a selective sweep can be detected by intra-specific studies only for a short period (< 0.5 N generations, where N is the effective human population size, i.e. approximately < 100 000–200 000 years ago) after fixation of the advantageous variant [40,41].

## Conclusion

In this study, we have examined the molecular evolution of the SRPX2 gene that causes brain disorders of the speech areas. One single, fixed amino acid change (R75K) occurred in the first sushi domain of SRPX2 after the human-chimpanzee split. Neither the primate analysis nor the population genetics separately demonstrated the existence of positive selection of SRPX2. Whether the single human-specific R75K mutation has sustained positive selection thus remains an open question. However the 3-D location of R75K right within an important functional domain of the SRPX2 protein, in the immediate vicinity of a pathogenic mutation, indicates that this evolutionary mutation may have functional importance. Because R75K occurred in a domain implicated in protein-protein interactions, it is possible that qualitative or quantitative changes in the interaction with one or several putative SRPX2 partners have been modified. Despite the present lack of any knowledge on either the actual function of SRPX2 or the proteins it interacts with, it is obvious that the comparison of the properties of the human and chimpanzee SRPX2 orthologues will help determine in the future if and how new functions were acquired. From this viewpoint, our study represents a first important step towards the analysis of the consequences of R75K on SRPX2 functioning and protein interactions.

## Methods

### DNA samples and sequencing

All human and nonhuman primate DNAs were extracted using standard procedures and according to the appropriate ethical committees and animal care rules, respectively. Polymerase chain reaction (PCR) fragments representing 9,908 kb from introns 3 and 4 were sequenced with specific primers in 12 human females from the major continental populations (two sub-Saharan Africans, four Europeans, three Asians, and three Magrebians). Each exon of the full-length coding region of SRPX2 was amplified by PCR and sequenced in the same 12 women and in six nonhuman primate species, using inter-specific consensus intronic primers. The nonhuman primates were of the hominidae, pongidae, hylobatidae, and cercopithecidae families: three chimpanzees (one male and two females; *Pan troglodytes*); one gorilla (*Gorilla gorilla*); one orangutan (*Pongo pygmaeus*), one gibbon (*Hylobates sp.*), one macaque (*Macaca mulatta*) and one baboon (*Papio sp.*). The primer sequences are listed in Additional file 4. Sequencing was performed using standard dye terminator chemistry. Sequences were analysed and polymorphisms detected using the Genalys 3.0 software [42]. Multiple sequence alignments were performed using the ClustalW program [43].

### Modeling analyses

Modeling of the first sushi module of human SRPX2 was undertaken based upon its closest homologue in the Protein Data Bank (PDB) [44,45]: the highest resolved X-ray structure of the first sushi module of complement receptor type 2 (PDB ID: 1LY2) [46], which was identified with a BLAST search [47,48]. The target and template shared 29% sequence identity from N- to C-terminal cysteines. The program Modeller release 8 version 1 [49] was used.

The alignment between the target and template sequences was based on an initial multiple sequence alignment of several individual human sushi-module sequences assigned to cluster-D (with the addition of template sequence CR2 module 1), using the program ProbCons
[50] to help position indels more appropriately (see Additional file 2 for target-template alignment). Cluster-D members, according to the recently published convention of sequence assignment [26], are characterized by the presence of a six-cysteine residue motif, forming three putative disulfide bridges. The target-template sequence alignment was subjected to further, minor manual editing guided by positioning of secondary structure elements more appropriately between the target and template sequences; secondary structure was predicted by the server PsiPred version 2.4 [51,52] for the target sequence, or identified by DSSP [53] for the template structure. The additional putative disulfide bridge present in the target (absent in its corresponding template) was restrained during model building. Twenty-five models were generated, and the one with the lowest objective function score [49] selected as the representative model. The loop KGGNY in the hypervariable loop region of SRPX2, for which CR2-template-derived restraints were largely absent was subsequently remodeled under SYBYL version 6.9 (Tripos Associates, St. Louis, MO, USA) after we conducted a loop search for that region against a high resolution subset of the PDB, in order to select the best root-mean square fitting matching loop conformation for that region.

Non-identical side-chain residues of the representative model were optimized using the side-chain replacement program SCWRL version 3 [54]. The model was then protonated and subjected to energy minimization using the Tripos forcefield [55] (20 steps steepest descent followed by 20 steps conjugate gradient) under SYBYL to remove clashes and bad geometries. The model structure was finally checked for valid stereochemistry using PRO-CHECK version 3.5.4 [56].

### Evolution data analyses

The phylogenetic analysis was performed using the phylogenomic analysis pipeline available in the FIGENIX [57] automated genomic annotation platform [58] with the human SRPX2 protein sequence [Genbank: NP_055282] as input and the NCBI nr or the Ensembl databases for BLAST searches (Additional file 1). The synonymous/non synonymous analyses were conducted in primates assuming the following unrooted tree topology: (Baboon, Macaque, (((Chimpanzee, Human), Gorilla), Orangutan), Gibbon)). Ancestral sequences were reconstructed with the pamp and codeml programs in the PAML package [59,60], using the parsimony method and the maximum likelihood method, respectively. From these data, the number of synonymous and non-synonymous substitutions and the Ka/Ks ratios were estimated using the DnaSP 4.0 package [61]. The codeml program from the PAML 3.15 packages was used to test for positive selection, using the branch models and the branch-site models, as previously described [62,63] (Additional file 3).

The DnaSP 4.0 program was also used for all population genetic analyses. Nucleotide diversity (π) and Watterson's θ were computed as described [64]. The neutral evolution hypothesis for the SRPX2 intronic polymorphisms was checked with the HKA (Hudson-Kreitman-Aguadé) test [39], using available data on worldwide polymorphisms in non-coding reference (neutral) autosomal and X-linked regions of the human genome (Table 2), as previously described [16].

### Competing interests

The author(s) declares that there are no competing interests.

### Authors' contributions

BR carried out the evolutionary genetics studies in nonhuman primates and in human; performed the statistical analyses; participated in the conception and the design of the study and in the drafting of the manuscript. DCS and PNB carried out the 3-D modeling analysis and participated in the drafting of the manuscript. REB and AB participated in the genetic study of nonhuman primates. PR and ARS participated in the sequence alignments and analyses. AL participated in the statistical analyses. PC participated in the design of the study and in the drafting of the manuscript. PP participated in the design of the study, in the phylogenetic analysis and in the drafting of the manuscript. PS conceived of the study and participated in its design and coordination and in the drafting of the manuscript. All authors read and approved the final manuscript.

### Additional material

#### Additional file 1

List of the ENSEMBL gene accession numbers. The sequences corresponding to the gene accession numbers were used to construct the phylogenetic tree depicted in Fig. 1.
Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2156-8-72-S1.doc]

#### Additional file 2

Target-template alignment used for modelling. This figure represents the output file for optimal sequence alignment between target sushi 1 of SRPX2 and template CR2 sushi module 1.
Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2156-8-72-S2.doc]

#### Additional file 3

Analysis of positive selection along the human lineage. The data represent the statistical analysis that was used to test for positive selection, using the branch models and the branch-site models.
Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2156-8-72-S3.doc]
Additional file 4

List of PCR primers. The 'Exons' PCR primers were used to amplify exons of the full-length coding region of SRPX2 in 12 human females and in six nonhuman primate species. The 'Introns' PCR primers were used to amplify human genomic fragments from introns 3 and 4 of SRPX2 in the same 12 women.

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2156-8-72-S4.doc]

Acknowledgements

We thank all the control individuals who participated in this study. Each participant gave informed consent prior to the study, according to the appropriate local ethical committee. We thank Ziheng Yang for his help with the use of the codeml program from the PAML package, Nicolas Galier for helpful discussion and comments and for critical reading of the manuscript, Donna Devine for help in editing the manuscript, and Gaby Doxiadis for the sending of nonhuman primate DNAs. Assistance from the "Genome Variation" core facilities (Marseille-Nice Génopôle) and from the Marseille Biological Resource Centre (BRC) was greatly appreciated. This study was supported by INSERM. BR is a recipient of a French MRT (Ministry of Research and Technology) PhD fellowship and PR has been a recipient of a LFCE (Ligue Française Contre l'Epilepsie) fellowship.

References

1. Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejarıwal A, Todd MA, Tanenbaum DM, Civello D, Lu F, Murphy B, Ferriera S, Wang G, Zheng X, White TJ, Sninsky JJ, Adams MD, Cargill M: **Inferring non-neutral evolution from human-chimp-mouse orthologous gene trios.** *Science* 2003, **302**:1960-1963.

2. Nielsen R, Bustamante C, Clark AG, Glanowski S, Sackton TB, Hubisz MJ, Fledel-Alon A, Tanenbaum DM, Civello D, White TJ, J Sninsky J, Adams MD, Cargill M: **A scan for positively selected genes in the genomes of humans and chimpanzees.** *PLoS Biol* 2005, **3**:e170.

3. Enard W, Khaitovich P, Klose J, Zollner S, Heissig F, Giavalisco P, Nieselt-Struwe K, Muchmore E, Varki A, Ravid R, Doxiadis GM, BonROP RE, Paabo S: **Intra- and interspecific variation in primate gene expression patterns.** *Science* 2002, **296**:340-343.

4. Caceres M, Lachuer J, Zapala MA, Redmond JC, Kudo L, Geschwind DH, Lockhart DJ, Preuss TM, Barlow C: **Elevated gene expression levels distinguish human from non-human primate brains.** *Proc Natl Acad Sci USA* 2003, **100**:13030-13035.

5. Khaitovich P, Muetzel B, She X, Lachmann M, Hellmann I, Dietzsch J, Steigele S, Do HH, Weiss G, Enard W, Heissig F, Arendt T, Nieselt-Struwe K, Eichler EE, Paabo S: **Regional patterns of gene expression in human and chimpanzee brains.** *Genome Res* 2004, **14**:1462-1473.

6. **Zhang J: Evolution of the human ASPM gene, a major determinant of brain size.** *Genetics* 2003, **165**:2063-2070.

7. Mekel-Bobrov N, Gilbert SL, Evans PD, Vallender EJ, Anderson JR, Hudson RR, Tishkoff SA, Lahn BT: **Ongoing adaptive evolution of ASPM, a brain size determinant in Homo sapiens.** *Science* 2005, **309**:1720-1722.

8. Evans PD, Anderson JR, Vallender EJ, Choi SS, Lahn BT: **Reconstructing the evolutionary history of microcephalin, a gene controlling human brain size.** *Hum Mol Genet* 2004, **13**:1139-1145.

9. Wang YQ, Su B: **Molecular evolution of microcephalin, a gene determining human brain size.** *Hum Mol Genet* 2004, **13**:1131-1137.

10. Evans PD, Gilbert SL, Mekel-Bobrov N, Vallender EJ, Anderson JR, Vaez-Azizi LM, Tishkoff SA, Hudson RR, Lahn BT: **Microcephalin, a genes regulating brain size, continues to evolve adaptively in humans.** *Science* 2005, **309**:1717-1720.

11. **Burki F, Kaessmann H:** **Birth and adaptive evolution of a hominoid gene that supports high neurotransmitter flux.** *Nat Genet* 2004, **36**:1061-1063.

12. Gilad Y, Rosenberg S, Przeworski M, Lancet D, Skorecki K: **Evidence for positive selection and population structure at the human MAO-A gene.** *Proc Natl Acad Sci USA* 2002, **99**:862-867.

13. Andres AM, Soldevila M, Navarro A, Kidd KK, Oliva B, Bertranpetit J: **Positive selection in MAOA gene is human exclusive: determination of the putative amino acid change selected in the human lineage.** *Hum Genet* 2004, **115**:377-386.

14. Dorus S, Anderson JR, Vallender EJ, Gilbert SL, Zhang L, Chemnick LG, Ryder OA, Li W, Lahn BT: **Sonic Hedgehog, a key development gene, experienced intensified molecular evolution in primates.** *Hum Mol Genet* 2006, **15**:2031-2037.

15. Enard W, Przeworski M, Fisher SF, Lai CS, Wiebe V, Kitano T, Monaco AP, Paabo S: **Molecular evolution of FOXP2, a gene involved in speech and language.** *Nature* 2002, **418**:869-872.

16. Zhang J, Webb DM, Podlaha O: **Accelerated protein evolution and origins of human-specific features: Foxp2 as an example.** *Genetics* 2002, **162**:1825-1835.

17. Fisher SE, Marcus GF: **The eloquent ape: genes, brains and the evolution of language.** *Nat Rev Genet* 2006, **7**:9-20.

18. Prabhakar S, Noonan JP, Paabo S, Rubin EM: **Accelerated evolution of conserved noncoding sequences in humans.** *Science* 2006, **314**:786.

19. Pollard KS, Salama SR, Lambert N, Lambot MA, Coppens S, Pedersen JS, Katzman S, King B, Onodera C, Siepel A, Kern AD, Dehay C, Igel H, Ares M Jr, Vanderhaeghen P, Haussler D: **An RNA gene expressed during cortical development evolved rapidly in humans.** *Nature* 2006, **443**:167-172.

20. Roll P, Rudolf G, Pereira S, Royer B, Scheffer IE, Massacrier A, Valenti MP, Roeckel-Trevisiol N, Jamali S, Beclin C, Seegmuller C, Metz-Lutz MN, Lemainque A, Delepine M, Caloustian C, de Saint Martin A, Brunea N, Depetris D, Mattei MG, Flori E, Robaglia-Schlupp A, Levy N, Neubauer BA, Ravid R, Marescaux C, Berkovic SF, Hirsch E, Lathrop M, Cau P, Szepetowski P: **SRPX2 mutations in disorders of language cortex and cognition.** *Hum Mol Genet* 2006, **15**:1195-1207.

21. **The Online Mendelian Inheritance in Man database** [http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=OMIM]

22. Jenkins HT, Mark L, Ball G, Persson J, Lindahl G, Uhrin D, Blom AM, Barlow PN: **Human C4b-binding protein – structural basis for interaction with streptococcal M protein, a major bacterial virulence factor.** *J Biol Chem* 2006, **281**:3690-3697.

23. **The Genome Browser Database** [http://genome.ucsc.edu/cgi-bin/hgGateway]

24. Kirkitadze MD, Barlow PN: **Structure and flexibility of the multiple domain proteins that regulate complement activation.** *Immunol Rev* 2001, **180**:146-161.

25. Oleszewski M, Gutwein P, von der Lieth W, Rauch U, Altevogt P: **Characterization of the L1-neurocan-binding site Implications for L1-L1 homophilic binding.** *J Biol Chem* 2000, **275**:34478-34485.

26. Soares DC, Gerloff DL, Syme NR, Coulson AF, Parkinson J, Barlow PN: **Large-scale modelling as a route to multiple surface comparisons of the CCP module family.** *Protein Eng Des Sel* 2005, **18**:379-388.

27. Soares DC, Barlow PN: **Structural Biology of the Complement System Edited by: Morikis D, Lambris JD. Boca Raton: CRC Press, Taylor & Francis Group; 2005:19-62.**

28. Richards A, Kemp EJ, Liszewski MK, Goodship JA, Lampe AK, Decorte R, Muslumanoglu MH, Kavukcu S, Filler G, Pirson Y, Wen LS, Atkinson JP, Goodship TH: **Mutations in human complement regulator, membrane cofactor protein (CD46), predispose to development of familial hemolytic uremic syndrome.** *Proc Natl Acad Sci* 2003, **100**:12966-12971.

29. Herbert AP, Uhrin D, Lyon M, Pangburn MK, Barlow PN: **Disease-associated sequence variations congregate in a polyanion-recognition patch on human factor H revealed in 3D structure.** *J Biol Chem* 2006, **281**:16512-16520.

30. Herbert A, Soares DC, Pangburn MK, Barlow PN: **Disease-associated sequence variants in factor H: a structural biology approach.** *Adv Exp Med Biol* 2006, **586**:313-327.

31. Wyckoff GJ, Wang W, Wu CI: **Rapid evolution of male reproductive genes in the descent of man.** *Nature* 2000, **403**:304-309.

32. Higuchi M, Kazazian HH Jr, Kasch L, Warren TC, McGinniss MJ, Phil-lips JA 3rd, Kasper C, Janco R, Antonarakis SE: **Molecular characterization of severe hemophilia A suggests that about half the mutations are not within the coding regions and splice