Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway

Disclosing complex mutational dynamics at a Y chromosome palindrome evolving through intra- and inter-chromosomal gene conversion

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1-P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10-5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10-8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

Human SRY Expression at the Sex-determining Period is Insufficient to Drive Testis Development in Mice

The sex-determining region of the Y chromosome, Sry/SRY, is an initiation factor for testis development in both humans and mice. Although the functional compatibility between murine SRY and human SRY was previously examined in transgenic mice, their equivalency remains inconclusive. Because molecular interaction and timeline of mammalian sex determination were mostly described in murine experiments, we generated a mouse model in which Sry was substituted with human SRY to verify the compatibility. The mouse model had the human SRY open reading frame at the locus of murine Sry exon 1-Sry(SRY) mice-and was generated using the CRISPR/Cas9 system. The reproductive system of the mice was analyzed. The expression of human SRY in the fetal gonadal ridge of Sry(SRY) mice was detected. The external and internal genitalia of adult Sry(SRY) mice were similar to those of wild-type females, without any significant difference in anogenital distance. Sry(SRY) mice obtained gonads, which were morphologically considered as ovaries. Histological analysis revealed that the cortical regions of gonads from adult Sry(SRY) mice contained few follicles. We successfully replaced genes on the Y chromosome with targeted genome editing using the CRISPR/Cas9 system. Because the Sry(SRY) XY mice did not develop testis, we concluded that human SRY was insufficient to drive testis development in mouse embryos. The difference in response elements and lack of glutamine-rich domains may have invalidated human SRY function in mice. Signal transduction between Sry/SRY expression and Sox9/SOX9 activation is possibly organized in a species-specific manner.

Genetic Disease and Therapy

Genetic diseases cause numerous complex and intractable pathologies. DNA sequences encoding each human's complexity and many disease risks are contained in the mitochondrial genome, nuclear genome, and microbial metagenome. Diagnosis of these diseases has unified around applications of next-generation DNA sequencing. However, translating specific genetic diagnoses into targeted genetic therapies remains a central goal. To date, genetic therapies have fallen into three broad categories: bulk replacement of affected genetic compartments with a new exogenous genome, nontargeted addition of exogenous genetic material to compensate for genetic errors, and most recently, direct correction of causative genetic alterations using gene editing. Generalized methods of diagnosis, therapy, and reagent delivery into each genetic compartment will accelerate the next generations of curative genetic therapies. We discuss the structure and variability of the mitochondrial, nuclear, and microbial metagenomic compartments, as well as the historical development and current practice of genetic diagnostics and gene therapies targeting each compartment.

Paternal gene pool of Malays in Southeast Asia and its applications for the early expansion of Austronesians

OBJECTIVES

The origin and differentiation of Austronesian populations and their languages have long fascinated linguists, archeologists, and geneticists. However, the founding process of Austronesians and when they separated from their close relatives, such as the Daic and Austro-Asiatic populations in the mainland of Asia, remain unclear. In this study, we explored the paternal origin of Malays in Southeast Asia and the early differentiation of Austronesians.

MATERIALS AND METHODS

We generated whole Y-chromosome sequences of 50 Malays and co-analyzed 200 sequences from other Austronesians and related populations. We generated a revised phylogenetic tree with time estimation.

RESULTS

We identified six founding paternal lineages among the studied Malays samples. These founding lineages showed a surprisingly coincident expansion age at 5000 to 6000 years ago. We also found numerous mostly close related samples of the founding lineages of Malays among populations from Mainland of Asia.

CONCLUSION

Our analyses provided a refined phylogenetic resolution for the dominant paternal lineages of Austronesians found by previous studies. We suggested that the co-expansion of numerous founding paternal lineages corresponds to the initial differentiation of the most recent common ancestor of modern Austronesians. The splitting time and divergence pattern in perspective of paternal Y-chromosome evidence are highly consistent with the previous theories of ethnologists, linguists, and archeologists.

Extinction of chromosomes due to specialization is a universal occurrence

The human X and Y chromosomes evolved from a pair of autosomes approximately 180 million years ago. Despite their shared evolutionary origin, extensive genetic decay has resulted in the human Y chromosome losing 97% of its ancestral genes while gene content and order remain highly conserved on the X chromosome. Five 'stratification' events, most likely inversions, reduced the Y chromosome's ability to recombine with the X chromosome across the majority of its length and subjected its genes to the erosive forces associated with reduced recombination. The remaining functional genes are ubiquitously expressed, functionally coherent, dosage-sensitive genes, or have evolved male-specific functionality. It is unknown, however, whether functional specialization is a degenerative phenomenon unique to sex chromosomes, or if it conveys a potential selective advantage aside from sexual antagonism. We examined the evolution of mammalian orthologs to determine if the selective forces that led to the degeneration of the Y chromosome are unique in the genome. The results of our study suggest these forces are not exclusive to the Y chromosome, and chromosomal degeneration may have occurred throughout our evolutionary history. The reduction of recombination could additionally result in rapid fixation through isolation of specialized functions resulting in a cost-benefit relationship during times of intense selective pressure.

Birth, expansion, and death of VCY-containing palindromes on the human Y chromosome

BACKGROUND

Large palindromes (inverted repeats) make up substantial proportions of mammalian sex chromosomes, often contain genes, and have high rates of structural variation arising via ectopic recombination. As a result, they underlie many genomic disorders. Maintenance of the palindromic structure by gene conversion between the arms has been documented, but over longer time periods, palindromes are remarkably labile. Mechanisms of origin and loss of palindromes have, however, received little attention.

RESULTS

Here, we use fiber-FISH, 10x Genomics Linked-Read sequencing, and breakpoint PCR sequencing to characterize the structural variation of the P8 palindrome on the human Y chromosome, which contains two copies of the VCY (Variable Charge Y) gene. We find a deletion of almost an entire arm of the palindrome, leading to death of the palindrome, a size increase by recruitment of adjacent sequence, and other complex changes including the formation of an entire new palindrome nearby. Together, these changes are found in ~ 1% of men, and we can assign likely molecular mechanisms to these mutational events. As a result, healthy men can have 1-4 copies of VCY.

CONCLUSIONS

Gross changes, especially duplications, in palindrome structure can be relatively frequent and facilitate the evolution of sex chromosomes in humans, and potentially also in other mammalian species.

Challenges and Approaches in Microbiome Research: From Fundamental to Applied

We face major agricultural challenges that remain a threat for global food security. Soil microbes harbor enormous potentials to provide sustainable and economically favorable solutions that could introduce novel approaches to improve agricultural practices and, hence, crop productivity. In this review we give an overview regarding the current state-of-the-art of microbiome research by discussing new technologies and approaches. We also provide insights into fundamental microbiome research that aim to provide a deeper understanding of the dynamics within microbial communities, as well as their interactions with different plant hosts and the environment. We aim to connect all these approaches with potential applications and reflect how we can use microbial communities in modern agricultural systems to realize a more customized and sustainable use of valuable resources (e.g., soil).

Toward a consensus on SNP and STR mutation rates on the human Y-chromosome

The mutation rate on the Y-chromosome matters for estimating the time-to-the-most-recent-common-ancestor (TMRCA, i.e. haplogroup age) in population genetics, as well as for forensic, medical, and genealogical studies. Large-scale sequencing efforts have produced several independent estimates of Y-SNP mutation rates. Genealogical, or pedigree, rates tend to be slightly faster than evolutionary rates obtained from ancient DNA or calibrations using dated (pre)historical events. It is, therefore, suggested to report TMRCAs using an envelope defined by the average aDNA-based rate and the average pedigree-based rate. The current estimate of the "envelope rate" is 0.75-0.89 substitutions per billion base pairs per year. The available Y-SNP mutation rates can be applied to high-coverage data from the entire X-degenerate region, but other datasets may demand recalibrated rates. While a consensus on Y-SNP rates is approaching, the debate on Y-STR rates has continued for two decades, because multiple genealogical rates were consistent with each other but three times faster than the single evolutionary estimate. Applying Y-SNP and Y-STR rates to the same haplogroups recently helped to clarify the issue. Genealogical and evolutionary STR rates typically provide lower and upper bounds of the "true" (SNP-based) age. The genealogical rate often-but not always-works well for haplogroups less than 7000 years old. The evolutionary rate, although calibrated using recent events, inflates ages of young haplogroups and deflates the age of the entire Y-chromosomal tree, but often provides reasonable estimates for intermediate ages (old haplogroups). Future rate estimates and accumulating case studies should further clarify the Y-SNP rates.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Systems-level chromosomal parameters represent a suprachromosomal basis for the non-random chromosomal arrangement in human interphase nuclei

Forty-six chromosome territories (CTs) are positioned uniquely in human interphase nuclei, wherein each of their positions can range from the centre of the nucleus to its periphery. A non-empirical basis for their non-random arrangement remains unreported. Here, we derive a suprachromosomal basis of that overall arrangement (which we refer to as a CT constellation), and report a hierarchical nature of the same. Using matrix algebra, we unify intrinsic chromosomal parameters (e.g., chromosomal length, gene density, the number of genes per chromosome), to derive an extrinsic effective gene density matrix, the hierarchy of which is dominated largely by extrinsic mathematical coupling of HSA19, followed by HSA17 (human chromosome 19 and 17, both preferentially interior CTs) with all CTs. We corroborate predicted constellations and effective gene density hierarchy with published reports from fluorescent in situ hybridization based microscopy and Hi-C techniques, and delineate analogous hierarchy in disparate vertebrates. Our theory accurately predicts CTs localised to the nuclear interior, which interestingly share conserved synteny with HSA19 and/or HSA17. Finally, the effective gene density hierarchy dictates how permutations among CT position represents the plasticity within its constellations, based on which we suggest that a differential mix of coding with noncoding genome modulates the same.

Y-Chromosome Structural Diversity in the Bonobo and Chimpanzee Lineages

The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ. We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

Large-scale recent expansion of European patrilineages shown by population resequencing

The proportion of Europeans descending from Neolithic farmers ∼ 10 thousand years ago (KYA) or Palaeolithic hunter-gatherers has been much debated. The male-specific region of the Y chromosome (MSY) has been widely applied to this question, but unbiased estimates of diversity and time depth have been lacking. Here we show that European patrilineages underwent a recent continent-wide expansion. Resequencing of 3.7 Mb of MSY DNA in 334 males, comprising 17 European and Middle Eastern populations, defines a phylogeny containing 5,996 single-nucleotide polymorphisms. Dating indicates that three major lineages (I1, R1a and R1b), accounting for 64% of our sample, have very recent coalescent times, ranging between 3.5 and 7.3 KYA. A continuous swathe of 13/17 populations share similar histories featuring a demographic expansion starting ∼ 2.1-4.2 KYA. Our results are compatible with ancient MSY DNA data, and contrast with data on mitochondrial DNA, indicating a widespread male-specific phenomenon that focuses interest on the social structure of Bronze Age Europe.

Genome sequence of a 45,000-year-old modern human from western Siberia

We present the high-quality genome sequence of a ∼45,000-year-old modern human male from Siberia. This individual derives from a population that lived before-or simultaneously with-the separation of the populations in western and eastern Eurasia and carries a similar amount of Neanderthal ancestry as present-day Eurasians. However, the genomic segments of Neanderthal ancestry are substantially longer than those observed in present-day individuals, indicating that Neanderthal gene flow into the ancestors of this individual occurred 7,000-13,000 years before he lived. We estimate an autosomal mutation rate of 0.4 × 10(-9) to 0.6 × 10(-9) per site per year, a Y chromosomal mutation rate of 0.7 × 10(-9) to 0.9 × 10(-9) per site per year based on the additional substitutions that have occurred in present-day non-Africans compared to this genome, and a mitochondrial mutation rate of 1.8 × 10(-8) to 3.2 × 10(-8) per site per year based on the age of the bone.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Endogenous retrovirus-mediated genomic variations in chimpanzees

Transposable elements (TEs) have played a significant role in the evolution of host genome by triggering genomic rearrangements. TEs have been studied in various research fields, ranging from population genomics to personalized medicines. Human-specific TEs and TEs existing in the human genome have been well studied. Unlike them, non-human primate-specific TEs remain shrouded in mystery. However, the study of TE-mediated genomic or genetic variations through comparative genomics is essential to understand mechanisms which TEs utilize to modify species-specific genome architecture and to cause species-specific diseases, Therefore, we have studied chimpanzee-specific TEs as well as human-specific TEs. At first, we identified human-specific HERV-K integrated into the human genome after the divergence of human and chimpanzee. Then, for a comparative study of HERV-Ks and non-human ERVs, we extracted chimpanzee-specific endogenous retroviruses (PtERVs) from the chimpanzee genome. We identified 256 chimpanzee-specific PtERVs and characterized them, focusing on their estimated evolutionary age, polymorphism level in chimpanzee populations, and potential impact on the difference between the human and chimpanzee genomes.

Longevity and plasticity of CFTR provide an argument for noncanonical SNP organization in hominid DNA

Like many other ancient genes, the cystic fibrosis transmembrane conductance regulator (CFTR) has survived for hundreds of millions of years. In this report, we consider whether such prodigious longevity of an individual gene – as opposed to an entire genome or species – should be considered surprising in the face of eons of relentless DNA replication errors, mutagenesis, and other causes of sequence polymorphism. The conventions that modern human SNP patterns result either from purifying selection or random (neutral) drift were not well supported, since extant models account rather poorly for the known plasticity and function (or the established SNP distributions) found in a multitude of genes such as CFTR. Instead, our analysis can be taken as a polemic indicating that SNPs in CFTR and many other mammalian genes may have been generated—and continue to accrue—in a fundamentally more organized manner than would otherwise have been expected. The resulting viewpoint contradicts earlier claims of ‘directional’ or ‘intelligent design-type’ SNP formation, and has important implications regarding the pace of DNA adaptation, the genesis of conserved non-coding DNA, and the extent to which eukaryotic SNP formation should be viewed as adaptive.

Evaluating the Y chromosomal timescale in human demographic and lineage dating

Y chromosome is a superb tool for inferring human evolution and recent demographic history from a paternal perspective. However, Y chromosomal substitution rates obtained using different modes of calibration vary considerably, and have produced disparate reconstructions of human history. Here, we discuss how substitution rate and date estimates are affected by the choice of different calibration points. We argue that most Y chromosomal substitution rates calculated to date have shortcomings, including a reliance on the ambiguous human-chimpanzee divergence time, insufficient sampling of deep-rooting pedigrees, and using inappropriate founding migrations, although the rates obtained from a single pedigree or calibrated with the peopling of the Americas seem plausible. We highlight the need for using more deep-rooting pedigrees and ancient genomes with reliable dates to improve the rate estimation.

Y chromosomes of 40% Chinese descend from three Neolithic super-grandfathers

Demographic change of human populations is one of the central questions for delving into the past of human beings. To identify major population expansions related to male lineages, we sequenced 78 East Asian Y chromosomes at 3.9 Mbp of the non-recombining region, discovered >4,000 new SNPs, and identified many new clades. The relative divergence dates can be estimated much more precisely using a molecular clock. We found that all the Paleolithic divergences were binary; however, three strong star-like Neolithic expansions at ∼6 kya (thousand years ago) (assuming a constant substitution rate of 1×10(-9)/bp/year) indicates that ∼40% of modern Chinese are patrilineal descendants of only three super-grandfathers at that time. This observation suggests that the main patrilineal expansion in China occurred in the Neolithic Era and might be related to the development of agriculture.

Chimpanzee-specific endogenous retrovirus generates genomic variations in the chimpanzee genome

Endogenous retroviruses (ERVs), eukaryotic transposable elements, exist as proviruses in vertebrates including primates and contribute to genomic changes during the evolution of their host genomes. Many studies about ERVs have focused on the elements residing in the human genome but only a few studies have focused on the elements which exist in non-human primate genomes. In this study, we identified 256 chimpanzee-specific endogenous retrovirus copies (PtERVs: Pan troglodyte endogenous retroviruses) from the chimpanzee reference genome sequence through comparative genomics. Among the chimpanzee-specific ERV copies, 121 were full-length chimpanzee-specific ERV elements while 110 were chimpanzee-specific solitary LTR copies. In addition, we found eight potential retrotransposition-competent full-length chimpanzee-specific ERV copies containing an intact env gene, and two of them were polymorphic in chimpanzee individuals. Through computational analysis and manual inspection, we found that some of the chimpanzee-specific ERVs have propagated via non-classical PtERV insertion (NCPI), and at least one of the PtERVs may have played a role in creating an alternative transcript of a chimpanzee gene. Based on our findings in this study, we state that the chimpanzee-specific ERV element is one of the sources of chimpanzee genomic variations, some of which might be related to the alternative transcripts in the chimpanzee population.

The 'extremely ancient' chromosome that isn't: a forensic bioinformatic investigation of Albert Perry's X-degenerate portion of the Y chromosome

Mendez and colleagues reported the identification of a Y chromosome haplotype (the A00 lineage) that lies at the basal position of the Y chromosome phylogenetic tree. Incorporating this haplotype, the authors estimated the time to the most recent common ancestor (TMRCA) for the Y tree to be 338,000 years ago (95% CI=237,000-581,000). Such an extraordinarily early estimate contradicts all previous estimates in the literature and is over a 100,000 years older than the earliest fossils of anatomically modern humans. This estimate raises two astonishing possibilities, either the novel Y chromosome was inherited after ancestral humans interbred with another species, or anatomically modern Homo sapiens emerged earlier than previously estimated and quickly became subdivided into genetically differentiated subpopulations. We demonstrate that the TMRCA estimate was reached through inadequate statistical and analytical methods, each of which contributed to its inflation. We show that the authors ignored previously inferred Y-specific rates of substitution, incorrectly derived the Y-specific substitution rate from autosomal mutation rates, and compared unequal lengths of the novel Y chromosome with the previously recognized basal lineage. Our analysis indicates that the A00 lineage was derived from all the other lineages 208,300 (95% CI=163,900-260,200) years ago.

Comparative analysis of mammalian Y chromosomes illuminates ancestral structure and lineage-specific evolution

Although more than thirty mammalian genomes have been sequenced to draft quality, very few of these include the Y chromosome. This has limited our understanding of the evolutionary dynamics of gene persistence and loss, our ability to identify conserved regulatory elements, as well our knowledge of the extent to which different types of selection act to maintain genes within this unique genomic environment. Here, we present the first MSY (male-specific region of the Y chromosome) sequences from two carnivores, the domestic dog and cat. By combining these with other available MSY data, our multiordinal comparison allows for the first accounting of levels of selection constraining the evolution of eutherian Y chromosomes. Despite gene gain and loss across the phylogeny, we show the eutherian ancestor retained a core set of 17 MSY genes, most being constrained by negative selection for nearly 100 million years. The X-degenerate and ampliconic gene classes are partitioned into distinct chromosomal domains in most mammals, but were radically restructured on the human lineage. We identified multiple conserved noncoding elements that potentially regulate eutherian MSY genes. The acquisition of novel ampliconic gene families was accompanied by signatures of positive selection and has differentially impacted the degeneration and expansion of MSY gene repertoires in different species.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.

Evidence that BRCA1- or BRCA2-associated cancers are not inevitable

Inheriting a BRCA1 or BRCA2 gene mutation can cause a deficiency in repairing complex DNA damage. This step leads to genomic instability and probably contributes to an inherited predisposition to breast and ovarian cancer. Complex DNA damage has been viewed as an integral part of DNA replication before cell division. It causes temporary replication blocks, replication fork collapse, chromosome breaks and sister chromatid exchanges (SCEs). Chemical modification of DNA may also occur spontaneously as a byproduct of normal processes. Pathways containing BRCA1 and BRCA2 gene products are essential to repair spontaneous complex DNA damage or to carry out SCEs if repair is not possible. This scenario creates a theoretical limit that effectively means there are spontaneous BRCA1/2-associated cancers that cannot be prevented or delayed. However, much evidence for high rates of spontaneous DNA mutation is based on measuring SCEs by using bromodeoxyuridine (BrdU). Here we find that the routine use of BrdU has probably led to overestimating spontaneous DNA damage and SCEs because BrdU is itself a mutagen. Evidence based on spontaneous chromosome abnormalities and epidemiologic data indicates strong effects from exogenous mutagens and does not support the inevitability of cancer in all BRCA1/2 mutation carriers. We therefore remove a theoretical argument that has limited efforts to develop chemoprevention strategies to delay or prevent cancers in BRCA1/2 mutation carriers.

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae)

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships. Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota. The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.

Factors affecting the relative abundance of nuclear copies of mitochondrial DNA (numts) in hominoids

Although nuclear copies of mitochondrial DNA (numts) can originate from any portion of the mitochondrial genome, evidence from humans suggests that more variable parts of the mitochondrial genome, such as the mitochondrial control region (MCR), are under-represented in the nucleus. This apparent deficit might arise from the erosion of sequence identity in numts originating from rapidly evolving mitochondrial sequences. However, the extent to which mitochondrial sequence properties impacts the number of numts detected in genomic surveys has not been evaluated. In order to address this question, we: (1) conducted exhaustive BLAST searches of MCR numts in three hominoid genomes; (2) assessed numt prevalence across the four MCR sub-domains (HV1, CCD, HV2, and MCR(F)); (3) estimated their insertion rates in great apes (Hominoidea); and (4) examined the relationship between mitochondrial DNA variability and numt prevalence in sequences originating from MCR and coding regions of the mitochondrial genome. Results indicate a marked deficit of numts from HV2 and MCR(F) MCR sub-domains in all three species. These MCR sub-domains exhibited the highest proportion of variable sites and the lowest number of detected numts per mitochondrial site. Variation in MCR insertion rate between lineages was also observed with a pronounced burst in recent integrations within chimpanzees and orangutans. A deficit of numts from HV2/MCR(F) was observed regardless of age, whereas HV1 is under-represented only in older numts (>25 million years). Finally, more variable mitochondrial genes also exhibit a lower identity with nuclear copies and because of this, appear to be under-represented in human numt databases.

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.

Construction of physical maps for the sex-specific regions of papaya sex chromosomes

Background

Papaya is a major fruit crop in tropical and subtropical regions worldwide. It is trioecious with three sex forms: male, female, and hermaphrodite. Sex determination is controlled by a pair of nascent sex chromosomes with two slightly different Y chromosomes, Y for male and Y^h for hermaphrodite. The sex chromosome genotypes are XY (male), XY^h (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Y^h chromosome region (HSY) is pericentromeric and heterochromatic. Physical mapping of HSY and its X counterpart is essential for sequencing these regions and uncovering the early events of sex chromosome evolution and to identify the sex determination genes for crop improvement.

Results

A reiterate chromosome walking strategy was applied to construct the two physical maps with three bacterial artificial chromosome (BAC) libraries. The HSY physical map consists of 68 overlapped BACs on the minimum tiling path, and covers all four HSY-specific Knobs. One gap remained in the region of Knob 1, the only knob structure shared between HSY and X, due to the lack of HSY-specific sequences. This gap was filled on the physical map of the HSY corresponding region in the X chromosome. The X physical map consists of 44 BACs on the minimum tiling path with one gap remaining in the middle, due to the nature of highly repetitive sequences. This gap was filled on the HSY physical map. The borders of the non-recombining HSY were defined genetically by fine mapping using 1460 F₂ individuals. The genetically defined HSY spanned approximately 8.5 Mb, whereas its X counterpart extended about 5.4 Mb including a 900 Kb region containing the Knob 1 shared by the HSY and X. The 8.5 Mb HSY corresponds to 4.5 Mb of its X counterpart, showing 4 Mb (89%) DNA sequence expansion.

Conclusion

The 89% increase of DNA sequence in HSY indicates rapid expansion of the Y^h chromosome after genetic recombination was suppressed 2–3 million years ago. The genetically defined borders coincide with the common BACs on the minimum tiling paths of HSY and X. The minimum tiling paths of HSY and its X counterpart are being used for sequencing these X and Y^h-specific regions.

Distribution of MGEs and their insertion sites in the Macaca mulatta genome

Mobile genetic elements (MGEs) are fragments of DNA that can move around within the genome through retrotransposition. These are responsible for various important events such as gene inactivation, transduction, regulation of gene expression and genome expansion. The present work involves the identification and study of the distribution of Alu and L1 retrotransposons in the genome of Macaca mulatta, an extensively used organism in biomedical studies. We also make comparisons with MGE distributions in other primate genomes and study the physicochemical properties of the local DNA structure around the transposon insertion site using ELAN. The present work also includes computational testing of the pre-insertion loci in order to detect unique features based on DNA structure, thermodynamic considerations and protein interaction measures. Although there is significant sequence divergence between the elements of M. mulatta and H. sapiens, their genome wide distribution is very similar; comparing the distribution of L1's in all available X chromosome sequences suggests a common mechanism behind the spread of MGE's in primate genomes.

Y-Chromosome variation in hominids: intraspecific variation is limited to the polygamous chimpanzee

Background

We have previously demonstrated that the Y-specific ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y) varied with respect to copy number and position among chimpanzees (Pan troglodytes). In comparison, seven Y-chromosomal lineages of the bonobo (Pan paniscus), the chimpanzee's closest living relative, showed no variation. We extend our earlier comparative investigation to include an analysis of the intraspecific variation of these genes in gorillas (Gorilla gorilla) and orangutans (Pongo pygmaeus), and examine the resulting patterns in the light of the species' markedly different social and mating behaviors.

Methodology/Principal Findings

Fluorescence in situ hybridization analysis (FISH) of DAZ and CDY in 12 Y-chromosomal lineages of western lowland gorilla (G. gorilla gorilla) and a single lineage of the eastern lowland gorilla (G. beringei graueri) showed no variation among lineages. Similar findings were noted for the 10 Y-chromosomal lineages examined in the Bornean orangutan (Pongo pygmaeus), and 11 Y-chromosomal lineages of the Sumatran orangutan (P. abelii). We validated the contrasting DAZ and CDY patterns using quantitative real-time polymerase chain reaction (qPCR) in chimpanzee and bonobo.

Conclusion/Significance

High intraspecific variation in copy number and position of the DAZ and CDY genes is seen only in the chimpanzee. We hypothesize that this is best explained by sperm competition that results in the variant DAZ and CDY haplotypes detected in this species. In contrast, bonobos, gorillas and orangutans—species that are not subject to sperm competition—showed no intraspecific variation in DAZ and CDY suggesting that monoandry in gorillas, and preferential female mate choice in bonobos and orangutans, probably permitted the fixation of a single Y variant in each taxon. These data support the notion that the evolutionary history of a primate Y chromosome is not simply encrypted in its DNA sequences, but is also shaped by the social and behavioral circumstances under which the specific species has evolved.

Evolutionary history of novel genes on the tammar wallaby Y chromosome: Implications for sex chromosome evolution

We report here the isolation and sequencing of 10 Y-specific tammar wallaby (Macropus eugenii) BAC clones, revealing five hitherto undescribed tammar wallaby Y genes (in addition to the five genes already described) and several pseudogenes. Some genes on the wallaby Y display testis-specific expression, but most have low widespread expression. All have partners on the tammar X, along with homologs on the human X. Nonsynonymous and synonymous substitution ratios for nine of the tammar XY gene pairs indicate that they are each under purifying selection. All 10 were also identified as being on the Y in Tasmanian devil (Sarcophilus harrisii; a distantly related Australian marsupial); however, seven have been lost from the human Y. Maximum likelihood phylogenetic analyses of the wallaby YX genes, with respective homologs from other vertebrate representatives, revealed that three marsupial Y genes (HCFC1X/Y, MECP2X/Y, and HUWE1X/Y) were members of the ancestral therian pseudoautosomal region (PAR) at the time of the marsupial/eutherian split; three XY pairs (SOX3/SRY, RBMX/Y, and ATRX/Y) were isolated from each other before the marsupial/eutherian split, and the remaining three (RPL10X/Y, PHF6X/Y, and UBA1/UBE1Y) have a more complex evolutionary history. Thus, the small marsupial Y chromosome is surprisingly rich in ancient genes that are retained in at least Australian marsupials and evolved from testis-brain expressed genes on the X.

Chimpanzees used for medical research shed light on the pathoetiology of leprosy

Leprosy is a chronic infectious disorder caused by Mycobacterium leprae, which mainly affects skin and peripheral nerves. It is classified as either paucibacillary or multibacillary based upon clinical manifestations and slit-skin smear results. It is speculated that leprosy develops after a long latency period following M. leprae infection. However, the actual time of infection and the duration of latency have never been proven in human patients. To date, four cases of spontaneous leprosy have been reported in chimpanzees who were caught in West Africa in infancy and used for medical research in the USA and Japan. One of these chimpanzees was extensively studied in Japan, and single-nucleotide polymorphism analysis for the M. leprae genome was conducted. This analysis revealed that the chimpanzee was infected with M. leprae during infancy in West Africa and the pathognomonic signs of leprosy appeared after at least 30 years of incubation. Analysis of leprosy in chimpanzees can contribute not only to medical research but also to the understanding of the pathoetiology of leprosy.

Locating protein-coding sequences under selection for additional, overlapping functions in 29 mammalian genomes

The degeneracy of the genetic code allows protein-coding DNA and RNA sequences to simultaneously encode additional, overlapping functional elements. A sequence in which both protein-coding and additional overlapping functions have evolved under purifying selection should show increased evolutionary conservation compared to typical protein-coding genes--especially at synonymous sites. In this study, we use genome alignments of 29 placental mammals to systematically locate short regions within human ORFs that show conspicuously low estimated rates of synonymous substitution across these species. The 29-species alignment provides statistical power to locate more than 10,000 such regions with resolution down to nine-codon windows, which are found within more than a quarter of all human protein-coding genes and contain ∼2% of their synonymous sites. We collect numerous lines of evidence that the observed synonymous constraint in these regions reflects selection on overlapping functional elements including splicing regulatory elements, dual-coding genes, RNA secondary structures, microRNA target sites, and developmental enhancers. Our results show that overlapping functional elements are common in mammalian genes, despite the vast genomic landscape.

A gene catalogue of the euchromatic male-specific region of the horse Y chromosome: comparison with human and other mammals

Studies of the Y chromosome in primates, rodents and carnivores provide compelling evidence that the male specific region of Y (MSY) contains functional genes, many of which have specialized roles in spermatogenesis and male-fertility. Little similarity, however, has been found between the gene content and sequence of MSY in different species. This hinders the discovery of species-specific male fertility genes and limits our understanding about MSY evolution in mammals. Here, a detailed MSY gene catalogue was developed for the horse – an odd-toed ungulate. Using direct cDNA selection from horse testis, and sequence analysis of Y-specific BAC clones, 37 horse MSY genes/transcripts were identified. The genes were mapped to the MSY BAC contig map, characterized for copy number, analyzed for transcriptional profiles by RT-PCR, examined for the presence of ORFs, and compared to other mammalian orthologs. We demonstrate that the horse MSY harbors 20 X-degenerate genes with known orthologs in other eutherian species. The remaining 17 genes are acquired or novel and have so far been identified only in the horse or donkey Y chromosomes. Notably, 3 transcripts were found in the heterochromatic part of the Y. We show that despite substantial differences between the sequence, gene content and organization of horse and other mammalian Y chromosomes, the functions of MSY genes are predominantly related to testis and spermatogenesis. Altogether, 10 multicopy genes with testis-specific expression were identified in the horse MSY, and considered likely candidate genes for stallion fertility. The findings establish an important foundation for the study of Y-linked genetic factors governing fertility in stallions, and improve our knowledge about the evolutionary processes that have shaped Y chromosomes in different mammalian lineages.

Identification of tammar wallaby SIRH12, derived from a marsupial-specific retrotransposition event

In humans and mice, there are 11 genes derived from sushi-ichi related retrotransposons, some of which are known to play essential roles in placental development. Interestingly, this family of retrotransposons was thought to exist only in eutherian mammals, indicating their significant contributions to the eutherian evolution, but at least one, PEG10, is conserved between marsupials and eutherians. Here we report a novel sushi-ichi retrotransposon-derived gene, SIRH12, in the tammar wallaby, an Australian marsupial species of the kangaroo family. SIRH12 encodes a protein highly homologous to the sushi-ichi retrotransposon Gag protein in the tammar wallaby, while SIRH12 in the South American short-tailed grey opossum is a pseudogene degenerated by accumulation of multiple nonsense mutations. This suggests that SIRH12 retrotransposition occurred only in the marsupial lineage but acquired and retained some as yet unidentified novel function, at least in the lineage of the tammar wallaby.

Expression of the human TSPY gene in the brains of transgenic mice suggests a potential role of this Y chromosome gene in neural functions

The testis specific protein Y-encoded (TSPY) is a member of TSPY/SET/NAP1 superfamily, encoded within the gonadoblastoma locus on the Y chromosome. TSPY shares a highly conserved SET/NAP-domain responsible for protein--protein interaction among TSPY/SET/NAP1 proteins. Accumulating data, so far, support the role of TSPY as the gonadoblastoma gene, involved in germ cell tumorigenesis. The X-chromosome homolog of TSPY, TSPX is expressed in various tissues at both fetal and adult stages, including the brain, and is capable of interacting with the multi-domain adapter protein CASK, thereby influencing the synaptic and transcriptional functions and developmental regulation of CASK in the brain and other neural tissues. Similar to TSPX, we demonstrated that TSPY could interact with CASK at its SET/NAP-domain in cultured cells. Transgenic mice harboring a human TSPY gene and flanking sequences showed specific expression of the human TSPY transgene in both testis and brain. The neural expression pattern of the human TSPY gene overlapped with those of the endogenous mouse Cask and Tspx gene. Similarly with TSPX, TSPY was co-localized with CASK in neuronal axon fibers in the brain, suggesting a potential role(s) of TSPY in development and/or physiology of the nervous system.

Major chimpanzee-specific structural changes in sperm development-associated genes

A comprehensive analysis of transcriptional structures of chimpanzee sperm development-associated genes is of significant interest for deeply understanding sperm development and male reproductive process. In this study, we sequenced 7,680 clones from a chimpanzee testis full-length cDNA library and obtained 1,933 nonredundant high-quality full-length cDNA sequences. Comparative analysis between human and chimpanzee showed that 78 sperm development-associated genes, most of which were yet uncharacterized, had undergone severe structural changes (mutations at the start/stop codons, INDELs, alternative splicing variations and fusion forms) on genomic and transcript levels throughout chimpanzee evolution. Specifically, among the 78 sperm development-associated genes, 39 including ODF2, UBC, and CD59 showed markedly chimpanzee-specific structural changes. Through dN/dS analysis, we found that 56 transcripts (including seven sperm development-associated genes) had values of greater than one when comparing human and chimpanzee DNA sequences, whereas the values were less than one when comparing humans and orangutans. Gene ontology annotation and expression profiling showed that the chimpanzee testis transcriptome was enriched with genes that are associated with chimpanzee male germ cell development. Taken together, our study provides the first comprehensive molecular evidence that many chimpanzee sperm development-associated genes had experienced severe structural changes over the course of evolution on genomic and transcript levels.

Evaluating the relationship between spermatogenic silencing of the X chromosome and evolution of the Y chromosome in chimpanzee and human

Chimpanzees and humans are genetically very similar, with the striking exception of their Y chromosomes, which have diverged tremendously. The male-specific region (MSY), representing the greater part of the Y chromosome, is inherited from father to son in a clonal fashion, with natural selection acting on the MSY as a unit. Positive selection might involve the performance of the MSY in spermatogenesis. Chimpanzees have a highly polygamous mating behavior, so that sperm competition is thought to provide a strong selective force acting on the Y chromosome in the chimpanzee lineage. In consequence of evolution of the heterologous sex chromosomes in mammals, meiotic sex chromosome inactivation (MSCI) results in a transcriptionally silenced XY body in male meiotic prophase, and subsequently also in postmeiotic repression of the sex chromosomes in haploid spermatids. This has evolved to a situation where MSCI has become a prerequisite for spermatogenesis. Here, by analysis of microarray testicular expression data representing a small number of male chimpanzees and men, we obtained information indicating that meiotic and postmeiotic X chromosome silencing might be more effective in chimpanzee than in human spermatogenesis. From this, we suggest that the remarkable reorganization of the chimpanzee Y chromosome, compared to the human Y chromosome, might have an impact on its meiotic interactions with the X chromosome and thereby on X chromosome silencing in spermatogenesis. Further studies will be required to address comparative functional aspects of MSCI in chimpanzee, human, and other placental mammals.

Y chromosomal variation tracks the evolution of mating systems in chimpanzee and bonobo

The male-specific regions of the Y chromosome (MSY) of the human and the chimpanzee (Pan troglodytes) are fully sequenced. The most striking difference is the dramatic rearrangement of large parts of their respective MSYs. These non-recombining regions include ampliconic gene families that are known to be important for male reproduction,and are consequently under significant selective pressure. However, whether the published Y-chromosomal pattern of ampliconic fertility genes is invariable within P. troglodytes is an open but fundamental question pertinent to discussions of the evolutionary fate of the Y chromosome in different primate mating systems. To solve this question we applied fluorescence in situ hybridisation (FISH) of testis-specific expressed ampliconic fertility genes to metaphase Y chromosomes of 17 chimpanzees derived from 11 wild-born males and 16 bonobos representing seven wild-born males. We show that of eleven P. troglodytes Y-chromosomal lines, ten Y-chromosomal variants were detected based on the number and arrangement of the ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y)-a so-far never-described variation of a species' Y chromosome. In marked contrast, no variation was evident among seven Y-chromosomal lines of the bonobo, P. paniscus, the chimpanzee's closest living relative. Although, loss of variation of the Y chromosome in the bonobo by a founder effect or genetic drift cannot be excluded, these contrasting patterns might be explained in the context of the species' markedly different social and mating behaviour. In chimpanzees, multiple males copulate with a receptive female during a short period of visible anogenital swelling, and this may place significant selection on fertility genes. In bonobos, however, female mate choice may make sperm competition redundant (leading to monomorphism of fertility genes), since ovulation in this species is concealed by the prolonged anogenital swelling, and because female bonobos can occupy high-ranking positions in the group and are thus able to determine mate choice more freely.

The human RPS4 paralogue on Yq11.223 encodes a structurally conserved ribosomal protein and is preferentially expressed during spermatogenesis

Background

The Y chromosome of mammals is particularly prone to accumulate genes related to male fertility. However, the high rate of molecular evolution on this chromosome predicts reduced power to the across-species comparative approach in identifying male-specific genes that are essential for sperm production in humans. We performed a comprehensive analysis of expression of Y-linked transcripts and their X homologues in several human tissues, and in biopsies of infertile patients, in an attempt to identify new testis-specific genes involved in human spermatogenesis.

Results

We present evidence that one of the primate-specific Y-linked ribosomal protein genes, RPS4Y2, has restricted expression in testis and prostate, in contrast with its X-linked homologue, which is ubiquitously expressed. Moreover, we have determined by highly specific quantitative real time PCR that RPS4Y2 is more highly expressed in testis biopsies containing germ cells. The in silico analysis of the promoter region of RPS4Y2 revealed several differences relative to RPS4Y1, the more widely expressed paralogue from which Y2 has originated through duplication. Finally, through comparative modelling we obtained the three dimensional models of the human S4 proteins, revealing a conserved structure. Interestingly, RPS4Y2 shows different inter-domain contacts and the potential to establish specific interactions.

Conclusions

These results suggest that one of the Y-linked copies of the ribosomal protein S4 is preferentially expressed during spermatogenesis and might be important for germ cell development. Even though RPS4Y2 has accumulated several amino acid changes following its duplication from RPS4Y1, approximately 35 million years ago, the evolution of the Y-encoded RPS4 proteins is structurally constrained. However, the exclusive expression pattern of RPS4Y2 and the novelties acquired at the C-terminus of the protein may indicate some degree of functional specialisation of this protein in spermatogenesis.

The (r)evolution of SINE versus LINE distributions in primate genomes: sex chromosomes are important

The densities of transposable elements (TEs) in the human genome display substantial variation both within individual chromosomes and among chromosome types (autosomes and the two sex chromosomes). Finding an explanation for this variability has been challenging, especially in light of genome landscapes unique to the sex chromosomes. Here, using a multiple regression framework, we investigate primate Alu and L1 densities shaped by regional genome features and location on a particular chromosome type. As a result of our analysis, first, we build statistical models explaining up to 79% and 44% of variation in Alu and L1 element density, respectively. Second, we analyze sex chromosome versus autosome TE densities corrected for regional genomic effects. We discover that sex-chromosome bias in Alu and L1 distributions not only persists after accounting for these effects, but even presents differences in patterns, confirming preferential Alu integration in the male germline, yet likely integration of L1s in both male and female germlines or in early embryogenesis. Additionally, our models reveal that local base composition (measured by GC content and density of L1 target sites) and natural selection (inferred via density of most conserved elements) are significant to predicting densities of L1s. Interestingly, measurements of local double-stranded breaks (a 13-mer associated with genome instability) strongly correlate with densities of Alu elements; little evidence was found for the role of recombination-driven deletion in driving TE distributions over evolutionary time. Thus, Alu and L1 densities have been influenced by the combination of distinct local genome landscapes and the unique evolutionary dynamics of sex chromosomes.

Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content

The human Y chromosome began to evolve from an autosome hundreds of millions of years ago, acquiring a sex-determining function and undergoing a series of inversions that suppressed crossing over with the X chromosome^1,2. Little is known about the Y chromosome’s recent evolution because only the human Y chromosome has been fully sequenced. Prevailing theories hold that Y chromosomes evolve by gene loss, the pace of which slows over time, eventually leading to a paucity of genes, and stasis^3,4. These theories have been buttressed by partial sequence data from newly emergent plant and animal Y chromosomes^5-8, but they have not been tested in older, highly evolved Y chromosomes like that of humans. We therefore finished sequencing the male-specific region of the Y chromosome (MSY) in our closest living relative, the chimpanzee, achieving levels of accuracy and completion previously reached for the human MSY. We then compared the MSYs of the two species and found that they differ radically in sequence structure and gene content, implying rapid evolution during the past 6 million years. The chimpanzee MSY harbors twice as many massive palindromes as the human MSY, yet it has lost large fractions of the MSY protein-coding genes and gene families present in the last common ancestor. We suggest that the extraordinary divergence of the chimpanzee and human MSYs was driven by four synergistic factors: the MSY’s prominent role in sperm production, genetic hitchhiking effects in the absence of meiotic crossing over, frequent ectopic recombination within the MSY, and species differences in mating behavior. While genetic decay may be the principal dynamic in the evolution of newly emergent Y chromosomes, wholesale renovation is the paramount theme in the ongoing evolution of chimpanzee, human, and perhaps other older MSYs.

Elucidating population histories using genomic DNA sequences

In 1993, Cliff Jolly suggested that rather than debating species definitions and classifications, energy would be better spent investigating multidimensional patterns of variation and gene flow among populations. Until now, however, genetic studies of wild primate populations have been limited to very small portions of the genome. Access to complete genome sequences of humans, chimpanzees, macaques, and other primates makes it possible to design studies surveying substantial amounts of DNA sequence variation at multiple genetic loci in representatives of closely related but distinct wild primate populations. Such data can be analyzed with new approaches that estimate not only when populations diverged but also the relative amounts and directions of subsequent gene flow. These analyses will reemphasize the difficulty of achieving consistent species and subspecies definitions by revealing the extent of variation in the amount and duration of gene flow accompanying population divergences.

Genomic analyses of sex chromosome evolution

Sex chromosomes have evolved multiple times in many taxa. The recent explosion in the availability of whole genome sequences from a variety of organisms makes it possible to investigate sex chromosome evolution within and across genomes. Comparative genomic studies have shown that quite distant species may share fundamental properties of sex chromosome evolution, while very similar species can evolve unique sex chromosome systems. Furthermore, within-species genomic analyses can illuminate chromosome-wide sequence and expression polymorphisms. Here, we explore recent advances in the study of vertebrate sex chromosomes achieved using genomic analyses.