Genome-wide recombination is essential for genome stability, evolution, and speciation. Mouse Tex11, an X-linked meiosis-specific gene, promotes meiotic recombination and chromosomal synapsis. Here, we report that TEX11 is mutated in infertile men with non-obstructive azoospermia and that an analogous mutation in the mouse impairs meiosis. Genetic screening of a large cohort of idiopathic infertile men reveals that TEX11 mutations, including frameshift and splicing acceptor site mutations, cause infertility in 1% of azoospermic men. Functional evaluation of three analogous human TEX11 missense mutations in transgenic mouse models identified one mutation (V748A) as a potential infertility allele and found two mutations non-causative. In the mouse model, an intronless autosomal Tex11 transgene functionally substitutes for the X-linked Tex11 gene, providing genetic evidence for the X-to-autosomal retrotransposition evolution phenomenon. Furthermore, we find that TEX11 protein levels modulate genome-wide recombination rates in both sexes. These studies indicate that TEX11 alleles affecting expression level or substituting single amino acids may contribute to variations in recombination rates between sexes and among individuals in humans.

Amplicons--large, nearly identical repeats in direct or inverted orientation--are abundant in the male-specific region of the human Y chromosome (MSY) and provide targets for intrachromosomal non-allelic homologous recombination (NAHR). Thus far, NAHR events resulting in deletions, duplications, inversions, or isodicentric chromosomes have been reported only for amplicon pairs located exclusively on the short arm (Yp) or the long arm (Yq). Here we report our finding of four men with Y chromosomes that evidently formed by intrachromosomal NAHR between inverted repeat pairs comprising one amplicon on Yp and one amplicon on Yq. In two men with spermatogenic failure, sister-chromatid crossing-over resulted in pseudoisoYp chromosome formation and loss of distal Yq. In two men with normal spermatogenesis, intrachromatid crossing-over generated pericentric inversions. These findings highlight the recombinogenic nature of the MSY, as intrachromosomal NAHR occurs for nearly all Y-chromosome amplicon pairs, even those located on opposing chromosome arms.

Although the mammalian X and Y chromosomes evolved from a single pair of autosomes, they are highly differentiated: the Y chromosome is dramatically smaller than the X and has lost most of its genes. The surviving genes are a specialized set with extraordinary evolutionary longevity. Most mammalian lineages have experienced delayed, or relatively recent, loss of at least one conserved Y-linked gene. An extreme example of this phenomenon is in the Japanese spiny rat, where the Y chromosome has disappeared altogether. In this species, many Y-linked genes were rescued by transposition to new genomic locations, but until our work presented here, this has been considered an isolated case.

Little is known about how human Y-Chromosome gene expression directly contributes to differences between XX (female) and XY (male) individuals in nonreproductive tissues. Here, we analyzed quantitative profiles of Y-Chromosome gene expression across 36 human tissues from hundreds of individuals. Although it is often said that Y-Chromosome genes are lowly expressed outside the testis, we report many instances of elevated Y-Chromosome gene expression in a nonreproductive tissue. A notable example is EIF1AY, which encodes eukaryotic translation initiation factor 1A Y-linked, together with its X-linked homolog EIF1AX Evolutionary loss of a Y-linked microRNA target site enabled up-regulation of EIF1AY, but not of EIF1AX, in the heart. Consequently, this essential translation initiation factor is nearly twice as abundant in male as in female heart tissue at the protein level. Divergence between the X and Y Chromosomes in regulatory sequence can therefore lead to tissue-specific Y-Chromosome-driven sex biases in expression of critical, dosage-sensitive regulatory genes.

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

It is widely believed that at least three nonoverlapping regions of the human Y chromosome-AZFa, AZFb, and AZFc ("azoospermia factors" a, b, and c)-are essential for normal spermatogenesis. These intervals are defined by interstitial Y-chromosome deletions that impair or extinguish spermatogenesis. Deletion breakpoints, mechanisms, and lengths, as well as inventories of affected genes, have been elucidated for deletions of AZFa and of AZFc but not for deletions of AZFb or of AZFb plus AZFc. We studied three deletions of AZFb and eight deletions of AZFb plus AZFc, as assayed by the STSs defining these intervals. Guided by Y-chromosome sequence, we localized breakpoints precisely and were able to sequence nine of the deletion junctions. Homologous recombination can explain seven of these deletions but not the remaining two. This fact and our discovery of breakpoint hotspots suggest that factors in addition to homology underlie these deletions. The deletions previously thought to define AZFb were found to extend from palindrome P5 to the proximal arm of palindrome P1, 1.5 Mb within AZFc. Thus, they do not define a genomic region separate from AZFc. We also found that the deletions of AZFb plus AZFc, as assayed by standard STSs heretofore available, in fact extend from P5 to the distal arm of P1 and spare distal AZFc. Both classes of deletions are massive: P5/proximal-P1 deletions encompass up to 6.2 Mb and remove 32 genes and transcripts; P5/distal-P1 deletions encompass up to 7.7 Mb and remove 42 genes and transcripts. To our knowledge, these are the largest of all human interstitial deletions for which deletion junctions and complete intervening sequence are available. The restriction of the associated phenotype to spermatogenic failure indicates the remarkable functional specialization of the affected regions of the Y chromosome.

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. Little is known about the recent evolution of the Y chromosome 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. These theories have been buttressed by partial sequence data from newly emergent plant and animal Y chromosomes, but they have not been tested in older, highly evolved Y chromosomes such as that of humans. Here we finished sequencing of 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. By comparing the MSYs of the two species we show that they differ radically in sequence structure and gene content, indicating rapid evolution during the past 6 million years. The chimpanzee MSY contains 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 prominent role of the MSY in sperm production, 'genetic hitchhiking' effects in the absence of meiotic crossing over, frequent ectopic recombination within the MSY, and species differences in mating behaviour. Although genetic decay may be the principal dynamic in the evolution of newly emergent Y chromosomes, wholesale renovation is the paramount theme in the continuing evolution of chimpanzee, human and perhaps other older MSYs.

Indian muntjac (Muntiacus muntjak vaginalis) has an extreme mammalian karyotype, with only six and seven chromosomes in the female and male, respectively. Chinese muntjac (Muntiacus reevesi) has a more typical mammalian karyotype, with 46 chromosomes in both sexes. Despite this disparity, the two muntjac species are morphologically similar and can even interbreed to produce viable (albeit sterile) offspring. Previous studies have suggested that a series of telocentric chromosome fusion events involving telomeric and/or satellite repeats led to the extant Indian muntjac karyotype.

Xcat mice display X-linked congenital cataracts and are a mouse model for the human X-linked cataract disease Nance Horan syndrome (NHS). The genetic defect in Xcat mice and NHS patients is not known. We isolated and sequenced a BAC contig representing a portion of the Xcat critical region. We combined our sequencing data with the most recent mouse sequence assemblies from both Celera and public databases. The sequence of the 2.2-Mb Xcat critical region was then analyzed for potential Xcat candidate genes. The coding regions of the seven known genes within this area (Rai2, Rbbp7, Ctps2, Calb3, Grpr, Reps2, and Syap1) were sequenced in Xcat mice and no mutations were detected. The expression of Rai2 was quantitatively identical in wild-type and Xcat mutant eyes. These results indicate that the Xcat mutation is within a novel, undiscovered gene.

Amplicons-large, highly identical segmental duplications-are a prominent feature of mammalian Y chromosomes. Although they encode genes essential for fertility, these amplicons differ vastly between species, and little is known about the selective constraints acting on them. Here, we develop computational tools to detect amplicon copy number with unprecedented accuracy from high-throughput sequencing data. We find that one-sixth (16.9%) of 1,216 males from the 1000 Genomes Project have at least one deleted or duplicated amplicon. However, each amplicon's reference copy number is scrupulously maintained among divergent branches of the Y chromosome phylogeny, including the ancient branch A00, indicating that the reference copy number is ancestral to all modern human Y chromosomes. Using phylogenetic analyses and simulations, we demonstrate that this pattern of variation is incompatible with neutral evolution and instead displays hallmarks of mutation-selection balance. We also observe cases of amplicon rescue, in which deleted amplicons are restored through subsequent duplications. These results indicate that, contrary to the lack of constraint suggested by the differences between species, natural selection has suppressed amplicon copy number variation in diverse human lineages.

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Genome-wide (GWAS) and copy number variant (CNV) association studies have reproducibly identified numerous risk alleles associated with bipolar disorder (BD), major depressive disorder (MDD), and schizophrenia (SCZ), but biological characterization of these alleles lags gene discovery, owing to the inaccessibility of live human brain cells and inadequate animal models for human psychiatric conditions. Human-derived induced pluripotent stem cells (iPSCs) provide a renewable cellular reagent that can be differentiated into living, disease-relevant cells and 3D brain organoids carrying the full complement of genetic variants present in the donor germline. Experimental studies of iPSC-derived cells allow functional characterization of risk alleles, establishment of causal relationships between genes and neurobiology, and screening for novel therapeutics. Here we report the creation and availability of an iPSC resource comprising clinical, genomic, and cellular data obtained from genetically isolated families with BD and related conditions. Results from the first 324 study participants, 61 of whom have validated pluripotent clones, show enrichment of rare single nucleotide variants and CNVs overlapping many known risk genes and pathogenic CNVs. This growing iPSC resource is available to scientists pursuing functional genomic studies of BD and related conditions.

Somatic cells of human males and females have 45 chromosomes in common, including the "active" X chromosome. In males the 46th chromosome is a Y; in females it is an "inactive" X (Xi). Through linear modeling of autosomal gene expression in cells from individuals with zero to three Xi and zero to four Y chromosomes, we found that Xi and Y impact autosomal expression broadly and with remarkably similar effects. Studying sex chromosome structural anomalies, promoters of Xi- and Y-responsive genes, and CRISPR inhibition, we traced part of this shared effect to homologous transcription factors-ZFX and ZFY-encoded by Chr X and Y. This demonstrates sex-shared mechanisms by which Xi and Y modulate autosomal expression. Combined with earlier analyses of sex-linked gene expression, our studies show that 21% of all genes expressed in lymphoblastoid cells or fibroblasts change expression significantly in response to Xi or Y chromosomes.

Spinal muscular atrophy (SMA) is one of the most common severe hereditary diseases of infancy and early childhood in North America, Europe, and Asia. SMA is usually caused by deletions of the survival motor neuron 1 (SMN1) gene. A closely related gene, SMN2, modifies the disease severity. SMA carriers have only 1 copy of SMN1 and are relatively common (1 in 30-50) in populations of European and Asian descent. SMN copy numbers and SMA carrier frequencies have not been reliably estimated in Malians and other sub-Saharan Africans.

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.

EndoC-βH1 is emerging as a critical human β cell model to study the genetic and environmental etiologies of β cell (dys)function and diabetes. Comprehensive knowledge of its molecular landscape is lacking, yet required, for effective use of this model. Here, we report chromosomal (spectral karyotyping), genetic (genotyping), epigenomic (ChIP-seq and ATAC-seq), chromatin interaction (Hi-C and Pol2 ChIA-PET), and transcriptomic (RNA-seq and miRNA-seq) maps of EndoC-βH1. Analyses of these maps define known (e.g., PDX1 and ISL1) and putative (e.g., PCSK1 and mir-375) β cell-specific transcriptional cis-regulatory networks and identify allelic effects on cis-regulatory element use. Importantly, comparison with maps generated in primary human islets and/or β cells indicates preservation of chromatin looping but also highlights chromosomal aberrations and fetal genomic signatures in EndoC-βH1. Together, these maps, and a web application we created for their exploration, provide important tools for the design of experiments to probe and manipulate the genetic programs governing β cell identity and (dys)function in diabetes.

The introduction of foreign DNA into cells and organisms has facilitated much of modern biological research, and it promises to become equally important in clinical practice. Locating sites of foreign DNA incorporation in mammalian genomes has proven burdensome, so the genomic location of most transgenes remains unknown. To address this challenge, we applied nanopore sequencing in search of the site of integration of Tg(Pou5f1-EGFP)2Mnn (also known as Oct4:EGFP), a widely used fluorescent reporter in mouse germ line research. Using this nanopore-based approach, we identified the site of Oct4:EGFP transgene integration near the telomere of Chromosome 9. This methodology simultaneously yielded an estimate of transgene copy number, provided direct evidence of transgene inversions, revealed contaminating E. coli genomic DNA within the transgene array, validated the integrity of neighboring genes, and enabled definitive genotyping. We suggest that such an approach provides a rapid, cost-effective method for identifying and analyzing transgene integration sites.

Chediak-Higashi Syndrome (CHS) is a lysosome-related organelle (LRO) disorder caused by biallelic mutations in the lysosomal trafficking regulator gene, LYST. The clinical features of CHS include oculocutaneous albinism, primary immunodeficiency, bleeding diathesis, risk for development of hemophagocyticlymphohistiocytosis,and progressive neurological problems. The pathophysiological mechanisms underlying this disease are unknown, so developing therapeutic options remains challenging. In this study,four induced pluripotent stem (iPSC) lines from unrelated CHS patients have been generated and successfully characterized for exploring the role of LYST in health and disease in diversecell types.

The meiotic prophase I to metaphase I (PI/MI) transition requires chromosome desynapsis and metaphase competence acquisition. However, control of these major meiotic events is poorly understood. Here, we identify an essential role for SKP1, a core subunit of the SKP1-Cullin-F-box (SCF) ubiquitin E3 ligase, in the PI/MI transition. SKP1 localizes to synapsed chromosome axes and evicts HORMAD proteins from these regions in meiotic spermatocytes. SKP1-deficient spermatocytes display premature desynapsis, precocious pachytene exit, loss of PLK1 and BUB1 at centromeres, but persistence of HORMAD, γH2AX, RPA2, and MLH1 in diplonema. Strikingly, SKP1-deficient spermatocytes show sharply reduced MPF activity and fail to enter MI despite treatment with okadaic acid. SKP1-deficient oocytes exhibit desynapsis, chromosome misalignment, and progressive postnatal loss. Therefore, SKP1 maintains synapsis in meiosis of both sexes. Furthermore, our results support a model where SKP1 functions as the long-sought intrinsic metaphase competence factor to orchestrate MI entry during male meiosis.

Glucocerebrosidase is a lysosomal hydrolase involved in the breakdown of glucosylceramide. Gaucher disease, a recessive lysosomal storage disorder, is caused by mutations in the gene GBA1 Dysfunctional glucocerebrosidase leads to accumulation of glucosylceramide and glycosylsphingosine in various cell types and organs. Mutations in GBA1 are also a common genetic risk factor for Parkinson disease and related synucleinopathies. In recent years, research on the pathophysiology of Gaucher disease, the molecular link between Gaucher and Parkinson disease, and novel therapeutics, have accelerated the need for relevant cell models with GBA1 mutations. Although induced pluripotent stem cells, primary rodent neurons, and transfected neuroblastoma cell lines have been used to study the effect of glucocerebrosidase deficiency on neuronal function, these models have limitations because of challenges in culturing and propagating the cells, low yield, and the introduction of exogenous mutant GBA1 To address some of these difficulties, we established a high yield, easy-to-culture mouse neuronal cell model with nearly complete glucocerebrosidase deficiency representative of Gaucher disease. We successfully immortalized cortical neurons from embryonic null allele gba(-/-) mice and the control littermate (gba(+/+)) by infecting differentiated primary cortical neurons in culture with an EF1α-SV40T lentivirus. Immortalized gba(-/-) neurons lack glucocerebrosidase protein and enzyme activity, and exhibit a dramatic increase in glucosylceramide and glucosylsphingosine accumulation, enlarged lysosomes, and an impaired ATP-dependent calcium-influx response; these phenotypical characteristics were absent in gba(+/+) neurons. This null allele gba(-/-) mouse neuronal model provides a much-needed tool to study the pathophysiology of Gaucher disease and to evaluate new therapies.