Analysis of mutation rates in the SMCY/SMCX genes shows that mammalian evolution is male driven

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

Sexually Dimorphic Neurosteroid Synthesis Regulates Neuronal Activity in the Murine Brain

Sex steroid hormones act on hypothalamic kisspeptin neurons to regulate reproductive neural circuits in the brain. Kisspeptin neurons start to express estrogen receptors in utero, suggesting steroid hormone action on these cells early during development. Whether neurosteroids are locally produced in the embryonic brain and impinge onto kisspeptin/reproductive neural circuitry is not known. To address this question, we analyzed aromatase expression, a key enzyme in estrogen synthesis, in male and female mouse embryos. We identified an aromatase neuronal network comprising ∼6000 neurons in the hypothalamus and amygdala. By birth, this network has become sexually dimorphic in a cluster of aromatase neurons in the arcuate nucleus adjacent to kisspeptin neurons. We demonstrate that male arcuate aromatase neurons convert testosterone to estrogen to regulate kisspeptin neuron activity. We provide spatiotemporal information on aromatase neuronal network development and highlight a novel mechanism whereby aromatase neurons regulate the activity of distinct neuronal populations expressing estrogen receptors.SIGNIFICANCE STATEMENT Sex steroid hormones, such as estradiol, are important regulators of neural circuits controlling reproductive physiology in the brain. Embryonic kisspeptin neurons in the hypothalamus express steroid hormone receptors, suggesting hormone action on these cells in utero Whether neurosteroids are locally produced in the brain and impinge onto reproductive neural circuitry is insufficiently understood. To address this question, we analyzed aromatase expression, a key enzyme in estradiol synthesis, in mouse embryos and identified a network comprising ∼6000 neurons in the brain. By birth, this network has become sexually dimorphic in a cluster of aromatase neurons in the arcuate nucleus adjacent to kisspeptin neurons. We demonstrate that male aromatase neurons convert testosterone to estradiol to regulate kisspeptin neuron activity.

Prenatal low-dose methyltestosterone, but not dihydrotestosterone, treatment induces penile formation in female mice and guinea pigs†

Genital tubercle has bisexual potential before sex differentiation. Females exposed to androgen during sex differentiation show masculinized external genitalia, but the effects of different androgens on tubular urethral and penile formation in females are mostly unknown. In this study, we compared the masculinization effects of commonly used androgens methyltestosterone, dihydrotestosterone, and testosterone on the induction of penile formation in females. Our results suggested that prenatal treatment with low doses of methyltestosterone, but not same doses of dihydrotestosterone or testosterone, could induce penile formation in female mice. The minimum dose of dihydrotestosterone and testosterone for inducing tubular urethral formation in female mice was, respectively, 50 and 20 times higher than that of methyltestosterone. In vivo methyltestosterone treatment induced more nuclear translocation of androgen receptors in genital tubercles of female mice, affected Wnt signaling gene expressions, and then led to similar patterns of cell proliferation and death in developing genital tubercles to those of control males. We further revealed that low-dose methyltestosterone, but not same dose of dihydrotestosterone or testosterone, treatment induced penile formation in female guinea pigs. Exposure of female mouse genital tubercle organ culture to methyltestosterone, dihydrotestosterone, or testosterone could induce nuclear translocation of androgen receptors, suggesting that the differential effect of the three androgens in vivo might be due to the hormonal profile in mother or fetus, rather than the local genital tissue. To understand the differential role of these androgens in masculinization process involved is fundamental to androgen replacement therapy for diseases related to external genital masculinization.

Mitotic antipairing of homologous and sex chromosomes via spatial restriction of two haploid sets

Mitotic recombination must be prevented to maintain genetic stability across daughter cells, but the underlying mechanism remains elusive. We report that mammalian cells impede homologous chromosome pairing during mitosis by keeping the two haploid chromosome sets apart, positioning them to either side of a meridional plane defined by the centrosomes. Chromosome oscillation analysis revealed collective genome behavior of noninteracting chromosome sets. Male translocation mice with a maternal-derived supernumerary chromosome display the tracer chromosome exclusively to the haploid set containing the X chromosome. This haploid set-based antipairing motif is shared by multiple cell types, is doubled in tetraploid cells, and is lost in carcinoma cells. The data provide a model of nuclear polarity through the antipairing of homologous chromosomes during mitosis.

Pairing homologous chromosomes is required for recombination. However, in nonmeiotic stages it can lead to detrimental consequences, such as allelic misregulation and genome instability, and is rare in human somatic cells. How mitotic recombination is prevented—and how genetic stability is maintained across daughter cells—is a fundamental, unanswered question. Here, we report that both human and mouse cells impede homologous chromosome pairing by keeping two haploid chromosome sets apart throughout mitosis. Four-dimensional analysis of chromosomes during cell division revealed that a haploid chromosome set resides on either side of a meridional plane, crossing two centrosomes. Simultaneous tracking of chromosome oscillation and the spindle axis, using fluorescent CENP-A and centrin1, respectively, demonstrates collective genome behavior/segregation of two haploid sets throughout mitosis. Using 3D chromosome imaging of a translocation mouse with a supernumerary chromosome, we found that this maternally derived chromosome is positioned by parental origin. These data, taken together, support the identity of haploid sets by parental origin. This haploid set-based antipairing motif is shared by multiple cell types, doubles in tetraploid cells, and is lost in a carcinoma cell line. The data support a mechanism of nuclear polarity that sequesters two haploid sets along a subcellular axis. This topological segregation of haploid sets revisits an old model/paradigm and provides implications for maintaining mitotic fidelity.

Permissiveness to form pluripotent stem cells may be an evolutionarily derived characteristic in Mus musculus

Mus musculus is the only known species from which embryonic stem cells (ESC) can be isolated under conditions requiring only leukemia inhibitory factor (LIF). Other species are non-permissive in LIF media, and form developmentally primed epiblast stem cells (EpiSC) similar to cells derived from post-implantation, egg cylinders. To evaluate whether non-permissiveness extends to induced pluripotent stem cells (iPSC), we derived iPSC from the eight founder strains of the mouse Collaborative Cross. Two strains, NOD/ShiLtJ and the WSB/EiJ, were non-permissive, consistent with the previous classification of NOD/ShiLtJ as non-permissive to ESC derivation. We determined non-permissiveness is recessive, and that non-permissive genomes do not compliment. We overcame iPSC non-permissiveness by using GSK3B and MEK inhibitors with serum, a technique we termed 2iS reprogramming. Although used for ESC derivation, GSK3B and MEK inhibitors have not been used during iPSC reprogramming because they inhibit survival of progenitor differentiated cells. iPSC derived in 2iS are more transcriptionally similar to ESC than EpiSC, indicating that 2iS reprogramming acts to overcome genetic background constraints. Finally, of species tested for ESC or iPSC derivation, only some M. musculus strains are permissive under LIF culture conditions suggesting that this is an evolutionarily derived characteristic in the M. musculus lineage.

Transcriptional activity of oestrogen receptors in the course of embryo development

Oestrogens are well-known proliferation and differentiation factors that play an essential role in the correct development of sex-related organs and behaviour in mammals. With the use of the ERE-Luc reporter mouse model, we show herein that throughout mouse development, oestrogen receptors (ERs) are active starting from day 12 post conception. Most interestingly, we show that prenatal luciferase expression in each organ is proportionally different in relation to the germ layer of the origin. The luciferase content is highest in ectoderm-derived organs (such as brain and skin) and is lowest in endoderm-derived organs (such as liver, lung, thymus and intestine). Consistent with the testosterone surge occurring in male mice at the end of pregnancy, in the first 2 days after birth, we observed a significant increase in the luciferase content in several organs, including the liver, bone, gonads and hindbrain. The results of the present study show a widespread transcriptional activity of ERs in developing embryos, pointing to the potential contribution of these receptors in the development of non-reproductive as well as reproductive organs. Consequently, the findings reported here might be relevant in explaining the significant differences in male and female physiopathology reported by a growing number of studies and may underline the necessity for more systematic analyses aimed at the identification of the prenatal effects of drugs interfering with ER signalling, such as aromatase inhibitors or endocrine disrupter chemicals.

Evolutionary Dynamics of the Gametologous CTNNB1 Gene on the Z and W Chromosomes of Snakes

Snakes exhibit genotypic sex determination with female heterogamety (ZZ males and ZW females), and the state of sex chromosome differentiation also varies among lineages. To investigate the evolutionary history of homologous genes located in the nonrecombining region of differentiated sex chromosomes in snakes, partial sequences of the gametologous CTNNB1 gene were analyzed for 12 species belonging to henophid (Cylindrophiidae, Xenopeltidae, and Pythonidae) and caenophid snakes (Viperidae, Elapidae, and Colubridae). Nonsynonymous/synonymous substitution ratios (Ka/Ks) in coding sequences were low (Ka/Ks < 1) between CTNNB1Z and CTNNB1W, suggesting that these 2 genes may have similar functional properties. However, frequencies of intron sequence substitutions and insertion–deletions were higher in CTNNB1Z than CTNNB1W, suggesting that Z-linked sequences evolved faster than W-linked sequences. Molecular phylogeny based on both intron and exon sequences showed the presence of 2 major clades: 1) Z-linked sequences of Caenophidia and 2) W-linked sequences of Caenophidia clustered with Z-linked sequences of Henophidia, which suggests that the sequence divergence between CTNNB1Z and CTNNB1W in Caenophidia may have occurred by the cessation of recombination after the split from Henophidia.

Vitamin D-deficiency and sex-specific dysregulation of placental inflammation

To investigate an immunomodulatory role for vitamin D in pregnancy we used mice raised on vitamin D-sufficient (SUFF), or -deficient (DEF) diets. At embryonic day 14, pregnant mice received intraperitoneal injection of lipopolysaccharide (LPS) or vehicle for 24h, with age-matched non-pregnant mice as controls. In non-pregnant mice, 6 serum analytes (IL-1β, IL-18, MDC/CCL22, MIP-1α/CCL3, EGF, IgA) were lower in DEF mice. In pregnant DEF mice only GH was higher. In non-pregnant mice LPS induced 28 analytes, with 5 (IL-18, IP-10/CXCL10, MCP-1/CCL2, MIP-1β/CCL4, MIP-3β/CCL19) being highest in DEF mice. In pregnant SUFF mice 16 serum analytes increased with LPS, and 6 of these (IP-10/CXCL10, MCP-1/CCL2, SAP, TIMP-1, VCAM-1, vWF) were higher and 1 (GCP-2/CXCL6) lower in DEF mice. Parallel analysis of placental mRNAs showed elevated mRNA for Il-6, Ccl2 and Cxcl10 in placentae from male and female fetuses in LPS-DEF mice. However, LPS-induced expression of Ifnγ, Tnfα, and Cxcl6 was only observed in female placentae from DEF mice. LPS-DEF mice also showed smaller litter sizes relative to control SUFF mice. Numbers of female fetuses per dam were significantly lower for DEF mice with or without LPS challenge. LPS had no effect on numbers of male fetuses from DEF mothers, but significantly decreased male fetuses from SUFF mothers. These data indicate that vitamin D is an important component of anti-inflammatory immune responses during pregnancy, with the placenta and fetal sex playing pivotal roles in this process.

Investigation of Histone Lysine-Specific Demethylase 5D KDM5D) Isoform Expression in Prostate Cancer Cell Lines: a System Approach

Background:

It is now well-demonstrated that histone demethylases play an important role in developmental controls, cell-fate decisions, and a variety of diseases such as cancer. Lysine-specific demethylase 5D (KDM5D) is a male-specific histone demethylase that specifically demethylates di- and tri-methyl H3K4 at the start site of active gene. In this light, the aim of this study was to investigate isoform/transcript-specific expression profiles of KDM5D in three prostate cancer cell lines, Du-145, LNCaP, and PC3.

Methods:

Real-time PCR analysis was performed to determine the expression levels of different KDM5D transcripts in the prostate cell lines. A gene regulatory network was established to analyze the gene expression profile.

Results:

Significantly different expression levels of both isoforms were found among the three cell lines. Interestingly, isoform I was expressed in three cell lines while isoform III did only in DU-145. The expression levels of both isoforms were higher in DU-145 when compared to other cell lines (P<0.0001). The observed expression profile was determined by using regulatory network analyses.

Conclusion:

The present study, for the first time, not only showed the expression profiles of KDM5D isoforms in prostate cancer cell lines but also evaluated the effects of the gene regulatory network on the expression profile of this gene.

Y-Single Nucleotide Polymorphisms Diversity in Chinese Indigenous Horse

In contrast to high genetic diversity of mitochondrial DNA (mtDNA), equine Y chromosome shows extremely low variability, implying limited patrilines in the domesticated horse. In this study, we applied direct sequencing and restriction fragment length polymorphism (RFLP) methods to investigate the polymorphisms of 33 Y chromosome specific loci in 304 Chinese indigenous horses from 13 breeds. Consequently, two Y-single nucleotide polymorphisms (SNPs) (Y-45701/997 and Y-50869) and one Y-indel (Y-45288) were identified. Of those, the Y-50869 (T>A) revealed the highest variation frequency (24.67%), whereas it was only 3.29% and 1.97% in Y-45288 (T/-) and Y-45701/997 (G>T) locus, respectively. These three mutations accounted for 27.96% of the total samples and identified five Y-SNP haplotypes, demonstrating genetic diversity of Y chromosome in Chinese horses. In addition, all the five Y-SNP haplotypes were shared by different breeds. Among 13 horse breeds analyzed, Balikun horse displayed the highest nucleotide diversity (π = 5.6×10(-4)) and haplotype diversity (h = 0.527), while Ningqiang horse showed the lowest nucleotide diversity (π = 0.00000) and haplotype diversity (h = 0.000). The results also revealed that Chinese horses had a different polymorphic pattern of Y chromosome from European and American horses. In conclusion, Chinese horses revealed genetic diversity of Y chromosome, however more efforts should be made to better understand the domestication and paternal origin of Chinese indigenous horses.

Murine arcuate nucleus kisspeptin neurons communicate with GnRH neurons in utero

Puberty is a transition period of reproductive development from juvenile stages to adulthood and depends upon the activity of gonadotropin-releasing hormone (GnRH) neurons. GnRH neurons are initially activated in utero but remain quiescent throughout the juvenile period. Premature reactivation of GnRH neurons results in precocious puberty in mice and humans, but the mechanisms underlying developmental control of GnRH neuron activity remain unknown. The neuropeptide kisspeptin, a potent activator of GnRH neurons that is implicated as a critical permissive signal triggering puberty and a major regulator of the adult female hypothalamus-pituitary-gonadal axis, is paradoxically produced by neurons in the developing brain well before puberty onset. Thus, the neural circuits controlling the timing of reproductive maturation remain elusive. Here, we delineate the underlying neural circuitry using conditional genetic transsynaptic tracing in female mouse embryos. We find that kisspeptin-producing neurons in the arcuate nucleus (ARC) already communicate with a specific subset of GnRH neurons in utero. We show that ARC kisspeptin neurons are upstream of GnRH neurons, and that GnRH neuron connectivity to ARC kisspeptin neurons does not depend on their spatial position in the brain. Furthermore, we demonstrate that the neural circuits between ARC kisspeptin and GnRH neurons are fully established and operative before birth. Finally, we find that most GnRH neurons express the kisspeptin receptor GPR54 upon circuit formation, suggesting that the signaling system implicated in gatekeeping puberty becomes operative in the embryo.

Rapid identification of targeted transgene integrations in ES cells by fluorescence detection

The generation of transgenic animals with a gain-of-function mutation is commonly achieved by procedures based on random DNA integration. The resulting transgenic founder lines are unique, not reproducible and have variable expression patterns. In contrast, the targeted integration of transgenes into a predetermined neutral genomic position solves most of the inadequacies of random integration methods. However, homologous recombination (HR) in mouse embryonic stem cells (ESCs) currently requires careful design of the targeting vector and a laborious procedure to identify clones with the correct insertion event. Here, we introduce a feasible strategy that employs a heterozygous double fluorescent reporter ESC line for simple identification of a knock-in HR event via detection of endogenous fluorescence expression. Following positive selection using antibiotics, the system offers a second selection step to identify targeted clones by the loss of one of two fluorescence reporters in lieu of the time consuming Southern blotting and PCR analysis routinely applied in conventional targeting experiments. Moreover, the method allows for the simple detection of chimerism (negating the need for appropriate coat colour combinations) and enables the early detection of germline transmission by fluorescence reporter expression in F1 neonates.

Sexually dimorphic placental development throughout gestation in the spiny mouse (Acomys cahirinus)

OBJECTIVES

It has been hypothesized that male fetuses down regulate placental growth during periods of accelerated fetal growth. We aimed to investigate this, and determine whether sexual dimorphism was apparent in the spiny mouse placenta. We hypothesized that expression of fetal growth promoters would be higher in placentas of males, whereas genes involved in placental structural development would be more highly expressed in placentas of females.

METHODS

Spiny mouse dams, a precocial rodent with an in utero endocrine milieu dissimilar from other rodents, but akin to humans, were sacrificed at gestational ages 15-37 (term = 39 days). Placentas were collected and processed for histology or qPCR analysis of selected genes (GCM1, MAP2K1, SLC2A1, NR3C1, IGF1, IGF1R).

RESULTS

Fetal and placental weights were similar for both sexes. Placentas of female fetuses had less spongy zone (P(SEX) < 0.0001), and more labyrinth (P(SEX) < 0.0001) than males. Early placenta and labyrinth expression of SLC2A1 was higher in males than females (P(SEX) < 0.05). Labyrinthine IGF1R remained constant until term in the female, compared with male where expression increased until term. Peak MAP2K1 expression occurred earlier in the male placenta than the female. Spongy zone SLC2A1 remained constant until term in the female, compared with male where expression increased until term.

CONCLUSIONS

The spiny mouse is a species that exhibits sexually dimorphic placental development. We suggest that these sex differences in placental gene expression and structure may underlie or compound the male vulnerability to a sub-optimal in utero environment.

Length difference between equine ZFX and ZFY genes and its application for molecular sex determination

PURPOSE

we analyzed the sex chromosome-encoding ZFX-ZFY genes and tested molecular sexing using the amplification patterns of intron 9 of ZFX-ZFY in the horse.

METHODS AND RESULTS

the amplification of the ZFX-ZFY produced two distinct patterns, reflecting sexual dimorphism based on a length difference between the X and Y chromosomes. The amplification products from foals showed two distinct bands: one was common to all foals and mares, indicating that this band was amplified from ZFX, while the other was specific to some foals, indicating that it was from ZFY. The result based on the PCR assay was identical to the results of amplification of the Y chromosome-specific SRY gene and those of investigations of the phenotypic gender in three different horse populations.

CONCLUSION

we suggest that this PCR strategy for determining sexes by comparing the amplification patterns of ZFX-ZFY genes is a convenient and precise method for discriminating sexes in horses.

Embryonic gonadotropin-releasing hormone signaling is necessary for maturation of the male reproductive axis

Gonadotropin-releasing hormone (GnRH) signaling regulates reproductive physiology in mammals. GnRH is released by a subset of hypothalamic neurons and binds to GnRH receptor (GnRHR) on gonadotropes in the anterior pituitary gland to control production and secretion of gonadotropins that in turn regulate the activity of the gonads. Central control of reproduction is well understood in adult animals, but GnRH signaling has also been implicated in the development of the reproductive axis. To investigate the role of GnRH signaling during development, we selectively ablated GnRHR-expressing cells in mice. This genetic strategy permitted us to identify an essential stage in male reproductive axis development, which depends on embryonic GnRH signaling. Our experiments revealed a striking dichotomy in the gonadotrope population of the fetal anterior pituitary gland. We show that luteinizing hormone-expressing gonadotropes, but not follicle-stimulating hormone-expressing gonadotropes, express the GnRHR at embryonic day 16.75. Furthermore, we demonstrate that an embryonic increase in luteinizing hormone secretion is needed to promote development of follicle-stimulating hormone-expressing gonadotropes, which might be mediated by paracrine interactions within the pituitary. Moreover, migration of GnRH neurons into the hypothalamus appeared normal with appropriate axonal connections to the median eminence, providing genetic evidence against autocrine regulation of GnRH neurons. Surprisingly, genetic ablation of GnRHR expressing cells significantly increased the number of GnRH neurons in the anterior hypothalamus, suggesting an unexpected role of GnRH signaling in establishing the size of the GnRH neuronal population. Our experiments define a functional role of embryonic GnRH signaling.

Y Chromosome microdeletion and altered sperm quality in human males with high concentration of seminal hexachlorocyclohexane (HCH)

Recent studies have shown Y chromosome microdeletions associated with male infertility. The factors responsible for Y chromosome microdeletions in spermatozoa remain unresolved. However, the environmental pollutants are known to damage DNA in differentiating and maturing germ cells in the male reproductive tract. Therefore, the aim of this study was to investigate the effects of seminal hexachlorocyclohexane (HCH) and its isomers, an environmental pollutant, in 50 fertile and 50 infertile males in relation to semen quality and the incidence of Y chromosome microdeletion in azoospermic factor (AZF) region. As compared to control, an increased HCH level and significantly decreased semen quality were observed in the infertile males. A positive significant association was found between sperm count with alpha-HCH and beta-HCH in the infertile males. A negative significant association was observed between sperm counts with gamma-HCH in asthenospermia patients and with beta-HCH and total HCH in oligo-asthenospermic patients. Out of 100 males studied, we found 10 patients with Yq deletion in AZFa and AZFc regions. Subdivision of infertile group revealed a deletion incidence of 61.5% in azoospermic patients, 11.1% in oligospermic patients and 16.6% in oligo-asthenospermic patients. The presence of Yq deletion in azoospermic patients with a significant mean difference of beta-HCH and total HCH in relation to reduced semen quality seem to corroborate with the mutagenic activity of HCH. The results of this study indicated the susceptibility of male germ line to mutagenic potential of HCH which is an acknowledged risk factor leading to spermatogenic failure.

Genetic interaction between Lrp6 and Wnt5a during mouse development

Lrp6 is generally described as a receptor required for signal transduction in the Wnt/beta-catenin pathway. Wnt5a, however, is a Wnt ligand that usually does not activate Wnt/beta-catenin but rather activates noncanonical Wnt signaling. We have previously shown that Lrp6 can inhibit noncanonical Wnt5a/Wnt11 signaling and that Lrp5/6 loss-of-function produces noncanonical gain-of function defects, which can be rescued by loss of Wnt5a. Here, we describe other phenotypes found in Wnt5a/Lrp6 compound mutant mice, including a worsening of individual Wnt5a or Lrp6 loss of function phenotypes. Lrp6 haploinsufficiency in a Wnt5a-/- background caused spina bifida and exacerbated posterior truncation. Wnt5a-/-Lrp6-/- embryos displayed presomitic mesoderm morphogenesis, somitogenesis, and neurogenesis defects, which are much more severe than in either of the single mutants. Interestingly these results reveal a further level of complexity in processes in which Wnt5a and LRP6 cooperate, or oppose each other, during mouse development.

Programmed genetic instability: a tumor-permissive mechanism for maintaining the evolvability of higher species through methylation-dependent mutation of DNA repair genes in the male germ line

Tumor suppressor genes are classified by their somatic behavior either as caretakers (CTs) that maintain DNA integrity or as gatekeepers (GKs) that regulate cell survival, but the germ line role of these disease-related gene subgroups may differ. To test this hypothesis, we have used genomic data mining to compare the features of human CTs (n = 38), GKs (n = 36), DNA repair genes (n = 165), apoptosis genes (n = 622), and their orthologs. This analysis reveals that repair genes are numerically less common than apoptosis genes in the genomes of multicellular organisms (P < 0.01), whereas CT orthologs are commoner than GK orthologs in unicellular organisms (P < 0.05). Gene targeting data show that CTs are less essential than GKs for survival of multicellular organisms (P < 0.0005) and that CT knockouts often permit offspring viability at the cost of male sterility. Patterns of human familial oncogenic mutations confirm that isolated CT loss is commoner than is isolated GK loss (P < 0.00001). In sexually reproducing species, CTs appear subject to less efficient purifying selection (i.e., higher Ka/Ks) than GKs (P = 0.000003); the faster evolution of CTs seems likely to be mediated by gene methylation and reduced transcription-coupled repair, based on differences in dinucleotide patterns (P = 0.001). These data suggest that germ line CT/repair gene function is relatively dispensable for survival, and imply that milder (e.g., epimutational) male prezygotic repair defects could enhance sperm variation—and hence environmental adaptation and speciation—while sparing fertility. We submit that CTs and repair genes are general targets for epigenetically initiated adaptive evolution, and propose a model in which human cancers arise in part as an evolutionarily programmed side effect of age- and damage-inducible genetic instability affecting both somatic and germ line lineages.

Characteristics, causes and evolutionary consequences of male-biased mutation

Mutation has traditionally been considered a random process, but this paradigm is challenged by recent evidence of divergence rate heterogeneity in different genomic regions. One facet of mutation rate variation is the propensity for genetic change to correlate with the number of germ cell divisions, reflecting the replication-dependent origin of many mutations. Haldane was the first to connect this association of replication and mutation to the difference in the number of cell divisions in oogenesis (low) and spermatogenesis (usually high), and the resulting sex difference in the rate of mutation. The concept of male-biased mutation has been thoroughly analysed in recent years using an evolutionary approach, in which sequence divergence of autosomes and/or sex chromosomes are compared to allow inference about the relative contribution of mothers and fathers in the accumulation of mutations. For instance, assuming that a neutral sequence is analysed, that rate heterogeneity owing to other factors is cancelled out by the investigation of many loci and that the effect of ancestral polymorphism is properly taken into account, the male-to-female mutation rate ratio, alpham, can be solved from the observed difference in rate of X and Y chromosome divergence. The male mutation bias is positively correlated with the relative excess of cell divisions in the male compared to the female germ line, as evidenced by a generation time effect: in mammals, alpham is estimated at approximately 4-6 in primates, approximately 3 in carnivores and approximately 2 in small rodents. Another life-history correlate is sexual selection: when there is intense sperm competition among males, increased sperm production will be associated with a larger number of mitotic cell divisions in spermatogenesis and hence an increase in alpham. Male-biased mutation has implications for important aspects of evolutionary biology such as mate choice in relation to mutation load, sexual selection and the maintenance of genetic diversity despite strong directional selection, the tendency for a disproportionate large role of the X (Z) chromosome in post-zygotic isolation, and the evolution of sex.

Substitution Rate Heterogeneity and the Male Mutation Bias

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (alpha(m)) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of alpha(m) and c. Alpha(m) is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of alpha(m) based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G15 = 37.6, p = 0.0010) as well as within W chromosome introns (G15 = 44.0, p = 0.0001). This heterogeneity also affected the estimates of alpha(m), which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F(13,14) = 2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of alpha(m).

Inefficient purifying selection: the mammalian Y chromosome in the rodent genus Mus

Two related genes with potentially similar functions, one on the Y chromosome and one on the X chromosome, were examined to determine if they evolved differently because of their chromosomal positions. Six hundred fifty-seven base pairs of coding sequence of Jarid1d (Smcy) on the Y chromosome and Jarid1c (Smcx) on the X chromosome were sequenced in 13 rodent taxa. An analysis of replacement and silent substitutions, using a counting method designed for samples with small evolutionary distances, showed a significant difference between the two genes. The different patterns of replacement and silent substitutions within Jarid1d and Jarid1c may be a result of evolutionary mechanisms that are particularly strong on the Y chromosome because of its unique properties. These findings are similar to results of previous studies of Y chromosomal genes in these and other mammalian taxa, suggesting that genes on the mammalian Y evolve in a chromosome-specific manner.

Fast accumulation of nonsynonymous mutations on the female-specific W chromosome in birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller's ratchet, result in a reduction in Ne, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (dN) to synonymous substitutions (dS). We tested this in an avian system by estimating the ratio between dN and dS in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher dN/dS ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the dN/dS ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d(N)/d(S) ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome.

Male-driven evolution in closely related species of the mouse genus Mus

Recently, other researchers have found that closely related primate species had a lower male-to-female mutation rate ratio (alpha) than distantly related species. To determine if this is a general phenomenon affecting other mammalian orders, eleven species or subspecies of the rodent genus Mus and two outgroup species were compared. Intron sequences from a gene in the nonrecombining region of the Y chromosome Jarid1d (Smcy) and its X chromosomal gametolog, Jarid1c (Smcx), were analyzed in a phylogenetic context. The male-to-female mutation rate ratio for all thirteen taxa is approximately 2.5, which is similar to previous estimates in more distantly related rodents. However, when branches with lengths of more than 2.5% were removed from the analysis, the male-to-female mutation rate ratio dropped to 0.9. Thus, in closely related rodents, as in closely related primates, the male-to-female mutation rate ratio is lower than expected.

High mutation rates in human and ape pseudoautosomal genes

It has been suggested that recombination may be mutagenic, which, if true, would inflate intraspecies diversity and interspecies silent divergence in regions of high recombination. Here, we test this hypothesis comparing human/orangutan genome-wide non-coding divergence (K) to that in the pseudoautosomal genes which were reported to recombine much more frequently than the rest of the genome. We demonstrate that, compared to the average human/orangutan non-coding divergence (K=3%), the substitution rate is significantly elevated in the introns of SHOX (K=5.7%), PPP2R3L (K=8.7%) and ASMT (K=6.5%) genes located in the human and orangutan Xp/Yp pseudoautosomal region (p-PAR), where recombination is over 20-fold higher than the genomic average. On the other hand, human/orangutan non-coding divergence at the Xp/Yp pseudoautosomal boundary (K=3.5%) and in the SYBL1 gene (K=2.7%), located in the human Xq/Yq pseudoautosomal region (q-PAR), where recombination is known to be less frequent than in p-PAR, was not significantly higher than the genome average. The data are consistent with the hypothesis that recombination may be mutagenic.

Male-biased transmission of deleterious mutations to the progeny in Arabidopsis thaliana

The extent and cause of male-biased mutation rates, the higher number of mutations in sperm than in eggs, is currently an active and controversial subject. Recent evidence indicates that this male (sperm) bias not only occurs in animals but also in plants. The higher mutation rate in plant sperm was inferred from rates of evolution of neutral DNA regions, and the results were confined to the mitochondria and chloroplasts of gymnosperms. However, the relative transmission rates of deleterious mutations, which have substantial evolutionary consequences, have rarely been studied. Here, an investigation is described by using the hermaphroditic self-compatible flowering plant Arabidopsis thaliana, in which we artificially increased the rate of mutation in pollen (i.e., sperm donor) and maternal (i.e., egg donor) parents, by using two kinds of UV irradiation in parallel and separate experiments, and assessed the deleterious effects on fitness of the F(2) generation. The results show that more deleterious induced mutations are transmitted to the progeny by a sperm than by an egg. These findings provide the first experimental evidence that more deleterious mutations are inherited from sperm than from an egg in any organism. Possible causes underlying this male bias are discussed.

Male-driven evolution of mitochondrial and chloroplastidial DNA sequences in plants

Although there is substantial evidence that, in animals, male-inherited neutral DNA evolves at a higher rate than female-inherited DNA, the relative evolutionary rate of male- versus female-inherited DNA has not been investigated in plants. We compared the substitution rates at neutral sites of maternally and paternally inherited organellar DNA in gymnosperms. The analysis provided substantial support for the presence of a higher evolutionary rate in both the mitochondrial and chloroplastidial DNA when the organelle was inherited paternally than when inherited maternally. These results suggest that, compared with eggs, sperm tend to carry a greater number of mutations in mitochondrial and chloroplastidial DNA. The existence of a male mutation bias in plants is remarkable because, unlike animals, the germ-lines are not separated from the somatic cells throughout an individual's lifetime. The data therefore suggest that even a brief period of male and female germ-line separation can cause gender-specific mutation rates. These results are the first to show that, at least in some species, germ-lines influence the number of mutations carried in the gametes. Possible causes of male mutation bias in plants are discussed.

Substitution rates in the X- and Y-linked genes of the plants, Silene latifolia and S. dioica

Theory predicts that selection should be less effective in the nonrecombining genes of Y-chromosomes, relative to the situation for genes on the other chromosomes, and this should lead to the accumulation of deleterious nonsynonymous substitutions. In addition, synonymous substitution rates may differ between X- and Y-linked genes because of the male-driven evolution effect and also because of actual differences in per-replication mutation rates between the sex chromosomes. Here, we report the first study of synonymous and nonsynonymous substitution rates on plant sex chromosomes. We sequenced two pairs of sex-linked genes, SlX1-SlY1 and SlX4-SlY4, from dioecious Silene latifolia and S. dioica, and their non-sex-linked homologues from nondioecious S. vulgaris and Lychnis flos-jovis, respectively. The rate of nonsynonymous substitutions in the SlY4 gene is significantly higher than that in the SlX4 gene. Silent substitution rates are also significantly higher in both Y-linked genes, compared with their X-linked homologues. The higher nonsynonymous substitution rate in the SlY4 gene is therefore likely to be caused by a mutation rate difference between the sex chromosomes. The difference in silent substitution rates between the SlX4 and SlY4 genes is too great to be explained solely by a higher per-generation mutation rate in males than females. It is thus probably caused by a difference in per-replication mutation rates between the sex chromosomes. This suggests that the local mutation rate can change in a relatively short evolutionary time.

Contrasting levels of nucleotide diversity on the avian Z and W sex chromosomes

Sex chromosomes may provide a context for studying the local effects of mutation rate on molecular evolution, since the two types of sex chromosomes are generally exposed to different mutational environments in male and female germ lines. Importantly, recent studies of some vertebrates have provided evidence for a higher mutation rate among males than among females. Thus, in birds, the Z chromosome, which spends two thirds of its time in the male germ line, is exposed to more mutations than the female-specific W chromosome. We show here that levels of nucleotide diversity are drastically higher on the avian Z chromosome than in paralogous sequences on the W chromosome. In fact, no intraspecific polymorphism whatsoever was seen in about 3.4 kb of CHD1W intron sequence from a total of >150 W chromosome copies of seven different bird species. In contrast, the amount of genetic variability in paralogous sequences on the Z chromosome was significant, with an average pairwise nucleotide diversity (d) of 0.0020 between CHD1Z introns and with 37 segregating sites in a total of 3.8 kb of Z sequence. The contrasting levels of genetic variability on the avian sex chromosomes are thus in a direction predicted from a male-biased mutation rate. However, although a low gene number, as well as some other factors, argues against background selection and/or selective sweeps shaping the genetic variability of the avian W chromosome, we cannot completely exclude selection as a contributor to the low levels of variation on the W chromosome.

Molecular evolution of the avian CHD1 genes on the Z and W sex chromosomes

Genes shared between the nonrecombining parts of the two types of sex chromosomes offer a potential means to study the molecular evolution of the same gene exposed to different genomic environments. We have analyzed the molecular evolution of the coding sequence of the first pair of genes found to be shared by the avian Z (present in both sexes) and W (female-specific) sex chromosomes, CHD1Z and CHD1W. We show here that these two genes evolve independently but are highly conserved at nucleotide as well as amino acid levels, thus not indicating a female-specific role of the CHD1W gene. From comparisons of sequence data from three avian lineages, the frequency of nonsynonymous substitutions (K(a)) was found to be higher for CHD1W (1.55 per 100 sites) than for CHD1Z (0.81), while the opposite was found for synonymous substitutions (K(s), 13.5 vs. 22.7). We argue that the lower effective population size and the absence of recombination on the W chromosome will generally imply that nonsynonymous substitutions accumulate faster on this chromosome than on the Z chromosome. The same should be true for the Y chromosome relative to the X chromosome in XY systems. Our data are compatible with a male-biased mutation rate, manifested by the faster rate of neutral evolution (synonymous substitutions) on the Z chromosome than on the female-specific W chromosome.

How important is DNA replication for mutagenesis?

Rates of mutation and substitution in mammals are generally greater in the germ lines of males. This is usually explained as resulting from the larger number of germ cell divisions during spermatogenesis compared with oogenesis, with the assumption made that mutations occur primarily during DNA replication. However, the rate of cell division is not the only difference between male and female germ lines, and mechanisms are known that can give rise to mutations independently of DNA replication. We investigate the possibility that there are other causes of male-biased mutation. First, we show that patterns of variation at approximately 5,200 short tandem repeat (STR) loci indicate a higher mutation rate in males. We estimate a ratio of male-to-female mutation rates of approximately 1.9. This is significantly greater than 1 and supports a greater rate of mutation in males, affecting the evolution of these loci. Second, we show that there are chromosome-specific patterns of nucleotide and dinucleotide composition in mammals that have been shaped by mutation at CpG dinucleotides. Comparable patterns occur in birds. In mammals, male germ lines are more methylated than female germ lines, and these patterns indicate that differential methylation has played a role in male-biased vertebrate evolution. However, estimates of male mutation bias obtained from both classes of mutation are substantially lower than estimates of cell division bias from anatomical data. This discrepancy, along with published data indicating slipped-strand mispairing arising at STR loci in nonreplicating DNA, suggests that a substantial percentage of mutation may occur in nonreplicating DNA.

Sex chromosomal transposable element accumulation and male-driven substitutional evolution in humans

We sequenced the genomic region containing the human Y-linked zinc finger gene (ZFY). Comparison of ZFY to the related region on the X chromosome (ZFX) and to autosomal sequences reveals a significant accumulation of transposable elements on the sex chromosomes. In addition, five times as many retroviruslike elements (RLEs) are present in the ZFY region as in the ZFX region. Thus, transposable elements accumulate more rapidly on the sex chromosomes, and the insertion of RLEs may occur more frequently in the male than in the female germ line. When the accumulation of substitutions in Alu elements was analyzed, it was found that the Alu elements at the Y-chromosomal locus diverged significantly faster than those at the X-chromosomal locus, whereas the divergence of autosomal Alu elements was intermediate. The male-to-female mutation rate ratio was estimated to be 2.5.

The human SWI-SNF complex protein p270 is an ARID family member with non-sequence-specific DNA binding activity

p270 is an integral member of human SWI-SNF complexes, first identified through its shared antigenic specificity with p300 and CREB binding protein. The deduced amino acid sequence of p270 reported here indicates that it is a member of an evolutionarily conserved family of proteins distinguished by the presence of a DNA binding motif termed ARID (AT-rich interactive domain). The ARID consensus and other structural features are common to both p270 and yeast SWI1, suggesting that p270 is a human counterpart of SWI1. The approximately 100-residue ARID sequence is present in a series of proteins strongly implicated in the regulation of cell growth, development, and tissue-specific gene expression. Although about a dozen ARID proteins can be identified from database searches, to date, only Bright (a regulator of B-cell-specific gene expression), dead ringer (a Drosophila melanogaster gene product required for normal development), and MRF-2 (which represses expression from the cytomegalovirus enhancer) have been analyzed directly in regard to their DNA binding properties. Each binds preferentially to AT-rich sites. In contrast, p270 shows no sequence preference in its DNA binding activity, thereby demonstrating that AT-rich binding is not an intrinsic property of ARID domains and that ARID family proteins may be involved in a wider range of DNA interactions.

A novel DNA-binding motif shares structural homology to DNA replication and repair nucleases and polymerases

A novel class of DNA-binding domains has been established from at least sixteen recently identified DNA-binding proteins. The three-dimensional structure of one of these domains, Mrf-2, has been solved using NMR methods. This structure is significantly different from known DNA-binding domain structures. The mechanism of DNA recognition by this motif has been suggested based on conserved residues, surface electrostatic potentials and chemical shift changes. This new DNA-binding motif shares structural homology with T4 RNase H, E. coli endonuclease III and Bacillus subtilis DNA polymerase I. The structural homology suggests a mechanism for substrate recognition by these enzymes.

Transposition of RhoA to the murine Y chromosome

In an effort to produce a more complete transcription map of the short (approximately 5 Mb) arm of the mouse Y chromosome, we have initiated exon trapping from Yp-derived YACs. Sequence analysis of the trapped products has identified exons of previously cloned mouse Y-located genes Zfy and SSty and potential exons homologous to the human Y-located Tspy gene family. In addition, a family of three Yp-located transcripts that show close homology to human RHOA (locus designation ARHA), a member of the Ras family of small GTPases, has been identified. To determine whether these Yp sequences had been transposed from an autosomal ancestor, we used this trapped product to isolate a full-length autosomal mouse RhoA cDNA that is 80% identical at the nucleotide level and 98% identical at the amino acid level to human RHOA and maps to mouse Chromosome 2 (locus designation ArhA). Sequence analysis indicates that the Y-linked copies have diverged from the autosomal form, with small deletions precluding maintenance of a significant open reading frame in all Yp copies. Yet RT-PCR analysis indicates that two of these pseudogenes, RhoAy1 and 3, are expressed in a testis-specific manner, in sharp contrast to the nearly ubiquitous expression pattern of the autosomal ancestor. The data indicate that the Y copies of RhoA have been transposed from an autosome, followed by subsequent duplication, sequence divergence, and acquisition of a testis-specific promoter/enhancer.

Cloning and mapping of bovine ZFX gene to the long arm of the X-chromosome (Xq34) and homologous mapping of ZFY gene to the distal region of the short arm of the bovine (Yp13), ovine (Yp12-p13), and caprine (Yp12-p13) Y chromosome

A part of mouse Zfy-2 sequence was synthesized and used to screen a genomic library of the spinous country-rat (Tokudaia osimensis spp., 2n = 45). An isolated clone had the C-terminal region of Zfy, which consisted of 1190 bp, encoded 336 amino acid residues, and harbored 11 out of 13 zinc finger motifs. With this as a probe, a bovine testis cDNA library was screened. Two ZFX clones were isolated and their sequences combined. The short sequence, lacking part of the 5' upstream region, was amplified by PCR or RT-PCR, cloned, and sequenced. A full-length ZFX was constructed by combining these three sequences. The bovine ZFX consisted of 5328 bp and encoded 800 amino acid residues, which contained 13 zinc finger motifs. ZFX was used as a probe for fluorescence in situ hybridization and was mapped to Xq34, different from its previously reported site at Xq21-q231. A SINE (short interspersed nuclear element) sequence consisting of 188 bp was found close to the end of the 3'-untranslated region of ZFX. The SINE sequence hybridized to all bovine chromosomes. ZFY is highly homologous with ZFX and, as a result, ZFY could be mapped simultaneously. ZFY was mapped to the distal region of the short arm of the Y Chromosome (Chr) (Yp13), contradicting the previously reported position Yq1. Ovine and caprine ZFY were also mapped with bovine ZFX. Both were mapped to the distal region of the short arm of the Y Chr (Yp12-p13). Ovine ZFX was mapped to a region close to the centromere of the X Chr (Xq13).