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.

Nonequilibrium Neutral Theory for Hitchhikers.

Selective sweep is a phenomenon of reduced variation at presumably neutrally evolving sites (hitchhikers) in the genome that is caused by the spread of a selected allele at a linked focal site, and is widely used to test for action of positive selection. Nonetheless, selective sweep may also provide an unprecedented opportunity for studying nonequilibrium properties of the neutral variation itself. We have demonstrated this possibility in relation to ancient selective sweep for modern human-specific changes and ongoing selective sweep for local population-specific changes.

Role of error catastrophe in transmission ability of virus.

The role played by error catastrophe is explicitly taken into account in a mathematical formulation to analyze COVID-19 data. The idea is to combine the mathematical genetics formalism of the error catastrophe of mutations in virus gene loci with the standard model of epidemics, which lacks the explicit incorporation of the effect of mutation on the spreading of viruses. We apply this formalism to the case of SARS-CoV-2 virus. We assume the universality of the error catastrophe in the process of analyzing the data. This means that some basic parameter to describe the error catastrophe is independent of which group (country or city) we deal with. Concretely, we analyze Omicron variant data from South Africa and then analyze cases from Japan using the same value of the basic parameter derived in the South Africa analysis. The excellent fit between the two sets of data, one from South Africa and the other from Japan, using the common values of genetic parameters, justifies our assumption of the universality of these parameters.

Evolutionary origin of peptidoglycan recognition proteins in vertebrate innate immune system.

BACKGROUND: Innate immunity is the ancient defense system of multicellular organisms against microbial infection. The basis of this first line of defense resides in the recognition of unique motifs conserved in microorganisms, and absent in the host. Peptidoglycans, structural components of bacterial cell walls, are recognized by Peptidoglycan Recognition Proteins (PGRPs). PGRPs are present in both vertebrates and invertebrates. Although some evidence for similarities and differences in function and structure between them has been found, their evolutionary history and phylogenetic relationship have remained unclear. Such studies have been severely hampered by the great extent of sequence divergence among vertebrate and invertebrate PGRPs. Here we investigate the birth and death processes of PGRPs to elucidate their origin and diversity. RESULTS: We found that (i) four rounds of gene duplication and a single domain duplication have generated the major variety of present vertebrate PGRPs, while in invertebrates more than ten times the number of duplications are required to explain the repertoire of present PGRPs, and (ii) the death of genes in vertebrates appears to be almost null whereas in invertebrates it is frequent. CONCLUSION: These results suggest that the emergence of new PGRP genes may have an impact on the availability of the repertoire and its function against pathogens. These striking differences in PGRP evolution of vertebrates and invertebrates should reflect the differences in the role of their innate immunity. Insights on the origin of PGRP genes will pave the way to understand the evolution of the interaction between host and pathogens and to lead to the development of new treatments for immune diseases that involve proteins related to the recognition of self and non-self.

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.

Frequent gene conversion events between the X and Y homologous chromosomal regions in primates.

BACKGROUND: Mammalian sex-chromosomes originated from a pair of autosomes. A step-wise cessation of recombination is necessary for the proper maintenance of sex-determination and, consequently, generates a four strata structure on the X chromosome. Each stratum shows a specific per-site nucleotide sequence difference (p-distance) between the X and Y chromosomes, depending on the time of recombination arrest. Stratum 4 covers the distal half of the human X chromosome short arm and the p-distance of the stratum is approximately 10%, on average. However, a 100-kb region, which includes KALX and VCX, in the middle of stratum 4 shows a significantly lower p-distance (1-5%), suggesting frequent sequence exchanges or gene conversions between the X and Y chromosomes in humans. To examine the evolutionary mechanism for this low p-distance region, sequences of a corresponding region including KALX/Y from seven species of non-human primates were analyzed. RESULTS: Phylogenetic analysis of this low p-distance region in humans and non-human primate species revealed that gene conversion like events have taken place at least ten times after the divergence of New World monkeys and Catarrhini (i.e., Old World monkeys and hominoids). A KALY-converted KALX allele in white-handed gibbons also suggests a possible recent gene conversion between the X and Y chromosomes. In these primate sequences, the proximal boundary of this low p-distance region is located in a LINE element shared between the X and Y chromosomes, suggesting the involvement of this element in frequent gene conversions. Together with a palindrome on the Y chromosome, a segmental palindrome structure on the X chromosome at the distal boundary near VCX, in humans and chimpanzees, may mediate frequent sequence exchanges between X and Y chromosomes. CONCLUSION: Gene conversion events between the X and Y homologous regions have been suggested, mainly in humans. Here, we found frequent gene conversions in the evolutionary course of primates. An insertion of a LINE element at the proximal end of the region may be a cause for these frequent conversions. This gene conversion in humans may also be one of the genetic causes of Kallmann syndrome.

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.

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.

Origin and evolution of processed pseudogenes that stabilize functional Makorin1 mRNAs in mice, primates and other mammals.

We investigate the origin and evolution of a mouse processed pseudogene, Makorin1-p1, whose transcripts stabilize functional Makorin1 mRNAs. It is shown that Makorin1-p1 originated almost immediately before the musculus and cervicolor species groups diverged from each other some 4 million years ago and that the Makorin1-p1 orthologs in various Mus species are transcribed. However, Mus caroli in the cervicolor species group expresses not only Makorin1-p1, but also another older Makorin1-derived processed pseudogene, demonstrating the rapid generation and turnover in subgenus Mus. Under this circumstance, transcribed processed pseudogenes (TPPs) of Makorin1 evolved in a strictly neutral fashion even with an enhanced substitution rate at CpG dinucleotide sites. Next, we extend our analyses to rats and other mammals. It is shown that although these species also possess their own Makorin1-derived TPPs, they occur rather infrequently in simian primates. Under this circumstance, it is hypothesized that already existing TPPs must be prevented from accumulating detrimental mutations by negative selection. This hypothesis is substantiated by the presence of two rather old TPPs, MKRNP1 and MKRN4, in humans and New World monkeys. The evolutionary rate and pattern of Makorin1-derived processed pseudogenes depend heavily on how frequently they are disseminated in the genome.

Evolutionary history of sex-linked mammalian amelogenin genes.

Amelogenin (AMEL) arose prior to the emergence of tetrapods and transposed into an intron of the Rho GTPase-activating protein 6 gene. In the mammalian lineage leading to eutherians, a pair of homologous autosomes with this nested gene structure fused with the then already differentiating sex chromosomes by suppressing homologous recombination. As sex-chromosomal differentiation extended to the fused region, a pair of homologous AMEL genes too differentiated from each other in two steps; first in the 5' region (the promoter region to transposon MER5 in intron 2) and second in the remaining 3' region. This resulted in gametologous AMELX and AMELY in the eutherian sex chromosomes. Although the early differentiation of the 5' region between AMELX and AMELY is consistent with the lowered expression level of AMELY, there is no indication for deterioration of AMELY at the amino acid level. Rather, both AMELX and AMELY in particular lineages might undergo positive selection, followed by negative selection to preserve established function. Based on patterns and levels of AMELX and AMELY polymorphisms in the human population, it is also argued that a recombination cold spot near AMELX might be related to the cause of the ancient pseudoautosomal boundary.

Lineage-specific loss of function of bitter taste receptor genes in humans and nonhuman primates.

Since the process of becoming dead genes or pseudogenes (pseudogenization) is irreversible and can occur rather rapidly under certain environmental circumstances, it is one plausible determinant for characterizing species specificity. To test this evolutionary hypothesis, we analyzed the tempo and mode of duplication and pseudogenization of bitter taste receptor (T2R) genes in humans as well as in 12 nonhuman primates. The results show that primates have accumulated more pseudogenes than mice after their separation from the common ancestor and that lineage-specific pseudogenization becomes more conspicuous in humans than in nonhuman primates. Although positive selection has operated on some amino acids in extracellular domains, functional constraints against T2R genes are more relaxed in primates than in mice and this trend has culminated in the rapid deterioration of the bitter-tasting capability in humans. Since T2R molecules play an important role in avoiding generally bitter toxic and harmful substances, substantial modification of the T2R gene repertoire is likely to reflect different responses to changes in the environment and to result from species-specific food preference during primate evolution.

Further evidence for recombination between mouse hemoglobin beta b1 and b2 genes based on the nucleotide sequences of intron, UTR, and intergenic spacer regions.

Nucleotide sequences of the intron regions and UTRs (Untranslated regions) of the hemoglobin beta adult genes, b1 and b2, and of the intergenic spacer region were determined for mouse strains representing the d, p, and w1 hemoglobin haplotypes defined by protein electrophoretic analyses. The hypothesis of recombination of the b1 and b2 genes between the d and w1 haplotypes previously reported in the cDNA nucleotide sequences was confirmed by neighbor-joining analyses of the intron regions and UTRs within the b1 and b2 genes, suggesting that all of the structures of hemoglobin beta adult genes support the hypothesis that the p haplotype was established by hybridization between d and w1 haplotype mice. The resultant recombinant of the p haplotype was found to have a d-like b1 gene and a w1-like b2 gene. In addition to the possible recombination, a break point was suggested around 2-3 kb downstream of the b1 gene within the intergenic spacer region, despite the absence of clear properties that could stimulate the recombination machinery. Some large insertions or deletions (indels) specific to the p or d haplotypes were located within the intergenic spacer region, in which the 1010-bp indel specific to the p haplotype was shared by all examined strains representing the p haplotype.

Evolutionary relationships of major histocompatibility complex class I genes in simian primates.

New World monkeys (NWMs) occupy a critical phylogenetic position in elucidating the evolutionary process of major histocompatibility complex (MHC) class I genes in primates. From three subfamilies of Aotinae, Cebinae, and Atelinae, the 5'-flanking regions of 18 class I genes are obtained and phylogenetically examined in terms of Alu/LINE insertion elements as well as the nucleotide substitutions. Two pairs of genes from Aotinae and Atelinae are clearly orthologous to human leukocyte antigen (HLA) -E and -F genes. Of the remaining 14 genes, 8 belong to the distinct group B, together with HLA-B and -C, to the exclusion of all other HLA class I genes. These NWM genes are classified into four groups, designated as NWM-B1, -B2, -B3, and -B4. Of these, NWM-B2 is orthologous to HLA-B/C. Also, orthologous relationships of NWM-B1, -B2, and -B3 exist among different families of Cebidae and Atelidae, which is in sharp contrast to the genus-specific gene organization within the subfamily Callitrichinae. The other six genes belong to the distinct group G. However, a clade of these NWM genes is almost equally related to HLA-A, -J, -G, and -K, and there is no evidence for their orthologous relationships to HLA-G. It is argued that class I genes in simian primates duplicated extensively in their common ancestral lineage and that subsequent evolution in descendant species has been facilitated mainly by independent loss of genes.

Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: implications for the origin and diversification of seagrasses in Japanese waters.

Seagrasses are composed of four families belonging to angiosperms and they are thought to become adaptive to aquatic life independently. Zosteraceae is one such family and because of the relatively high species diversity around Japan and Korea coast areas, the family might have arisen therefrom. To elucidate the origin and evolution of Zosteraceae which consists of three genera, Phyllospadix, Zostera, and Heterozostera, 2.8 kb nucleotide sequences of rbcL and matK genes in the chloroplast genome were examined for various species, including cosmopolitan Z. marina and endemic Z. caulescens. The phylogenetic analysis reveals the following three features. First, based on the synonymous nucleotide substitution rate of the rice chloroplast genome, we estimated the divergence times between Zosteraceae and its closest relative, Potamogetonaceae, and between different genera, Zostera and Phyllospadix, as approximately 100 million years (myr) and 36 myr, respectively, suggesting that Zosteraceae emerged somewhere in the period from 36 myr ago to 100 myr ago. Second, two subgenera of Zostera, Zostera and Zosterella, exhibit their reciprocal monophyly and appear to have differentiated from each other approximately 33 myr ago. However, the third genus Heterozostera branched off only 5 myr ago from the stem lineage leading to Zosterella and this seems too recent in comparison with the ancient divergence of the two subgenera. Third, we estimated the most recent common ancestor of subgenus Zostera as 6 myr. In Z. marina four haplotypes were found in the sample and have diversified in the past 1.5 myr. One haplotype is shared by both sides of the Japan Archipelago and its closely related haplotypes occur also in eastern Pacific Ocean. Based on these phylogeographic analyses, we propose a provisional age related classification of Zosteraceae to argue the origin and evolution.

Use of Cumulative Poisson Probability Distribution as an Estimator of the Recombination Rate in an Expanding Population: Example of the Macaca fascicularis Major Histocompatibility Complex.

We describe a method to estimate the rate of recombination per generation from the genotypes of a large individual sample of an expanding population, for which the founding event is dated. The approach is illustrated with an application to estimating the major histocompatibility complex (MHC) recombination rate in the Mauritian macaque population. We genotyped 750 macaques by means of 17 microsatellites across the MHC region and reconstructed the seven most frequent haplotypes assumed to represent the founding haplotypes (H(rec(0))) as well as the 31% recombinant haplotypes (H(rec(h))) resulting from a variable number "h" of recombinations between the founding haplotypes. The relative frequencies of the various classes of haplotypes (H(rec(0)) and H(rec(h))) follow a Poisson distribution. By using a maximum likelihood method, we calculated the mean of the Poisson distribution that best fits the data. By dividing this mean by the number of generations (50-100) from the date of the population founding, we deduced that rate of recombination in the MHC is approximately 0.004 to 0.008 in the Mauritian macaque population. When the founding date of the population is precisely known, our method presents a useful alternative to the coalescent method.

Genomic structure and evolution of multigene families: "flowers" on the human genome.

We report the results of an extensive investigation of genomic structures in the human genome, with a particular focus on relatively large repeats (>50 kb) in adjacent chromosomal regions. We named such structures "Flowers" because the pattern observed on dot plots resembles a flower. We detected a total of 291 Flowers in the human genome. They were predominantly located in euchromatic regions. Flowers are gene-rich compared to the average gene density of the genome. Genes involved in systems receiving environmental information, such as immunity and detoxification, were overrepresented in Flowers. Within a Flower, the mean number of duplication units was approximately four. The maximum and minimum identities between homologs in a Flower showed different distributions; the maximum identity was often concentrated to 100% identity, while the minimum identity was evenly distributed in the range of 78% to 100%. Using a gene conversion detection test, we found frequent and/or recent gene conversion events within the tested Flowers. Interestingly, many of those converted regions contained protein-coding genes. Computer simulation studies suggest that one role of such frequent gene conversions is the elongation of the life span of gene families in a Flower by the resurrection of pseudogenes.