Sex chromosomes evolve distinctive types of chromatin from a pair of ancestral autosomes that are usually euchromatic. In Drosophila, the dosage-compensated X becomes enriched for hyperactive chromatin in males (mediated by H4K16ac), while the Y chromosome acquires silencing heterochromatin (enriched for H3K9me2/3). Drosophila autosomes are typically mostly euchromatic but the small dot chromosome has evolved a heterochromatin-like milieu (enriched for H3K9me2/3) that permits the normal expression of dot-linked genes, but which is different from typical pericentric heterochromatin. In Drosophila busckii, the dot chromosomes have fused to the ancestral sex chromosomes, creating a pair of 'neo-sex' chromosomes. Here we collect genomic, transcriptomic and epigenomic data from D. busckii, to investigate the evolutionary trajectory of sex chromosomes from a largely heterochromatic ancestor. We show that the neo-sex chromosomes formed <1 million years ago, but nearly 60% of neo-Y linked genes have already become non-functional. Expression levels are generally lower for the neo-Y alleles relative to their neo-X homologs, and the silencing heterochromatin mark H3K9me2, but not H3K9me3, is significantly enriched on silenced neo-Y genes. Despite rampant neo-Y degeneration, we find that the neo-X is deficient for the canonical histone modification mark of dosage compensation (H4K16ac), relative to autosomes or the compensated ancestral X chromosome, possibly reflecting constraints imposed on evolving hyperactive chromatin in an originally heterochromatic environment. Yet, neo-X genes are transcriptionally more active in males, relative to females, suggesting the evolution of incipient dosage compensation on the neo-X. Our data show that Y degeneration proceeds quickly after sex chromosomes become established through genomic and epigenetic changes, and are consistent with the idea that the evolution of sex-linked chromatin is influenced by its ancestral configuration.

Dosage compensation has been thought to be a ubiquitous property of sex chromosomes that are represented differently in males and females. The expression of most X-borne genes is equalized between XX females and XY males in therian mammals (marsupials and "placentals") by inactivating one X chromosome in female somatic cells. However, compensation seems not to be strictly required to equalize the expression of most Z-borne genes between ZZ male and ZW female birds. Whether dosage compensation operates in the third mammal lineage, the egg-laying monotremes, is of considerable interest, since the platypus has a complex sex chromosome system in which five X and five Y chromosomes share considerable genetic homology with the chicken ZW sex chromosome pair, but not with therian XY chromosomes. The assignment of genes to four platypus X chromosomes allowed us to examine X dosage compensation in this unique species. Quantitative PCR showed a range of compensation, but SNP analysis of several X-borne genes showed that both alleles are transcribed in a heterozygous female. Transcription of 14 BACs representing 19 X-borne genes was examined by RNA-FISH in female and male fibroblasts. An autosomal control gene was expressed from both alleles in nearly all nuclei, and four pseudoautosomal BACs were usually expressed from both alleles in male as well as female nuclei, showing that their Y loci are active. However, nine X-specific BACs were usually transcribed from only one allele. This suggests that while some genes on the platypus X are not dosage compensated, other genes do show some form of compensation via stochastic transcriptional inhibition, perhaps representing an ancestral system that evolved to be more tightly controlled in placental mammals such as human and mouse.

With a draft genome-sequence assembly for the chimpanzee available, it is now possible to perform genome-wide analyses to identify, at a submicroscopic level, structural rearrangements that have occurred between chimpanzees and humans. The goal of this study was to investigate chromosomal regions that are inverted between the chimpanzee and human genomes. Using the net alignments for the builds of the human and chimpanzee genome assemblies, we identified a total of 1,576 putative regions of inverted orientation, covering more than 154 mega-bases of DNA. The DNA segments are distributed throughout the genome and range from 23 base pairs to 62 mega-bases in length. For the 66 inversions more than 25 kilobases (kb) in length, 75% were flanked on one or both sides by (often unrelated) segmental duplications. Using PCR and fluorescence in situ hybridization we experimentally validated 23 of 27 (85%) semi-randomly chosen regions; the largest novel inversion confirmed was 4.3 mega-bases at human Chromosome 7p14. Gorilla was used as an out-group to assign ancestral status to the variants. All experimentally validated inversion regions were then assayed against a panel of human samples and three of the 23 (13%) regions were found to be polymorphic in the human genome. These polymorphic inversions include 730 kb (at 7p22), 13 kb (at 7q11), and 1 kb (at 16q24) fragments with a 5%, 30%, and 48% minor allele frequency, respectively. Our results suggest that inversions are an important source of variation in primate genome evolution. The finding of at least three novel inversion polymorphisms in humans indicates this type of structural variation may be a more common feature of our genome than previously realized.

Human recombination rates vary along the chromosomes as well as between the two sexes. There is growing evidence that epigenetic factors may have an important influence on recombination rates, as well as on crossover position. Using both public database analysis and wet-bench approaches, we revisited the relationship between increased rates of meiotic recombination and genome imprinting. We constructed metric linkage disequilibrium (LD) maps for all human chromosomal regions known to contain one or more imprinted genes. We show that imprinted regions contain significantly more LD units (LDU) and have significantly more haplotype blocks of smaller sizes than flanking nonimprinted regions. There is also an excess of hot-spots of recombination at imprinted regions, and this is likely to do with the presence of imprinted genes, per se. These findings indicate that imprinted chromosomal regions are historical "hot-spots" of recombination. We also demonstrate, by direct segregation analysis at the 11p15.5 imprinted region, that there is remarkable agreement between sites of meiotic recombination and steps in LD maps. Although the increase in LDU/Megabase at imprinted regions is not associated with any significant enrichment for any particular sequence class, major sequence determinants of recombination rates seem to differ between imprinted and control regions. Interestingly, fine-mapping of recombination events within the most male meiosis-specific recombination hot-spot of Chromosome 11p15.5 indicates that many events may occur within or directly adjacent to regions that are differentially methylated in somatic cells. Taken together, these findings support the involvement of a combination of specific DNA sequences and epigenetic factors as major determinants of hot-spots of recombination at imprinted chromosomal regions.

We investigated eight families with a novel subtype of congenital generalized lipodystrophy (CGL4) of whom five members had died from sudden cardiac death during their teenage years. ECG studies revealed features of long-QT syndrome, bradycardia, as well as supraventricular and ventricular tachycardias. Further symptoms comprised myopathy with muscle rippling, skeletal as well as smooth-muscle hypertrophy, leading to impaired gastrointestinal motility and hypertrophic pyloric stenosis in some children. Additionally, we found impaired bone formation with osteopenia, osteoporosis, and atlanto-axial instability. Homozygosity mapping located the gene within 2 Mbp on chromosome 17. Prioritization of 74 candidate genes with GeneDistiller for high expression in muscle and adipocytes suggested PTRF-CAVIN (Polymerase I and transcript release factor/Cavin) as the most probable candidate leading to the detection of homozygous mutations (c.160delG, c.362dupT). PTRF-CAVIN is essential for caveolae biogenesis. These cholesterol-rich plasmalemmal vesicles are involved in signal-transduction and vesicular trafficking and reside primarily on adipocytes, myocytes, and osteoblasts. Absence of PTRF-CAVIN did not influence abundance of its binding partner caveolin-1 and caveolin-3. In patient fibroblasts, however, caveolin-1 failed to localize toward the cell surface and electron microscopy revealed reduction of caveolae to less than 3%. Transfection of full-length PTRF-CAVIN reestablished the presence of caveolae. The loss of caveolae was confirmed by Atomic Force Microscopy (AFM) in combination with fluorescent imaging. PTRF-CAVIN deficiency thus presents the phenotypic spectrum caused by a quintessential lack of functional caveolae.

During meiotic prophase I, homologous chromosomes pair, synapse and recombine in a tightly regulated process that ensures the generation of genetically variable haploid gametes. Although the mechanisms underlying meiotic cell division have been well studied in model species, our understanding of the dynamics of meiotic prophase I in non-traditional model mammals remains in its infancy. Here, we reveal key meiotic features in previously uncharacterised marsupial species (the tammar wallaby and the fat-tailed dunnart), plus the fat-tailed mouse opossum, with a focus on sex chromosome pairing strategies, recombination and meiotic telomere homeostasis. We uncovered differences between phylogroups with important functional and evolutionary implications. First, sex chromosomes, which lack a pseudo-autosomal region in marsupials, had species specific pairing and silencing strategies, with implications for sex chromosome evolution. Second, we detected two waves of γH2AX accumulation during prophase I. The first wave was accompanied by low γH2AX levels on autosomes, which correlated with the low recombination rates that distinguish marsupials from eutherian mammals. In the second wave, γH2AX was restricted to sex chromosomes in all three species, which correlated with transcription from the X in tammar wallaby. This suggests non-canonical functions of γH2AX on meiotic sex chromosomes. Finally, we uncover evidence for telomere elongation in primary spermatocytes of the fat-tailed dunnart, a unique strategy within mammals. Our results provide new insights into meiotic progression and telomere homeostasis in marsupials, highlighting the importance of capturing the diversity of meiotic strategies within mammals.

Chromosome segregation errors during meiosis result in the formation of aneuploid gametes and are the leading cause of pregnancy loss and birth defects in humans. Proper chromosome segregation requires pairwise associations of maternal and paternal homologous chromosomes. Chiasmata, which are the cytological manifestations of crossovers (COs), provide a physical link that holds the homologs together as a pair, facilitating their orientation on the spindle at meiosis I. Although CO-promoting activities ensure a balanced number and position of COs, their identity and mechanism of action in mammals remain understudied. Previous work in yeast and Arabidopsis has shown that Zip2 and Shoc1 are ortholog proteins with an important role in promoting the formation of COs. Our work is the first study in mammals showing the in vivo and in vitro function of mouse and human SHOC1. We show that purified recombinant human SHOC1, an XPF/MUS81 family member, preferentially binds branched DNA molecules but apparently lacks in vitro endonuclease activity, despite its conserved ERCC4-(HhH)2 core structure. Cytological observations suggest that initial steps of recombination are normal in a majority of spermatocytes from SHOC1 hypomorphic mice. However, late stages of recombination appear abnormal, as chromosomal localization of MLH1 is reduced. In agreement, chiasma formation is reduced, and cells arrest at metaphase I with a few lagging chromosomes and subsequent apoptosis. This analysis of SHOC1-deficient mice and the selective localization of SHOC1 to a subset of recombination sites show that SHOC1 acts at key mid-stage steps of the CO formation process. The formation of chromosome axial elements and homologous pairing are apparently normal, but synapsis is altered with SYCP1 frequently failing to extend the full length of the chromosome axes. Finally, we describe that SHOC1 interacts with TEX11, another protein important for the formation of COs, connecting SHOC1 to chromosome axis and structure.

In many mammals, including humans and mice, the zinc finger histone methyltransferase PRDM9 performs the first step in meiotic recombination by specifying the locations of hotspots, the sites of genetic recombination. PRDM9 binds to DNA at hotspots through its zinc finger domain and activates recombination by trimethylating histone H3K4 on adjacent nucleosomes through its PR/SET domain. Recently, the isolated PR/SET domain of PRDM9 was shown capable of also trimethylating H3K36 in vitro, raising the question of whether this reaction occurs in vivo during meiosis, and if so, what its function might be. Here, we show that full-length PRDM9 does trimethylate H3K36 in vivo in mouse spermatocytes. Levels of H3K4me3 and H3K36me3 are highly correlated at hotspots, but mutually exclusive elsewhere. In vitro, we find that although PRDM9 trimethylates H3K36 much more slowly than it does H3K4, PRDM9 is capable of placing both marks on the same histone molecules. In accord with these results, we also show that PRDM9 can trimethylate both K4 and K36 on the same nucleosomes in vivo, but the ratio of K4me3/K36me3 is much higher for the pair of nucleosomes adjacent to the PRDM9 binding site compared to the next pair further away. Importantly, H3K4me3/H3K36me3-double-positive nucleosomes occur only in regions of recombination: hotspots and the pseudoautosomal (PAR) region of the sex chromosomes. These double-positive nucleosomes are dramatically reduced when PRDM9 is absent, showing that this signature is PRDM9-dependent at hotspots; the residual double-positive nucleosomes most likely come from the PRDM9-independent PAR. These results, together with the fact that PRDM9 is the only known mammalian histone methyltransferase with both H3K4 and H3K36 trimethylation activity, suggest that trimethylation of H3K36 plays an important role in the recombination process. Given the known requirement of H3K36me3 for double strand break repair by homologous recombination in somatic cells, we suggest that it may play the same role in meiosis.

Cells closely coordinate cell division with chromosome replication and segregation; however, the mechanisms responsible for this coordination still remain largely unknown. Here, we analyzed the spatial arrangement and temporal dynamics of the 9.1 Mb circular chromosome in the rod-shaped cells of Myxococcus xanthus. For chromosome segregation, M. xanthus uses a parABS system, which is essential, and lack of ParB results in chromosome segregation defects as well as cell divisions over nucleoids and the formation of anucleate cells. From the determination of the dynamic subcellular location of six genetic loci, we conclude that in newborn cells ori, as monitored following the ParB/parS complex, and ter regions are localized in the subpolar regions of the old and new cell pole, respectively and each separated from the nearest pole by approximately 1 µm. The bulk of the chromosome is arranged between the two subpolar regions, thus leaving the two large subpolar regions devoid of DNA. Upon replication, one ori region remains in the original subpolar region while the second copy segregates unidirectionally to the opposite subpolar region followed by the rest of the chromosome. In parallel, the ter region of the mother chromosome relocates, most likely passively, to midcell, where it is replicated. Consequently, after completion of replication and segregation, the two chromosomes show an ori-ter-ter-ori arrangement with mirror symmetry about a transverse axis at midcell. Upon completion of segregation of the ParB/parS complex, ParA localizes in large patches in the DNA-free subpolar regions. Using an Ssb-YFP fusion as a proxy for replisome localization, we observed that the two replisomes track independently of each other from a subpolar region towards ter. We conclude that M. xanthus chromosome arrangement and dynamics combine features from previously described systems with new features leading to a novel spatiotemporal arrangement pattern.

Malaria parasites elude eradication attempts both within the human host and across nations. At the individual level, parasites evade the host immune responses through antigenic variation. At the global level, parasites escape drug pressure through single nucleotide variants and gene copy amplification events conferring drug resistance. Despite their importance to global health, the rates at which these genomic alterations emerge have not been determined. We studied the complete genomes of different Plasmodium falciparum clones that had been propagated asexually over one year in the presence and absence of drug pressure. A combination of whole-genome microarray analysis and next-generation deep resequencing (totaling 14 terabases) revealed a stable core genome with only 38 novel single nucleotide variants appearing in seventeen evolved clones (avg. 5.4 per clone). In clones exposed to atovaquone, we found cytochrome b mutations as well as an amplification event encompassing the P. falciparum multidrug resistance associated protein (mrp1) on chromosome 1. We observed 18 large-scale (>1 kb on average) deletions of telomere-proximal regions encoding multigene families, involved in immune evasion (9.5×10(-6) structural variants per base pair per generation). Six of these deletions were associated with chromosomal crossovers generated during mitosis. We found only minor differences in rates between genetically distinct strains and between parasites cultured in the presence or absence of drug. Using these derived mutation rates for P. falciparum (1.0-9.7×10(-9) mutations per base pair per generation), we can now model the frequency at which drug or immune resistance alleles will emerge under a well-defined set of assumptions. Further, the detection of mitotic recombination events in var gene families illustrates how multigene families can arise and change over time in P. falciparum. These results will help improve our understanding of how P. falciparum evolves to evade control efforts within both the individual hosts and large populations.

The repeated rDNA array gives rise to the nucleolus, an organelle that is central to cellular processes as varied as stress response, cell cycle regulation, RNA modification, cell metabolism, and genome stability. The rDNA array is also responsible for the production of more than 70% of all cellular RNAs (the ribosomal RNAs). The rRNAs are produced from two sets of loci: the 5S rDNA array resides exclusively on human chromosome 1 while the 45S rDNA arrays reside on the short arm of five human acrocentric chromosomes. These critical genome elements have remained unassembled and have been excluded from all Hi-C analyses to date. Here we built the first high resolution map of 5S and 45S rDNA array contacts with the rest of the genome combining over 15 billion Hi-C reads from several experiments. The data enabled sufficiently high coverage to map rDNA-genome interactions with 1MB resolution and identify rDNA-gene contacts. The map showed that the 5S and 45S arrays display preferential contact at common sites along the genome but are not themselves sufficiently close to yield 5S-45S Hi-C contacts. Ribosomal DNA contacts are enriched in segments of closed, repressed, and late replicating chromatin, as well as CTCF binding sites. Finally, we identified functional categories whose dispersed genes coalesced in proximity to the rDNA arrays or instead avoided proximity with the rDNA arrays. The observations further our understanding of the spatial localization of rDNA arrays and their contribution to the architecture of the cell nucleus.

Canine hip dysplasia is a common, non-congenital, complex and hereditary disorder. It can inflict severe pain via secondary osteoarthritis and lead to euthanasia. An analogous disorder exists in humans. The genetic background of hip dysplasia in both species has remained ambiguous despite rigorous studies. We aimed to investigate the genetic causes of this disorder in one of the high-risk breeds, the German Shepherd. We performed genetic analyses with carefully phenotyped case-control cohorts comprising 525 German Shepherds. In our genome-wide association studies we identified four suggestive loci on chromosomes 1 and 9. Targeted resequencing of the two loci on chromosome 9 from 24 affected and 24 control German Shepherds revealed deletions of variable sizes in a putative enhancer element of the NOG gene. NOG encodes for noggin, a well-described bone morphogenetic protein inhibitor affecting multiple developmental processes, including joint development. The deletion was associated with the healthy controls and mildly dysplastic dogs suggesting a protective role against canine hip dysplasia. Two enhancer variants displayed a decreased activity in a dual luciferase reporter assay. Our study identifies novel loci and candidate genes for canine hip dysplasia, with potential regulatory variants in the NOG gene. Further research is warranted to elucidate how the identified variants affect the expression of noggin in canine hips, and what the potential effects of the other identified loci are.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Children affected by Specific Language Impairment (SLI) fail to acquire age appropriate language skills despite adequate intelligence and opportunity. SLI is highly heritable, but the understanding of underlying genetic mechanisms has proved challenging. In this study, we use molecular genetic techniques to investigate an admixed isolated founder population from the Robinson Crusoe Island (Chile), who are affected by a high incidence of SLI, increasing the power to discover contributory genetic factors. We utilize exome sequencing in selected individuals from this population to identify eight coding variants that are of putative significance. We then apply association analyses across the wider population to highlight a single rare coding variant (rs144169475, Minor Allele Frequency of 4.1% in admixed South American populations) in the NFXL1 gene that confers a nonsynonymous change (N150K) and is significantly associated with language impairment in the Robinson Crusoe population (p = 2.04 × 10-4, 8 variants tested). Subsequent sequencing of NFXL1 in 117 UK SLI cases identified four individuals with heterozygous variants predicted to be of functional consequence. We conclude that coding variants within NFXL1 confer an increased risk of SLI within a complex genetic model.

The Drosophila Nonspecific Lethal (NSL) complex is a major transcriptional regulator of housekeeping genes. It contains at least seven subunits that are conserved in the human KANSL complex: Nsl1/Wah (KANSL1), Dgt1/Nsl2 (KANSL2), Rcd1/Nsl3 (KANSL3), Rcd5 (MCRS1), MBD-R2 (PHF20), Wds (WDR5) and Mof (MOF/KAT8). Previous studies have shown that Dgt1, Rcd1 and Rcd5 are implicated in centrosome maintenance. Here, we analyzed the mitotic phenotypes caused by RNAi-mediated depletion of Rcd1, Rcd5, MBD-R2 or Wds in greater detail. Depletion of any of these proteins in Drosophila S2 cells led to defects in chromosome segregation. Consistent with these findings, Rcd1, Rcd5 and MBD-R2 RNAi cells showed reduced levels of both Cid/CENP-A and the kinetochore component Ndc80. In addition, RNAi against any of the four genes negatively affected centriole duplication. In Wds-depleted cells, the mitotic phenotypes were similar but milder than those observed in Rcd1-, Rcd5- or MBD-R2-deficient cells. RT-qPCR experiments and interrogation of published datasets revealed that transcription of many genes encoding centromere/kinetochore proteins (e.g., cid, Mis12 and Nnf1b), or involved in centriole duplication (e.g., Sas-6, Sas-4 and asl) is substantially reduced in Rcd1, Rcd5 and MBD-R2 RNAi cells, and to a lesser extent in wds RNAi cells. During mitosis, both Rcd1-GFP and Rcd5-GFP accumulate at the centrosomes and the telophase midbody, MBD-R2-GFP is enriched only at the chromosomes, while Wds-GFP accumulates at the centrosomes, the kinetochores, the midbody, and on a specific chromosome region. Collectively, our results suggest that the mitotic phenotypes caused by Rcd1, Rcd5, MBD-R2 or Wds depletion are primarily due to reduced transcription of genes involved in kinetochore assembly and centriole duplication. The differences in the subcellular localizations of the NSL components may reflect direct mitotic functions that are difficult to detect at the phenotypic level, because they are masked by the transcription-dependent deficiency of kinetochore and centriolar proteins.

Unisexual taxa are commonly considered short-lived as the absence of meiotic recombination is supposed to accumulate deleterious mutations and hinder the creation of genetic diversity. However, the gynogenetic gibel carp (Carassius gibelio) with high genetic diversity and wide ecological distribution has outlived its predicted extinction time of a strict unisexual reproduction population. Unlike other unisexual vertebrates, males associated with supernumerary microchromosomes have been observed in gibel carp, which provides a unique system to explore the rationales underlying male occurrence in unisexual lineage and evolution of unisexual reproduction. Here, we identified a massively expanded satellite DNA cluster on microchromosomes of hexaploid gibel carp via comparing with the ancestral tetraploid crucian carp (Carassius auratus). Based on the satellite cluster, we developed a method for single chromosomal fluorescence microdissection and isolated three male-specific microchromosomes in a male metaphase cell. Genomic anatomy revealed that these male-specific microchromosomes contained homologous sequences of autosomes and abundant repetitive elements. Significantly, several potential male-specific genes with transcriptional activity were identified, among which four and five genes displayed male-specific and male-biased expression in gonads, respectively, during the developmental period of sex determination. Therefore, the male-specific microchromosomes resembling common features of sex chromosomes may be the main driving force for male occurrence in gynogenetic gibel carp, which sheds new light on the evolution of unisexual reproduction.