Different Host-Endogenous Retrovirus Relationships between Mammals and Birds Reflected in Genome-Wide Evolutionary Interaction Patterns.

Mammals and birds differ largely in their average endogenous retrovirus loads, namely the proportion of endogenous retrovirus in the genome. The host-endogenous retrovirus relationships, including conflict and co-option, have been hypothesized among the causes of this difference. However, there has not been studies about the genomic evolutionary signal of constant host-endogenous retrovirus interactions in a long-term scale and how such interactions could lead to the endogenous retrovirus load difference. Through a phylogeny-controlled correlation analysis on ∼5,000 genes between the dN/dS ratio of each gene and the load of endogenous retrovirus in 12 mammals and 21 birds, separately, we detected genes that may have evolved in association with endogenous retrovirus loads. Birds have a higher proportion of genes with strong correlation between dN/dS and the endogenous retrovirus load than mammals. Strong evidence of association is found between the dN/dS of the coding gene for leucine-rich repeat-containing protein 23 and endogenous retrovirus load in birds. Gene set enrichment analysis shows that gene silencing rather than immunity and DNA recombination may have a larger contribution to the association between dN/dS and the endogenous retrovirus load for both mammals and birds. The above results together showing different evolutionary patterns between bird and mammal genes can partially explain the apparently lower endogenous retrovirus loads of birds, while gene silencing may be a universal mechanism that plays a remarkable role in the evolutionary interaction between the host and endogenous retrovirus. In summary, our study presents signals that the host genes might have driven or responded to endogenous retrovirus load changes in long-term evolution.

Japanese wolves are most closely related to dogs and share DNA with East Eurasian dogs.

Although the domestic dog's origin is still unclear, this lineage is believed to have been domesticated from an extinct population of gray wolves, which is expected to be more closely related to dogs than to other populations of gray wolves. Here, we sequence the whole genomes of nine Japanese wolves (7.5-100x: Edo to Meiji periods) and 11 modern Japanese dogs and analyze them together with those from other populations of dogs and wolves. A phylogenomic tree shows that, among the gray wolves, Japanese wolves are closest to the dog, suggesting that the ancestor of dogs is closely related to the ancestor of the Japanese wolf. Based on phylogenetic and geographic relationships, the dog lineage has most likely originated in East Asia, where it diverged from a common ancestor with the Japanese wolf. Since East Eurasian dogs possess Japanese wolf ancestry, we estimate an introgression event from the ancestor of the Japanese wolf to the ancestor of the East Eurasian dog that occurred before the dog's arrival in the Japanese archipelago.

Population dynamics in the Japanese Archipelago since the Pleistocene revealed by the complete mitochondrial genome sequences.

The Japanese Archipelago is widely covered with acidic soil made of volcanic ash, an environment which is detrimental to the preservation of ancient biomolecules. More than 10,000 Palaeolithic and Neolithic sites have been discovered nationwide, but few skeletal remains exist and preservation of DNA is poor. Despite these challenging circumstances, we succeeded in obtaining a complete mitogenome (mitochondrial genome) sequence from Palaeolithic human remains. We also obtained those of Neolithic (the hunting-gathering Jomon and the farming Yayoi cultures) remains, and over 2,000 present-day Japanese. The Palaeolithic mitogenome sequence was not found to be a direct ancestor of any of Jomon, Yayoi, and present-day Japanese people. However, it was an ancestral type of haplogroup M, a basal group of the haplogroup M. Therefore, our results indicate continuity in the maternal gene pool from the Palaeolithic to present-day Japanese. We also found that a vast increase of population size happened and has continued since the Yayoi period, characterized with paddy rice farming. It means that the cultural transition, i.e. rice agriculture, had significant impact on the demographic history of Japanese population.

Evolutionary History of the Risk of SNPs for Diffuse-Type Gastric Cancer in the Japanese Population.

A genome wide association study reported that the T allele of rs2294008 in a cancer-related gene, PSCA, is a risk allele for diffuse-type gastric cancer. This allele has the highest frequency (0.63) in Japanese in Tokyo (JPT) among 26 populations in the 1000 Genomes Project database. FST ≈ 0.26 at this single nucleotide polymorphism is one of the highest between JPT and the genetically close Han Chinese in Beijing (CHB). To understand the evolutionary history of the alleles in PSCA, we addressed: (i) whether the C non-risk allele at rs2294008 is under positive selection, and (ii) why the mainland Japanese population has a higher T allele frequency than other populations. We found that haplotypes harboring the C allele are composed of two subhaplotypes. We detected that positive selection on both subhaplotypes has occurred in the East Asian lineage. However, the selection on one of the subhaplotypes in JPT seems to have been relaxed or ceased after divergence from the continental population; this may have caused the elevation of T allele frequency. Based on simulations under the dual structure model (a specific demography for the Japanese) and phylogenetic analysis with ancient DNA, the T allele at rs2294008 might have had high frequency in the Jomon people (one of the ancestral populations of the modern Japanese); this may explain the high T allele frequency in the extant Japanese.

Do Genes Associated with Dyslexia of Chinese Characters Evolve Neutrally?

Dyslexia, or reading disability, is found to have a genetic basis, and several related genes have been reported. We investigated whether natural selection has acted on single nucleotide polymorphisms (SNPs) that were reported to be associated with risk/non-risk for the reading disability of Chinese characters. We applied recently developed 2D SFS-based statistics to SNP data of East Asian populations to examine whether there is any sign of selective sweep. While neutrality was not rejected for most SNPs, significant signs of selection were detected for two linkage disequilibrium (LD) regions containing the reported SNPs of GNPTAB and DCDC2. Furthermore, we searched for a selection target site among the SNPs in these LD regions, because a causal site is not necessarily a reported SNP but could instead be a tightly linked site. In both LD regions, we found candidate target sites, which may have an effect on expression regulation and have been selected, although which genes these SNPs affect remains unknown. Because most people were not engaged in reading until recently, it is unlikely that there has been selective pressure on reading ability itself. Consistent with this, our results suggest a possibility of genetic hitchhiking, whereby alleles of the reported SNPs may have increased in frequency together with the selected target, which could have functions for other genes and traits apart from reading ability.

Characterization of TRPA channels in the starfish Patiria pectinifera: involvement of thermally activated TRPA1 in thermotaxis in marine planktonic larvae.

The vast majority of marine invertebrates spend their larval period as pelagic plankton and are exposed to various environmental cues. Here we investigated the thermotaxis behaviors of the bipinnaria larvae of the starfish, Patiria pectinifera, in association with TRPA ion channels that serve as thermal receptors in various animal species. Using a newly developed thermotaxis assay system, we observed that P. pectinifera larvae displayed positive thermotaxis toward high temperatures, including toward temperatures high enough to cause death. In parallel, we identified two TRPA genes, termed PpTRPA1 and PpTRPA basal, from this species. We examined the phylogenetic position, spatial expression, and channel properties of each PpTRPA. Our results revealed the following: (1) The two genes diverged early in animal evolution; (2) PpTRPA1 and PpTRPA basal are expressed in the ciliary band and posterior digestive tract of the larval body, respectively; and (3) PpTRPA1 is activated by heat stimulation as well as by known TRPA1 agonists. Moreover, knockdown and rescue experiments demonstrated that PpTRPA1 is involved in positive thermotaxis in P. pectinifera larvae. This is the first report to reveal that TRPA1 channels regulate the behavioral response of a marine invertebrate to temperature changes during its planktonic larval period.

Acropora digitifera Encodes the Largest Known Family of Fluorescent Proteins that Has Persisted during the Evolution of Acropora Species.

Fluorescent proteins (FPs) are well known and broadly used as bio-imaging markers in molecular biology research. Many FP genes were cloned from anthozoan species and it was suggested that multi-copies of these genes are present in their genomes. However, the full complement of FP genes in any single coral species remained unidentified. In this study, we analyzed the FP genes in two stony coral species. FP cDNA sequences from Acropora digitifera and Acropora tenuis revealed the presence of a multi-gene family with an unexpectedly large number of genes, separated into short-/middle-wavelength emission (S/MWE), middle-/long-wavelength emission (M/LWE), and chromoprotein (CP) clades. FP gene copy numbers in the genomes of four A. digitifera colonies were estimated as 16-22 in the S/MWE, 3-6 in the M/LWE, and 8-12 in the CP clades, and, in total, 35, 31, 33, and 33 FP gene copies per individual shown by quantitative PCR. To the best of our knowledge, these are the largest sets of FP genes per genome. The fluorescent light produced by recombinant protein products encoded by the newly isolated genes explained the fluorescent range of live A. digitifera, suggesting that the high copy multi-FP gene family generates coral fluorescence. The functionally diverse multi-FP gene family must have existed in the ancestor of Acropora species, as suggested by molecular phylogenetic and evolutionary analyses. The persistence of a diverse function and high copy number multi-FP gene family may indicate the biological importance of diverse fluorescence emission and light absorption in Acropora species.

mtDNA diversity of the Zapotec in Mexico suggests a population decline long before the first contact with Europeans.

The New World is the last continent colonized by anatomically modern humans, Homo sapiens. The first migrants entered the New World from Asia through Beringia. It is suggested that there were three streams of Asian gene flow, one major and two additional minor gene flows. The first major migrants took a Pacific coastal route and began spreading to the American continent before the opening of the ice-free corridor. We investigated the diversity of full-length mitochondrial DNA genomes of the Zapotec population, residing in the Mesoamerican region, and reconstructed their demographic history using Bayesian Skyline Plots. We estimated the initial date of gene flow into the New World by Zapotec ancestors at around 17 000­19 000 years ago,which is highly concordant with previous studies. We also show a population decline after the initial expansion. This decline started 4000 years ago, long before European contact with Native Americans. This indicates that other factors including climatec hange should be considered to explain the observed demographic pattern.

Complete mitogenome analysis of indigenous populations in Mexico: its relevance for the origin of Mesoamericans.

Mesoamerica has an important role in the expansion of Paleoamericans as the route to South America. In this study, we determined complete mitogenome sequences of 113 unrelated individuals from two indigenous populations of Mesoamerica, Mazahua and Zapotec. All newly sequenced mitogenomes could be classified into haplogroups A2, B2, C1 and D1, but one sequence in Mazahua was D4h3a, a subclade of haplogroup D4. This haplogroup has been mostly found in South America along the Pacific coast. Haplogroup X2a was not found in either population. Genetic similarity obtained using phylogenetic tree construction and principal component analysis showed that these two populations are distantly related to each other. Actually, the Mazahua and the Zapotec shared no sequences (haplotypes) in common, while each also showed a number of unique subclades. Surprisingly, Zapotec formed a cluster with indigenous populations living in an area from central Mesoamerica to Central America. By contrast, the Mazahua formed a group with indigenous populations living in external areas, including southwestern North America and South America. This intriguing genetic relationship suggests the presence of two paleo-Mesoamerican groups, invoking a scenario in which one group had expanded into South America and the other resided in Mesoamerica.

An autosomal analysis gives no genetic evidence for complex speciation of humans and chimpanzees.

There have been conflicting arguments as to what happened in the human-chimpanzee speciation event. Patterson et al. (2006, Genetic evidence for complex speciation of humans and chimpanzees. Nature 441:1103-1108) proposed a hypothesis that the human-chimpanzee speciation event involved a complicated demographic process: that is, the ancestral lineages of humans and chimpanzees experienced temporal isolation followed by a hybridization event. This hypothesis stemmed from two major observations: a wide range of human-chimpanzee nucleotide divergence across the autosomal genome and very low divergence in the X chromosome. In contrast, Innan and Watanabe (2006, The effect of gene flow on the coalescent time in the human-chimpanzee ancestral population. Mol Biol Evol. 23:1040-1047) demonstrated that the null model of instantaneous speciation fits the genome-wide divergence data for the two species better than alternative models involving partial isolation and migration. To reconcile these two conflicting reports, we first reexamined the analysis of autosomal data by Patterson et al. (2006). By providing a theoretical framework for their analysis, we demonstrated that their observation is what is theoretically expected under the null model of instantaneous speciation with a large ancestral population. Our analysis indicated that the observed wide range of autosomal divergence is simply due to the coalescent process in the large ancestral population of the two species. To further verify this, we developed a maximum likelihood function to detect evidence of hybridization in genome-wide divergence data. Again, the null model with no hybridization best fits the data. We conclude that the simplest speciation model with instantaneous split adequately describes the human-chimpanzee speciation event, and there is no strong reason to involve complicated factors in explaining the autosomal data.

Comparative genetics of the poly-Q tract of ataxin-1 and its binding protein PQBP-1.

Human PQBP-1 is known to interact with triplet repeat disease gene products such as ataxin and huntingtin through their poly-glutamine (poly-Q) tracts. The poly-Q tracts show extensive variation in both the number and the configuration of repeats among species. A surface plasmon resonance assay showed clear interaction between human PQBP-1 and Q(11), representative of the poly-Q tract of the ataxin-1 of Old World monkeys. No response was observed using Q(2)PQ(2)P(4)Q(2), representative of the poly-Q tract of the ataxin-1 of New World monkeys. This implies that the interaction of human PQBP-1 with ataxin-1 is limited to humans and closely related species. Comparison of the human and mouse PQBP-1 sequences showed an elevated amino acid substitution rate in the polar amino acid-rich domain of PQBP-1 that is responsible for binding to poly-Q tracts. This could have been advantageous to the new biological function of human PQBP-1 through poly-Q tracts.

Elevated evolutionary rate in genes with homopolymeric amino acid repeats constituting nondisordered structure.

Homopolymeric amino acid repeats are tandem repeats of single amino acids. About 650 genes are known to have repeats of this kind comprising seven residues or more in the human genome. According to the evolutionary conservativeness, we classified the repeats into three categories: those whose length is conserved among mammals (CM), those whose length differs among nonprimate mammals but is conserved among primates (CP), and those whose length differs among primates (VP). The frequency of each repeat, especially Ala, Leu, Pro, and Glu repeats, varies greatly in each category. The 3D structure of homopolymeric amino acid repeats is considered to be intrinsically disordered. As expected, a large proportion of the repeats had a disordered structure, and nearly half of the repeats were predicted as completely disordered. However, a number of the repeats predicted to have nondisordered structure: 13% and 25% of the repeats for categories CM and VP, respectively. Comparison of the substitution rates showed a higher Ka/Ks ratio for the genes with not disordered repeats than the genes with disordered repeats. These results indicate that amino acid substitution rates have been elevated in the genes with nondisordered repeats.

Potential of fish opsin gene duplications to evolve new adaptive functions.

The duplication of four cone-opsin gene families is heavily involved in visual adaptation in bony fish. We found that two gene families for the middle-wave range of the vision spectrum have, on average, older duplications followed by accelerated amino acid substitution, in comparison with the other two families that define the boundaries. This could be due to the difference in the potential to evolve new functions; middle-wave genes should have greater contribution to adaptive vision evolution through gene duplication.

Quantifying dominance and deleterious effect on human disease genes.

Human genes responsible for inherited diseases are important for the understanding of human disease. We investigated the degree of polymorphism and divergence in the human disease genes to elucidate the effect of natural selection on human disease genes. In particular, the effect of disease dominance was incorporated into the analysis. Both dominant disease genes (DDG) and recessive disease genes (RDG) had a higher mutation rate per site and encoded longer proteins than the nondisease genes, which exposed the disease genes to a faster flux of new mutations. Using an unbiased polymorphism dataset, we found that, proportionally, RDG harbor more nonsynonymous polymorphisms compared with DDG. We estimated the selection intensity on the disease genes using polymorphism and divergence data and determined whether the different patterns of polymorphism and divergence between DDG and RDG could be explained by the difference in only dominance. Even after the dominance effect was considered, the selection intensity on RDG was significantly different from DDG, suggesting that the deleterious effect of the dominant and recessive disease mutations are fundamentally different.

Adaptive evolution in humans revealed by the negative correlation between the polymorphism and fixation phases of evolution.

The selective forces acting on amino acid substitutions may be different in the two phases of molecular evolution: polymorphism and fixation. Negative selection and genetic drift may dominate the first phase, whereas positive selection may become much more significant in the second phase. However, the conventional dichotomy of synonymous vs. nonsynonymous changes does not offer the resolution needed to study the dynamics of these two phases. Following previously published methods, we separated amino acid changes into 75 elementary types (1-bp substitution between their respective codons). The likelihood of each type of amino acid change becoming polymorphic (PI, which stands for "polymorphic index"), relative to synonymous changes, can then be calculated. Similarly, the likelihood of fixation (FI, for "fixation index"), conditional on common polymorphisms, is also calculated. Using Perlegen and HapMap data on human polymorphisms and the chimpanzee sequences as the outgroup, we compared the evolutionary dynamics of the 75 elementary changes in the two phases. We found a strong "L-shaped" negative correlation (P < 0.001) between FI and PI. Only those changes with low PIs show FI > 1, which is often a signature of adaptive evolution. These patterns suggest that negative and positive selection operate more effectively on the same set of amino acid changes and that approximately 10-13% of amino acid substitutions between humans and chimpanzee may be adaptive.