Molecular mechanisms and topological consequences of drastic chromosomal rearrangements of muntjac deer

Muntjac deer have experienced drastic karyotype changes during their speciation, making it an ideal model for studying mechanisms and functional consequences of mammalian chromosome evolution. Here we generated chromosome-level genomes for Hydropotes inermis (2n = 70), Muntiacus reevesi (2n = 46), female and male M. crinifrons (2n = 8/9) and a contig-level genome for M. gongshanensis (2n = 8/9). These high-quality genomes combined with Hi-C data allowed us to reveal the evolution of 3D chromatin architectures during mammalian chromosome evolution. We find that the chromosome fusion events of muntjac species did not alter the A/B compartment structure and topologically associated domains near the fusion sites, but new chromatin interactions were gradually established across the fusion sites. The recently borne neo-Y chromosome of M. crinifrons, which underwent male-specific inversions, has dramatically restructured chromatin compartments, recapitulating the early evolution of canonical mammalian Y chromosomes. We also reveal that a complex structure containing unique centromeric satellite, truncated telomeric and palindrome repeats might have mediated muntjacs’ recurrent chromosome fusions. These results provide insights into the recurrent chromosome tandem fusion in muntjacs, early evolution of mammalian sex chromosomes, and reveal how chromosome rearrangements can reshape the 3D chromatin regulatory conformations during species evolution.

MRE11 and UBR5 Co-Operate to Suppress RNF168-Mediated Fusion of Dysfunctional Telomeres

TRF2 is part of the shelterin complex that hides telomeric DNA ends and prevents the activation of the cNHEJ pathway that can lead to chromosomal fusion. TRF2, however, also actively suppresses the cNHEJ pathway by recruiting two proteins, MRE11 and UBR5. MRE11 binds BRCC3, which in turn deubiquitinates γH2AX deposited at exposed telomeric DNA ends and limits RNF168 recruitment to the telomere. UBR5, in contrast directly ubiquitinates and destroys RNF168. The loss of telomeric RNF168 in turn blocks the subsequent recruitment of 53BP1 and prevents the cNHEJ-mediated fusion of chromosomes with exposed telomeric DNA ends. Although MRE11 and UBR5 are both involved in the control of telomeric RNF168 levels and the chromosome fusion process, their relative contributions have not been directly addressed. To do so we genetically suppressed MRE11 and UBR5 alone or in combination in glioma cell lines which we previously showed contained dysfunctional telomeres that were dependent on TRF2 for suppression of telomeric fusion and monitored the effects on events associated with telomere fusion. We here show that while suppression of either MRE11 or UBR5 alone had minimal effects on RNF168 telomeric accumulation, 53BP1 recruitment, and telomeric fusion, their combined suppression led to significant increases in RNF168 and 53BP1 telomeric recruitment and telomeric fusion and eventually cell death, all of which were reversible by suppression of RNF168 itself. These results show that MRE11 and UBR5 co-operate to suppress fusion at dysfunctional telomeres.

A high-quality chromosome-scale assembly of the centipedegrass [Eremochloa ophiuroides (Munro) Hack.] genome provides insights into chromosomal structural evolution and prostrate growth habit

Centipedegrass [Eremochloa ophiuroides (Munro) Hack.], a member of the Panicoideae subfamily, is one of the most important warm-season turfgrasses originating from China. This grass has an extremely developed prostrate growth habit and has been widely used in transitional and warm climatic regions. To better understand the genetic basis of important biological characteristics, such as prostrate growth and seed yield, in warm-season turfgrasses, we present a high-quality reference genome for centipedegrass and use PacBio, BioNano, and Hi-C technologies to anchor the 867.43 Mb genome assembly into nine pseudochromosomes, with a scaffold N50 of 86.05 Mb and 36,572 annotated genes. Centipedegrass was most closely related to sorghum and diverged from their common ancestor ~16.8 Mya. We detected a novel chromosome reshuffling event in centipedegrass, namely, the nest chromosome fusion event in which fusion of chromosomes 8 and 10 of sorghum into chromosome 3 of centipedegrass likely occurred after the divergence of centipedegrass from sorghum. The typical prostrate growth trait in centipedegrass may be linked to the expansion of candidate PROSTRATE GROWTH 1 (PROG1) genes on chromosome 2. Two orthologous genes of OsPROG1, EoPROG1, and EoPROG2, were confirmed to increase the stem number and decrease the stem angle in Arabidopsis. Collectively, our assembled reference genome of centipedegrass offers new knowledge and resources to dissect the genome evolution of Panicoideae and accelerate genome-assisted breeding and improvement of plant architecture in turf plants.

Epigenetic dynamics of centromeres and neocentromeres in Cryptococcus deuterogattii

Deletion of native centromeres in the human fungal pathogen Cryptococcus deuterogattii leads to neocentromere formation. Native centromeres span truncated transposable elements, while neocentromeres do not and instead span actively expressed genes. To explore the epigenetic organization of neocentromeres, we analyzed the distribution of the heterochromatic histone modification H3K9me2, 5mC DNA methylation and the euchromatin mark H3K4me2. Native centromeres are enriched for both H3K9me2 and 5mC DNA methylation marks and are devoid of H3K4me2, while neocentromeres do not exhibit any of these features. Neocentromeres in cen10Δ mutants are unstable and chromosome-chromosome fusions occur. After chromosome fusion, the neocentromere is inactivated and the native centromere of the chromosome fusion partner remains as the sole, active centromere. In the present study, the active centromere of a fused chromosome was deleted to investigate if epigenetic memory promoted the re-activation of the inactive neocentromere. Our results show that the inactive neocentromere is not re-activated and instead a novel neocentromere forms directly adjacent to the deleted centromere of the fused chromosome. To study the impact of transcription on centromere stability, the actively expressed URA5 gene was introduced into the CENP-A bound regions of a native centromere. The introduction of the URA5 gene led to a loss of CENP-A from the native centromere, and a neocentromere formed adjacent to the native centromere location. Remarkably, the inactive, native centromere remained enriched for heterochromatin, yet the integrated gene was expressed and devoid of H3K9me2. A cumulative analysis of multiple CENP-A distribution profiles revealed centromere drift in C. deuterogattii, a previously unreported phenomenon in fungi. The CENP-A-binding shifted within the ORF-free regions and showed a possible association with a truncated transposable element. Taken together, our findings reveal that neocentromeres in C. deuterogattii are highly unstable and are not marked with an epigenetic memory, distinguishing them from native centromeres.

Multiple recent sex chromosome fusions in Drosophila virilis associated with elevated satellite DNA abundance

Repetitive satellite DNA is highly variable both within and between species, and is often located near centromeres. However, the abundance or array length of satellite DNA may be constrained or have maximum limits. Drosophila virilis contains among the highest relative satellite abundances, with almost half its genome composed of three related 7 bp satellites. We discovered a strain of D. virilis that has 15% more pericentromeric satellite DNA compared to other strains, and also underwent two independent centromere-to-centromere sex chromosome fusion events. These fusions are presumably caused by DNA breakage near the pericentromeric satellites followed by repair using similar repetitive regions of nonhomologous chromosomes. We hypothesized that excess satellite DNA might increase the risk of DNA breaks and genome instability when stressed, which would be consistent with the apparent high rate of fusions we found in this strain. To directly quantify DNA breakage levels between strains with different satellite DNA abundances, we performed the comet assay after feeding flies gemcitabine and administering low-dose gamma radiation. We found a positive correlation between the rate of DNA breakage and satellite DNA abundance. This was further supported by a significant decrease in DNA breakage in an otherwise genetically identical substrain that lost the chromosome fusion and several megabases of satellite DNA. We find that the centromere-to-centromere fusions resulted in up to a 21% nondisjunction rate between the X and Y chromosomes in males, adding a fitness cost. Finally, we propose a model consistent with our data that implicates genome instability as a critical evolutionary constraint to satellite abundance.

Diploid chromosome-scale assembly of the Muscadinia rotundifolia genome supports chromosome fusion and disease resistance gene expansion during Vitis and Muscadinia divergence

Muscadinia rotundifolia, the muscadine grape, has been cultivated for centuries in the southeastern United States. M. rotundifolia is resistant to many of the pathogens that detrimentally affect Vitis vinifera, the grape species commonly used for winemaking. For this reason, M. rotundifolia is a valuable genetic resource for breeding. Single-molecule real-time reads were combined with optical maps to reconstruct the two haplotypes of each of the 20 M. rotundifolia cv. Trayshed chromosomes. The completeness and accuracy of the assembly were confirmed using a high-density linkage map of M. rotundifolia. Protein-coding genes were annotated using an integrated and comprehensive approach. This included using Full-length cDNA sequencing (Iso-Seq) to improve gene structure and hypothetical spliced variant predictions. Our data strongly support that Muscadinia chromosomes 7 and 20 are fused in Vitis and pinpoint the location of the fusion in Cabernet Sauvignon and PN40024 chromosome 7. Disease-related gene numbers in Trayshed and Cabernet Sauvignon were similar, but their clustering locations were different. A dramatic expansion of the Toll/Interleukin-1 Receptor-like Nucleotide-Binding Site Leucine-Rich Repeat (TIR-NBS-LRR) class was detected on Trayshed chromosome 12 at the Resistance to Uncinula necator 1 (RUN1)/ Resistance to Plasmopara viticola 1 (RPV1) locus, which confers strong dominant resistance to powdery and downy mildews. A genome browser for Trayshed, its annotation, and an associated Blast tool are available at .www.grapegenomics.com

Cytogenetics of Four Species of the Green Clade Aplastodiscus Lutz, 1950 (Anura: Cophomantinae): New Insights into the Chromosomal Evolution of the Genus

The tree frog <i>Aplastodiscus</i> is a Neotropical taxon that encompasses 15 species in the Atlantic forest biome, with one isolated species in the Central Brazilian Cerrado. To date, only 8 species have been karyotyped, showing high levels of diploid number variation, which allowed clustering species in chromosome number groups: 2n = 24 (<i>Aplastodiscus perviridis</i> group), 2n = 22 (<i>Aplastodiscus albofrenatus</i> group), 2n = 20, and 2n = 18 (both within <i>Aplastodiscus albosignatus</i> group). This study aims to report karyotypic information on 4 species from the last 2 groups using classical and molecular cytogenetic techniques and hypothesize chromosomal evolutionary trends within the species groups. <i>Aplastodiscus weygoldti</i> showed 2n = 22; Ag-NOR and FISH 18S rDNA signals were located in the interstitial region of the short arms of chromosome pair 6. <i>Aplastodiscus cavicola, Aplastodiscus</i> sp. 4, and <i>Aplastodiscus</i> sp. 6 showed 2n = 18; Ag-NOR and FISH 18S rDNA bands were located in the terminal region of the long arm of chromosome pair 9. Our results support multiple and independent chromosome fusion events within <i>Aplastodiscus</i>, including a new chromosome fission event<i>.</i> Ag-NOR and FISH 18S rDNA patterns were restricted to the small chromosome pairs, similar to the other species within this genus, and confirm overall chromosome morphology conservation among the genera of Cophomantinae.

The effects of chromosome fusions on genetic diversity and evolutionary turnover of functional loci consistently depends on chromosome size

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with the chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 21 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. Following fusion, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.