Triangulation supports agricultural spread of the Transeurasian languages.

The origin and early dispersal of speakers of Transeurasian languages-that is, Japanese, Korean, Tungusic, Mongolic and Turkic-is among the most disputed issues of Eurasian population history1-3. A key problem is the relationship between linguistic dispersals, agricultural expansions and population movements4,5. Here we address this question by 'triangulating' genetics, archaeology and linguistics in a unified perspective. We report wide-ranging datasets from these disciplines, including a comprehensive Transeurasian agropastoral and basic vocabulary; an archaeological database of 255 Neolithic-Bronze Age sites from Northeast Asia; and a collection of ancient genomes from Korea, the Ryukyu islands and early cereal farmers in Japan, complementing previously published genomes from East Asia. Challenging the traditional 'pastoralist hypothesis'6-8, we show that the common ancestry and primary dispersals of Transeurasian languages can be traced back to the first farmers moving across Northeast Asia from the Early Neolithic onwards, but that this shared heritage has been masked by extensive cultural interaction since the Bronze Age. As well as marking considerable progress in the three individual disciplines, by combining their converging evidence we show that the early spread of Transeurasian speakers was driven by agriculture.

Infectious diseases may have arrested the southward advance of microblades in Upper Palaeolithic East Asia.

An unsolved archaeological puzzle of the East Asian Upper Palaeolithic is why the southward expansion of an innovative lithic technology represented by microblades stalled at the Qinling-Huaihe Line. It has been suggested that the southward migration of foragers with microblades stopped there, which is consistent with ancient DNA studies showing that populations to the north and south of this line had differentiated genetically by 19 000 years ago. Many infectious pathogens are believed to have been associated with hominins since the Palaeolithic, and zoonotic pathogens in particular are prevalent at lower latitudes, which may have produced a disease barrier. We propose a mathematical model to argue that mortality due to infectious diseases may have arrested the wave-of-advance of the technologically advantaged foragers from the north.

An analysis of the demographic history of the risk allele R4810K in RNF213 of moyamoya disease.

BACKGROUND: Ring finger protein 213 (RNF213) is a susceptibility gene of moyamoya disease (MMD). A previous case-control study and a family analysis demonstrated a strong association of the East Asian-specific variant, R4810K (rs112735431), with MMD. Our aim is to uncover evolutionary history of R4810K in East Asian populations. METHODS: The RNF213 locus of 24 MMD patients in Japan were sequenced using targeted-capture sequencing. Based on the sequence data, we conducted population genetic analysis and estimated the age of R4810K using coalescent simulation. RESULTS: The diversity of the RNF213 gene was higher in Africans than non-Africans, which can be explained by bottleneck effect of the out-of-Africa migration. Coalescent simulation showed that the risk variant was born in East Asia 14,500-5100 years ago and came to the Japanese archipelago afterward, probably in the period when the known migration based on archaeological evidences occurred. CONCLUSIONS: Although clinical data show that the symptoms varies, all sequences harboring the risk allele are almost identical with a small number of exceptions, suggesting the MMD phenotypes are unaffected by the variants of this gene and rather would be more affected by environmental factors.

Confounding effects of microbiome on the susceptibility of TNFSF15 to Crohn's disease in the Ryukyu Islands.

Crohn's disease (CD) involves chronic inflammation in the gastrointestinal tract due to dysregulation of the host immune response to the gut microbiome. Even though the host-microbiome interactions are likely contributors to the development of CD, a few studies have detected genetic variants that change bacterial compositions and increase CD risk. We focus on one of the well-replicated susceptible genes, tumor necrosis factor superfamily member 15 (TNFSF15), and apply statistical analyses for personal profiles of genotypes and salivary microbiota collected from CD cases and controls in the Ryukyu Islands, southernmost islands of the Japanese archipelago. Our association test confirmed the susceptibility of TNFSF15 in the Ryukyu Islands. We found that the recessive model was supported to fit the observed genotype frequency of risk alleles slightly better than the additive model, defining the genetic effect on CD if a pair of the chromosomes in an individual consists of all risk alleles. The combined analysis of haplotypes and salivary microbiome from a small set of samples showed a significant association of the genetic effect with the increase of Prevotella, which led to a significant increase of CD risk. However, the genetic effect on CD disappeared if the abundance of Prevotella was low, suggesting the genetic contribution to CD is conditionally independent given a fixed amount of Prevotella. Although our statistical power is limited due to the small sample size, these results support an idea that the genetic susceptibility of TNFSF15 to CD may be confounded, in part, by the increase of Prevotella.

The allele frequency of ALDH2*Glu504Lys and ADH1B*Arg47His for the Ryukyu islanders and their history of expansion among East Asians.

OBJECTIVES: A cline of frequencies of the derived allele of the ALDH2 gene, which causes a deficiency of an enzyme and "facial flushing" in humans who drink alcohol, has been known among the people of the Japanese archipelago. This cline is conventionally explained by admixture with immigrants from the Asian continent occurring during the Yayoi period. Previous studies lack sufficient data from the peripheral regions of the indigenous Jomon people, and those data the ADH1B gene that is involved in the Class I ADH gene cluster and contains another variant leading to a functional change. METHODS: We focused on the southwestern-most people from the Ryukyu Islands (n = 218) and those from northern Kyushu (n = 21) where the Yayoi immigrants likely arrived. We investigated both the Class I ADH and ALDH2 loci, as well as neutral genetic markers. RESULTS: In the Ryukyu Islands, the frequencies of the ancestral alleles in both loci were always higher than those in mainland Japan, while the frequencies of ADH1B were less than those of the derived allele. A haplotype block was not observed in ALDH2 but was in Class I ADH. DISCUSSION: Our data suggest that the derived allele of ALDH2 came with the Yayoi immigrants from the Asian continent to the Japanese archipelago. However, the derived allele of ADH1B is unlikely to be related to the Yayoi migration. Therefore, we postulate that the expansion of the derived allele of ADHIB in East Asia could be traced back to the last glacial period.

Model-based verification of hypotheses on the origin of modern Japanese revisited by Bayesian inference based on genome-wide SNP data.

Various hypotheses for the peopling of the Japanese archipelago have been proposed, which can be classified into three models: transformation, replacement, and hybridization. In recent years, one of the hybridization models ("dual-structure model") has been widely accepted. According to this model, Neolithic hunter-gatherers known as Jomon, who are assumed to have originated in southeast Asia and lived in the Japanese archipelago greater than 10,000 years ago, admixed with an agricultural people known as Yayoi, whom were migrants from the East Asian continent 2,000-3,000 years ago. Meanwhile, some anthropologists propose that rather, morphological differences between the Jomon and Yayoi people can be explained by microevolution following the lifestyle change. To resolve this controversy, we compared three demographic models by approximate Bayesian computation using genome-wide single nucleotide polymorphism (gwSNP) data from the Ainu people who are thought to be direct descendants of indigenous Jomon. If we assume Chinese people sampled in Beijing from HapMap have the same ancestry as Yayoi, then the hybridization model is predicted to be between 29 and 63 times more likely than the replacement and transformation models, respectively. Furthermore, our data provide strong support for a model in which the Jomon lineages had population structure diversified in local areas before the admixture event. Initial divergence between the Jomon and Yayoi ancestries was dated to late Pleistocene, followed by the divergence of Jomon lineages at early Holocene. These results suggest gwSNP data provides a detailed picture of the complex hybridization model for Japanese population history.

Meiotic recombination counteracts male-biased mutation (male-driven evolution).

Meiotic recombination is believed to produce greater genetic variation despite the fact that deoxyribonucleic acid (DNA)-replication errors are a major source of mutations. In some vertebrates, mutation rates are higher in males than in females, which developed the theory of male-driven evolution (male-biased mutation). However, there is little molecular evidence regarding the relationships between meiotic recombination and male-biased mutation. Here we tested the theory using the frog Rana rugosa, which has both XX/XY- and ZZ/ZW-type sex-determining systems within the species. The male-to-female mutation-rate ratio (α) was calculated from homologous sequences on the X/Y or Z/W sex chromosomes, which supported male-driven evolution. Surprisingly, each α value was notably higher in the XX/XY-type group than in the ZZ/ZW-type group, although α should have similar values within a species. Interestingly, meiotic recombination between homologous chromosomes did not occur except at terminal regions in males of this species. Then, by subdividing α into two new factors, a replication-based male-to-female mutation-rate ratio (ß) and a meiotic recombination-based XX-to-XY/ZZ-to-ZW mutation-rate ratio (γ), we constructed a formula describing the relationship among a nucleotide-substitution rate and the two factors, ß and γ. Intriguingly, the ß- and γ-values were larger and smaller than 1, respectively, indicating that meiotic recombination might reduce male-biased mutations.

Genome-wide SNP analysis reveals population structure and demographic history of the ryukyu islanders in the southern part of the Japanese archipelago.

The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.

Natural allelic variations of xenobiotic-metabolizing enzymes affect sexual dimorphism in Oryzias latipes.

Sexual dimorphisms, which are phenotypic differences between males and females, are driven by sexual selection. Interestingly, sexually selected traits show geographical variations within species despite strong directional selective pressures. This paradox has eluded many evolutionary biologists for some time, and several models have been proposed (e.g. 'indicator model' and 'trade-off model'). However, disentangling which of these theories explains empirical patterns remains difficult, because genetic polymorphisms that cause variation in sexual differences are still unknown. In this study, we show that polymorphisms in cytochrome P450 (CYP) 1B1, which encodes a xenobiotic-metabolizing enzyme, are associated with geographical differences in sexual dimorphism in the anal fin morphology of medaka fish (Oryzias latipes). Biochemical assays and genetic cross experiments show that high- and low-activity CYP1B1 alleles enhanced and declined sex differences in anal fin shapes, respectively. Behavioural and phylogenetic analyses suggest maintenance of the high-activity allele by sexual selection, whereas the low-activity allele possibly has experienced positive selection due to by-product effects of CYP1B1 in inferred ancestral populations. The present data can elucidate evolutionary mechanisms behind genetic variations in sexual dimorphism and indicate trade-off interactions between two distinct mechanisms acting on the two alleles with pleiotropic effects of xenobiotic-metabolizing enzymes.

Human genetic research, race, ethnicity and the labeling of populations: recommendations based on an interdisciplinary workshop in Japan.

BACKGROUND: A challenge in human genome research is how to describe the populations being studied. The use of improper and/or imprecise terms has the potential to both generate and reinforce prejudices and to diminish the clinical value of the research. The issue of population descriptors has not attracted enough academic attention outside North America and Europe. In January 2012, we held a two-day workshop, the first of its kind in Japan, to engage in interdisciplinary dialogue between scholars in the humanities, social sciences, medical sciences, and genetics to begin an ongoing discussion of the social and ethical issues associated with population descriptors. DISCUSSION: Through the interdisciplinary dialogue, we confirmed that the issue of race, ethnicity and genetic research has not been extensively discussed in certain Asian communities and other regions. We have found, for example, the continued use of the problematic term, "Mongoloid" or continental terms such as "European," "African," and "Asian," as population descriptors in genetic studies. We, therefore, introduce guidelines for reporting human genetic studies aimed at scientists and researchers in these regions. CONCLUSION: We need to anticipate the various potential social and ethical problems entailed in population descriptors. Scientists have a social responsibility to convey their research findings outside of their communities as accurately as possible, and to consider how the public may perceive and respond to the descriptors that appear in research papers and media articles.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago.

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

Generation of chimpanzee induced pluripotent stem cell lines for cross-species comparisons.

As humans' closest living relatives, chimpanzees offer valuable insights into human evolution. However, technical and ethical limitations hinder investigations into the molecular and cellular foundations that distinguish chimpanzee and human traits. Recently, induced pluripotent stem cells (iPSCs) have emerged as a novel model for functional comparative studies and provided a non-invasive alternative for studying embryonic phenomena. In this study, we generated five new chimpanzee iPSC lines from peripheral blood cells and skin fibroblasts with SeV vectors carrying four reprogramming factors (human OCT3/4, SOX2, KLF4, and L-MYC) and characterized their pluripotency and differentiation potential. We also examined the expression of a human-specific non-coding RNA, HSTR1, which is predicted to be involved in human brain development. Our results show that the chimpanzee iPSCs possess pluripotent characteristics and can differentiate into various cell lineages. Moreover, we found that HSTR1 is expressed in human iPSCs and their neural derivatives but not in chimpanzee counterparts, supporting its possible role in human-specific brain development. As iPSCs are inherently variable due to genetic and epigenetic differences in donor cells or reprogramming procedures, it is essential to expand the number of chimpanzee iPSC lines to comprehensively capture the molecular and cellular properties representative of chimpanzees. Hence, our cells provide a valuable resource for investigating the function and regulation of human-specific transcripts such as HSTR1 and for understanding human evolution more generally.

Generation of monoclonal antibodies against the Galß1-4Gal epitope: a key tool in studies of species-specific glycans expressed in fish, amphibians and birds.

Whereas the Galß1-4Gal epitope is rarely found in mammalian glycans, it has been found in glycans of various species of non-mammalian vertebrates, such as fish, amphibians and birds. Although glycans containing Galß1-4Gal in these vertebrates were detected by precise structural analysis of the glycans using mass spectrometry and/or NMR spectrometry, there are no convenient methods to detect Galß1-4Gal from various samples. To analyze systematically the distribution of Galß1-4Gal in nature, we generated mouse monoclonal antibodies (mAbs) specific for Galß1-4Gal using extracts of medaka eggs as an immunogen. Four mAbs (two immunoglobulin (Ig)Ms and two IgG1s) were obtained by enzyme-linked immunosorbent assay-based screening. The specificities of these mAbs were evaluated by frontal affinity chromatography using 142 kinds of 2-aminopyridine (PA)-derivatized oligosaccharides. While all mAbs interacted with (Galß1-4Gal)-containing oligosaccharides at their non-reducing termini with dissociation constants (K(d)) ranging from 1.0 x 10⁻5 to 2.8 x 10⁻4 M, no apparent interaction was observed with any other glycans. The number of branches containing Galß1-4Gal on N-glycans did not significantly affect K(d) of mAbs of IgG1 subclasses, but those of IgM mAbs were decreased by ∼1 order of magnitude, in increments of the number of branches present. Using the mAbs, we established that Galß1-4Gal is also expressed on glycoproteins in various tissues from the African clawed frog. Immunohistochemical staining of medaka sections revealed that Galß1-4Gal epitopes were expressed in the endothelium, epithelium and epidermis, which directly contact the external environment or invading organisms. Thus, these mAbs are useful for systematically investigating the species-specific expression of glycans, which may act as a barrier against infection.

Crohn's disease risk alleles on the NOD2 locus have been maintained by natural selection on standing variation.

Risk alleles for complex diseases are widely spread throughout human populations. However, little is known about the geographic distribution and frequencies of risk alleles, which may contribute to differences in disease susceptibility and prevalence among populations. Here, we focus on Crohn's disease (CD) as a model for the evolutionary study of complex disease alleles. Recent genome-wide association studies and classical linkage analyses have identified more than 70 susceptible genomic regions for CD in Europeans, but only a few have been confirmed in non-European populations. Our analysis of eight European-specific susceptibility genes using HapMap data shows that at the NOD2 locus the CD-risk alleles are linked with a haplotype specific to CEU at a frequency that is significantly higher compared with the entire genome. We subsequently examined nine global populations and found that the CD-risk alleles spread through hitchhiking with a high-frequency haplotype (H1) exclusive to Europeans. To examine the neutrality of NOD2, we performed phylogenetic network analyses, coalescent simulation, protein structural prediction, characterization of mutation patterns, and estimations of population growth and time to most recent common ancestor (TMRCA). We found that while H1 was significantly prevalent in European populations, the H1 TMRCA predated human migration out of Africa. H1 is likely to have undergone negative selection because 1) the root of H1 genealogy is defined by a preexisting amino acid substitution that causes serious conformational changes to the NOD2 protein, 2) the haplotype has almost become extinct in Africa, and 3) the haplotype has not been affected by the recent European expansion reflected in the other haplotypes. Nevertheless, H1 has survived in European populations, suggesting that the haplotype is advantageous to this group. We propose that several CD-risk alleles, which destabilize and disrupt the NOD2 protein, have been maintained by natural selection on standing variation because the deleterious haplotype of NOD2 is advantageous in diploid individuals due to heterozygote advantage and/or intergenic interactions.

Endocrine, inflammatory and immune responses and individual differences in acute hypobaric hypoxia in lowlanders.

When lowlanders are exposed to environments inducing hypobaric hypoxia (HH) such as high mountains, hemodynamic changes occur to maintain oxygen levels in the body. However, changes to other physiological functions under such conditions have yet to be clarified. This study investigated changes in endocrine, inflammatory and immune parameters and individual differences during acute HH exposure using a climatic chamber (75 min of exposure to conditions mimicking 3500 m) in healthy lowlanders. Aldosterone and cortisol were significantly decreased and interleukin (IL)-6, IL-8 and white blood cell (WBC) counts were significantly increased after HH. Lower peripheral oxygen saturation (SpO2) was associated with higher IL-6 and WBC counts, and higher IL-8 was associated with higher cortisol. These findings suggest that endocrine, inflammatory and immune responses are evoked even with a short 75-min exposure to HH and individuals with lower SpO2 seemed to show more pronounced responses. Our results provide basic data for understanding the physiological responses and interactions of homeostatic systems during acute HH.

A common variation in EDAR is a genetic determinant of shovel-shaped incisors.

Shovel shape of upper incisors is a common characteristic in Asian and Native American populations but is rare or absent in African and European populations. Like other common dental traits, genetic polymorphisms involved in the tooth shoveling have not yet been clarified. In ectodysplasin A receptor (EDAR), where dysfunctional mutations cause hypohidrotic ectodermal dysplasia, there is a nonsynonymous-derived variant, 1540C (rs3827760), that has a geographic distribution similar to that of the tooth shoveling. This allele has been recently reported to be associated with Asian-specific hair thickness. We aimed to clarify whether EDAR 1540C is also associated with dental morphology. For this purpose, we measured crown diameters and tooth-shoveling grades and analyzed the correlations between the dental traits and EDAR genotypes in two Japanese populations, inhabitants around Tokyo and in Sakishima Islands. The number of EDAR 1540C alleles in an individual was strongly correlated with the tooth-shoveling grade (p = 7.7 x 10(-10)). The effect of the allele was additive and explained 18.9% of the total variance in the shoveling grade, which corresponds to about one-fourth of the heritability of the trait reported previously. For data reduction of individual-level metric data, we applied a principal-component analysis, which yielded PC1-4, corresponding to four patterns of tooth size; this result implies that multiple factors are involved in dental morphology. The 1540C allele also significantly affected PC1 (p = 4.9 x 10(-3)), which denotes overall tooth size, and PC2 (p = 2.6 x 10(-3)), which denotes the ratio of mesiodistal diameter to buccolingual diameter.

Effects of an Asian-specific nonsynonymous EDAR variant on multiple dental traits.

Dental morphology is highly diverse among individuals and between human populations. Although it is thought that genetic factors mainly determine common dental variations, only a few such genetic factors have been identified. One study demonstrated that a nonsynonymous single-nucleotide polymorphism (370V/A, rs3827760) in the ectodysplasin A receptor gene (EDAR) is associated with shoveling and double-shoveling grades of upper first incisors and tooth crown size. Here, we examined the association of EDAR 370V/A with several dental characters in Korean and Japanese subjects. A meta-analysis that combined analyses of Korean and Japanese subjects revealed that the Asian-specific 370A allele is associated with an increase in the grades of shoveling and double shoveling, as previously found. We also showed a highly significant association between EDAR 370V/A genotype and crown size, especially mesiodistal diameters of anterior teeth. Moreover, we found that the 370A allele was associated with the presence of hypoconulids of lower second molars. These results indicated that the EDAR polymorphism is responsible, in part, for the Sinodonty and Sundadonty dichotomy in Asian populations, and clearly demonstrated that the EDAR polymorphism has pleiotropic effects on tooth morphology. As the 370A allele is known to be a most likely target of positive selection in Asian populations, some phenotypes associated with the variant may be 'hitchhiking phenotypes', while others may be actual targets of selection.

The history of human populations in the Japanese Archipelago inferred from genome-wide SNP data with a special reference to the Ainu and the Ryukyuan populations.

The Japanese Archipelago stretches over 4000 km from north to south, and is the homeland of the three human populations; the Ainu, the Mainland Japanese and the Ryukyuan. The archeological evidence of human residence on this Archipelago goes back to >30 000 years, and various migration routes and root populations have been proposed. Here, we determined close to one million single-nucleotide polymorphisms (SNPs) for the Ainu and the Ryukyuan, and compared these with existing data sets. This is the first report of these genome-wide SNP data. Major findings are: (1) Recent admixture with the Mainland Japanese was observed for more than one third of the Ainu individuals from principal component analysis and frappe analyses; (2) The Ainu population seems to have experienced admixture with another population, and a combination of two types of admixtures is the unique characteristics of this population; (3) The Ainu and the Ryukyuan are tightly clustered with 100% bootstrap probability followed by the Mainland Japanese in the phylogenetic trees of East Eurasian populations. These results clearly support the dual structure model on the Japanese Archipelago populations, though the origins of the Jomon and the Yayoi people still remain to be solved.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine.

The transition to agriculture occurred relatively late in Eastern Europe, leading researchers to debate whether it was a gradual, interactive process or a colonisation event. In the forest and forest-steppe regions of Ukraine, farming appeared during the fifth millennium BCE, associated with the Cucuteni-Trypillia cultural complex (CTCC, ~ 5000-3000 BCE). Across Europe, the Neolithisation process was highly variable across space and over time. Here, we investigate the population dynamics of early agriculturalists from the eastern forest-steppe region based on the analyses of 20 ancient genomes from the site of Verteba Cave (3935-825 cal BCE). Results reveal that the CTCC individuals' ancestry is related to both western hunter-gatherers and Near Eastern farmers, has no local ancestry associated with Ukrainian Neolithic hunter-gatherers and has steppe ancestry. An Early Bronze Age individual has an ancestry profile related to the Yamnaya expansions but with 20% of ancestry related to the other Trypillian individuals, which suggests admixture between the Trypillians and the incoming populations carrying steppe-related ancestry. A Late Bronze Age individual dated to 980-825 cal BCE has a genetic profile indicating affinity to Beaker-related populations, detected close to 1000 years after the end of the Bell Beaker phenomenon during the third millennium BCE.

Gene conversion and purifying selection shape nucleotide variation in gibbon L/M opsin genes.

BACKGROUND: Routine trichromatic color vision is a characteristic feature of catarrhines (humans, apes and Old World monkeys). This is enabled by L and M opsin genes arrayed on the X chromosome and an autosomal S opsin gene. In non-human catarrhines, genetic variation affecting the color vision phenotype is reported to be absent or rare in both L and M opsin genes, despite the suggestion that gene conversion has homogenized the two genes. However, nucleotide variation of both introns and exons among catarrhines has only been examined in detail for the L opsin gene of humans and chimpanzees. In the present study, we examined the nucleotide variation of gibbon (Catarrhini, Hylobatidae) L and M opsin genes. Specifically, we focused on the 3.6~3.9-kb region that encompasses the centrally located exon 3 through exon 5, which encode the amino acid sites functional for the spectral tuning of the genes. RESULTS: Among 152 individuals representing three genera (Hylobates, Nomascus and Symphalangus), all had both L and M opsin genes and no L/M hybrid genes. Among 94 individuals subjected to the detailed DNA sequencing, the nucleotide divergence between L and M opsin genes in the exons was significantly higher than the divergence in introns in each species. The ratio of the inter-LM divergence to the intra-L/M polymorphism was significantly lower in the introns than that in synonymous sites. When we reconstructed the phylogenetic tree using the exon sequences, the L/M gene duplication was placed in the common ancestor of catarrhines, whereas when intron sequences were used, the gene duplications appeared multiple times in different species. Using the GENECONV program, we also detected that tracts of gene conversions between L and M opsin genes occurred mostly within the intron regions. CONCLUSIONS: These results indicate the historical accumulation of gene conversions between L and M opsin genes in the introns in gibbons. Our study provides further support for the homogenizing role of gene conversion between the L and M opsin genes and for the purifying selection against such homogenization in the central exons to maintain the spectral difference between L and M opsins in non-human catarrhines.