<title>Background</title> <p>Gene duplication has long been considered as one of the main mechanisms of the adaptive evolution and as an important source of the genetic novelty [##UREF##0##1##]. Differential duplications and deletions of chromosomal regions including coding genes provide a powerful source for the evolution of species-specific biological differences [##REF##12368483##2##]. Compared to other mammals, the genomes of primates show an enrichment of large segmental duplications with high levels (&gt;90%) of sequence identity [##REF##16770338##3##]. In the human genome particularly pronounced expansions of the copy number have been reported for genes involved in the structure and function of the brain [##REF##15252450##4##]. In comparison of human and its closest relative chimpanzee, large duplications contribute considerably (2.7%; [##REF##16136132##5##]) to the overall divergence compared to single base pair substitutions (1.2–1.5%; [##REF##12368483##2##,##REF##12672966##6##, ####REF##12799463##7##, ##REF##11992255##8##, ##REF##16136131##9##, ##REF##15164055##10##, ##REF##17075697##11##, ##REF##11948215##12####11948215##12##]). In addition to providing the substrate for non-allelic homologous recombination mediating genomic disorders (reviewed by [##REF##11818139##13##]), the duplication architecture of a genome may also influence normal phenotypic variation. It has been estimated that ~20% of segmental duplications are polymorphic within human and chimpanzee populations contributing to intraspecies diversity [##REF##16136132##5##,##REF##15918152##14##]. Despite the fact that segmental duplications cover a substantial fraction of the great apes genomes, the experimental data on the divergence and detailed evolutionary dynamics of duplicated gene regions is still limited. Sequence comparison of duplicated genes in sister-species would assist in understanding the mechanisms behind primate evolution and in associating the genetic divergence with phenotypic diversification.</p> <p>One of the genomic regions that has evolved through several gene duplication events in primate lineage is the Luteinizing Hormone Beta (<italic>LHB</italic>)/Chorionic Gonadotropin Beta (<italic>CGB</italic>) gene cluster locating in human at 19q13.32. The <italic>LHB/CGB </italic>genes have an essential role in reproduction: placentally expressed HCG hormone contributes to the implantation process of the embryo during the early stages of pregnancy, pituitary expressed luteinizing hormone promotes the ovulation and luteinization of follicles and stimulates the steroidogenesis. In human, the cluster consists of seven highly homologous genes: an ancestral <italic>LHB </italic>and six duplicated <italic>CGB </italic>genes [##REF##9582510##15##]. The data from other primates supports the hypothesis of several sequential duplication events increasing gradually the number of <italic>CGB </italic>genes among primates from one (New World monkeys: the owl monkey, <italic>Aotus trivirgatus </italic>and the dusky titi monkey, <italic>Callicebus moloch</italic>) to six in human (Figure ##FIG##0##1A##). The mapped copy number of the <italic>CGB </italic>gene among Old World monkeys varies: three in rhesus macaque (<italic>Macaca mulatta</italic>), five in guereza monkey (<italic>Colobus guereza</italic>) and dusky leaf monkey (<italic>Presbytis obscura</italic>), four in orangutan (<italic>Pongo pygmaeus</italic>) [##REF##11861891##16##]. It has been suggested that <italic>CGB </italic>gene first arose in the common ancestor of the anthropoid primates after diverging from tarsiers [##REF##11861891##16##].</p> <p>We have chosen <italic>LHB/CGB </italic>genomic region as a model to study the evolution of recent primate duplications. Although the chimpanzee genome has been sequenced, there are still large gaps and uncertainties concerning the segmental duplications regions, including the <italic>LHB/CGB </italic>region. To obtain a high-quality DNA sequence, we constructed and sequenced a shot-gun library, and hereby report the complete sequence of the entire <italic>LHB/CGB </italic>cluster in the common chimpanzee. We are addressing the following aspects regarding to the evolution of duplicated genes in closely related species: (i) in-depth comparison of the human and chimpanzee <italic>LHB/CGB </italic>genome clusters; (ii) variation in substitution rates; (iii) genic and intergenic divergences; (iv) impact of intraspecies gene conversion in phylogeny, divergence and transversion/transition ratios; (v) evidence of natural selection. To our knowledge, this is the first detailed report of parallel independent duplication events initiated within a duplicated genome cluster and leading to structural divergence in the two sister-species, human and common chimpanzee.</p>
<title>Methods</title> <title>Bacterial Artificial Chromosome (BAC) screening and shotgun library construction</title> <p>BAC library of common chimpanzee (<italic>Pan troglodytes</italic>) RPCI-43 was obtained from BACPAC Resource Center at the Children's Hospital Oakland Research Institute (Oakland, CA). In order to identify BAC clones containing <italic>LHB/CGB </italic>genome cluster we used recommended protocols and performed hybridization screening with a PCR-product containing <italic>LHB</italic>-specific sequence amplified from chimpanzee genomic DNA. Probe DNA was labeled with [<sup>32</sup>P]dCTP by random primer extension using DecaLabel™ DNA Labeling Kit (MBI Fermentas, Vilnius, Lithuania). BAC DNA was isolated using NucleoBond<sup>® </sup>BAC 100 plasmid purification kit (Macherey-Nagel GmbH &amp; Co. KG, Düren, Germany). Two overlapping BAC clones (68P2 and 109B10) containing <italic>LHB/CGB </italic>genome cluster were sheared by nebulization to approximately 5 kb long fragments and used for shotgun library construction with TOPO<sup>® </sup>Shotgun Subcloning Kit (Invitrogen, Carlsbad, CA) according to manufacturers' instructions.</p> <title>DNA sequencing and data analysis</title> <p>Plasmid DNA was purified with NucleoSpin<sup>®</sup>-Plasmid kit (Macherey-Nagel GmbH &amp; Co. KG) and sequenced on ABI 3730 × l sequencer using BigDye<sup>® </sup>Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). Plasmid ends were sequenced using M13F and M13R primers, additional primers for primer walking were designed with the web-based version of the Primer3 software [##UREF##3##84##]. Sequencing primers are available in Additional file ##SUPPL##5##6##. <italic>LHB/CGB </italic>genome cluster from two chimpanzee BAC clones were sequenced with an average redundancy of 7×, which was sufficient for assembly. Sequences were assembled using ContigExpress program from Vector NTI Suite 9 (Invitrogen) and the chimpanzee sequence was compared to human <italic>LHB/CGB </italic>genome cluster (GenBank: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>). The full Chimpanzee <italic>LHB/CGB </italic>cluster sequence has been deposited to GenBank (accession number <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>). Sequence alignments were performed and homologies determined by the web-based ClustalW [##UREF##4##85##] and Stretcher implemented in the EMBOSS package [##UREF##5##86##]. The aligned sequences of the major transcripts of the chimp and human <italic>LHB/CGB </italic>genes are given in Additional file ##SUPPL##6##7##. Substitution and indel divergences were calculated as the percentage of the number of substitutions and the number of nucleotides in indels divided with the total number of aligned nucleotides in the specific genomic region. Phylogenetic trees were constructed by MEGA3.1 [##REF##15260895##87##] using Kimura's two parameter model to infer the neighbor-joining and the branch-and-bound algorithms to find maximum parsimony trees with 1000 replications for bootstrapping.</p> <p>For coding regions maximum likelihood method [##REF##7968486##39##] was used to estimate non-synonymous/synonymous rate ratio ω (= d<sub>n</sub>/d<sub>s</sub>) by CODEML implemented in PAML package version 4 [##REF##9367129##88##,##REF##17483113##89##]. Codon frequencies were estimated from the dataset using the F3 × 4 option, other settings were as default. The simplest model M0 or one-ratio model was used to estimate the ω (an average over all the sites). As for <italic>CGB1 </italic>alternative reading frames have been predicted and no functional protein has been characterized so far. We defined <italic>CGB1 </italic>mRNA sequence and subsequently the predicted reading frame as supported by the published experimental data [##REF##1371113##18##,##REF##8954017##90##]. In parallel, the number of non-synonymous substitutions per non-synonymous site (d<sub>n</sub>) and synonymous substitutions per synonymous site (d<sub>s</sub>) were estimated using an alternative method – the Li93 method [##REF##8433381##40##,##REF##8487630##41##]. The significance of the difference between d<sub>n </sub>and d<sub>s </sub>was examined by a two-tailed Z-test [##UREF##1##42##] using MEGA3.1. To address which segments of the genes are evolving more rapidly, we performed CRANN analysis [##REF##12468084##47##,##REF##15593409##48##] using sliding and overlapping windows and for the results in visual documentation of rate heterogeneities of d<sub>n </sub>and d<sub>s</sub>, Window size was set on 20 and shift size on 10 codons using the Li93 method.</p> <p>Repetitive elements were detected by the REPEATMASKER program [##UREF##6##91##].</p>
<title>Results and Discussion</title> <title>Human and chimpanzee <italic>LHB/CGB </italic>genome clusters differ considerably in size</title> <p>The total length of the sequenced chimpanzee (<italic>Ch</italic>) <italic>LHB/CGB </italic>genomic region obtained from two overlapping BAC clones was 43,945 bp. It encompasses 1,029 bp of the flanking (centromeric side) <italic>RUVBL2 </italic>gene, 39,876 bp of the entire <italic>LHB/CGB </italic>cluster and 2,986 bp of the flanking (telomeric side) <italic>NTF5 </italic>gene (Figure ##FIG##0##1B##; Genbank submission: <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>). Compared to the human (<italic>Hu) LHB/CGB </italic>region (Genbank: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>), the sequence of the <italic>ChLHB/CGB </italic>cluster is 5,289 bp shorter. The sequence characteristics of the <italic>LHB/CGB </italic>genomic regions are similar in these two species: extremely high GC-nucleotide content (57% compared to average 41% for <italic>Hu </italic>and <italic>Ch </italic>[##REF##15164055##10##]), high fraction of CpG islands (<italic>Hu </italic>6.6%, <italic>Ch </italic>6.1% compared to estimated 1–3.5% for <italic>Hu </italic>and <italic>Ch </italic>[##REF##12672966##6##,##REF##11992255##8##]) and repetitive sequences (<italic>Hu </italic>26.9%, <italic>Ch </italic>25.15%), especially SINEs (<italic>Hu </italic>23.23%, <italic>Ch </italic>21.81%) (Additional file ##SUPPL##0##1##). High repeat content is also characteristic to several other duplicated regions, such as <italic>MHC </italic>class I region [##REF##12799463##7##] and <italic>Apolipoprotein CI </italic>genomic segment [##REF##10795830##17##].</p> <p>As expected, there is a considerable similarity between the genomic organization of human and chimpanzee <italic>LHB/CGB </italic>clusters (Figure ##FIG##0##1B##). We identified two highly identical, apparently orthologous segments within the cluster: <italic>RUVBL2/LHB/intergenic region A </italic>(Ch 8,084 bp, Hu 7,973 bp; 96% sequence identity) and the region spanning from <italic>CGB1 </italic>to <italic>NTF5 </italic>(Ch 29,136 bp, Hu 28,568 bp; 94.8% sequence identity). However, a large species-specific structural rearrangement was localized between the <italic>intergenic region A </italic>and <italic>CGB1 </italic>gene, resulting in discordant size of human (45,165 bp) and chimpanzee (39,876 bp) clusters as well as species-specific number of duplicated gene copies, seven for human (1 <italic>LHB </italic>+ 6 <italic>CGB </italic>genes) and six (1 <italic>LHB </italic>+ 5 <italic>CGB</italic>) for chimpanzee. In human the rearranged region (12,700 bp) harbors one HCG beta coding <italic>CGB </italic>gene and one <italic>CGB1/2</italic>-like gene (<italic>CGB2</italic>) recognized by a specific promoter-segment [##REF##1371113##18##,##REF##7534731##19##], while in chimp (rearranged region 6,725 bp) only a <italic>CGB1/2</italic>-like gene (<italic>CGB1B</italic>) is present in an inverted orientation compared to human. In addition, <italic>ChLHB/CGB </italic>cluster lacks the whole <italic>intergenic region C' </italic>and has a considerably shorter inverted <italic>intergenic region B' </italic>(Figure ##FIG##0##1B##).</p> <title>Evidence for independent duplication events within <italic>LHB/CGB </italic>gene cluster for two closely related primate species</title> <p>We considered alternative scenarios that may have led to structural differences in <italic>LHB/CGB </italic>clusters in two sister-species. The scenario of species-specific parallel duplications was supported by several lines of evidence (Figure ##FIG##0##1B##, Figure ##FIG##1##2##). First, it was the most parsimonious solution requiring the smallest number of rearrangement events. Assuming that the ancestral <italic>Hu-ChLHB/CGB </italic>cluster consisted of the present-day highly identical segments (Figure ##FIG##0##1B##; see above), only one evolutionary event would explain the current structure of <italic>ChLHB/CGB </italic>cluster – a direct duplication of <italic>CGB1 </italic>and most of <italic>intergenic region B </italic>(excluding <italic>psNTF6G'</italic>) giving arise to the segment <italic>CGB1B/Ch-B'</italic>. In humans, two events would have lead to present <italic>HuLHB/CGB </italic>(not in order) – an inverted duplication of the entire region from <italic>CGB1 </italic>to <italic>CGB5 </italic>gene (segment from <italic>CGB </italic>to <italic>CGB2</italic>) and a direct duplication and translocation of <italic>region C </italic>creating <italic>C' </italic>next to <italic>CGB1 </italic>gene (Figure ##FIG##0##1B##). These parallel duplications in two sister-species might have been initiated by non-allelic homologous recombinations between multiple <italic>Alu </italic>SINE sequences (<italic>Alu</italic>Sx, <italic>Alu</italic>Sp, <italic>Alu</italic>Sq) locating at the junctions of both human- and chimp-specific duplications. Consistently, the phylogenetic trees that were constructed using different fragments of intergenic regions show that <italic>Ch-B </italic>and <italic>Ch-B'</italic>, <italic>Hu-B </italic>and <italic>Hu-B' </italic>as well as <italic>Hu-C </italic>and <italic>Hu-C' </italic>clearly cluster together supporting their paralogous status (Figure ##FIG##1##2B–F##). The <italic>Hu-D </italic>and <italic>Ch-D </italic>region form a well-supported clade (Figure ##FIG##1##2B–D, G##) giving evidence that these are orthologous segments. The phylogenetic relationship supports the ancestral status of <italic>region C </italic>(present in both species) compared to <italic>Hu-C' </italic>(human-specific duplication). Also the sequence similarity calculations are consistent with the scenario of independent species-specific duplications: the identity of <italic>Ch-B </italic>and <italic>Ch-B' </italic>regions is 98.1%, and of <italic>ChCGB1 </italic>and <italic>ChCGB1B </italic>98.7%. The latter exceeds the estimates for any orthologous genes between human and chimp in <italic>LHB/CGB </italic>cluster (from 98.2% in <italic>LHB </italic>to 97.4% in <italic>CGB5 </italic>and <italic>CGB8</italic>). However, the possibility of more intensive gene conversion between these intergenic segments resulting in higher intraspecies similarity cannot be excluded. It is well known that past intraspecies gene conversion events might be reflected by tree phylogenies and could lead to erroneous conclusions [##REF##11861891##16##,##REF##16216919##20##]. The footprint of gene conversion is also reflected on the phylogenic tree of human and chimpanzee <italic>LHB/CGB </italic>genes (Figure ##FIG##1##2A##). Instead of two separate clades for orthologous <italic>CGB5 </italic>and <italic>CGB8</italic>, the genes within one species cluster together.</p> <p>Alternative scenarios leading to discordant gene number in <italic>LHB/CGB </italic>gene clusters in these two species are less supported. A minimum of three rearrangement events (Figure ##FIG##0##1B##) would have been required for the chimpanzee-specific deletion: loss of <italic>CGB </italic>and <italic>psNTF6A </italic>gene accompanied with the inversion of <italic>CGB2 </italic>and <italic>intergenic region B' </italic>(giving rise to <italic>CGB1B </italic>and chimp <italic>B'</italic>), and a separate deletion of <italic>region C' </italic>in chimp. Also the scenario of human-specific duplication (Figure ##FIG##0##1B##) would have required at least three events: an inversion of <italic>CGB1B </italic>and <italic>intergenic region Ch-B' </italic>(giving rise to part of <italic>Hu-B' </italic>and <italic>CGB2</italic>); either a direct duplication and translocation of <italic>CGB8 </italic>gene along with <italic>psNTF6A' </italic>or an inverted duplication and translocation of <italic>CGB5 </italic>along with <italic>psNTF6G' </italic>(resulting in <italic>CGB </italic>and <italic>psNTF6A</italic>), and a direct duplication and translocation of <italic>intergenic region C </italic>creating <italic>Hu-C' </italic>next to <italic>CGB1 </italic>gene.</p> <p>A number of gene families have been characterized where the gene number differs between human and chimpanzee due to species-specific indels [##REF##12799463##7##,##REF##15246530##21##]. To our knowledge, this is the first report where parallel independent duplications arisen within the same region in human and chimpanzee genomes give the best explanation for the observed structural differences between two sister-species. However, there are examples of independent duplications among primates resulting in convergent functions. A more recent duplication of <italic>X-linked opsin </italic>gene in New World howler monkeys (<italic>Alouatta seniculus </italic>and <italic>Alouatta caraya</italic>) compared to Old World primates, has lead to full trichromacy [##REF##9893841##22##, ####REF##9893843##23##, ##REF##10413401##24####10413401##24##] and also there are independently arisen functionally close genes within the <italic>Growth Hormone/Somatomammotropin </italic>genome cluster in New World monkeys and Old World monkeys/hominoid lineages [##REF##16216919##20##,##REF##12200231##25##].</p> <title>Comparative nucleotide divergence profiling of the orthologous regions reveals non-uniform substitution rates across the <italic>LHB/CGB </italic>region</title> <p>We generated a comparative nucleotide divergence profile of substitutions and indels across the orthologous genomic regions of <italic>LHB/CGB </italic>cluster (Figure ##FIG##0##1B##, Figure ##FIG##2##3##). Human (36,541 bp) and chimpanzee (37,220 bp) aligned genomic sequences were analyzed using a non-overlapping sliding window of 500 bp. Several studies have suggested that the majority of the genomic divergence between human and chimpanzee comprises of indels (3.0–11.9%) compared to contribution of nucleotide substitutions (1.2–1.5%) [##REF##12368483##2##,##REF##12672966##6##, ####REF##12799463##7##, ##REF##11992255##8##, ##REF##16136131##9##, ##REF##15164055##10##, ##REF##17075697##11####17075697##11##,##REF##11170892##26##]. In the duplicated <italic>LHB/CGB </italic>region indels (2.7%) and substitutions (2.3%) contributed approximately equally to the total divergence (5%) between the two species. In total, 61 indels (mean 16; range 1–637 bp) were identified, 26 as human and 35 as chimpanzee deletions. The size distribution of these indels was consistent with previous reports [##REF##12368483##2##,##REF##12672966##6##,##REF##16136131##9##, ####REF##15164055##10##, ##REF##17075697##11####17075697##11##] revealing an excess of short indels: 44% involved a single basepair, 77% 1–5 bp and 96.7% &lt;100 bp (Additional files ##SUPPL##0##1## and ##SUPPL##1##2##). Two large indels of &gt;100 bp co-locate with repetitive elements: the 128 bp long indel in <italic>intergenic region A </italic>is located in a simple-repeat rich region; the 637 bp indel in <italic>intergenic region B </italic>is flanked and composed of SINE sequences (Figure ##FIG##0##1B##, Figure ##FIG##2##3##). The latter indel could be defined as a recent human-specific sequence loss since the duplicated human-specific <italic>intergenic region B' </italic>lacks this deletion. High proportion (&gt;50%) of indels identified between human and chimpanzee has been shown to contain repetitive elements [##REF##12672966##6##,##REF##12799463##7##]. The regional content of repeats was also correlated with the non-uniform distribution of the nucleotide identity across the cluster (Figure ##FIG##2##3##; Additional file ##SUPPL##0##1##) consistent with the data that repetitive sequences, e.g. <italic>Alu </italic>repeats have a higher rate of base substitutions [##REF##2175877##27##].</p> <title>Transition to transversion ratio in duplicated regions differs from the estimations for unique genomic sequences: possible role of gene conversion</title> <p>We investigated the distribution of nucleotide substitutions within orthologous regions of <italic>LHB/CGB </italic>cluster in more detail (Figure ##FIG##3##4##). Transitions (C⇔T, A⇔G) and transversions (T, C⇔A, G) were found to contribute 62% and 38% of the total substitutions in five orthologous <italic>LHB/CGB </italic>genes, respectively (Figure ##FIG##3##4A##). The corresponding estimates for the whole orthologous region (including genes) were 57% for transitions and 43% for transversions (Figure ##FIG##3##4B##). Notably, the contribution of transitions is ~10% lower than reported in previous studies comparing human and chimpanzee genomic regions (68.87% – 70.3%) [##REF##12799463##7##,##REF##11992255##8##]. The most frequent transversions are G⇔C substitutions, contributing 16% of all substitutions in genes and 18% across entire orthologous region, exceeding previous estimations by twofold (9.14%, 9%) [##REF##12799463##7##,##REF##11992255##8##]. Consistently, a large excess of G⇔C transversional pairs as compared to other substitutions has been reported for human <italic>HSP70 </italic>and mouse <italic>Hsp70 </italic>orthologous duplicate genes [##REF##15084682##28##]. The high proportion of transversions in the <italic>LHB/CGB </italic>cluster can be explained by biased gene conversion [##REF##16251463##29##] that leads to a high GC content (57% in human and chimp <italic>LHB/CGB </italic>region) and thus increases the probability of G⇔C substitutions due to the altered base composition [##REF##15084682##28##,##REF##17418442##30##,##REF##12801726##31##].</p> <p>Among the transitions, we observed an excess of C⇔T substitutions in <italic>LHB/CGB </italic>genes (38%) versus the whole region (28%). It is generally accepted that a high proportion of transitions are C to T substitutions in CpG dinucleotides, exhibiting about 10 times higher mutation rate than the genomic average [##REF##16136131##9##,##REF##10978293##32##]. A higher GC content (64% vs 57%) and presence of CpG islands could explain an excess of C⇔T substitutions in <italic>LHB/CGB </italic>genes compared to intergenic regions.</p> <title>Sequence divergence between human and chimpanzee duplicated <italic>LHB/CGB </italic>genes is higher than estimates for single copy genes</title> <p>We studied the human and the chimpanzee orthologous genes <italic>LHB</italic>, <italic>CGB1</italic>, <italic>CGB5</italic>, <italic>CGB8 </italic>and <italic>CGB7 </italic>for nucleotide divergence. The nucleotide divergence ranged from 1.8% for <italic>LHB </italic>to 2.59% for <italic>CGB5 </italic>and <italic>CGB8 </italic>genomic sequences (Figure ##FIG##4##5##; Additional file ##SUPPL##2##3##). Despite that the coding regions are most conserved among the species, the exonic divergence rates (mean 1.39%; range 1 – 1.88%) exceeded many times the previous estimates for single-copy regions [##REF##11992255##8##, ####REF##16136131##9##, ##REF##15164055##10####15164055##10##,##REF##11948215##12##,##REF##11170892##26##,##REF##12826612##33##]. In non-coding regions, the average nucleotide divergence for promoter regions was as high as 3.22% (range 1 – 5.1%), for introns 2.62% (range 2.04 – 3.24%) and for 5' UTR 2.54% (range 0 – 3.83%). The comparative estimates for single-copy genes are much lower: in promoters 0.75–0.88%, in exons 0.51–1.09%, in introns 1.03–1.47% and in 5' UTR 1.00–1.41% [##REF##15164055##10##, ####REF##17075697##11##, ##REF##11948215##12####11948215##12##,##REF##11170892##26##,##REF##12826612##33##, ####REF##12766228##34##, ##REF##16136134##35##, ##REF##16432233##36####16432233##36##]. It has been suggested that interaction of selection and gene conversion contributes to a higher divergence and diversity in multigene families compared to single-copy genes [##REF##3855539##37##,##REF##1862096##38##]. As any <italic>de novo </italic>mutation has a potential to be spread by gene conversion from the original locus to other gene copies, every duplicate could accumulate substitutions arisen in neighboring genes.</p> <p>Notably, when the polymorphism data from re-sequencing studies of the human (n = 95 individuals) [##REF##16251463##29##] and the chimpanzee (n = 11 individuals) [Hallast et al, unpublished] <italic>LHB/CGB </italic>genes were incorporated into calculations, the sequence divergence in genes dropped from 1.8% to 1.26% in <italic>LHB</italic>, from 2.59% to 1.84% in <italic>CGB5</italic>, from 2.59% to 1.9% in <italic>CGB8 </italic>and from 2.18% to 1.57% in <italic>CGB7 </italic>(Figure ##FIG##4##5##; Additional file ##SUPPL##2##3##). Still, a higher divergence compared to the published data for single-copy genes remained.</p> <p>Most importantly, our data indicates that the divergence estimates between human and chimpanzee might be substantially lower than reported when the intraspecies variation is taken into account.</p> <title>Duplicated, highly homologous <italic>LHB/CGB </italic>genes evolve under different selective constraints</title> <p>A number of non-synonymous changes were identified in human-chimpanzee comparison: two in <italic>LHB</italic>, two in <italic>CGB8</italic>, three in <italic>CGB5 </italic>and <italic>CGB7 </italic>and four in <italic>CGB1 </italic>(Additional file ##SUPPL##2##3##). None of the differences found in exonic regions in chimp altered the ORF nor created a preliminary stop-codon compared to human genes. The maximum likelihood (ML) method was used to estimate ω, non-synonymous (d<sub>n</sub>)/synonymous (d<sub>s</sub>) rate ratio for the orthologous genes (n = 5) in the human and chimpanzee <italic>LHB/CGB </italic>gene cluster (Table ##TAB##0##1##) [##REF##7968486##39##]. For four of the compared genes (<italic>LHB, CGB5, CGB8, CGB7</italic>) the parameter was ω &lt; 1, indicating purifying selection. An especially low ω estimate was obtained for <italic>LHB </italic>(ω = 0.087; 6 synonymous, 2 non-synonymous changes; amino acid divergence 1.42%) and <italic>CGB8 </italic>(ω = 0.099; 5 synonymous, 2 non-synonymous changes; amino acid divergence 1.21%). Consistently, the d<sub>n</sub>/d<sub>s</sub>ratio calculated by an alternative Li93 method [##REF##8433381##40##,##REF##8487630##41##] was statistically significant (d<sub>n</sub>/d<sub>s </sub>= 0.146; p = 0.049 two-tailed Z-test [##UREF##1##42##]) for <italic>LHB </italic>and reached borderline significance for the <italic>CGB8 </italic>(d<sub>n</sub>/d<sub>s </sub>= 0.162; p = 0.078) (Additional file ##SUPPL##3##4##). Indeed, the <italic>LHB </italic>gene is highly conserved throughout vertebrate evolution in association with its function – it is coding for the receptor-binding beta subunit of luteinizing hormone (LH) critical for successful reproduction [##REF##9582510##15##]. The functional constraint of the human <italic>CGB8 </italic>may be associated with its major role in contributing to the <italic>HCG beta </italic>subunit mRNA transcript pool compared to other <italic>CGB </italic>genes [##REF##16123088##43##].</p> <p><italic>CGB1 </italic>gene (ω = 2.658, 1 synonymous, 4 non-synonymous changes; amino acid divergence 3.03%) stood out as the only locus in the gene cluster with estimated ω &gt;1, which would be consistent with positive or adaptive evolution. It has been suggested that <italic>CGB1 </italic>has arisen in the common ancestor of African great apes through a duplication event accompanied by an insertion of a novel putative promoter, 5'UTR and exon 1. So far, the detection of <italic>CGB1 </italic>gene has been unsuccessful in orangutan <italic>Pongo pygmaeus </italic>by using the PCR approach [##REF##17055150##44##] and in macaque <italic>Macaca mulatta </italic>(Genbank: <ext-link ext-link-type="gen" xlink:href="AC202849">AC202849</ext-link>) by <italic>in silico </italic>search of the current genome assembly. It has been shown that in human the contributions of <italic>CGB1 </italic>and its duplicate human-specific <italic>CGB2 </italic>to the summarized expression of the six <italic>CGB </italic>genes in placenta is much lower (1/1000 to 1/10000) compared to their gene dosage (two genes out of six total) [##REF##16123088##43##,##REF##18048458##45##]. However, in testis the proportional contribution of <italic>CGB1/2 </italic>to the total <italic>CGB </italic>transcript pool is as high as 1/3 [##REF##18048458##45##], which may indicate a possible role of these genes in male reproductive tract. Indeed, a recent study has shown that HCG alpha and HCG beta free subunits are produced in high amounts in the prostate and testes and are subsequently observed in seminal plasma [##REF##17097221##46##].</p> <p>In order to address which parts of the studied genes exhibit signals of evolving under positive selection, we used CRANN analysis calculating d<sub>n </sub>and d<sub>s </sub>values for sliding and overlapping windows along individual genes (Additional file ##SUPPL##4##5##) [##REF##12468084##47##,##REF##15593409##48##]. In case of <italic>HCG beta </italic>subunit coding genes (<italic>CGB5</italic>, <italic>CGB8 </italic>and <italic>CGB7</italic>), the patterns of nucleotide differences between the sister species were similar – across the protein the synonymous substitutions exceeded the non-synonymous ones and were concentrated in the N- and C-terminus of the protein. In contrast, in <italic>CGB1 </italic>the non-synonymous substitutions dominated and were distributed in the signal peptide (amino acids 1–20) and the centre of the protein.</p> <title>Impact of gene duplications in primate divergence</title> <p>Primate-specific gene duplications have involved loci regulating immunity (e.g. <italic>MHC, beta-defensin</italic>, <italic>CD33rSiglec </italic>gene clusters), reproduction (e.g. <italic>LHB/CGB</italic>, <italic>GH/CSH</italic>, <italic>PRAME genes; Y-chromosomal gene families</italic>), development and adaptation (e.g. <italic>Beta Globin, Opsin</italic>, <italic>Rh blood group</italic>, <italic>Class 1 ADH</italic>, <italic>PRDM </italic>and <italic>FAM90A </italic>gene families), and brain functions (e.g. <italic>NAIP, ROCK1, USP10 and MGC8902 </italic>genes) [##REF##15252450##4##,##REF##11861891##16##,##REF##16216919##20##,##REF##9893841##22##,##REF##17666543##49##, ####REF##12679545##50##, ##REF##12734011##51##, ##REF##15331780##52##, ##REF##16159394##53##, ##REF##1908094##54##, ##REF##8105088##55##, ##REF##17204375##56##, ##REF##17916234##57##, ##REF##17684299##58##, ##REF##16946073##59####16946073##59##]. It has been suggested that these duplication events may have been facilitated by non-allelic homologous recombination between <italic>Alu </italic>sequences, expanded into millions of copies all over primate genomes [##REF##1908094##54##,##REF##14505274##60##]. There are several examples of independent duplication events in the New World monkey (NWM) and Old World monkey (OWM)/hominoid lineages as well as in distinct primate species. For example, the OWM and apes have three <italic>Opsin </italic>genes and are trichromats due to gene duplication at the base of the OWM lineage. In NWM, the situation is more variable: most species exhibit two <italic>Opsin </italic>genes, but in the howler monkey an additional gene duplication has led to full trichromacy [##REF##9893843##23##,##REF##10413401##24##]. In <italic>Growth Hormone </italic>gene cluster (five to eight gene copies) some of the duplicate genes in the OWM/hominoid lineages have acquired a novel function and code for Chorionic Somatomammotropin (<italic>CSH </italic>genes) involved in the glucose metabolism of the fetus and the mother. However, the CSH genes are missing from the genomes of NW monkeys [##REF##2744760##61##]. Further species-specific duplications of <italic>GH/CSH </italic>genes have been reported for gibbon, macaque, chimpanzee and human [##REF##15848116##62##,##REF##8404617##63##]. In human-chimpanzee comparison, only three <italic>GH/CSH </italic>genes are clearly orthologous [##REF##15246530##21##]. Other gene clusters with independent gene duplications in human, apes and macaque lineages are <italic>MHC </italic>and testes-expressed <italic>PRAME </italic>genes [##REF##16159394##53##,##REF##17431167##64##]. The ancestral <italic>MHC-B </italic>duplicated into <italic>MHC-B </italic>and <italic>MHC-C </italic>in hominoids. While human <italic>MHC-C </italic>orthologs are found in African apes and orangutans, they are not present in macaque or any other OW monkeys [##REF##8648109##65##,##REF##1560209##66##]. Species-specific evolutionary scenarios have also been reported for <italic>LHB/CGB </italic>genes [##REF##11861891##16##]. In addition to the structural differences between human and chimpanzee <italic>LHB/CGB genes </italic>reported in this study, an expansion of <italic>CGB </italic>genes up to 50 gene copies has been shown in gorilla (<italic>Gorilla gorilla</italic>) [##REF##15252450##4##,##REF##17666543##49##]. Despite the high number of structural differences that has been shown between human and chimpanzee genomes by <italic>in silico </italic>whole-genome analysis (reviewed in [##REF##17597777##67##]), the experimental data for copy number differences between these sister-species has been reported only for a few gene clusters (e.g. <italic>GH/CSH, LHB/CGB, PRAME, MGC8902, CXYorf1, KGF, CD33rSiglec, NANOG</italic>) [##REF##11861891##16##,##REF##16216919##20##,##REF##15331780##52##,##REF##16159394##53##,##REF##16946073##59##,##REF##10655549##68##, ####REF##9326632##69##, ##REF##16469101##70####16469101##70##].</p> <p>In addition to creating structural divergence among the species, duplications provide also the bases for diversification of gene functions. For primate-specific gene duplications, there is evidence of variability in evolutionary rates among the gene copies within and among the species, and of different selective constraints acting on different members of the gene clusters, such as <italic>MHC, beta-globin</italic>, <italic>GH/CSH, PRAME</italic>, <italic>Rh blood group</italic>, <italic>CD33rSiglec </italic>and <italic>beta-defensin </italic>genes [##REF##3412472##71##, ####REF##15486692##72##, ##REF##16843621##73####16843621##73##]. For example, in human-chimpanzee comparison, the chimp <italic>MHC </italic>class I loci <italic>A, B </italic>and <italic>C </italic>are characterized by lower intra-species allelic variation compared to human, providing evidence that ancestral chimpanzee populations may have experienced a selective sweep [##REF##12186979##74##,##REF##15758205##75##]. In marmoset <italic>LHB/CGB </italic>region, a switch of functions has happened between the ancestral <italic>LHB </italic>and the derived <italic>CGB </italic>gene. Although <italic>LHB </italic>and <italic>CGB </italic>genes are both present at the genomic level, only <italic>CGB </italic>gene is expressed in the pituitary and placenta tissues. Thus, Chorionic Gonadotropin Hormone is the only gonadotropin carrying also the luteinizing function fulfilled by Luteinizing Hormone in mammals [##REF##7492691##76##,##REF##14765996##77##].</p> <p>Duplicated genes tend to evolve <italic>in consort </italic>facilitated by active inter-locus gene conversion increasing and preserving sequence similarity among the gene copies [##REF##3855539##37##,##REF##1862096##38##,##REF##17846636##78##]. Concerted evolution within species may lead to erroneous phylogenetic trees and to the overestimation of interspecies divergence dates [##REF##16216919##20##,##REF##10331249##79##,##REF##17009125##80##]. Incorporation of gene conversion data into the equations for calculating interspecies divergence may be still a challenge requiring detailed knowledge of a particular genomic region.</p> <p>So far, only a limited number of reports has been published that focus on detailed variation patterns of duplicated gene families in primates (<italic>MHC</italic>, <italic>Globin</italic>, <italic>GH/CSH </italic>and <italic>LHB/CGB </italic>genes) [##REF##12799463##7##,##REF##16251463##29##,##REF##16702430##81##, ####REF##17223997##82##, ##UREF##2##83####2##83##]. However, the common observation is that compared to single-copy segments, duplicated regions tend to exhibit higher interspecies diversity that could be explained by relaxed selective pressures and/or gene conversions spreading mutations. Thus, when calculating the divergence in duplicated primate-specific regions, the inclusion or exclusion of intraspecies variation data into equations may have a substantial impact on the divergence estimates.</p>
<title>Results and Discussion</title> <title>Human and chimpanzee <italic>LHB/CGB </italic>genome clusters differ considerably in size</title> <p>The total length of the sequenced chimpanzee (<italic>Ch</italic>) <italic>LHB/CGB </italic>genomic region obtained from two overlapping BAC clones was 43,945 bp. It encompasses 1,029 bp of the flanking (centromeric side) <italic>RUVBL2 </italic>gene, 39,876 bp of the entire <italic>LHB/CGB </italic>cluster and 2,986 bp of the flanking (telomeric side) <italic>NTF5 </italic>gene (Figure ##FIG##0##1B##; Genbank submission: <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>). Compared to the human (<italic>Hu) LHB/CGB </italic>region (Genbank: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>), the sequence of the <italic>ChLHB/CGB </italic>cluster is 5,289 bp shorter. The sequence characteristics of the <italic>LHB/CGB </italic>genomic regions are similar in these two species: extremely high GC-nucleotide content (57% compared to average 41% for <italic>Hu </italic>and <italic>Ch </italic>[##REF##15164055##10##]), high fraction of CpG islands (<italic>Hu </italic>6.6%, <italic>Ch </italic>6.1% compared to estimated 1–3.5% for <italic>Hu </italic>and <italic>Ch </italic>[##REF##12672966##6##,##REF##11992255##8##]) and repetitive sequences (<italic>Hu </italic>26.9%, <italic>Ch </italic>25.15%), especially SINEs (<italic>Hu </italic>23.23%, <italic>Ch </italic>21.81%) (Additional file ##SUPPL##0##1##). High repeat content is also characteristic to several other duplicated regions, such as <italic>MHC </italic>class I region [##REF##12799463##7##] and <italic>Apolipoprotein CI </italic>genomic segment [##REF##10795830##17##].</p> <p>As expected, there is a considerable similarity between the genomic organization of human and chimpanzee <italic>LHB/CGB </italic>clusters (Figure ##FIG##0##1B##). We identified two highly identical, apparently orthologous segments within the cluster: <italic>RUVBL2/LHB/intergenic region A </italic>(Ch 8,084 bp, Hu 7,973 bp; 96% sequence identity) and the region spanning from <italic>CGB1 </italic>to <italic>NTF5 </italic>(Ch 29,136 bp, Hu 28,568 bp; 94.8% sequence identity). However, a large species-specific structural rearrangement was localized between the <italic>intergenic region A </italic>and <italic>CGB1 </italic>gene, resulting in discordant size of human (45,165 bp) and chimpanzee (39,876 bp) clusters as well as species-specific number of duplicated gene copies, seven for human (1 <italic>LHB </italic>+ 6 <italic>CGB </italic>genes) and six (1 <italic>LHB </italic>+ 5 <italic>CGB</italic>) for chimpanzee. In human the rearranged region (12,700 bp) harbors one HCG beta coding <italic>CGB </italic>gene and one <italic>CGB1/2</italic>-like gene (<italic>CGB2</italic>) recognized by a specific promoter-segment [##REF##1371113##18##,##REF##7534731##19##], while in chimp (rearranged region 6,725 bp) only a <italic>CGB1/2</italic>-like gene (<italic>CGB1B</italic>) is present in an inverted orientation compared to human. In addition, <italic>ChLHB/CGB </italic>cluster lacks the whole <italic>intergenic region C' </italic>and has a considerably shorter inverted <italic>intergenic region B' </italic>(Figure ##FIG##0##1B##).</p> <title>Evidence for independent duplication events within <italic>LHB/CGB </italic>gene cluster for two closely related primate species</title> <p>We considered alternative scenarios that may have led to structural differences in <italic>LHB/CGB </italic>clusters in two sister-species. The scenario of species-specific parallel duplications was supported by several lines of evidence (Figure ##FIG##0##1B##, Figure ##FIG##1##2##). First, it was the most parsimonious solution requiring the smallest number of rearrangement events. Assuming that the ancestral <italic>Hu-ChLHB/CGB </italic>cluster consisted of the present-day highly identical segments (Figure ##FIG##0##1B##; see above), only one evolutionary event would explain the current structure of <italic>ChLHB/CGB </italic>cluster – a direct duplication of <italic>CGB1 </italic>and most of <italic>intergenic region B </italic>(excluding <italic>psNTF6G'</italic>) giving arise to the segment <italic>CGB1B/Ch-B'</italic>. In humans, two events would have lead to present <italic>HuLHB/CGB </italic>(not in order) – an inverted duplication of the entire region from <italic>CGB1 </italic>to <italic>CGB5 </italic>gene (segment from <italic>CGB </italic>to <italic>CGB2</italic>) and a direct duplication and translocation of <italic>region C </italic>creating <italic>C' </italic>next to <italic>CGB1 </italic>gene (Figure ##FIG##0##1B##). These parallel duplications in two sister-species might have been initiated by non-allelic homologous recombinations between multiple <italic>Alu </italic>SINE sequences (<italic>Alu</italic>Sx, <italic>Alu</italic>Sp, <italic>Alu</italic>Sq) locating at the junctions of both human- and chimp-specific duplications. Consistently, the phylogenetic trees that were constructed using different fragments of intergenic regions show that <italic>Ch-B </italic>and <italic>Ch-B'</italic>, <italic>Hu-B </italic>and <italic>Hu-B' </italic>as well as <italic>Hu-C </italic>and <italic>Hu-C' </italic>clearly cluster together supporting their paralogous status (Figure ##FIG##1##2B–F##). The <italic>Hu-D </italic>and <italic>Ch-D </italic>region form a well-supported clade (Figure ##FIG##1##2B–D, G##) giving evidence that these are orthologous segments. The phylogenetic relationship supports the ancestral status of <italic>region C </italic>(present in both species) compared to <italic>Hu-C' </italic>(human-specific duplication). Also the sequence similarity calculations are consistent with the scenario of independent species-specific duplications: the identity of <italic>Ch-B </italic>and <italic>Ch-B' </italic>regions is 98.1%, and of <italic>ChCGB1 </italic>and <italic>ChCGB1B </italic>98.7%. The latter exceeds the estimates for any orthologous genes between human and chimp in <italic>LHB/CGB </italic>cluster (from 98.2% in <italic>LHB </italic>to 97.4% in <italic>CGB5 </italic>and <italic>CGB8</italic>). However, the possibility of more intensive gene conversion between these intergenic segments resulting in higher intraspecies similarity cannot be excluded. It is well known that past intraspecies gene conversion events might be reflected by tree phylogenies and could lead to erroneous conclusions [##REF##11861891##16##,##REF##16216919##20##]. The footprint of gene conversion is also reflected on the phylogenic tree of human and chimpanzee <italic>LHB/CGB </italic>genes (Figure ##FIG##1##2A##). Instead of two separate clades for orthologous <italic>CGB5 </italic>and <italic>CGB8</italic>, the genes within one species cluster together.</p> <p>Alternative scenarios leading to discordant gene number in <italic>LHB/CGB </italic>gene clusters in these two species are less supported. A minimum of three rearrangement events (Figure ##FIG##0##1B##) would have been required for the chimpanzee-specific deletion: loss of <italic>CGB </italic>and <italic>psNTF6A </italic>gene accompanied with the inversion of <italic>CGB2 </italic>and <italic>intergenic region B' </italic>(giving rise to <italic>CGB1B </italic>and chimp <italic>B'</italic>), and a separate deletion of <italic>region C' </italic>in chimp. Also the scenario of human-specific duplication (Figure ##FIG##0##1B##) would have required at least three events: an inversion of <italic>CGB1B </italic>and <italic>intergenic region Ch-B' </italic>(giving rise to part of <italic>Hu-B' </italic>and <italic>CGB2</italic>); either a direct duplication and translocation of <italic>CGB8 </italic>gene along with <italic>psNTF6A' </italic>or an inverted duplication and translocation of <italic>CGB5 </italic>along with <italic>psNTF6G' </italic>(resulting in <italic>CGB </italic>and <italic>psNTF6A</italic>), and a direct duplication and translocation of <italic>intergenic region C </italic>creating <italic>Hu-C' </italic>next to <italic>CGB1 </italic>gene.</p> <p>A number of gene families have been characterized where the gene number differs between human and chimpanzee due to species-specific indels [##REF##12799463##7##,##REF##15246530##21##]. To our knowledge, this is the first report where parallel independent duplications arisen within the same region in human and chimpanzee genomes give the best explanation for the observed structural differences between two sister-species. However, there are examples of independent duplications among primates resulting in convergent functions. A more recent duplication of <italic>X-linked opsin </italic>gene in New World howler monkeys (<italic>Alouatta seniculus </italic>and <italic>Alouatta caraya</italic>) compared to Old World primates, has lead to full trichromacy [##REF##9893841##22##, ####REF##9893843##23##, ##REF##10413401##24####10413401##24##] and also there are independently arisen functionally close genes within the <italic>Growth Hormone/Somatomammotropin </italic>genome cluster in New World monkeys and Old World monkeys/hominoid lineages [##REF##16216919##20##,##REF##12200231##25##].</p> <title>Comparative nucleotide divergence profiling of the orthologous regions reveals non-uniform substitution rates across the <italic>LHB/CGB </italic>region</title> <p>We generated a comparative nucleotide divergence profile of substitutions and indels across the orthologous genomic regions of <italic>LHB/CGB </italic>cluster (Figure ##FIG##0##1B##, Figure ##FIG##2##3##). Human (36,541 bp) and chimpanzee (37,220 bp) aligned genomic sequences were analyzed using a non-overlapping sliding window of 500 bp. Several studies have suggested that the majority of the genomic divergence between human and chimpanzee comprises of indels (3.0–11.9%) compared to contribution of nucleotide substitutions (1.2–1.5%) [##REF##12368483##2##,##REF##12672966##6##, ####REF##12799463##7##, ##REF##11992255##8##, ##REF##16136131##9##, ##REF##15164055##10##, ##REF##17075697##11####17075697##11##,##REF##11170892##26##]. In the duplicated <italic>LHB/CGB </italic>region indels (2.7%) and substitutions (2.3%) contributed approximately equally to the total divergence (5%) between the two species. In total, 61 indels (mean 16; range 1–637 bp) were identified, 26 as human and 35 as chimpanzee deletions. The size distribution of these indels was consistent with previous reports [##REF##12368483##2##,##REF##12672966##6##,##REF##16136131##9##, ####REF##15164055##10##, ##REF##17075697##11####17075697##11##] revealing an excess of short indels: 44% involved a single basepair, 77% 1–5 bp and 96.7% &lt;100 bp (Additional files ##SUPPL##0##1## and ##SUPPL##1##2##). Two large indels of &gt;100 bp co-locate with repetitive elements: the 128 bp long indel in <italic>intergenic region A </italic>is located in a simple-repeat rich region; the 637 bp indel in <italic>intergenic region B </italic>is flanked and composed of SINE sequences (Figure ##FIG##0##1B##, Figure ##FIG##2##3##). The latter indel could be defined as a recent human-specific sequence loss since the duplicated human-specific <italic>intergenic region B' </italic>lacks this deletion. High proportion (&gt;50%) of indels identified between human and chimpanzee has been shown to contain repetitive elements [##REF##12672966##6##,##REF##12799463##7##]. The regional content of repeats was also correlated with the non-uniform distribution of the nucleotide identity across the cluster (Figure ##FIG##2##3##; Additional file ##SUPPL##0##1##) consistent with the data that repetitive sequences, e.g. <italic>Alu </italic>repeats have a higher rate of base substitutions [##REF##2175877##27##].</p> <title>Transition to transversion ratio in duplicated regions differs from the estimations for unique genomic sequences: possible role of gene conversion</title> <p>We investigated the distribution of nucleotide substitutions within orthologous regions of <italic>LHB/CGB </italic>cluster in more detail (Figure ##FIG##3##4##). Transitions (C⇔T, A⇔G) and transversions (T, C⇔A, G) were found to contribute 62% and 38% of the total substitutions in five orthologous <italic>LHB/CGB </italic>genes, respectively (Figure ##FIG##3##4A##). The corresponding estimates for the whole orthologous region (including genes) were 57% for transitions and 43% for transversions (Figure ##FIG##3##4B##). Notably, the contribution of transitions is ~10% lower than reported in previous studies comparing human and chimpanzee genomic regions (68.87% – 70.3%) [##REF##12799463##7##,##REF##11992255##8##]. The most frequent transversions are G⇔C substitutions, contributing 16% of all substitutions in genes and 18% across entire orthologous region, exceeding previous estimations by twofold (9.14%, 9%) [##REF##12799463##7##,##REF##11992255##8##]. Consistently, a large excess of G⇔C transversional pairs as compared to other substitutions has been reported for human <italic>HSP70 </italic>and mouse <italic>Hsp70 </italic>orthologous duplicate genes [##REF##15084682##28##]. The high proportion of transversions in the <italic>LHB/CGB </italic>cluster can be explained by biased gene conversion [##REF##16251463##29##] that leads to a high GC content (57% in human and chimp <italic>LHB/CGB </italic>region) and thus increases the probability of G⇔C substitutions due to the altered base composition [##REF##15084682##28##,##REF##17418442##30##,##REF##12801726##31##].</p> <p>Among the transitions, we observed an excess of C⇔T substitutions in <italic>LHB/CGB </italic>genes (38%) versus the whole region (28%). It is generally accepted that a high proportion of transitions are C to T substitutions in CpG dinucleotides, exhibiting about 10 times higher mutation rate than the genomic average [##REF##16136131##9##,##REF##10978293##32##]. A higher GC content (64% vs 57%) and presence of CpG islands could explain an excess of C⇔T substitutions in <italic>LHB/CGB </italic>genes compared to intergenic regions.</p> <title>Sequence divergence between human and chimpanzee duplicated <italic>LHB/CGB </italic>genes is higher than estimates for single copy genes</title> <p>We studied the human and the chimpanzee orthologous genes <italic>LHB</italic>, <italic>CGB1</italic>, <italic>CGB5</italic>, <italic>CGB8 </italic>and <italic>CGB7 </italic>for nucleotide divergence. The nucleotide divergence ranged from 1.8% for <italic>LHB </italic>to 2.59% for <italic>CGB5 </italic>and <italic>CGB8 </italic>genomic sequences (Figure ##FIG##4##5##; Additional file ##SUPPL##2##3##). Despite that the coding regions are most conserved among the species, the exonic divergence rates (mean 1.39%; range 1 – 1.88%) exceeded many times the previous estimates for single-copy regions [##REF##11992255##8##, ####REF##16136131##9##, ##REF##15164055##10####15164055##10##,##REF##11948215##12##,##REF##11170892##26##,##REF##12826612##33##]. In non-coding regions, the average nucleotide divergence for promoter regions was as high as 3.22% (range 1 – 5.1%), for introns 2.62% (range 2.04 – 3.24%) and for 5' UTR 2.54% (range 0 – 3.83%). The comparative estimates for single-copy genes are much lower: in promoters 0.75–0.88%, in exons 0.51–1.09%, in introns 1.03–1.47% and in 5' UTR 1.00–1.41% [##REF##15164055##10##, ####REF##17075697##11##, ##REF##11948215##12####11948215##12##,##REF##11170892##26##,##REF##12826612##33##, ####REF##12766228##34##, ##REF##16136134##35##, ##REF##16432233##36####16432233##36##]. It has been suggested that interaction of selection and gene conversion contributes to a higher divergence and diversity in multigene families compared to single-copy genes [##REF##3855539##37##,##REF##1862096##38##]. As any <italic>de novo </italic>mutation has a potential to be spread by gene conversion from the original locus to other gene copies, every duplicate could accumulate substitutions arisen in neighboring genes.</p> <p>Notably, when the polymorphism data from re-sequencing studies of the human (n = 95 individuals) [##REF##16251463##29##] and the chimpanzee (n = 11 individuals) [Hallast et al, unpublished] <italic>LHB/CGB </italic>genes were incorporated into calculations, the sequence divergence in genes dropped from 1.8% to 1.26% in <italic>LHB</italic>, from 2.59% to 1.84% in <italic>CGB5</italic>, from 2.59% to 1.9% in <italic>CGB8 </italic>and from 2.18% to 1.57% in <italic>CGB7 </italic>(Figure ##FIG##4##5##; Additional file ##SUPPL##2##3##). Still, a higher divergence compared to the published data for single-copy genes remained.</p> <p>Most importantly, our data indicates that the divergence estimates between human and chimpanzee might be substantially lower than reported when the intraspecies variation is taken into account.</p> <title>Duplicated, highly homologous <italic>LHB/CGB </italic>genes evolve under different selective constraints</title> <p>A number of non-synonymous changes were identified in human-chimpanzee comparison: two in <italic>LHB</italic>, two in <italic>CGB8</italic>, three in <italic>CGB5 </italic>and <italic>CGB7 </italic>and four in <italic>CGB1 </italic>(Additional file ##SUPPL##2##3##). None of the differences found in exonic regions in chimp altered the ORF nor created a preliminary stop-codon compared to human genes. The maximum likelihood (ML) method was used to estimate ω, non-synonymous (d<sub>n</sub>)/synonymous (d<sub>s</sub>) rate ratio for the orthologous genes (n = 5) in the human and chimpanzee <italic>LHB/CGB </italic>gene cluster (Table ##TAB##0##1##) [##REF##7968486##39##]. For four of the compared genes (<italic>LHB, CGB5, CGB8, CGB7</italic>) the parameter was ω &lt; 1, indicating purifying selection. An especially low ω estimate was obtained for <italic>LHB </italic>(ω = 0.087; 6 synonymous, 2 non-synonymous changes; amino acid divergence 1.42%) and <italic>CGB8 </italic>(ω = 0.099; 5 synonymous, 2 non-synonymous changes; amino acid divergence 1.21%). Consistently, the d<sub>n</sub>/d<sub>s</sub>ratio calculated by an alternative Li93 method [##REF##8433381##40##,##REF##8487630##41##] was statistically significant (d<sub>n</sub>/d<sub>s </sub>= 0.146; p = 0.049 two-tailed Z-test [##UREF##1##42##]) for <italic>LHB </italic>and reached borderline significance for the <italic>CGB8 </italic>(d<sub>n</sub>/d<sub>s </sub>= 0.162; p = 0.078) (Additional file ##SUPPL##3##4##). Indeed, the <italic>LHB </italic>gene is highly conserved throughout vertebrate evolution in association with its function – it is coding for the receptor-binding beta subunit of luteinizing hormone (LH) critical for successful reproduction [##REF##9582510##15##]. The functional constraint of the human <italic>CGB8 </italic>may be associated with its major role in contributing to the <italic>HCG beta </italic>subunit mRNA transcript pool compared to other <italic>CGB </italic>genes [##REF##16123088##43##].</p> <p><italic>CGB1 </italic>gene (ω = 2.658, 1 synonymous, 4 non-synonymous changes; amino acid divergence 3.03%) stood out as the only locus in the gene cluster with estimated ω &gt;1, which would be consistent with positive or adaptive evolution. It has been suggested that <italic>CGB1 </italic>has arisen in the common ancestor of African great apes through a duplication event accompanied by an insertion of a novel putative promoter, 5'UTR and exon 1. So far, the detection of <italic>CGB1 </italic>gene has been unsuccessful in orangutan <italic>Pongo pygmaeus </italic>by using the PCR approach [##REF##17055150##44##] and in macaque <italic>Macaca mulatta </italic>(Genbank: <ext-link ext-link-type="gen" xlink:href="AC202849">AC202849</ext-link>) by <italic>in silico </italic>search of the current genome assembly. It has been shown that in human the contributions of <italic>CGB1 </italic>and its duplicate human-specific <italic>CGB2 </italic>to the summarized expression of the six <italic>CGB </italic>genes in placenta is much lower (1/1000 to 1/10000) compared to their gene dosage (two genes out of six total) [##REF##16123088##43##,##REF##18048458##45##]. However, in testis the proportional contribution of <italic>CGB1/2 </italic>to the total <italic>CGB </italic>transcript pool is as high as 1/3 [##REF##18048458##45##], which may indicate a possible role of these genes in male reproductive tract. Indeed, a recent study has shown that HCG alpha and HCG beta free subunits are produced in high amounts in the prostate and testes and are subsequently observed in seminal plasma [##REF##17097221##46##].</p> <p>In order to address which parts of the studied genes exhibit signals of evolving under positive selection, we used CRANN analysis calculating d<sub>n </sub>and d<sub>s </sub>values for sliding and overlapping windows along individual genes (Additional file ##SUPPL##4##5##) [##REF##12468084##47##,##REF##15593409##48##]. In case of <italic>HCG beta </italic>subunit coding genes (<italic>CGB5</italic>, <italic>CGB8 </italic>and <italic>CGB7</italic>), the patterns of nucleotide differences between the sister species were similar – across the protein the synonymous substitutions exceeded the non-synonymous ones and were concentrated in the N- and C-terminus of the protein. In contrast, in <italic>CGB1 </italic>the non-synonymous substitutions dominated and were distributed in the signal peptide (amino acids 1–20) and the centre of the protein.</p> <title>Impact of gene duplications in primate divergence</title> <p>Primate-specific gene duplications have involved loci regulating immunity (e.g. <italic>MHC, beta-defensin</italic>, <italic>CD33rSiglec </italic>gene clusters), reproduction (e.g. <italic>LHB/CGB</italic>, <italic>GH/CSH</italic>, <italic>PRAME genes; Y-chromosomal gene families</italic>), development and adaptation (e.g. <italic>Beta Globin, Opsin</italic>, <italic>Rh blood group</italic>, <italic>Class 1 ADH</italic>, <italic>PRDM </italic>and <italic>FAM90A </italic>gene families), and brain functions (e.g. <italic>NAIP, ROCK1, USP10 and MGC8902 </italic>genes) [##REF##15252450##4##,##REF##11861891##16##,##REF##16216919##20##,##REF##9893841##22##,##REF##17666543##49##, ####REF##12679545##50##, ##REF##12734011##51##, ##REF##15331780##52##, ##REF##16159394##53##, ##REF##1908094##54##, ##REF##8105088##55##, ##REF##17204375##56##, ##REF##17916234##57##, ##REF##17684299##58##, ##REF##16946073##59####16946073##59##]. It has been suggested that these duplication events may have been facilitated by non-allelic homologous recombination between <italic>Alu </italic>sequences, expanded into millions of copies all over primate genomes [##REF##1908094##54##,##REF##14505274##60##]. There are several examples of independent duplication events in the New World monkey (NWM) and Old World monkey (OWM)/hominoid lineages as well as in distinct primate species. For example, the OWM and apes have three <italic>Opsin </italic>genes and are trichromats due to gene duplication at the base of the OWM lineage. In NWM, the situation is more variable: most species exhibit two <italic>Opsin </italic>genes, but in the howler monkey an additional gene duplication has led to full trichromacy [##REF##9893843##23##,##REF##10413401##24##]. In <italic>Growth Hormone </italic>gene cluster (five to eight gene copies) some of the duplicate genes in the OWM/hominoid lineages have acquired a novel function and code for Chorionic Somatomammotropin (<italic>CSH </italic>genes) involved in the glucose metabolism of the fetus and the mother. However, the CSH genes are missing from the genomes of NW monkeys [##REF##2744760##61##]. Further species-specific duplications of <italic>GH/CSH </italic>genes have been reported for gibbon, macaque, chimpanzee and human [##REF##15848116##62##,##REF##8404617##63##]. In human-chimpanzee comparison, only three <italic>GH/CSH </italic>genes are clearly orthologous [##REF##15246530##21##]. Other gene clusters with independent gene duplications in human, apes and macaque lineages are <italic>MHC </italic>and testes-expressed <italic>PRAME </italic>genes [##REF##16159394##53##,##REF##17431167##64##]. The ancestral <italic>MHC-B </italic>duplicated into <italic>MHC-B </italic>and <italic>MHC-C </italic>in hominoids. While human <italic>MHC-C </italic>orthologs are found in African apes and orangutans, they are not present in macaque or any other OW monkeys [##REF##8648109##65##,##REF##1560209##66##]. Species-specific evolutionary scenarios have also been reported for <italic>LHB/CGB </italic>genes [##REF##11861891##16##]. In addition to the structural differences between human and chimpanzee <italic>LHB/CGB genes </italic>reported in this study, an expansion of <italic>CGB </italic>genes up to 50 gene copies has been shown in gorilla (<italic>Gorilla gorilla</italic>) [##REF##15252450##4##,##REF##17666543##49##]. Despite the high number of structural differences that has been shown between human and chimpanzee genomes by <italic>in silico </italic>whole-genome analysis (reviewed in [##REF##17597777##67##]), the experimental data for copy number differences between these sister-species has been reported only for a few gene clusters (e.g. <italic>GH/CSH, LHB/CGB, PRAME, MGC8902, CXYorf1, KGF, CD33rSiglec, NANOG</italic>) [##REF##11861891##16##,##REF##16216919##20##,##REF##15331780##52##,##REF##16159394##53##,##REF##16946073##59##,##REF##10655549##68##, ####REF##9326632##69##, ##REF##16469101##70####16469101##70##].</p> <p>In addition to creating structural divergence among the species, duplications provide also the bases for diversification of gene functions. For primate-specific gene duplications, there is evidence of variability in evolutionary rates among the gene copies within and among the species, and of different selective constraints acting on different members of the gene clusters, such as <italic>MHC, beta-globin</italic>, <italic>GH/CSH, PRAME</italic>, <italic>Rh blood group</italic>, <italic>CD33rSiglec </italic>and <italic>beta-defensin </italic>genes [##REF##3412472##71##, ####REF##15486692##72##, ##REF##16843621##73####16843621##73##]. For example, in human-chimpanzee comparison, the chimp <italic>MHC </italic>class I loci <italic>A, B </italic>and <italic>C </italic>are characterized by lower intra-species allelic variation compared to human, providing evidence that ancestral chimpanzee populations may have experienced a selective sweep [##REF##12186979##74##,##REF##15758205##75##]. In marmoset <italic>LHB/CGB </italic>region, a switch of functions has happened between the ancestral <italic>LHB </italic>and the derived <italic>CGB </italic>gene. Although <italic>LHB </italic>and <italic>CGB </italic>genes are both present at the genomic level, only <italic>CGB </italic>gene is expressed in the pituitary and placenta tissues. Thus, Chorionic Gonadotropin Hormone is the only gonadotropin carrying also the luteinizing function fulfilled by Luteinizing Hormone in mammals [##REF##7492691##76##,##REF##14765996##77##].</p> <p>Duplicated genes tend to evolve <italic>in consort </italic>facilitated by active inter-locus gene conversion increasing and preserving sequence similarity among the gene copies [##REF##3855539##37##,##REF##1862096##38##,##REF##17846636##78##]. Concerted evolution within species may lead to erroneous phylogenetic trees and to the overestimation of interspecies divergence dates [##REF##16216919##20##,##REF##10331249##79##,##REF##17009125##80##]. Incorporation of gene conversion data into the equations for calculating interspecies divergence may be still a challenge requiring detailed knowledge of a particular genomic region.</p> <p>So far, only a limited number of reports has been published that focus on detailed variation patterns of duplicated gene families in primates (<italic>MHC</italic>, <italic>Globin</italic>, <italic>GH/CSH </italic>and <italic>LHB/CGB </italic>genes) [##REF##12799463##7##,##REF##16251463##29##,##REF##16702430##81##, ####REF##17223997##82##, ##UREF##2##83####2##83##]. However, the common observation is that compared to single-copy segments, duplicated regions tend to exhibit higher interspecies diversity that could be explained by relaxed selective pressures and/or gene conversions spreading mutations. Thus, when calculating the divergence in duplicated primate-specific regions, the inclusion or exclusion of intraspecies variation data into equations may have a substantial impact on the divergence estimates.</p>
<title>Conclusion</title> <p>We compared the human and chimpanzee duplicated <italic>LHB/CGB </italic>genome clusters and hereby present the detailed evidence for parallel independent duplication events in the two sister-species resulting in discordant number of <italic>CGB </italic>genes (6 in human, 5 in chimpanzee). To our knowledge, this is the first detailed report of parallel duplications in these sister-species leading to structural divergence. The evolutionary fate of duplicated genes is shaped by the interaction of gene conversion and selection. In <italic>LHB/CGB </italic>gene cluster, active gene conversion may have contributed to higher interspecies sequence divergence (both genic and intergenic sequences) and altered transition/transversion ratio compared to the single-copy loci. This higher divergence remained when intraspecies variation was taken into account. However, the drop in divergence estimates after incorporating the intraspecies variation data promotes to reanalyze previously studied loci, where the human-chimpanzee divergence may be substantially lower than initially calculated. Despite the high sequence homology among <italic>LHB/CGB </italic>genes (85–99%), there are signs of functional differentiation among the gene copies. To reconstruct the full evolutionary history of <italic>LHB/CGB </italic>gene cluster, further studies are required comprising high-quality sequence data from several primate species.</p>
<title>Background</title> <p>Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (<italic>LHB</italic>)/Chorionic Gonadotropin Beta (<italic>CGB</italic>) gene cluster as a model. The placental <italic>CGB </italic>genes that are essential for implantation have evolved from an ancestral pituitary <italic>LHB </italic>gene by duplications in the primate lineage.</p> <title>Results</title> <p>We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) <italic>LHB/CGB </italic>regions and hereby present evidence for structural variation resulting in discordant number of <italic>CGB </italic>genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of <italic>LHB/CGB </italic>cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the <italic>LHB/CGB </italic>gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among <italic>LHB/CGB </italic>genes, there are signs of functional differentiation among the gene copies. Estimates for d<sub>n</sub>/d<sub>s </sub>rate ratio suggested a purifying selection on <italic>LHB </italic>and <italic>CGB8</italic>, and a positive evolution of <italic>CGB1</italic>.</p> <title>Conclusion</title> <p>If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.</p>
<title>Abbreviations</title> <p><italic>ADH</italic>: alcohol dehydrogenase; <italic>CD33rSiglec</italic>: CD33-related Siglecs; <italic>CGB</italic>: chorionic gonadotropin beta subunit gene; FAM90A: family with sequence similarity 90; GH/CSH: growth hormone/somatomammotropin; HCG: human chorionic gonadotropin; KGF: keratinocyte growth factor; LH: luteinizing hormone: LHB: luteinizing hormone beta subnit gene; MHC: major histocompatibility complex; NAIP: neuronal apoptosis inhibitory protein; NTF5: neurotrophin 5 gene; PRAME gene family: preferentially expressed antigen of melanoma; psNTF6: neurotrophin 6 pseudogene(s); ROCK1: Rho-dependent protein kinase; RUVBL2: RuvB-like 2, homologue of the bacterial RuvB gene; USP10: ubiquitin-specific protease.</p> <title>Authors' contributions</title> <p>PH, JS, AP and ML designed the research; PH performed the research; JS, AP contributed to the analytic tools; PH and ML analyzed the data and wrote the paper. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Tõnu Margus, Siim Sõber, Tarmo Annilo, Pekka Ellonen, Mari Kaunisto, Maija Wessman and Verneri Anttila for discussions and advice, and Kärt Tomberg for editing the English language. M.L. is a Wellcome Trust International Senior Research Fellow (grants no. 070191/Z/03/Z) in Biomedical Science in Central Europe and a HHMI International Scholar (grant #55005617). Additionally, the study has been supported by the Estonian Ministry of Education and Science core grant no. 0182721s06 and the Estonian Science Foundation grant no. 5796 (M.L., P.H.), as well as personal scholarships from the Centre for International Mobility (CIMO), Kristjan Jaak Stipend Program and World Federation of Scientists' (P.H.), the Center of Excellence grant of Complex Disease Genetics of the Academy of Finland and the Sigrid Juselius Foundation (J.S., A.P).</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Evolution of the <italic>Gonadotropin Hormone Beta </italic>(<italic>CGB</italic>) genes</bold>. Duplication of ancestral <italic>Luteinizing Hormone Beta </italic>(<italic>LHB</italic>) gene in primate lineage has given rise to a novel gene, <italic>CGB</italic>. (A) A simplified schematic presentation of the evolution of <italic>LHB/CGB </italic>genes in primates [##REF##9582510##15##,##REF##11861891##16##]. (B) Comparative structure of the human (GenBank reference: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>) and the chimpanzee <italic>LHB/CGB </italic>cluster (this study, Genbank: accession number <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>) drawn to an approximate scale. Coding genes are depicted as wide empty arrows in the direction of transcription on the sense strand. <italic>A </italic>– <italic>E </italic>indicate the <italic>intergenic regions</italic>. <italic>B' </italic>and <italic>C' </italic>denote putative duplications of the <italic>intergenic regions B </italic>and <italic>C</italic>. Identical color and pattern codes refer to the DNA segments within the cluster with highly similar sequences, the direction of the DNA sequence is indicated on the sense strand. Sequence identity within the cluster: between coding genes 85–99%; <italic>intergenic regions A </italic>and <italic>E </italic>81%; <italic>C </italic>and <italic>C' </italic>96%; <italic>B</italic>, <italic>B' </italic>and <italic>D </italic>ranging 81–98%.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Neighbor-joining trees based on genic (A) and intergenic (B-G) regions within <italic>LHB/CGB </italic>gene cluster</bold>. (A) A phylogenetic tree of the full sequences of <italic>LHB/CGB </italic>gene from <italic>Homo sapiens</italic>, <italic>Pan troglodytes</italic>, <italic>Gorilla gorilla </italic>and <italic>Pongo pygmeaus</italic>. (*) denotes sequences from [##REF##17055150##44##]. Phylogenetic analysis of intergenic regions was conducted with segments without (B-D) and with (E-G) covering <italic>NTF6 </italic>pseudogenes. The homologous segments used for each respective phylogenetic analysis are indicated with a circle on a consensus structure of the intergenic regions in <italic>LHB/CGB </italic>cluster (boxed; from Figure 1B). The nomenclature of the intergenic regions is as on Figure 1B. Bootstrap support values are shown at the nodes (1000 bootstrap replications). Abbreviations: hu – human, ch – chimpanzee, gor – gorilla, orang – orangutan.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Divergence profile between orthologous regions of the human and the chimpanzee <italic>LHB/CGB </italic>clusters</bold>. In total the compared region covered 37,220 bp (Ch)/36,541 (Hu), including 8,084 bp (Ch)/7,973 bp (Hu) from <italic>RUVBL2 </italic>gene to the end of <italic>intergenic region A </italic>and 29,136 bp (Ch)/28,568 bp (Hu) from <italic>CGB1 </italic>to <italic>NTF5 </italic>gene. The species-specific large duplications (human 12,700 bp, chimp 6,725 bp) have been excluded from the comparison. The percents of nucleotide substitutions and indels are calculated per 500 bp non-overlapping windows. Grey arrows indicate the locations of coding genes drawn to an approximate scale. <italic>A </italic>– <italic>E </italic>denote intergenic regions from Figure 1B. Intergenic repeat fraction includes <italic>SINEs, LINEs</italic>, satellites, simple repeats and low complexity DNA sequences within each intergenic region.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Profile of nucleotide substitutions in human and chimpanzee orthologous <italic>LHB/CGB </italic>genes</bold>. Grouping of nucleotide substitutions: (A) Nucleotide substitutions in orthologous <italic>LHB</italic>, <italic>CGB1, CGB5, CGB8 </italic>and <italic>CGB7 </italic>genes (in total 6,878 bp; GC-nucleotide content 64%; 161 substitutions). (B) Nucleotide substitutions in the whole orthologous region of the <italic>LHB/CGB </italic>genome cluster (in total 36,211 bp, GC-nucleotide content 57%, 835 substitutions). Percents for all substitution types are shown with summarized information for transversions and transitions.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Nucleotide divergence (%) between human and chimpanzee orthologous <italic>LHB </italic>and <italic>CGB </italic>genes calculated for promoter, 5'UTR, genic, mRNA, exonic and intronic regions.</bold> Divergence was estimated by using the human (GenBank: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>) and chimpanzee (this study, Genbank: <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>) reference sequences alone or by incorporating the diversity data obtained from re-sequencing for one or both species into the calculations. The re-sequencing data for human (n = 95) originated from the published study[##REF##16251463##29##] and for chimpanzee (n = 11) from the unpublished dataset of the authors.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Maximum likelihood estimation of ω (= d<sub>n</sub>/d<sub>s</sub>) values by PAML analysis and amino acid divergence in human and chimpanzee orthologous genes.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Gene</th><th align="left">n<sup>a</sup></th><th align="left">ω<sup>b</sup></th><th align="left">AA %<sup>c</sup></th></tr></thead><tbody><tr><td align="left">LHB</td><td align="left">141</td><td align="left">0.087</td><td align="left">1.42</td></tr><tr><td align="left"><italic>CGB1</italic></td><td align="left">132</td><td align="left">2.658</td><td align="left">3.03</td></tr><tr><td align="left"><italic>CGB5</italic></td><td align="left">165</td><td align="left">0.180</td><td align="left">1.82</td></tr><tr><td align="left"><italic>CGB8</italic></td><td align="left">165</td><td align="left">0.099</td><td align="left">1.21</td></tr><tr><td align="left"><italic>CGB7</italic></td><td align="left">165</td><td align="left">0.378</td><td align="left">1.82</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Table for sequence parameters. Sequence parameters for the intergenic regions and for the whole <italic>LHB/CGB </italic>cluster in human and chimpanzee.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Distribution of species-specific indels. For simplification we defined all identified gaps in sequence alignments of human and chimpanzee <italic>LHB/CGB </italic>region orthologous segments (Figure ##FIG##0##1B##) as deletions in one of the species. The figure shows the number of species-specific gaps (Y-axes) relative to their length in base pairs (X-axis) and the contribution of each deletion class (1–5; 6–10; 11–15; 16–20; 20–25; &gt;100 bp) to the total length of species-specific gapped sequence.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Table for nucleotide sequence divergence. Nucleotide sequence divergence in orthologous human and chimpanzee <italic>LHB/CGB </italic>genes. Divergence was estimated by using the human (GenBank: <ext-link ext-link-type="gen" xlink:href="NG_000019">NG_000019</ext-link>) and chimpanzee (this study, Genbank: <ext-link ext-link-type="gen" xlink:href="EU000308">EU000308</ext-link>) reference sequences alone or by incorporating the re-sequencing data for one or both species (for human n = 95 [##REF##16251463##29##]; for chimpanzee, n = 11, unpublished data of the authors) into the calculations.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Table of estimated d<sub>n </sub>and d<sub>s </sub>values by Li93 method. Application of Li93 method [##REF##8433381##40##,##REF##8487630##41##] for estimating non-synonymous (d<sub>n</sub>) and synonymous (d<sub>s</sub>) substitutions, and amino acid divergence in human and chimpanzee orthologous genes and for testing significance in the deviation of d<sub>n</sub>/d<sub>s </sub>ratio from expectation under neutrality.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Results of the CRANN analyses of human and chimpanzee orthologous genes (A) <italic>LHB</italic>, (B) <italic>CGB5</italic>, (C) <italic>CGB8</italic>, (D) <italic>CGB7</italic>, (E) <italic>CGB1</italic>. Results of moving window analysis carried out with CRANN [##REF##12468084##47##,##REF##15593409##48##]. X-axis shows the successive windows of 20 codon sites (window size: 20 codons, shift size: 10 codons). As the number of substitutions calculated in each moving window for human and chimpanzee orthologous genes was low, the d<sub>n </sub>and d<sub>s </sub>values were mostly zero and thus the d<sub>n</sub>/d<sub>s </sub>ratio was not shown.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Primer sequences. Primers for sequencing of subcloned BACs 68P2 and 109B10 originating from common chimpanzee (<italic>Pan troglodytes</italic>) BAC library RPCI-43 (BACPAC Resource Center at the Children's Hospital Oakland Research Institute; Oakland, CA).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Alignments of the chimp and human <italic>LHB/CGB </italic>transcripts. The aligned sequences of the major transcripts of the chimp and human <italic>LHB/CGB </italic>genes. The ATG of each gene has been indicated in red and underlined font. Translation STOP codons have been boxed.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a </sup>number of codons in the sequence</p><p><sup>b </sup>non-synonymous/synonymous rate ratio, averaged over sites (d<sub>n</sub>/d<sub>s</sub>).</p><p><sup>c </sup>Divergence at amino acid level (%).</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-195-1"/> <graphic xlink:href="1471-2148-8-195-2"/> <graphic xlink:href="1471-2148-8-195-3"/> <graphic xlink:href="1471-2148-8-195-4"/> <graphic xlink:href="1471-2148-8-195-5"/>
<media xlink:href="1471-2148-8-195-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-195-S7.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Ohno"], "given-names": ["S"], "source": ["Evolution by gene duplication"], "year": ["1970"], "publisher-name": ["New York , Springer"]}, {"surname": ["Nei", "Kumar"], "given-names": ["M", "S"], "source": ["Molecular Evolution and Phylogenetics"], "year": ["2000"], "publisher-name": ["New York , OXFORD University press"]}, {"surname": ["Sedman", "Padhukasahasram", "Kelgo", "Laan"], "given-names": ["L", "B", "P", "M"], "article-title": ["Complex signatures of locus-specific selective pressures and gene conversion on Human Growth Hormone/Chorionic Somatomammotropin genes"], "source": ["Hum Mutat"], "year": ["2008"]}, {"article-title": ["Primer3 software"], "ext-link": ["http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi"]}, {"article-title": ["ClustalW"], "ext-link": ["http://www.ebi.ac.uk/clustalw/"]}, {"article-title": ["EMBOSS package"], "ext-link": ["http://emboss.sourceforge.net/"]}, {"article-title": ["REPEATMASKER program"], "ext-link": ["http://www.repeatmasker.org/"]}]
BMC Evol Biol. 2008 Jul 7; 8:195