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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Characterization of natural killer cells in tamarins: a technical basis for studies of innate immunity.

Natural killer (NK) cells are capable of regulating viral infection without major histocompatibility complex restriction. Hepatitis C is caused by chronic infection with hepatitis C virus (HCV), and impaired activity of NK cells may contribute to the control of the disease progression, although the involvement of NK cells in vivo remains to be proven. GB virus B (GBV-B), which is genetically most closely related to HCV, induces acute and chronic hepatitis upon experimental infection of tamarins. This non-human primate model seems likely to be useful for unveiling the roles of NK cells in vivo. Here we characterized the biological phenotypes of NK cells in tamarins and found that depletion of the CD16+ subset in vivo by administration of a monoclonal antibody significantly reduced the number and activity of NK cells.

Identification of Plasmodium malariae, a human malaria parasite, in imported chimpanzees.

It is widely believed that human malaria parasites infect only man as a natural host. However, earlier morphological observations suggest that great apes are likely to be natural reservoirs as well. To identify malaria parasites in great apes, we screened 60 chimpanzees imported into Japan. Using the sequences of small subunit rRNA and the mitochondrial genome, we identified infection of Plasmodium malariae, a human malaria parasite, in two chimpanzees that were imported about thirty years ago. The chimpanzees have been asymptomatic to the present. In Japan, indigenous malaria disappeared more than fifty years ago; and thus, it is most likely inferred that the chimpanzees were infected in Africa, and P. malariae isolates were brought into Japan from Africa with their hosts, suggesting persistence of parasites at low level for thirty years. Such a long term latent infection is a unique feature of P. malariae infection in humans. To our knowledge, this is the first to report P. malariae infection in chimpanzees and a human malaria parasite from nonhuman primates imported to a nonendemic country.

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.

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.

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.

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.

Human-specific expression of Siglec-6 in the placenta.

CD33-related-Siglecs are lectins on immune cells that recognize sialic acids via extracellular domains, and deliver negative signals via cytosolic tyrosine-based regulatory motifs. We report that while Siglec-6/OB-BP1 (which can also bind leptin) is expressed on immune cells of both humans and the closely related great apes, placental trophoblast expression is human-specific, with little or no expression in ape placentae. Human-specific transcription factor recognition site changes in the Siglec-6 promoter region can help explain the human-specific expression. Human placenta also expresses natural ligands for Siglec-6 (a mixture of glycoproteins carrying cognate sialylated targets), in areas adjacent to Siglec-6 expression. Ligands were also found in uterine endometrium and on cell lines of trophoblastic or endometrial origin. Thus, Siglec-6 was recruited to placental expression during human evolution, presumably to interact with sialylated ligands for specific negative signaling functions and/or to regulate leptin availability. The control of human labor is poorly understood, but involves multiple cues, including placental signaling. Human birthing is also prolonged in comparison to that in our closest evolutionary relatives, the great apes. We found that Siglec-6 levels are generally low in placentae from elective surgical deliveries without known labor and the highest following completion of labor. We therefore speculate that the negative signaling potential of Siglec-6 was recruited to human-specific placental expression, to slow the tempo of the human birth process. The leptin-binding ability of Siglec-6 is also consistent with this hypothesis, as leptin-deficient mice have increased parturition times.

An Alu retrotransposition-mediated deletion of CHD7 in a patient with CHARGE syndrome.

CHD7 mutations account for about 60-65% among more than 200 CHARGE syndrome cases. When rare whole gene deletion cases associated with chromosomal abnormalities are excluded, all mutations of CHD7 reported to date have been point mutations and small deletions and insertions, rather than exonic deletions. To test whether exonic deletions represent a common pathogenic mechanism, we assessed exon copy number by using a recently developed method, the multiplex PCR/liquid chromatography assay (MP/LC). Multiple exons were amplified using unlabeled primers, then separated by ion-pair reversed-phase high-performance liquid chromatography, and quantitated by fluorescence detection using a post-column intercalation dye under the premise that the relative peak intensities for each target directly reflect exon copy number. By using MP/LC, we identified one CHARGE syndrome patient who had a de novo deletion encompassing exons 8-12 among 13 classic CHARGE patients in whom screening by denaturing high-performance liquid chromatography (DHPLC) failed to identify point mutations and small insertions/deletions in CHD7. This is the first CHARGE patient who was documented to have exonic deletion of CHD7. The deletion closely recapitulated the Alu-mediated inactivation of the human CMP-N-acetylneuraminic acid hydroxylase gene (CMP-Neu5Ac hydroxylase), which is regarded as a novel molecular mechanism in the evolution from non-human primates to humans. As demonstrated in this study, MP/LC is a promising method for characterizing exonic deletions, which are largely left unexamined in most routine mutation analysis.