Coevolution of Siglec-11 and Siglec-16 via gene conversion in primates.

BACKGROUND: Siglecs-11 and -16 are members of the sialic acid recognizing Ig-like lectin family, and expressed in same cells. Siglec-11 functions as an inhibitory receptor, whereas Siglec-16 exhibits activating properties. In humans, SIGLEC11 and SIGLEC16 gene sequences are extremely similar in the region encoding the extracellular domain due to gene conversions. Human SIGLEC11 was converted by the nonfunctional SIGLEC16P allele, and the converted SIGLEC11 allele became fixed in humans, possibly because it provides novel neuroprotective functions in brain microglia. However, the detailed evolutionary history of SIGLEC11 and SIGLEC16 in other primates remains unclear. RESULTS: We analyzed SIGLEC11 and SIGLEC16 gene sequences of multiple primate species, and examined glycan binding profiles of these Siglecs. The phylogenetic tree demonstrated that gene conversions between SIGLEC11 and SIGLEC16 occurred in the region including the exon encoding the sialic acid binding domain in every primate examined. Functional assays showed that glycan binding preference is similar between Siglec-11 and Siglec-16 in all analyzed hominid species. Taken together with the fact that Siglec-11 and Siglec-16 are expressed in the same cells, Siglec-11 and Siglec-16 are regarded as paired receptors that have maintained similar ligand binding preferences via gene conversions. Relaxed functional constraints were detected on the SIGLEC11 and SIGLEC16 exons that underwent gene conversions, possibly contributing to the evolutionary acceptance of repeated gene conversions. The frequency of nonfunctional SIGLEC16P alleles is much higher than that of SIGLEC16 alleles in every human population. CONCLUSIONS: Our findings indicate that Siglec-11 and Siglec-16 have been maintained as paired receptors by repeated gene conversions under relaxed functional constraints in the primate lineage. The high prevalence of the nonfunctional SIGLEC16P allele and the fixation of the converted SIGLEC11 imply that the loss of Siglec-16 and the gain of Siglec-11 in microglia might have been favored during the evolution of human lineage.

Specific inactivation of two immunomodulatory SIGLEC genes during human evolution.

Sialic acid-recognizing Ig-like lectins (Siglecs) are signaling receptors that modulate immune responses, and are targeted for interactions by certain pathogens. We describe two primate Siglecs that were rendered nonfunctional by single genetic events during hominin evolution after our common ancestor with the chimpanzee. SIGLEC13 was deleted by an Alu-mediated recombination event, and a single base pair deletion disrupted the ORF of SIGLEC17. Siglec-13 is expressed on chimpanzee monocytes, innate immune cells that react to bacteria. The human SIGLEC17P pseudogene mRNA is still expressed at high levels in human natural killer cells, which bridge innate and adaptive immune responses. As both resulting pseudogenes are homozygous in all human populations, we resurrected the originally encoded proteins and examined their functions. Chimpanzee Siglec-13 and the resurrected human Siglec-17 recruit a signaling adapter and bind sialic acids. Expression of either Siglec in innate immune cells alters inflammatory cytokine secretion in response to Toll-like receptor-4 stimulation. Both Siglecs can also be engaged by two potentially lethal sialylated bacterial pathogens of newborns and infants, agents with a potential impact on reproductive fitness. Neanderthal and Denisovan genomes show human-like sequences at both loci, corroborating estimates that the initial pseudogenization events occurred in the common ancestral population of these hominins. Both loci also show limited polymorphic diversity, suggesting selection forces predating the origin of modern humans. Taken together, these data suggest that genetic elimination of Siglec-13 and/or Siglec-17 represents signatures of infectious and/or other inflammatory selective processes contributing to population restrictions during hominin origins.

Evolution of siglec-11 and siglec-16 genes in hominins.

We previously reported a human-specific gene conversion of SIGLEC11 by an adjacent paralogous pseudogene (SIGLEC16P), generating a uniquely human form of the Siglec-11 protein, which is expressed in the human brain. Here, we show that Siglec-11 is expressed exclusively in microglia in all human brains studied-a finding of potential relevance to brain evolution, as microglia modulate neuronal survival, and Siglec-11 recruits SHP-1, a tyrosine phosphatase that modulates microglial biology. Following the recent finding of a functional SIGLEC16 allele in human populations, further analysis of the human SIGLEC11 and SIGLEC16/P sequences revealed an unusual series of gene conversion events between two loci. Two tandem and likely simultaneous gene conversions occurred from SIGLEC16P to SIGLEC11 with a potentially deleterious intervening short segment happening to be excluded. One of the conversion events also changed the 5' untranslated sequence, altering predicted transcription factor binding sites. Both of the gene conversions have been dated to ~1-1.2 Ma, after the emergence of the genus Homo, but prior to the emergence of the common ancestor of Denisovans and modern humans about 800,000 years ago, thus suggesting involvement in later stages of hominin brain evolution. In keeping with this, recombinant soluble Siglec-11 binds ligands in the human brain. We also address a second-round more recent gene conversion from SIGLEC11 to SIGLEC16, with the latter showing an allele frequency of ~0.1-0.3 in a worldwide population study. Initial pseudogenization of SIGLEC16 was estimated to occur at least 3 Ma, which thus preceded the gene conversion of SIGLEC11 by SIGLEC16P. As gene conversion usually disrupts the converted gene, the fact that ORFs of hSIGLEC11 and hSIGLEC16 have been maintained after an unusual series of very complex gene conversion events suggests that these events may have been subject to hominin-specific selection forces.

Power of neutrality tests for detecting natural selection.

Detection of natural selection is one of the main interests in population genetics. Thus, many tests have been developed for detecting natural selection using genomic data. Although it is recognized that the utility of tests depends on several evolutionary factors, such as the timing of selection, strength of selection, frequency of selected alleles, demographic events, and initial frequency of selected allele when selection started acting (softness of selection), the relationships between such evolutionary factors and the power of tests are not yet entirely clear. In this study, we investigated the power of 4 tests: Tajiama's D, Fay and Wu's H, relative extended haplotype homozygosity (rEHH), and integrated haplotype score (iHS), under ranges of evolutionary parameters and demographic models to quantitatively expand the understanding of approaches for detecting selection. The results show that each test detects selection within a limited parameter range, and there are still wide ranges of parameters for which none of these tests work effectively. In addition, the parameter space in which each test shows the highest power overlaps the empirical results of previous research. These results indicate that our present perspective of adaptation is limited to only a part of actual adaptation.

Lower promoter activity of the ST8SIA2 gene has been favored in evolving human collective brains.

ST8SIA2 is an important molecule regulating expression of the phenotype involved in schizophrenia. Lowered promoter activity of the ST8SIA2 gene is considered to be protective against schizophrenia by conferring tolerance to psychosocial stress. Here, we examined the promoter-type composition of anatomically modern humans (AMHs) and archaic humans (AHs; Neanderthals and Denisovans), and compared the promoter activity at the population level (population promoter activity; PPA) between them. In AMHs, the TCT-type, showing the second lowest promoter activity, was most prevalent in the ancestral population of non-Africans. However, the detection of only the CGT-type from AH samples and recombination tracts in AH sequences showed that the CGT- and TGT-types, exhibiting the two highest promoter activities, were common in AH populations. Furthermore, interspecies gene flow occurred into AMHs from AHs and into Denisovans from Neanderthals, influencing promoter-type compositions independently in both AMHs and AHs. The difference of promoter-type composition makes PPA unique in each population. East and Southeast Asian populations show the lowest PPA. This results from the selective increase of the CGC-type, showing the lowest promoter activity, in these populations. Every non-African population shows significantly lower PPA than African populations, resulting from the TCT-type having the highest prevalence in the ancestral population of non-Africans. In addition, PPA reduction is also found among subpopulations within Africa via a slight increase of the TCT-type. These findings indicate a trend toward lower PPA in the spread of AMHs, interpreted as a continuous adaptation to psychosocial stress arising in migration. This trend is considered as genetic tuning for the evolution of collective brains. The inferred promoter-type composition of AHs differed markedly from that of AMHs, resulting in higher PPA in AHs than in AMHs. This suggests that the trend toward lower PPA is a unique feature in AMH spread.

Genetic basis of human brain evolution: accelerating along the primate speedway.

Using novel variations of traditional methods, report in the December 29(th) issue of Cell that diverse genes involved in neural biology (particularly those critical in development) show higher rates of protein evolution in primates than in rodents-particularly in the lineage leading to humans.

Geographical origin of Plasmodium vivax in the Republic of Korea: haplotype network analysis based on the parasite's mitochondrial genome.

BACKGROUND: The Republic of Korea (South Korea) is one of the countries where vivax malaria had been successfully eradicated by the late 1970s. However, re-emergence of vivax malaria in South Korea was reported in 1993. Several epidemiological studies and some genetic studies using antigenic molecules of Plasmodium vivax in the country have been reported, but the evolutionary history of P. vivax has not been fully understood. In this study, the origin of the South Korean P. vivax population was estimated by molecular phylogeographic analysis. METHODS: A haplotype network analysis based on P. vivax mitochondrial (mt) DNA sequences was conducted on 11 P. vivax isolates from South Korea and another 282 P. vivax isolates collected worldwide. RESULTS: The network analysis of P. vivax mtDNA sequences showed that the coexistence of two different groups (A and B) in South Korea. Groups A and B were identical or close to two different populations in southern China. CONCLUSIONS: Although the direct introduction of the two P. vivax populations in South Korea were thought to have been from North Korea, the results of this analysis suggest the genealogical origin to be the two different populations in southern China.

Two-dimensional site frequency spectrum for detecting, classifying and dating incomplete selective sweeps.

The two-dimensional site frequency spectrum (2D SFS) was investigated to describe the intra-allelic variability (IAV) maintained within a derived allele (D) group that has undergone an incomplete selective sweep against an ancestral allele group. We observed that recombination certainly muddles the ancestral relationships of allelic lineages between the two allele groups; however, the 2D SFS reveals intriguing signatures of recombination as well as the genealogical structure of the D group, particularly the size of a mutation and the time to the most recent common ancestor (TMRCA). Coalescent simulations were performed to achieve powerful and robust 2D SFS-based statistics with special reference to accurate evaluation of IAV, significance of recombination effects, and distinction between hard and soft selective sweeps. These studies were extended to a case wherein an incomplete selective sweep is no longer in progress and ceased in the recent past. The 2D SFS-based method was applied to 100 intronic linkage disequilibrium regions randomly chosen from the East Asian population of modern humans to examine the P value distributions of the summary statistics under the null hypothesis of neutrality in a nonequilibrium demographic model. We argue that about 96% of intronic variants are non-adaptive with a 10% false discovery rate. Furthermore, this method was applied to six genomic regions in Eurasian populations that were claimed to have experienced recent selective sweeps. We found that two of these genomic regions did not have significant signals of selective sweeps, but the remaining four had undergone hard and soft sweeps and were dated, in terms of TMRCA, after the major out-of-Africa dispersal of modern humans.

Big bang in the evolution of extant malaria parasites.

Malaria parasites (genus Plasmodium) infect all classes of terrestrial vertebrates and display host specificity in their infections. It is therefore assumed that malaria parasites coevolved intimately with their hosts. Here, we propose a novel scenario of malaria parasite-host coevolution. A phylogenetic tree constructed using the malaria parasite mitochondrial genome reveals that the extant primate, rodent, bird, and reptile parasite lineages rapidly diverged from a common ancestor during an evolutionary short time period. This rapid diversification occurred long after the establishment of the primate, rodent, bird, and reptile host lineages, which implies that host-switch events contributed to the rapid diversification of extant malaria parasite lineages. Interestingly, the rapid diversification coincides with the radiation of the mammalian genera, suggesting that adaptive radiation to new mammalian hosts triggered the rapid diversification of extant malaria parasite lineages.

Genetic diversity and population structure of Plasmodium falciparum in the Philippines.

BACKGROUND: In the Philippines, malaria morbidity and mortality have decreased since the 1990s by effective malaria control. Several epidemiological surveys have been performed in the country, but the characteristics of the Plasmodium falciparum populations are not yet fully understood. In this study, the genetic structure of P. falciparum populations in the Philippines was examined. METHODS: Population genetic analyses based on polymorphisms of 10 microsatellite loci of the parasite were conducted on 92 isolates from three provinces (Kalinga, Palawan, and Davao del Norte) with different malaria endemicity. RESULTS: The levels of genetic diversity and the effective population sizes of P. falciparum in the Philippines were similar to those reported in the mainland of Southeast Asia or South America. In the low malaria transmission area (Kalinga), there was a low level of genetic diversity and a strong linkage disequilibrium (LD) when the single-clone haplotype (SCH) was used in the multilocus LD analysis, while in the high malaria transmission areas (Palawan and Davao del Norte), there was a high level of genetic diversity and a weak LD when SCH was used in the multilocus LD analysis. On the other hand, when the unique haplotypes were used in the multilocus LD analysis, no significant LD was observed in the Kalinga and the Palawan populations. The Kalinga and the Palawan populations were, therefore, estimated to have an epidemic population structure. The three populations were moderately differentiated from each other. CONCLUSION: In each area, the level of genetic diversity correlates with the local malaria endemicity. These findings confirm that population genetic analyses using microsatellite loci are a useful tool for evaluating malaria endemicity.

A macroscopic examination of M. biceps femoris and M. gluteus maximus in the orangutan.

The musculature of the hip and thigh in the orangutan has been described previously. Anatomically, there are various descriptions among primates in those structures, in particular, the relationship between M. biceps femoris and M. gluteus maximus, their derivatives, and the muscle segment. However, a detailed innervation system to this ischiofemoral part has not been described, thus there is still uncertainty as to with which muscle it is associated. In this analysis, we examined the gross anatomy of the hip and thigh muscles of the orangutan and chimpanzee, including their innervation. Also, a comparison was made with documented data of other primates. As a result of these observations, it was found that the ischiofemoral part in the orangutan is innervated by the same sciatic nerve branch (the common peroneal nerve) as the long head of M. biceps femoris, but not by the same nerve as M. gluteus maximus. Therefore, the ischiofemoral part is appropriately considered as a part of the long head of M. biceps femoris. It appears that this morphologic feature is an adjustment to the arboreal life of the orangutan. The development of the flexor complex of the thigh is necessary for this arboreal adaptation, resulting in a unique musculature of M. biceps femoris in the orangutan.

A new inference method for detecting an ongoing selective sweep.

A simple method was developed to detect signatures of ongoing selective sweeps in single nucleotide polymorphism (SNP) data. Based largely on the traditional site frequency spectrum (SFS), the method additionally incorporates linkage disequilibrium (LD) between pairs of SNP sites and uniquely represents both SFS and LD information as hierarchical "barcodes." This barcode representation allows the identification of a hitchhiking genomic region surrounding a putative target site of positive selection, or a core site. Sweep signals at linked neutral sites are then measured by the proportion (Fc) of derived alleles within the hitchhiking region that are linked in the derived allele group defined at the core site. In measuring Fc or intra-allelic variability in an informative way, certain conditions for derived allele frequencies are required, as illustrated with the human ST8SIA2 locus. Coalescent simulators with and without positive selection are used to assess the false-positive and false-negative rates of the Fc statistic. To demonstrate its power, the method was further applied to the LCT, OCA2, EDAR, SLC24A5 and ASPM loci, which are known to have undergone positive selection in human populations. Overall, the method is powerful and can be used to identify core sites responsible for ongoing selective sweeps.

Positive selection on schizophrenia-associated ST8SIA2 gene in post-glacial Asia.

A number of loci are associated with highly heritable schizophrenia and the prevalence of this mental illness has had considerable negative fitness effects on human populations. Here we focused on one particular schizophrenia-associated gene that encodes a sialyltransferase (ST8SIA2) and is expressed preferentially in the brain with the level being largely determined by three SNPs in the promoter region. It is suggested that the expression level of the ST8SIA2 gene is a genetic determinant of schizophrenia risk, and we found that a geographically differentiated non-risk SNP type (CGC-type) has significantly reduced promoter activity. A newly developed method for detecting ongoing positive selection was applied to the ST8SIA2 genomic region with the identification of an unambiguous sweep signal in a rather restricted region of 18 kb length surrounding the promoter. We also found that while the CGC-type emerged in anatomically modern humans in Africa over 100 thousand years ago, it has increased its frequency in Asia only during the past 20-30 thousand years. These findings support that the positive selection is driven by psychosocial stress due to changing social environments since around the last glacial maximum, and raise a possibility that schizophrenia extensively emerged during the Upper Paleolithic and Neolithic era.

Recent independent evolution of msp1 polymorphism in Plasmodium vivax and related simian malaria parasites.

The Plasmodium MSP-1 is a promising malaria vaccine candidate. However, the highly polymorphic nature of the MSP-1 gene (msp1) presents a potential obstacle for effective vaccine development. To investigate the evolutionary history of msp1 polymorphism in P. vivax, we construct phylogenetic trees of msp1 from P. vivax and related monkey malaria parasite species. All P. vivax msp1 alleles cluster in the P. vivax lineage and are not distributed among other species. Similarly, all P. cynomolgi msp1 alleles cluster in the P. cynomolgi lineage. This suggests that, in contrast to presumed ancient origin of P. falciparum msp1 polymorphism, the origin of P. vivax msp1 polymorphism is relatively recent. We observed positive selection in the P. vivax lineage but not in P. cynomolgi. Also, positive selection acts on different regions of msp1 in P. vivax and P. falciparum. This study shows that the evolutionary history of msp1 differs greatly among parasite lineages.

Fixation of the human-specific CMP-N-acetylneuraminic acid hydroxylase pseudogene and implications of haplotype diversity for human evolution.

The human CMP-N-acetylneuraminic acid hydroxylase gene (CMAH) suffered deletion of an exon that encodes an active center for the enzyme approximately 3.2 million years ago (MYA). We analyzed a 7.3-kb intronic region of 132 CMAH genes to explore the fixation process of this pseudogene and the demographic implication of its haplotype diversity. Fifty-six variable sites were sorted into 18 different haplotypes with significant linkage disequilibrium. Despite the rather low nucleotide diversity, the most recent common ancestor at CMAH dates to 2.9 MYA. This deep genealogy follows shortly after the original exon deletion, indicating that the deletion has fixed in the population, although whether this fixation was facilitated by natural selection remains to be resolved. Remarkable features are exceptionally long persistence of two lineages and the confinement of one lineage in Africa, implying that some African local populations were in relative isolation while others were directly involved in multiple African exoduses of the genus Homo. Importantly, haplotypes found in Eurasia suggest interbreeding between then-contemporaneous human species. Although population structure within Africa complicates the interpretation of phylogeographic information of haplotypes, the data support a single origin of modern humans, but not with complete replacement of archaic inhabitants by modern humans.

Discovery of Siglec-14, a novel sialic acid receptor undergoing concerted evolution with Siglec-5 in primates.

Immune receptors that show high mutual sequence similarity and have antagonizing signaling properties are called paired receptors, and are believed to fine-tune immune responses. Siglecs are sialic acid-recognizing receptors of the immunoglobulin (Ig) superfamily expressed on immune cells. Human Siglec-5, encoded by SIGLEC5 gene, has four extracellular Ig-like domains and a cytosolic inhibitory motif. We discovered human Siglec-14 with three Ig-like domains, encoded by the SIGLEC14 gene, adjacent to SIGLEC5. Human Siglec-14 has almost complete sequence identity with human Siglec-5 at the first two Ig-like domains, shows a glycan binding preference similar to that of human Siglec-5, and associates with the activating adapter protein DAP12. Thus, Siglec-14 and Siglec-5 appear to be the first glycan binding paired receptors. Near-complete sequence identity of the amino-terminal part of human Siglec-14 and Siglec-5 indicates partial gene conversion between SIGLEC14 and SIGLEC5. Remarkably, SIGLEC14 and SIGLEC5 in other primates also show evidence of gene conversions within each lineage. Evidently, balancing the interactions between Siglec-14, Siglec-5 and their common ligand(s) had selective advantage during the course of evolution. The "essential arginine" critical for sialic acid recognition in both Siglec-14 and Siglec-5 is present in humans but mutated in almost all great ape alleles.

The time course of lymph drainage from the peritoneal cavity in beagle dogs.

Lymph drainage routes from the abdominal and pelvic cavities in beagle dogs were observed serially by following the time course of India ink administered intraperitoneally. Four systems of lymph drainage routes from the peritoneal cavity were observed in this study. The earliest drainage returned to the cranial mediastinal lymph nodes via the sternal lymph vessels; subsequently, the sternal lymph nodes located along the internal thoracic artery became involved. Then, a drainage route via the lymph vessel along the left vagus nerve was observed. The final drainage route flowed into the lateral lymph vessel through the thoracic duct located on the vertebra. These results show that India ink is absorbed from the peritoneal cavity, and that the lymph drainage first flows mainly towards the cranial mediastinal lymph nodes through the ventral lymphatic channels. Our serial observations suggest that, over time, the lymph drainage routes changed from the ventral abdominal to the dorsal thoracic lymphatic channels in the thorax.

Emergence of infectious malignant thrombocytopenia in Japanese macaques (Macaca fuscata) by SRV-4 after transmission to a novel host.

We discovered a lethal hemorrhagic syndrome arising from severe thrombocytopenia in Japanese macaques kept at the Primate Research Institute, Kyoto University. Extensive investigation identified that simian retrovirus type 4 (SRV-4) was the causative agent of the disease. SRV-4 had previously been isolated only from cynomolgus macaques in which it is usually asymptomatic. We consider that the SRV-4 crossed the so-called species barrier between cynomolgus and Japanese macaques, leading to extremely severe acute symptoms in the latter. Infectious agents that cross the species barrier occasionally amplify in virulence, which is not observed in the original hosts. In such cases, the new hosts are usually distantly related to the original hosts. However, Japanese macaques are closely related to cynomolgus macaques, and can even hybridize when given the opportunity. This lethal outbreak of a novel pathogen in Japanese macaques highlights the need to modify our expectations about virulence with regards crossing species barriers.

Age of the last common ancestor of extant Plasmodium parasite lineages.

Parasites of the genus Plasmodium infect all classes of amniotes (mammals, birds and reptiles) and display host specificity in their infections. It is therefore generally believed that Plasmodium parasites co-evolved intimately with their hosts. Here, we report that based on an evolutionary analysis using 22 genes in the nuclear genome, extant lineages of Plasmodium parasites originated roughly in the Oligocene epoch after the emergence of their hosts. This timing on the age of the common ancestor of extant Plasmodium parasites suggest the importance of host switches and lends support to the evolutionary scenario of a "malaria big bang" that was proposed based on the evolutionary analysis using the mitochondrial genome.

Improved capacity of a monkey-tropic HIV-1 derivative to replicate in cynomolgus monkeys with minimal modifications.

Human immunodeficiency virus type 1 (HIV-1) hardly replicates in Old World monkeys. Recently, a mutant HIV-1 clone, NL-DT5R, in which a small part of gag and the entire vif gene are replaced with SIVmac239-derived ones, was shown to be able to replicate in pigtail monkeys but not in rhesus monkeys (RM). In the present study, we found that a modified monkey-tropic HIV-1 (HIV-1mt), MN4-5S, acquired the ability to replicate efficiently in cynomolgus monkeys as compared with the NL-DT5R, while neither NL-DT5R nor MN4-5S replicated in RM cells. These results suggest that multiple determinants may be involved in the restriction of HIV-1 replication in macaques, depending on the species of macaques. The new HIV-1mt clone will be useful for studying molecular mechanisms by which anti-viral host factors regulate HIV-1 replication in macaques.