The history of the Y chromosome in man

Studies of the Y chromosome over the past few decades have opened a window into the history of our species, through the reconstruction and exploitation of a patrilineal (Y-genealogical) tree based on several hundred single-nucleotide variants (SNVs). A new study validates, refines and extends this tree by incorporating >65,000 Y-linked variants identified in 1,244 men representing worldwide diversity

Copy number variation arising from gene conversion on the human Y chromosome

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0–3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fbre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identifed at least 20 mutational events, and investigation of fanking paralogous sequence variants showed that the mutations involved fanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0–2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur

High divergence in primate-specific duplicated regions: Human and chimpanzee Chorionic Gonadotropin Beta genes

<p>Abstract</p> <p>Background</p> <p>Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (<it>LHB</it>)/Chorionic Gonadotropin Beta (<it>CGB</it>) gene cluster as a model. The placental <it>CGB </it>genes that are essential for implantation have evolved from an ancestral pituitary <it>LHB </it>gene by duplications in the primate lineage.</p> <p>Results</p> <p>We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) <it>LHB/CGB </it>regions and hereby present evidence for structural variation resulting in discordant number of <it>CGB </it>genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of <it>LHB/CGB </it>cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the <it>LHB/CGB </it>gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among <it>LHB/CGB </it>genes, there are signs of functional differentiation among the gene copies. Estimates for d_n/d_srate ratio suggested a purifying selection on <it>LHB </it>and <it>CGB8</it>, and a positive evolution of <it>CGB1</it>.</p> <p>Conclusion</p> <p>If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.</p

Recombination dynamics of a human Y-chromosomal palindrome:rapid GC-biased gene conversion, multi-kilobase conversion tracts, and rare inversions

The male-specific region of the human Y chromosome (MSY) includes eight large inverted repeats (palindromes) in which arm-to-arm similarity exceeds 99.9%, due to gene conversion activity. Here, we studied one of these palindromes, P6, in order to illuminate the dynamics of the gene conversion process. We genotyped ten paralogous sequence variants (PSVs) within the arms of P6 in 378 Y chromosomes whose evolutionary relationships within the SNP-defined Y phylogeny are known. This allowed the identification of 146 historical gene conversion events involving individual PSVs, occurring at a rate of 2.9-8.4×10(-4) events per generation. A consideration of the nature of nucleotide change and the ancestral state of each PSV showed that the conversion process was significantly biased towards the fixation of G or C nucleotides (GC-biased), and also towards the ancestral state. Determination of haplotypes by long-PCR allowed likely co-conversion of PSVs to be identified, and suggested that conversion tract lengths are large, with a mean of 2068 bp, and a maximum in excess of 9 kb. Despite the frequent formation of recombination intermediates implied by the rapid observed gene conversion activity, resolution via crossover is rare: only three inversions within P6 were detected in the sample. An analysis of chimpanzee and gorilla P6 orthologs showed that the ancestral state bias has existed in all three species, and comparison of human and chimpanzee sequences with the gorilla outgroup confirmed that GC bias of the conversion process has apparently been active in both the human and chimpanzee lineages

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequenceas a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale-multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8-15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals speciesspecific branch length variation that may reflect differences in long-term generation times

Analysis of SNP profiles in patients with major depressive disorder

The present study focused on 91 single-nucleotide polymorphisms (SNPs) in 21 candidate genes to find associations with major depressive disorder (MDD). In total, 160 healthy controls and 177 patients with MDD were studied. We applied arrayed primer extension (APEX) based genotyping technology followed by association and haplotype analysis. SNPs in CCKAR, DRD1, DRD2, and HTR2C genes showed nom- inally significant associations with MDD. None of these associations remained significant after adjustment for multiple testing. Haplotype analysis revealed CCKAR haplotypes to be associated with MDD (global p=0.004). More precisely, we found the GAGT haplotype to be associated with increased risk for MDD (OR 7.42, 95% CI 2.13–25.85, p=0.002). This haplotype effect remained significant after Bonferroni correction (p=0.04 after Bonferroni’s adjustment). Altogether we were able to find some nominal associations, but due to small sample size these results should be taken as exploratory. However, the effect of GAGT haplotype on the CCKAR gene may be considered as increasing the risk for MDD

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

Abstract: 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

Analysis of SNP profiles in patients with major depressive disorder

The present study focused on 91 single-nucleotide polymorphisms (SNPs) in 21 candidate genes to find associations with major depressive disorder (MDD). In total, 160 healthy controls and 177 patients with MDD were studied. We applied arrayed primer extension (APEX) based genotyping technology followed by association and haplotype analysis. SNPs in CCKAR, DRD1, DRD2, and HTR2C genes showed nominally significant associations with MDD. None of these associations remained significant after adjustment for multiple testing. Haplotype analysis revealed CCKAR haplotypes to be associated with MDD (global p=0.004). More precisely, we found the GAGT haplotype to be associated with increased risk for MDD (OR 7.42, 95% CI 2.13–25.85, p=0.002). This haplotype effect remained significant after Bonferroni correction (p=0.04 after Bonferroni's adjustment). Altogether we were able to find some nominal associations, but due to small sample size these results should be taken as exploratory. However, the effect of GAGT haplotype on the CCKAR gene may be considered as increasing the risk for MDD

Polymorphisms in wolframin (WFS1) gene are possibly related to increased risk for mood disorders

Wolfram syndrome gene (WFS1) has been suggested to have a role in the susceptibility for mood disorders. A 26-fold increased risk for psychiatric disorders in WFS1 mutation carriers has been sug- gested. In this study we tested the hypothesis that the WFS1 gene is related to the risk for mood disorders. We analysed 28 single-nucleotide polymorphisms (SNPs) of the WFS1 gene in 224 unrelated patients with major depressive disorder and bipolar disorder and in 160 healthy control subjects. Patients were further stratified according to their comorbidity with anxiety disorders. We applied arrayed primer extension (APEX)-based genotyping technology followed by association and haplotype analysis. Five SNPs in the WFS1 gene were associated with major depressive disorder, and three SNPs with bipolar disorder. Haplotype analysis revealed a common GTA haplotype, formed by SNPs 684C/ G, 1185C/T and 1832G/A, conferring risk for affective disorders. Specifically, for major depression the GTA haplotype has an OR of 1.59 (p=0.01) and for bipolar disorder an OR of 1.89 (p=0.03). These results support the hypothesis that the WFS1 gene is involved in the genetic predisposition for mood disorders

High level of complexity and global diversity of the 3q29 locus revealed by optical mapping and long-read sequencing.

BACKGROUND: High sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood. METHODS: We employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints. RESULTS: In total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes. CONCLUSIONS: We have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one proband with del3q29S. These results should provide a better understanding of the underlying genetic architecture that contributes to the etiology of del3q29S and dup3q29S