Spatial and temporal diversity of positive selection on shared haplotypes at the PSCA locus among worldwide human populations.

Selection on standing genetic variation is important for rapid local genetic adaptation when the environment changes. We report that, for the prostate stem cell antigen (PSCA) gene, different populations have different target haplotypes, even though haplotypes are shared among populations. The C-C-A haplotype, whereby the first C is located at rs2294008 of PSCA and is a low risk allele for gastric cancer, has become a target of positive selection in Asia. Conversely, the C-A-G haplotype carrying the same C allele has become a selection target mainly in Africa. However, Asian and African share both haplotypes, consistent with the haplotype divergence time (170 kya) prior to the out-of-Africa dispersal. The frequency of C-C-A/C-A-G is 0.344/0.278 in Asia and 0.209/0.416 in Africa. Two-dimensional site frequency spectrum analysis revealed that the extent of intra-allelic variability of the target haplotype is extremely small in each local population, suggesting that C-C-A or C-A-G is under ongoing hard sweeps in local populations. From the time to the most recent common ancestor (TMRCA) of selected haplotypes, the onset times of positive selection were recent (3-55 kya), concurrently with population subdivision from a common ancestor. Additionally, estimated selection coefficients from ABC analysis were up to ~3%, similar to those at other loci under recent positive selection. Phylogeny of local populations and TMRCA of selected haplotypes revealed that spatial and temporal switching of positive selection targets is a unique and novel feature of ongoing selection at PSCA. This switching may reflect the potential of rapid adaptability to distinct environments.

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.

Development of a novel monoclonal antibody that binds to most HLA-A allomorphs in a conformation-dependent yet peptide-promiscuous fashion.

Specificity analyses of peptide binding to human leukocyte antigen (HLA)-A molecules have been hampered due to a lack of proper monoclonal antibodies (mAbs) for certain allomorphs, such as the prevalent HLA-A1 for Caucasians and HLA-A11 for Asians. We developed a mAb that recognizes a conformational epitope common to most HLA-A allomorphs. The mAb, named A-1, does not discriminate peptides by amino acid sequences, making it suitable for measuring peptide binding. A stabilization assay using TAP-deficient cell lines and A-1 was developed to investigate the specificity of peptide binding to HLA-A molecules. Regarding the evolution of HLA-A genes, the A-1 epitope has been conserved among most HLA-A allomorphs but was lost when the HLA-A gene diversified into the HLA-A*32, HLA-A*31, and HLA-A*33 lineages together with HLA-A*29 after bifurcating from the HLA-A*25 and HLA-A*26 branchs. The establishment of A-1 is expected to help researchers investigate the peptide repertoire and develop computational tools to identify cognate peptides. Since no HLA-A locus-specific mAb has been available, A-1 will also be useful for analyzing the locus-specific regulation of the HLA gene expression.

Functional Evolution of Avian RIG-I-Like Receptors.

RIG-I-like receptors (retinoic acid-inducible gene-I-like receptors, or RLRs) are family of pattern-recognition receptors for RNA viruses, consisting of three members: retinoic acid-inducible gene I (RIG-I), melanoma differentiation-associated gene 5 (MDA5) and laboratory of genetics and physiology 2 (LGP2). To understand the role of RLRs in bird evolution, we performed molecular evolutionary analyses on the coding genes of avian RLRs using filtered predicted coding sequences from 62 bird species. Among the three RLRs, conservation score and dN/dS (ratio of nonsynonymous substitution rate over synonymous substitution rate) analyses indicate that avian MDA5 has the highest conservation level in the helicase domain but a lower level in the caspase recruitment domains (CARDs) region, which differs from mammals; LGP2, as a whole gene, has a lower conservation level than RIG-I or MDA5. We found evidence of positive selection across all bird lineages in RIG-I and MDA5 but only on the stem lineage of Galliformes in LGP2, which could be related to the loss of RIG-I in Galliformes. Analyses also suggest that selection relaxation may have occurred in LGP2 during the middle of bird evolution and the CARDs region of MDA5 contains many positively selected sites, which might explain its conservation level. Spearman's correlation test indicates that species-to-ancestor dN/dS of RIG-I shows a negative correlation with endogenous retroviral abundance in bird genomes, suggesting the possibility of interaction between immunity and endogenous retroviruses during bird evolution.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II.

BACKGROUND: In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs. RESULTS: We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago's brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens. CONCLUSIONS: This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

A pre-metazoan origin of the CRK gene family and co-opted signaling network.

CRK and CRKL adapter proteins play essential roles in development and cancer through their SRC homology 2 and 3 (SH2 and SH3) domains. To gain insight into the origin of their shared functions, we have investigated their evolutionary history. We propose a term, crk/crkl ancestral (crka), for orthologs in invertebrates before the divergence of CRK and CRKL in the vertebrate ancestor. We have isolated two orthologs expressed in the choanoflagellate Monosiga brevicollis, a unicellular relative to the metazoans. Consistent with its highly-conserved three-dimensional structure, the SH2 domain of M. brevicollis crka1 can bind to the mammalian CRK/CRKL SH2 binding consensus phospho-YxxP, and to the SRC substrate/focal adhesion protein BCAR1 (p130CAS) in the presence of activated SRC. These results demonstrate an ancient origin of the CRK/CRKL SH2-target recognition specificity. Although BCAR1 orthologs exist only in metazoans as identified by an N-terminal SH3 domain, YxxP motifs, and a C-terminal FAT-like domain, some pre-metazoan transmembrane proteins include several YxxP repeats in their cytosolic region, suggesting that they are remotely related to the BCAR1 substrate domain. Since the tyrosine kinase SRC also has a pre-metazoan origin, co-option of BCAR1-related sequences may have rewired the crka-dependent network to mediate adhesion signals in the metazoan ancestor.

Nonsynonymous substitution rate heterogeneity in the peptide-binding region among different HLA-DRB1 lineages in humans.

An extraordinary diversity of amino acid sequences in the peptide-binding region (PBR) of human leukocyte antigen [HLA; human major histocompatibility complex (MHC)] molecules has been maintained by balancing selection. The process of accumulation of amino acid diversity in the PBR for six HLA genes (HLA-A, B, C, DRB1, DQB1, and DPB1) shows that the number of amino acid substitutions in the PBR among alleles does not linearly correlate with the divergence time of alleles at the six HLA loci. At these loci, some pairs of alleles show significantly less nonsynonymous substitutions at the PBR than expected from the divergence time. The same phenomenon was observed not only in the HLA but also in the rat MHC. To identify the cause for this, DRB1 sequences, a representative case of a typical nonlinear pattern of substitutions, were examined. When the amino acid substitutions in the PBR were placed with maximum parsimony on a maximum likelihood tree based on the non-PBR substitutions, heterogeneous rates of nonsynonymous substitutions in the PBR were observed on several branches. A computer simulation supported the hypothesis that allelic pairs with low PBR substitution rates were responsible for the stagnation of accumulation of PBR nonsynonymous substitutions. From these observations, we conclude that the nonsynonymous substitution rate at the PBR sites is not constant among the allelic lineages. The deceleration of the rate may be caused by the coexistence of certain pathogens for a substantially long time during HLA evolution.

Evolutionary history of the cancer immunity antigen MAGE gene family.

The evolutionary mode of a multi-gene family can change over time, depending on the functional differentiation and local genomic environment of family members. In this study, we demonstrate such a change in the melanoma antigen (MAGE) gene family on the mammalian X chromosome. The MAGE gene family is composed of ten subfamilies that can be categorized into two types. Type I genes are of relatively recent origin, and they encode epitopes for human leukocyte antigen (HLA) in cancer cells. Type II genes are relatively ancient and some of their products are known to be involved in apoptosis or cell proliferation. The evolutionary history of the MAGE gene family can be divided into four phases. In phase I, a single-copy state of an ancestral gene and the evolutionarily conserved mode had lasted until the emergence of eutherian mammals. In phase II, eight subfamily ancestors, with the exception for MAGE-C and MAGE-D subfamilies, were formed via retrotransposition independently. This would coincide with a transposition burst of LINE elements at the eutherian radiation. However, MAGE-C was generated by gene duplication of MAGE-A. Phase III is characterized by extensive gene duplication within each subfamily and in particular the formation of palindromes in the MAGE-A subfamily, which occurred in an ancestor of the Catarrhini. Phase IV is characterized by the decay of a palindrome in most Catarrhini, with the exception of humans. Although the palindrome is truncated by frequent deletions in apes and Old World monkeys, it is retained in humans. Here, we argue that this human-specific retention stems from negative selection acting on MAGE-A genes encoding epitopes of cancer cells, which preserves their ability to bind to highly divergent HLA molecules. These findings are interpreted with consideration of the biological factors shaping recent human MAGE-A genes.

Identification of CTL epitopes in hepatitis C virus by a genome-wide computational scanning and a rational design of peptide vaccine.

Developing a peptide-based vaccine for the highly variable hepatitis C virus (HCV) remains a challenging task. Variant viruses not only escape antigen presentation but also persist in a patient as quasi-species. Such variants are often antagonistic to the responding T cell repertoire. To overcome these problems, we herein propose a cocktail vaccine consisting of a few epitope peptides, which make it possible to outpace the emergence of variant viruses. To design such a vaccine, we developed a way to identify HLA-A*2402-binding peptides efficiently by means of the computational scanning of the whole genome of the pathogen. Most of the predicted peptides exhibited strong binding to the HLA-A*2402 molecule, while also inducing CD8 T cell responses from the patients' peripheral blood mononuclear cells (PBMCs). Peptide-induced T cells were capable of lysing HCV-expressing HepG2 cells which process antigens endogenously. The amount of HCV core antigen in the patients' livers suggested that the lytic activity of the peptide-induced T cells was clearly in a range suitable for therapeutic use. If T cells were activated under optimal conditions by high density peptides, then they tended to be relatively tolerant of single amino acid variations for cytolysis. Finally, an analysis of the viral population isolated in Japan suggested no obvious changes due to immune evasion in the viral genome even in a host population highly biased toward HLA-A*2402.

The dominance of alleles controlling self-incompatibility in Brassica pollen is regulated at the RNA level.

Self-incompatibility (SI) in Brassica is controlled sporophytically by the multiallelic S-locus. The SI phenotype of pollen in an S-heterozygote is determined by the relationship between the two S-haplotypes it carries, and dominant/recessive relationships often are observed between the two S-haplotypes. The S-locus protein 11 (SP11, also known as the S-locus cysteine-rich protein) gene has been cloned from many pollen-dominant S-haplotypes (class I) and shown to encode the pollen S-determinant. However, SP11 from pollen-recessive S-haplotypes (class II) has never been identified by homology-based cloning strategies, and how the dominant/recessive interactions between the two classes occur was not known. We report here the identification and molecular characterization of SP11s from six class II S-haplotypes of B. rapa and B. oleracea. Phylogenetic analysis revealed that the class II SP11s form a distinct group separated from class I SP11s. The promoter sequences and expression patterns of SP11s also were different between the two classes. The mRNA of class II SP11, which was detected predominantly in the anther tapetum in homozygotes, was not detected in the heterozygotes of class I and class II S-haplotypes, suggesting that the dominant/recessive relationships of pollen are regulated at the mRNA level of SP11s.