Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems.

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

An endogenous retrovirus presumed to have been endogenized or relocated recently in a marsupial, the red-necked wallaby.

An albino infant wallaby was born to a mother with wild-type body color. PCR and sequencing analyses of TYR (encoding tyrosinase, which is essential for melanin biosynthesis) of this albino wallaby revealed a 7.1-kb-long DNA fragment inserted in the first exon. Since the fragment carried long terminal repeats, we assumed it to be a copy of an endogenous retrovirus, which we named walb. We cloned other walb copies residing in the genomes of this species and of another wallaby species. The copies exhibited length variation, and the longest copy (>8.0 kb) contained open reading frames whose deduced amino acid sequences were well aligned with those of gag, pol, and env of retroviruses. It is unknown through which of the following likely processes the walb copy was inserted into TYR: endogenization (infection of a germline cell by an exogenous virus), reinfection (infection by a virus produced from a previously endogenized provirus), or retrotransposition (intracellular relocation of a provirus). In any case, the insertion into TYR is considered to have been a recent event on an evolutionary timescale because albino mutant alleles generally do not persist for long because of their deleterious effects in wild circumstances.

Complex intragene deletion leads to oculocutaneous albinism in tanuki (Japanese raccoon dog).

Tanuki (Nyctereutes procyonoides viverrinus), or Japanese raccoon dog, is a canine native to Japan. Tanuki with complete oculocutaneous albinism are relatively frequent in mountainous areas of mainland Japan. Tyrosinase, which is encoded by the TYR gene, is an enzyme essential for the biosynthesis of melanin pigment. We examined the structure and nucleotide sequence of TYR in an albino tanuki and found that the third exon was removed due to a deletion of approximately 11 kb. In addition, two nonsynonymous nucleotide substitutions were found in the fifth exon. These mutations are possible causes of the albino phenotype; however, the order of occurrence is unclear. Even if the 11-kb deletion was not the first of these mutations, it is considered to cause a total loss of the tyrosinase function because the third exon carries codons for one of the two copper-binding sites of tyrosinase and these sites are essential for the enzyme function. Intriguingly, the deletion was not a simple removal of an 11-kb segment: an internal portion was retained as a segment in the reverse orientation. We propose possible formation processes for this mutation that involve multiple DNA scission events, or an inversion followed by a deletion.

The evolutionary process of mammalian sex determination genes focusing on marsupial SRYs.

BACKGROUND: Maleness in mammals is genetically determined by the Y chromosome. On the Y chromosome SRY is known as the mammalian male-determining gene. Both placental mammals (Eutheria) and marsupial mammals (Metatheria) have SRY genes. However, only eutherian SRY genes have been empirically examined by functional analyses, and the involvement of marsupial SRY in male gonad development remains speculative. RESULTS: In order to demonstrate that the marsupial SRY gene is similar to the eutherian SRY gene in function, we first examined the sequence differences between marsupial and eutherian SRY genes. Then, using a parsimony method, we identify 7 marsupial-specific ancestral substitutions, 13 eutherian-specific ancestral substitutions, and 4 substitutions that occurred at the stem lineage of therian SRY genes. A literature search and molecular dynamics computational simulations support that the lineage-specific ancestral substitutions might be involved with the functional differentiation between marsupial and eutherian SRY genes. To address the function of the marsupial SRY gene in male determination, we performed luciferase assays on the testis enhancer of Sox9 core (TESCO) using the marsupial SRY. The functional assay shows that marsupial SRY gene can weakly up-regulate the luciferase expression via TESCO. CONCLUSIONS: Despite the sequence differences between the marsupial and eutherian SRY genes, our functional assay indicates that the marsupial SRY gene regulates SOX9 as a transcription factor in a similar way to the eutherian SRY gene. Our results suggest that SRY genes obtained the function of male determination in the common ancestor of Theria (placental mammals and marsupials). This suggests that the marsupial SRY gene has a function in male determination, but additional experiments are needed to be conclusive.

The Reliability and Stability of an Inferred Phylogenetic Tree from Empirical Data.

The reliability of a phylogenetic tree obtained from empirical data is usually measured by the bootstrap probability (Pb) of interior branches of the tree. If the bootstrap probability is high for most branches, the tree is considered to be reliable. If some interior branches show relatively low bootstrap probabilities, we are not sure that the inferred tree is really reliable. Here, we propose another quantity measuring the reliability of the tree called the stability of a subtree. This quantity refers to the probability of obtaining a subtree (Ps) of an inferred tree obtained. We then show that if the tree is to be reliable, both Pb and Ps must be high. We also show that Ps is given by a bootstrap probability of the subtree with the closest outgroup sequence, and computer program RESTA for computing the Pb and Ps values will be presented.

FGFR2 mutation in 46,XY sex reversal with craniosynostosis.

Patients with 46,XY gonadal dysgenesis (GD) exhibit genital anomalies, which range from hypospadias to complete male-to-female sex reversal. However, a molecular diagnosis is made in only 30% of cases. Heterozygous mutations in the human FGFR2 gene cause various craniosynostosis syndromes including Crouzon and Pfeiffer, but testicular defects were not reported. Here, we describe a patient whose features we would suggest represent a new FGFR2-related syndrome, craniosynostosis with XY male-to-female sex reversal or CSR. The craniosynostosis patient was chromosomally XY, but presented as a phenotypic female due to complete GD. DNA sequencing identified the FGFR2c heterozygous missense mutation, c.1025G>C (p.Cys342Ser). Substitution of Cys342 by Ser or other amino acids (Arg/Phe/Try/Tyr) has been previously reported in Crouzon and Pfeiffer syndrome. We show that the 'knock-in' Crouzon mouse model Fgfr2c(C342Y/C342Y) carrying a Cys342Tyr substitution displays XY gonadal sex reversal with variable expressivity. We also show that despite FGFR2c-Cys342Tyr being widely considered a gain-of-function mutation, Cys342Tyr substitution in the gonad leads to loss of function, as demonstrated by sex reversal in Fgfr2c(C342Y/-) mice carrying the knock-in allele on a null background. The rarity of our patient suggests the influence of modifier genes which exacerbated the testicular phenotype. Indeed, patient whole exome analysis revealed several potential modifiers expressed in Sertoli cells at the time of testis determination in mice. In summary, this study identifies the first FGFR2 mutation in a 46,XY GD patient. We conclude that, in certain rare genetic contexts, maintaining normal levels of FGFR2 signaling is important for human testis determination.

Sex chromosome cycle as a mechanism of stable sex determination.

Recent advances in DNA sequencing technology have enabled the precise decoding of genomes in non-model organisms, providing a basis for unraveling the patterns and mechanisms of sex chromosome evolution. Studies of different species have yielded conflicting results regarding the traditional theory that sex chromosomes evolve from autosomes via the accumulation of deleterious mutations and degeneration of the Y (or W) chromosome. The concept of the 'sex chromosome cycle,' emerging from this context, posits that at any stage of the cycle (i.e., differentiation, degeneration, or loss), sex chromosome turnover can occur while maintaining stable sex determination. Thus, understanding the mechanisms that drive both the persistence and turnover of sex chromosomes at each stage of the cycle is crucial. In this review, we integrate recent findings on the mechanisms underlying maintenance and turnover, with a special focus on several organisms having unique sex chromosomes. Our review suggests that the diversity of sex chromosomes in the maintenance of stable sex determination is underappreciated and emphasizes the need for more research on the sex chromosome cycle.

Genomic differentiation within East Asian Helicobacter pylori.

The East Asian region, including China, Japan and Korea, accounts for half of gastric cancer deaths. However, different areas have contrasting gastric cancer incidences and the population structure of Helicobacter pylori in this ethnically diverse region is yet unknown. We aimed to investigate genomic differences in H. pylori between these areas to identify sequence polymorphisms associated with increased cancer risk. We analysed 381 H. pylori genomes collected from different areas of the three countries using phylogenetic and population genetic tools to characterize population differentiation. The functional consequences of SNPs with a highest fixation index (Fst) between subpopulations were examined by mapping amino acid changes on 3D protein structure, solved or modelled. Overall, 329/381 genomes belonged to the previously identified hspEAsia population indicating that import of bacteria from other regions of the world has been uncommon. Seven subregional clusters were found within hspEAsia, related to subpopulations with various ethnicities, geographies and gastric cancer risks. Subpopulation-specific amino acid changes were found in multidrug exporters (hefC), transporters (frpB-4), outer membrane proteins (hopI) and several genes involved in host interaction, such as a catalase site, involved in H2O2 entrance, and a flagellin site mimicking host glycosylation. Several of the top hits, including frpB-4, hefC, alpB/hopB and hofC, have been found to be differentiated within the Americas in previous studies, indicating that a handful of genes may be key to local geographic adaptation. H. pylori within East Asia are not homogeneous but have become differentiated geographically at multiple loci that might have facilitated adaptation to local conditions and hosts. This has important implications for further evaluation of these changes in relation to the varying gastric cancer incidence between geographical areas in this region.

Comparative genomics of Glandirana rugosa using unsupervised AI reveals a high CG frequency.

The Japanese wrinkled frog (Glandirana rugosa) is unique in having both XX-XY and ZZ-ZW types of sex chromosomes within the species. The genome sequencing and comparative genomics with other frogs should be important to understand mechanisms of turnover of sex chromosomes within one species or during a short period. In this study, we analyzed the newly sequenced genome of G. rugosa using a batch-learning self-organizing map which is unsupervised artificial intelligence for oligonucleotide compositions. To clarify genome characteristics of G. rugosa, we compared its short oligonucleotide compositions in all 1-Mb genomic fragments with those of other six frog species (Pyxicephalus adspersus, Rhinella marina, Spea multiplicata, Leptobrachium leishanense, Xenopus laevis, and Xenopus tropicalis). In G. rugosa, we found an Mb-level large size of repeat sequences having a high identity with the W chromosome of the African bullfrog (P. adspersus). Our study concluded that G. rugosa has unique genome characteristics with a high CG frequency, and its genome is assumed to heterochromatinize a large size of genome via methylataion of CG.

The Reproducibility of an Inferred Tree and the Diploidization of Gene Segregation after Genome Duplication.

We previously introduced a numerical quantity called the stability (Ps) of an inferred tree and showed that for the tree to be reliable this stability as well as the reliability of the tree, which is usually computed as the bootstrap probability (Pb), must be high. However, if genome duplication occurs in a species, a gene family of the genome also duplicates, and for this reason alone some Ps values can be high in a tree of the duplicated gene families. In addition, the topology of the duplicated gene family can be similar to that of the original gene family if such gene families are identifiable. After genome duplication, however, the gene families are often partially deleted or partially duplicated, and the duplicated gene family may not show the same topology as that of the original family. It is therefore necessary to compute the similarity of the topologies of the duplicated and the original gene families. In this paper, we introduce another quantity called the reproducibility (Pr) for measuring the similarity of the two gene families. To show how to compute the Pr values, we first compute the Pb and Ps values for each of the MHC class II α and ß chain gene families, which were apparently generated by genome duplication. We then compute the Pr values for the MHC class II α and ß chain gene families. The Pr values for the α and ß chain gene families are now low, and this suggests that the diploidization of gene segregation has occurred after the genome duplication. Currently higher animals, defined as animals with complex phenotypic characters, generally have a higher genome size, and this increase in genome size appears to have been caused by genome duplication and diploidization of gene segregation after genome duplication.

Heterogeneity of synonymous substitution rates in the Xenopus frog genome.

With the increasing availability of high quality genomic data, there is opportunity to deeply explore the genealogical relationships of different gene loci between closely related species. In this study, we utilized genomes of Xenopus laevis (XLA, a tetraploid species with (L) and (S) sub-genomes) and X. tropicalis (XTR, a diploid species) to investigate whether synonymous substitution rates among orthologous or homoeologous genes displayed any heterogeneity. From over 1500 orthologous/homoeologous genes collected, we calculated proportion of synonymous substitutions between genomes/sub-genomes (k) and found variation within and between chromosomes. Within most chromosomes, we identified higher k with distance from the centromere, likely attributed to higher substitution rates and recombination in these regions. Using maximum likelihood methods, we identified further evidence supporting rate heterogeneity, and estimated species divergence times and ancestral population sizes. Estimated species divergence times (XLA.L-XLA.S: ~25.5 mya; XLA-XTR: ~33.0 mya) were slightly younger compared to a past study, attributed to consideration of population size in our study. Meanwhile, we found very large estimated population size in the ancestral populations of the two species (NA = 2.55 x 106). Local hybridization and population structure, which have not yet been well elucidated in frogs, may be a contributing factor to these possible large population sizes.

Evolutionary Medicine of Retroviruses in the Human Genome.

Humans are infected with many viruses, and the immune system mostly removes viruses and the infected cells. However, certain viruses have entered the human genome. Of the human genome, ∼45% is composed of transposable elements (long interspersed nuclear elements [LINEs], short interspersed nuclear elements [SINEs] and transposons) and 5-8% is derived from viral sequences with similarity to infectious retroviruses. If integration of retrovirus occurs in a germline, the integrated viral sequences are heritable. Accumulation of viral sequences has created the current human genome. This article summarizes recent studies of retroviruses in humans and bridges clinical fields and evolutionary genetics. First, we report the repertories of human-infective retroviruses. Second, we review endogenous retroviruses in the human genome and diseases associated with endogenous retroviruses. Third, we discuss the biological functions of endogenous retroviruses and propose the concept of accelerated human evolution via viruses. Finally, we present perspectives of virology in the field of evolutionary medicine.

In vivo function and evolution of the eutherian-specific pluripotency marker UTF1.

Embryogenesis in placental mammals is sustained by exquisite interplay between the embryo proper and placenta. UTF1 is a developmentally regulated gene expressed in both cell lineages. Here, we analyzed the consequence of loss of the UTF1 gene during mouse development. We found that homozygous UTF1 mutant newborn mice were significantly smaller than wild-type or heterozygous mutant mice, suggesting that placental insufficiency caused by the loss of UTF1 expression in extra-embryonic ectodermal cells at least in part contributed to this phenotype. We also found that the effects of loss of UTF1 expression in embryonic stem cells on their pluripotency were very subtle. Genome structure and sequence comparisons revealed that the UTF1 gene exists only in placental mammals. Our analyses of a family of genes with homology to UTF1 revealed a possible mechanism by which placental mammals have evolved the UTF1 genes.

No evidence for a second evolutionary stratum during the early evolution of mammalian sex chromosomes.

Mammalian sex chromosomes originated from a pair of autosomes, and homologous genes on the sex chromosomes (gametologs) differentiated through recombination arrest between the chromosomes. It was hypothesized that this differentiation in eutherians took place in a stepwise fashion and left a footprint on the X chromosome termed "evolutionary strata." The evolutionary stratum hypothesis claims that strata 1 and 2 (which correspond to the first two steps of chromosomal differentiation) were generated in the stem lineage of Theria or before the divergence between eutherians and marsupials. However, this prediction relied solely on the molecular clock hypothesis between pairs of human gametologs, and molecular evolution of marsupial sex chromosomal genes has not yet been investigated. In this study, we analyzed the following 7 pairs of marsupial gametologs, together with their eutherian orthologs that reside in stratum 1 or 2: SOX3/SRY, RBMX/Y, RPS4X/Y, HSFX/Y, XKRX/Y, SMCX/Y (KDM5C/D, JARID1C/D), and UBE1X/Y (UBA1/UBA1Y). Phylogenetic analyses and estimated divergence time of these gametologs reveal that they all differentiated at the same time in the therian ancestor. We have also provided strong evidence for gene conversion that occurred in the 3' region of the eutherian stratum 2 genes (SMCX/Y and UBE1X/Y). The results of the present study show that (1) there is no compelling evidence for the second stratum in the stem lineage of Theria; (2) gene conversion, which may have occurred between SMCX/Y and UBE1X/Y in the eutherian lineage, potentially accounts for their apparently lower degree of overall divergence.

Evolution of genomic structures on Mammalian sex chromosomes.

Throughout mammalian evolution, recombination between the two sex chromosomes was suppressed in a stepwise manner. It is thought that the suppression of recombination led to an accumulation of deleterious mutations and frequent genomic rearrangements on the Y chromosome. In this article, we review three evolutionary aspects related to genomic rearrangements and structures, such as inverted repeats (IRs) and palindromes (PDs), on the mammalian sex chromosomes. First, we describe the stepwise manner in which recombination between the X and Y chromosomes was suppressed in placental mammals and discuss a genomic rearrangement that might have led to the formation of present pseudoautosomal boundaries (PAB). Second, we describe ectopic gene conversion between the X and Y chromosomes, and propose possible molecular causes. Third, we focus on the evolutionary mode and timing of PD formation on the X and Y chromosomes. The sequence of the chimpanzee Y chromosome was recently published by two groups. Both groups suggest that rapid evolution of genomic structure occurred on the Y chromosome. Our re-analysis of the sequences confirmed the species-specific mode of human and chimpanzee Y chromosomal evolution. Finally, we present a general outlook regarding the rapid evolution of mammalian sex chromosomes.

Identification of an enhancer in the Ad4BP/SF-1 gene specific for fetal Leydig cells.

Adrenal 4 binding protein/steroidogenic factor 1 (Ad4BP/SF-1) (Nr5a1) is a nuclear receptor essential for reproductive tissue development and endocrine regulation. This factor is expressed in steroidogenic tissues (e.g. adrenal glands and gonads), and expression of this factor is tightly regulated in a tissue and cell type-specific manner. Our previous studies have identified tissue and cell type-specific enhancers in the introns of the Ad4BP/SF-1 gene in fetal adrenal glands, ventromedial hypothalamus, and pituitary gonadotrope. Characterization of the enhancers had provided new insights into tissue and cell development. However, these studies have failed to identify any gonad-specific enhancer. Here, we identified a fetal Leydig cell-specific enhancer in the upstream region of the mouse Ad4BP/SF-1 gene using transgenic mouse assays. Alignment of the upstream regions among vertebrate animal species demonstrated that the enhancer consisted of three conserved regions, whereby the most highly conserved region contained an Ad4BP/SF-1 binding sequence and an E-box. Mutation of each sequence abolished the enhancer activity and led to a loss of reporter gene expression. These results suggested that Ad4BP/SF-1 gene expression in the fetal Leydig cell is regulated by a yet unidentified E-box binding protein(s) and by an autoregulatory loop formed by Ad4BP/SF-1. Although fetal Leydig cells have been thought to play crucial roles for masculinization of various fetal tissues through androgen production, other functions have remained elusive. Our identification of a fetal Leydig cell-specific enhancer in the Ad4BP/SF-1 gene would be a powerful tool to address these gaps in the knowledge base.

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.

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.