Inbreeding drives maize centromere evolution

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

High-quality Gossypium hirsutum and Gossypium barbadense genome assemblies reveal the landscape and evolution of centromeres

Centromere positioning and organization are crucial for genome evolution; however, research on centromere biology is largely influenced by the quality of available genome assemblies. Here, we combined Oxford Nanopore and Pacific Biosciences technologies to de novo assemble two high-quality reference genomes for Gossypium hirsutum (TM-1) and Gossypium barbadense (3-79). Compared with previously published reference genomes, our assemblies show substantial improvements, with the contig N50 improved by 4.6-fold and 5.6-fold, respectively, and thus represent the most complete cotton genomes to date. These high-quality reference genomes enable us to characterize 14 and 5 complete centromeric regions for G. hirsutum and G. barbadense, respectively. Our data revealed that the centromeres of allotetraploid cotton are occupied by members of the centromeric repeat for maize (CRM) and Tekay long terminal repeat families, and the CRM family reshapes the centromere structure of the At subgenome after polyploidization. These two intertwined families have driven the convergent evolution of centromeres between the two subgenomes, ensuring centromere function and genome stability. In addition, the repositioning and high sequence divergence of centromeres between G. hirsutum and G. barbadense have contributed to speciation and centromere diversity. This study sheds light on centromere evolution in a significant crop and provides an alternative approach for exploring the evolution of polyploid plants.

Plasticity of parental CENH3 incorporation into the centromeres in wheat × barley F1 hybrids

Incorporating the centromere-specific histone H3 protein CENH3 into the centromeric nucleosomes is indispensable for accurate centromere function and balanced chromosome segregation in most eukaryotes, including higher plants. In the cell nuclei of interspecific hybrids, divergent centromeric DNAs cohabit and lead the corresponding parental chromosomes through the mitotic and meiotic cell divisions. Depending on the transmission of the parental chromosomes carrying the CENH3-encoding genes, CENH3 proteins from one or both parents may be present in these hybrids. The incorporation of parental CENH3 proteins into the divergent centromeres and their role in the chromosome elimination process in interspecific hybrids is still poorly understood. Here, we produced wheat × barley F1 hybrids that carried different combinations of barley chromosomes with genes encoding for either one (αCENH3) or both barley CENH3 protein variants (α- and βCENH3). We generated specific antibodies distinguishing between the wheat CENH3 proteins and barley αCENH3 and applied them together with FISH probes to detect the precise pattern of parental CENH3 deposition into the wheat and barley centromeric nucleosomes. Analysis of somatic and meiotic nuclei of the wheat × barley hybrids revealed the plasticity of the maternal (wheat) CENH3 proteins to become incorporated into the paternal (barley) centromeric nucleosomes. However, no evidence for paternal CENH3 plasticity was detected in this study. The significance of the unilateral centromere plasticity and possible patterns of CENH3 incorporation into centromeres in interspecific hybrids are discussed.

Topological data analysis expands the genotype to phenotype map for 3D maize root system architecture

A central goal of biology is to understand how genetic variation produces phenotypic variation, which has been described as a genotype to phenotype (G to P) map. The plant form is continuously shaped by intrinsic developmental and extrinsic environmental inputs, and therefore plant phenomes are highly multivariate and require comprehensive approaches to fully quantify. Yet a common assumption in plant phenotyping efforts is that a few pre-selected measurements can adequately describe the relevant phenome space. Our poor understanding of the genetic basis of root system architecture is at least partially a result of this incongruence. Root systems are complex 3D structures that are most often studied as 2D representations measured with relatively simple univariate traits. In prior work, we showed that persistent homology, a topological data analysis method that does not pre-suppose the salient features of the data, could expand the phenotypic trait space and identify new G to P relations from a commonly used 2D root phenotyping platform. Here we extend the work to entire 3D root system architectures of maize seedlings from a mapping population that was designed to understand the genetic basis of maize-nitrogen relations. Using a panel of 84 univariate traits, persistent homology methods developed for 3D branching, and multivariate vectors of the collective trait space, we found that each method captures distinct information about root system variation as evidenced by the majority of non-overlapping QTL, and hence that root phenotypic trait space is not easily exhausted. The work offers a data-driven method for assessing 3D root structure and highlights the importance of non-canonical phenotypes for more accurate representations of the G to P map.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

Pan-centromere reveals widespread centromere repositioning of soybean genomes

Centromere repositioning refers to a de novo centromere formation at another chromosomal position without sequence rearrangement. This phenomenon was frequently encountered in both mammalian and plant species and has been implicated in genome evolution and speciation. To understand the dynamic of centromeres on soybean genome, we performed the pan-centromere analysis using CENH3-ChIP-seq data from 27 soybean accessions, including 3 wild soybeans, 9 landraces, and 15 cultivars. Building upon the previous discovery of three centromere satellites in soybean, we have identified two additional centromere satellites that specifically associate with chromosome 1. These satellites reveal significant rearrangements in the centromere structures of chromosome 1 across different accessions, consequently impacting the localization of CENH3. By comparative analysis, we reported a high frequency of centromere repositioning on 14 out of 20 chromosomes. Most newly emerging centromeres formed in close proximity to the native centromeres and some newly emerging centromeres were apparently shared in distantly related accessions, suggesting their emergence is independent. Furthermore, we crossed two accessions with mismatched centromeres to investigate how centromere positions would be influenced in hybrid genetic backgrounds. We found that a significant proportion of centromeres in the S9 generation undergo changes in size and position compared to their parental counterparts. Centromeres preferred to locate at satellites to maintain a stable state, highlighting a significant role of centromere satellites in centromere organization. Taken together, these results revealed extensive centromere repositioning in soybean genome and highlighted how important centromere satellites are in constraining centromere positions and supporting centromere function.

Centromere Plasticity With Evolutionary Conservation and Divergence Uncovered by Wheat 10+ Genomes

Centromeres (CEN) are the chromosomal regions that play a crucial role in maintaining genomic stability. The underlying highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of CEN in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat CEN. We found that most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core CEN. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned CEN. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize CEN. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

Centromere repositioning and shifts in wheat evolution

The centromere is the region of a chromosome that directs its separation and plays an important role in cell division and reproduction of organisms. Elucidating the dynamics of centromeres is an alternative strategy for exploring the evolution of wheat. Here, we comprehensively analyzed centromeres from the de novo-assembled common wheat cultivar Aikang58 (AK58), Chinese Spring (CS), and all sequenced diploid and tetraploid ancestors by chromatin immunoprecipitation sequencing, whole-genome bisulfite sequencing, RNA sequencing, assay for transposase-accessible chromatin using sequencing, and comparative genomics. We found that centromere-associated sequences were concentrated during tetraploidization and hexaploidization. Centromeric repeats of wheat (CRWs) have undergone expansion during wheat evolution, with strong interweaving between the A and B subgenomes post tetraploidization. We found that CENH3 prefers to bind with younger CRWs, as directly supported by immunocolocalization on two chromosomes (1A and 2A) of wild emmer wheat with dicentromeric regions, only one of which bound with CENH3. In a comparison of AK58 with CS, obvious centromere repositioning was detected on chromosomes 1B, 3D, and 4D. The active centromeres showed a unique combination of lower CG but higher CHH and CHG methylation levels. We also found that centromeric chromatin was more open than pericentromeric chromatin, with higher levels of gene expression but lower gene density. Frequent introgression between tetraploid and hexaploid wheat also had a strong influence on centromere position on the same chromosome. This study also showed that active wheat centromeres were genetically and epigenetically determined.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

Centromeres: From chromosome biology to biotechnology applications and synthetic genomes in plants

Centromeres are the genomic regions that organize and regulate chromosome behaviours during cell cycle, and their variations are associated with genome instability, karyotype evolution and speciation in eukaryotes. The highly repetitive and epigenetic nature of centromeres were documented during the past half century. With the aid of rapid expansion in genomic biotechnology tools, the complete sequence and structural organization of several plant and human centromeres were revealed recently. Here, we systematically summarize the current knowledge of centromere biology with regard to the DNA compositions and the histone H3 variant (CENH3)-dependent centromere establishment and identity. We discuss the roles of centromere to ensure cell division and to maintain the three-dimensional (3D) genomic architecture in different species. We further highlight the potential applications of manipulating centromeres to generate haploids or to induce polyploids offspring in plant for breeding programs, and of targeting centromeres with CRISPR/Cas for chromosome engineering and speciation. Finally, we also assess the challenges and strategies for de novo design and synthesis of centromeres in plant artificial chromosomes. The biotechnology applications of plant centromeres will be of great potential for the genetic improvement of crops and precise synthetic breeding in the future.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

The genetics and epigenetics of satellite centromeres

Centromeres, the chromosomal loci where spindle fibers attach during cell division to segregate chromosomes, are typically found within satellite arrays in plants and animals. Satellite arrays have been difficult to analyze because they comprise megabases of tandem head-to-tail highly repeated DNA sequences. Much evidence suggests that centromeres are epigenetically defined by the location of nucleosomes containing the centromere-specific histone H3 variant cenH3, independently of the DNA sequences where they are located; however, the reason that cenH3 nucleosomes are generally found on rapidly evolving satellite arrays has remained unclear. Recently, long-read sequencing technology has clarified the structures of satellite arrays and sparked rethinking of how they evolve, and new experiments and analyses have helped bring both understanding and further speculation about the role these highly repeated sequences play in centromere identification.

Characterizing mobile element insertions in 5675 genomes

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact

Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

Diverse mechanisms of centromere specification

The centromere performs a universally conserved function, to accurately partition genetic information upon cell division. Yet, centromeres are among the most rapidly evolving regions of the genome and are bound by a varying assortment of centromere-binding factors that are themselves highly divergent at the protein-sequence level. A common thread in most species is the dependence on the centromere-specific histone variant CENP-A for the specification of the centromere site. However, CENP-A is not universally required in all species or cell types, making the identification of a general mechanism for centromere specification challenging. In this review, we examine our current understanding of the mechanisms of centromere specification in CENP-A-dependent and independent systems, focusing primarily on recent work.

Oligo-FISH barcode in beans: a new chromosome identification system

An Oligo-FISH barcode system was developed for two model legumes, allowing the identification of all cowpea and common bean chromosomes in a single FISH experiment, and revealing new chromosome rearrangements. The FISH barcode system emerges as an effective tool to understand the chromosome evolution of economically important legumes and their related species. Current status on plant cytogenetic and cytogenomic research has allowed the selection and design of oligo-specific probes to individually identify each chromosome of the karyotype in a target species. Here, we developed the first chromosome identification system for legumes based on oligo-FISH barcode probes. We selected conserved genomic regions between Vigna unguiculata (Vu, cowpea) and Phaseolus vulgaris (Pv, common bean) (diverged ~ 9.7-15 Mya), using cowpea as a reference, to produce a unique barcode pattern for each species. We combined our oligo-FISH barcode pattern with a set of previously developed FISH probes based on BACs and ribosomal DNA sequences. In addition, we integrated our FISH maps with genome sequence data. Based on this integrated analysis, we confirmed two translocation events (involving chromosomes 1, 5, and 8; and chromosomes 2 and 3) between both species. The application of the oligo-based probes allowed us to demonstrate the participation of chromosome 5 in the translocation complex for the first time. Additionally, we detailed a pericentric inversion on chromosome 4 and identified a new paracentric inversion on chromosome 10. We also detected centromere repositioning associated with chromosomes 2, 3, 5, 7, and 9, confirming previous results for chromosomes 2 and 3. This first barcode system for legumes can be applied for karyotyping other Phaseolinae species, especially non-model, orphan crop species lacking genomic assemblies and cytogenetic maps, expanding our understanding of the chromosome evolution and genome organization of this economically important legume group.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

De novo centromere formation on chromosome fragments with an inactive centromere in maize (Zea mays)

The B chromosome of maize undergoes nondisjunction at the second pollen mitosis as part of its accumulation mechanism. Previous work identified 9-Bic-1 (9-B inactivated centromere-1), which comprises an epigenetically silenced B chromosome centromere that was translocated to the short arm of chromosome 9(9S). This chromosome is stable in isolation, but when normal B chromosomes are added to the genotype, it will attempt to undergo nondisjunction during the second pollen mitosis and usually fractures the chromosome in 9S. These broken chromosomes allow a test of whether the inactive centromere is reactivated or whether a de novo centromere is formed elsewhere on the chromosome to allow recovery of fragments. Breakpoint determination on the B chromosome and chromosome 9 showed that mini chromosome B1104 has the same breakpoint as 9-Bic-1 in the B centromere region and includes a portion of 9S. CENH3 binding was found on the B centromere region and on 9S, suggesting both centromere reactivation and de novo centromere formation. Another mini chromosome, B496, showed evidence of rearrangement, but it also only showed evidence for a de novo centromere. Other mini chromosome fragments recovered were directly derived from the B chromosome with breakpoints concentrated near the centromeric knob region, which suggests that the B chromosome is broken at a low frequency due to the failure of the sister chromatids to separate at the second pollen mitosis. Our results indicate that both reactivation and de novo centromere formation could occur on fragments derived from the progenitor possessing an inactive centromere.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

Comparative analyses of gibbon centromeres reveal dynamic genus specific shifts in repeat composition

Centromeres are functionally conserved chromosomal loci essential for proper chromosome segregation during cell division, yet they show high sequence diversity across species. Despite their variation, a near universal feature of centromeres is the presence of repetitive sequences, such as DNA satellites and transposable elements (TEs). Because of their rapidly evolving karyotypes, gibbons represent a compelling model to investigate divergence of functional centromere sequences across short evolutionary timescales. In this study, we use ChIP-seq, RNA-seq, and fluorescence in situ hybridization to comprehensively investigate the centromeric repeat content of the four extant gibbon genera (Hoolock, Hylobates, Nomascus, and Siamang). In all gibbon genera, we find that CENP-A nucleosomes and the DNA-proteins that interface with the inner kinetochore preferentially bind retroelements of broad classes rather than satellite DNA. A previously identified gibbon-specific composite retrotransposon, LAVA, known to be expanded within the centromere regions of one gibbon genus (Hoolock), displays centromere- and species-specific sequence differences, potentially as a result of its co-option to a centromeric function. When dissecting centromere satellite composition, we discovered the presence of the retroelement-derived macrosatellite SST1 in multiple centromeres of Hoolock, whereas alpha-satellites represent the predominate satellite in the other genera, further suggesting an independent evolutionary trajectory for Hoolock centromeres. Finally, using de novo assembly of centromere sequences, we determined that transcripts originating from gibbon centromeres recapitulate the species-specific TE composition. Combined, our data reveal dynamic shifts in the repeat content that define gibbon centromeres and coincide with the extensive karyotypic diversity within this lineage.

Sequence, Chromatin and Evolution of Satellite DNA

Satellite DNA consists of abundant tandem repeats that play important roles in cellular processes, including chromosome segregation, genome organization and chromosome end protection. Most satellite DNA repeat units are either of nucleosomal length or 5–10 bp long and occupy centromeric, pericentromeric or telomeric regions. Due to high repetitiveness, satellite DNA sequences have largely been absent from genome assemblies. Although few conserved satellite-specific sequence motifs have been identified, DNA curvature, dyad symmetries and inverted repeats are features of various satellite DNAs in several organisms. Satellite DNA sequences are either embedded in highly compact gene-poor heterochromatin or specialized chromatin that is distinct from euchromatin. Nevertheless, some satellite DNAs are transcribed into non-coding RNAs that may play important roles in satellite DNA function. Intriguingly, satellite DNAs are among the most rapidly evolving genomic elements, such that a large fraction is species-specific in most organisms. Here we describe the different classes of satellite DNA sequences, their satellite-specific chromatin features, and how these features may contribute to satellite DNA biology and evolution. We also discuss how the evolution of functional satellite DNA classes may contribute to speciation in plants and animals.

Population and subspecies diversity at mouse centromere satellites

BACKGROUND

Mammalian centromeres are satellite-rich chromatin domains that execute conserved roles in kinetochore assembly and chromosome segregation. Centromere satellites evolve rapidly between species, but little is known about population-level diversity across these loci.

RESULTS

We developed a k-mer based method to quantify centromere copy number and sequence variation from whole genome sequencing data. We applied this method to diverse inbred and wild house mouse (Mus musculus) genomes to profile diversity across the core centromere (minor) satellite and the pericentromeric (major) satellite repeat. We show that minor satellite copy number varies more than 10-fold among inbred mouse strains, whereas major satellite copy numbers span a 3-fold range. In contrast to widely held assumptions about the homogeneity of mouse centromere repeats, we uncover marked satellite sequence heterogeneity within single genomes, with diversity levels across the minor satellite exceeding those at the major satellite. Analyses in wild-caught mice implicate subspecies and population origin as significant determinants of variation in satellite copy number and satellite heterogeneity. Intriguingly, we also find that wild-caught mice harbor dramatically reduced minor satellite copy number and elevated satellite sequence heterogeneity compared to inbred strains, suggesting that inbreeding may reshape centromere architecture in pronounced ways.

CONCLUSION

Taken together, our results highlight the power of k-mer based approaches for probing variation across repetitive regions, provide an initial portrait of centromere variation across Mus musculus, and lay the groundwork for future functional studies on the consequences of natural genetic variation at these essential chromatin domains.

Maize centromeric chromatin scales with changes in genome size

Centromeres are defined by the location of Centromeric Histone H3 (CENP-A/CENH3) which interacts with DNA to define the locations and sizes of functional centromeres. An analysis of 26 maize genomes including 110 fully assembled centromeric regions revealed positive relationships between centromere size and genome size. These effects are independent of variation in the amounts of the major centromeric satellite sequence CentC. We also backcrossed known centromeres into two different lines with larger genomes and observed consistent increases in functional centromere sizes for multiple centromeres. Although changes in centromere size involve changes in bound CENH3, we could not mimic the effect by overexpressing CENH3 by threefold. Literature from other fields demonstrate that changes in genome size affect protein levels, organelle size and cell size. Our data demonstrate that centromere size is among these scalable features, and that multiple limiting factors together contribute to a stable centromere size equilibrium.

Grouping of genomic markers in populations with family structure

BACKGROUND

Linkage and linkage disequilibrium (LD) between genome regions cause dependencies among genomic markers. Due to family stratification in populations with non-random mating in livestock or crop, the standard measures of population LD such as [Formula: see text] may be biased. Grouping of markers according to their interdependence needs to account for the actual population structure in order to allow proper inference in genome-based evaluations.

RESULTS

Given a matrix reflecting the strength of association between markers, groups are built successively using a greedy algorithm; largest groups are built at first. As an option, a representative marker is selected for each group. We provide an implementation of the grouping approach as a new function to the R package hscovar. This package enables the calculation of the theoretical covariance between biallelic markers for half- or full-sib families and the derivation of representative markers. In case studies, we have shown that the number of groups comprising dependent markers was smaller and representative SNPs were spread more uniformly over the investigated chromosome region when the family stratification was respected compared to a population-LD approach. In a simulation study, we observed that sensitivity and specificity of a genome-based association study improved if selection of representative markers took family structure into account.

CONCLUSIONS

Chromosome segments which frequently recombine in the underlying population can be identified from the matrix of pairwise dependence between markers. Representative markers can be exploited, for instance, for dimension reduction prior to a genome-based association study or the grouping structure itself can be employed in a grouped penalization approach.

Multiple wheat genomes reveal global variation in modern breeding

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome¹, and the lack of genome-assembly data for multiple wheat lines^2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses^4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm1⁶, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.

Rapid Birth or Death of Centromeres on Fragmented Chromosomes in Maize

Comparative genomics has revealed common occurrences in karyotype evolution such as chromosomal end-to-end fusions and insertions of one chromosome into another near the centromere, as well as many cases of de novo centromeres that generate positional polymorphisms. However, how rearrangements such as dicentrics and acentrics persist without being destroyed or lost remains unclear. Here, we sought experimental evidence for the frequency and timeframe for inactivation and de novo formation of centromeres in maize (Zea mays). The pollen from plants with supernumerary B chromosomes was gamma-irradiated and then applied to normal maize silks of a line without B chromosomes. In ∼8,000 first-generation seedlings, we found many B-A translocations, centromere expansions, and ring chromosomes. We also found many dicentric chromosomes, but a fraction of these show only a single primary constriction, which suggests inactivation of one centromere. Chromosomal fragments were found without canonical centromere sequences, revealing de novo centromere formation over unique sequences; these were validated by immunolocalization with Thr133-phosphorylated histone H2A, a marker of active centromeres, and chromatin immunoprecipitation-sequencing with the CENH3 antibody. These results illustrate the regular occurrence of centromere birth and death after chromosomal rearrangement during a narrow window of one to potentially only a few cell cycles for the rearranged chromosomes to be recognized in this experimental regime.

Extraordinary Sequence Diversity and Promiscuity of Centromeric Satellites in the Legume Tribe Fabeae

Satellite repeats are major sequence constituents of centromeres in many plant and animal species. Within a species, a single family of satellite sequences typically occupies centromeres of all chromosomes and is absent from other parts of the genome. Due to their common origin, sequence similarities exist among the centromere-specific satellites in related species. Here, we report a remarkably different pattern of centromere evolution in the plant tribe Fabeae, which includes genera Pisum, Lathyrus, Vicia, and Lens. By immunoprecipitation of centromeric chromatin with CENH3 antibodies, we identified and characterized a large and diverse set of 64 families of centromeric satellites in 14 species. These families differed in their nucleotide sequence, monomer length (33-2,979 bp), and abundance in individual species. Most families were species-specific, and most species possessed multiple (2-12) satellites in their centromeres. Some of the repeats that were shared by several species exhibited promiscuous patterns of centromere association, being located within CENH3 chromatin in some species, but apart from the centromeres in others. Moreover, FISH experiments revealed that the same family could assume centromeric and noncentromeric positions even within a single species. Taken together, these findings suggest that Fabeae centromeres are not shaped by the coevolution of a single centromeric satellite with its interacting CENH3 proteins, as proposed by the centromere drive model. This conclusion is also supported by the absence of pervasive adaptive evolution of CENH3 sequences retrieved from Fabeae species.

Linkage Disequilibrium-Based Inference of Genome Homology and Chromosomal Rearrangements Between Species

The aim of this study was to analyze the genomic homology between cattle (Bos taurus) and buffaloes (Bubalus bubalis) and to propose a rearrangement of the buffalo genome through linkage disequilibrium analyses of buffalo SNP markers referenced in the cattle genome assembly and also compare it to the buffalo genome assembly. A panel of bovine SNPs (single nucleotide polymorphisms) was used for hierarchical, non-hierarchical and admixture cluster analyses. Thus, the linkage disequilibrium information between markers of a specific panel of buffalo was used to infer chromosomal rearrangement. Haplotype diversity and imputation accuracy of the submetacentric chromosomes were also analyzed. The genomic homology between the species enabled us to use the bovine genome assembly to recreate a buffalo genomic reference by rearranging the submetacentric chromosomes. The centromere of the submetacentric chromosomes exhibited high linkage disequilibrium and low haplotype diversity. It allowed hypothesizing about chromosome evolution. It indicated that buffalo submetacentric chromosomes are a centric fusion of ancestral acrocentric chromosomes. The chronology of fusions was also suggested. Moreover, a linear regression between buffalo and cattle rearranged assembly and the imputation accuracy indicated that the rearrangement of the chromosomes was adequate. When using the bovine reference genome assembly, the rearrangement of the buffalo submetacentric chromosomes could be done by SNP BTA (chromosome of Bos taurus) calculations: shorter BTA (shorter arm of buffalo chromosome) was given as [(shorter BTA length - SNP position in shorter BTA)] and larger BTA length as [shorter BTA length + (larger BTA length - SNP position in larger BTA)]. Finally, the proposed linkage disequilibrium-based method can be applied to elucidate other chromosomal rearrangement events in other species with the possibility of better understanding the evolutionary relationship between their genomes.

QTL mapping of spontaneous haploid genome doubling using genotyping-by-sequencing in maize (Zea mays L.)

A major QTL for SHGD was identified on chromosome 5 with stable expression across environments. The introgression this QTL can overcome the need of colchicine in DH lines development. Genome doubling of haploids is one of the major constraints of large-scale doubled haploid (DH) technology. Improving spontaneous haploid genome doubling (SHGD) is an alternative to overcome this limitation. In this study, we aimed to construct a high-density linkage map based on genotyping by sequencing of single nucleotide polymorphism, to detect QTL and QTL by environment (Q by E) interactions affecting SHGD and to identify the best trait for mapping and selection of haploid male fertility (HMF). To this end, a biparental population of 220 F_2:3 families was developed from a cross between A427 (high HMF) and CR1Ht (moderate HMF) to be used as donor. A high-density linkage map was constructed containing 4171 SNP markers distributed over 10 chromosomes with an average distance between adjacent markers of 0.51 cM. QTL mapping for haploid fertile anther emergence, pollen production, tassel size, and HMF, identified 27 QTL across three environments, and Q by E interactions were significant. A major QTL was identified on chromosome 5. This QTL explained over 45% of the observed variance for all traits across all environments. The introgression of this major QTL, using marker-assisted backcrossing, has great potential to overcome the need of using colchicine in DH line development.

Epigenetics and genome stability

Maintaining genome stability is essential to an organism's health and survival. Breakdown of the mechanisms protecting the genome and the resulting genome instability are an important aspect of the aging process and have been linked to diseases such as cancer. Thus, a large network of interconnected pathways is responsible for ensuring genome integrity in the face of the continuous challenges that induce DNA damage. While these pathways are diverse, epigenetic mechanisms play a central role in many of them. DNA modifications, histone variants and modifications, chromatin structure, and non-coding RNAs all carry out a variety of functions to ensure that genome stability is maintained. Epigenetic mechanisms ensure the functions of centromeres and telomeres that are essential for genome stability. Epigenetic mechanisms also protect the genome from the invasion by transposable elements and contribute to various DNA repair pathways. In this review, we highlight the integral role of epigenetic mechanisms in the maintenance of genome stability and draw attention to issues in need of further study.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha

Marchantia polymorpha has recently become a prime model for cellular, evo-devo, synthetic biological, and evolutionary investigations. We present a pseudomolecule-scale assembly of the M. polymorpha genome, making comparative genome structure analysis and classical genetic mapping approaches feasible. We anchored 88% of the M. polymorpha draft genome to a high-density linkage map resulting in eight pseudomolecules. We found that the overall genome structure of M. polymorpha is in some respects different from that of the model moss Physcomitrella patens. Specifically, genome collinearity between the two bryophyte genomes and vascular plants is limited, suggesting extensive rearrangements since divergence. Furthermore, recombination rates are greatest in the middle of the chromosome arms in M. polymorpha like in most vascular plant genomes, which is in contrast with P. patens where recombination rates are evenly distributed along the chromosomes. Nevertheless, some other properties of the genome are shared with P. patens. As in P. patens, DNA methylation in M. polymorpha is spread evenly along the chromosomes, which is in stark contrast with the angiosperm model Arabidopsis thaliana, where DNA methylation is strongly enriched at the centromeres. Nevertheless, DNA methylation and recombination rate are anticorrelated in all three species. Finally, M. polymorpha and P. patens centromeres are of similar structure and marked by high abundance of retroelements unlike in vascular plants. Taken together, the highly contiguous genome assembly we present opens unexplored avenues for M. polymorpha research by linking the physical and genetic maps, making novel genomic and genetic analyses, including map-based cloning, feasible.

Genome Evolution in Arabideae Was Marked by Frequent Centromere Repositioning

Centromere position may change despite conserved chromosomal collinearity. Centromere repositioning and evolutionary new centromeres (ENCs) were frequently encountered during vertebrate genome evolution but only rarely observed in plants. The largest crucifer tribe, Arabideae (∼550 species; Brassicaceae, the mustard family), diversified into several well-defined subclades in the virtual absence of chromosome number variation. Bacterial artificial chromosome-based comparative chromosome painting uncovered a constancy of genome structures among 10 analyzed genomes representing seven Arabideae subclades classified as four genera: Arabis, Aubrieta, Draba, and Pseudoturritis Interestingly, the intra-tribal diversification was marked by a high frequency of ENCs on five of the eight homoeologous chromosomes in the crown-group genera, but not in the most ancestral Pseudoturritis genome. From the 32 documented ENCs, at least 26 originated independently, including 4 ENCs recurrently formed at the same position in not closely related species. While chromosomal localization of ENCs does not reflect the phylogenetic position of the Arabideae subclades, centromere seeding was usually confined to long chromosome arms, transforming acrocentric chromosomes to (sub)metacentric chromosomes. Centromere repositioning is proposed as the key mechanism differentiating overall conserved homoeologous chromosomes across the crown-group Arabideae subclades. The evolutionary significance of centromere repositioning is discussed in the context of possible adaptive effects on recombination and epigenetic regulation of gene expression.

Artificial generation of centromeres and kinetochores to understand their structure and function

The centromere is an essential genomic region that provides the surface to form the kinetochore, which binds to the spindle microtubes to mediate chromosome segregation during mitosis and meiosis. Centromeres of most organisms possess highly repetitive sequences, making it difficult to study these loci. However, an unusual centromere called a "neocentromere," which does not contain repetitive sequences, was discovered in a patient and can be generated experimentally. Recent advances in genome biology techniques allow us to analyze centromeric chromatin using neocentromeres. In addition to neocentromeres, artificial kinetochores have been generated on non-centromeric loci, using protein tethering systems. These are powerful tools to understand the mechanism of the centromere specification and kinetochore assembly. In this review, we introduce recent studies utilizing the neocentromeres and artificial kinetochores and discuss current problems in centromere biology.

What makes a centromere?

Centromeres are the eukaryotic chromosomal sites at which the kinetochore forms and attaches to spindle microtubules to orchestrate chromosomal segregation in mitosis and meiosis. Although centromeres are essential for cell division, their sequences are not conserved and evolve rapidly. Centromeres vary dramatically in size and organization. Here we categorize their diversity and explore the evolutionary forces shaping them. Nearly all centromeres favor AT-rich DNA that is gene-free and transcribed at a very low level. Repair of frequent centromere-proximal breaks probably contributes to their rapid sequence evolution. Point centromeres are only ~125 bp and are specified by common protein-binding motifs, whereas short regional centromeres are 1-5 kb, typically have unique sequences, and may have pericentromeric repeats adapted to facilitate centromere clustering. Transposon-rich centromeres are often ~100-300 kb and are favored by RNAi machinery that silences transposons, by suppression of meiotic crossovers at centromeres, and by the ability of some transposons to target centromeres. Megabase-length satellite centromeres arise in plants and animals with asymmetric female meiosis that creates centromere competition, and favors satellite monomers one or two nucleosomes in length that position and stabilize centromeric nucleosomes. Holocentromeres encompass the length of a chromosome and may differ dramatically between mitosis and meiosis. We propose a model in which low level transcription of centromeres facilitates the formation of non-B DNA that specifies centromeres and promotes loading of centromeric nucleosomes.

Engineered minichromosomes in plants

Artificial chromosome platforms are described in plants. Because the function of centromeres is largely epigenetic, attempts to produce artificial chromosomes with plant centromere DNA have failed. The removal of the centromeric sequences from the cell strips off the centromeric histone that is the apparent biochemical marker of centromere activity. Thus, engineered minichromosomes have been produced by telomere mediated chromosomal truncation. The introduction of telomere repeats will cleave the chromosome at the site of insertion and attach the accompanying transgenes in the process. Such truncation events have been documented in maize, Arabidopsis, barley, rice, Brassica and wheat. Truncation of the nonvital supernumerary B chromosome of maize is a favorite target but engineered minichromosomes derived from the normal A chromosomes have also been recovered. Transmission through mitosis of small chromosomes is apparently normal but there is loss during meiosis. Potential solutions to address this issue are discussed. With procedures now well established to produce the foundation for artificial chromosomes in plants, current efforts are directed at building them up to specification using gene stacking methods and editing techniques.

Genetic and epigenetic effects on centromere establishment

Endogenous chromosomes contain centromeres to direct equal chromosomal segregation in mitosis and meiosis. The location and function of existing centromeres is usually maintained through cell cycles and generations. Recent studies have investigated how the centromere-specific histone H3 variant CENP-A is assembled and replenished after DNA replication to epigenetically propagate the centromere identity. However, existing centromeres occasionally become inactivated, with or without change in underlying DNA sequences, or lost after chromosomal rearrangements, resulting in acentric chromosomes. New centromeres, known as neocentromeres, may form on ectopic, non-centromeric chromosomal regions to rescue acentric chromosomes from being lost, or form dicentric chromosomes if the original centromere is still active. In addition, de novo centromeres can form after chromatinization of purified DNA that is exogenously introduced into cells. Here, we review the phenomena of naturally occurring and experimentally induced new centromeres and summarize the genetic (DNA sequence) and epigenetic features of these new centromeres. We compare the characteristics of new and native centromeres to understand whether there are different requirements for centromere establishment and propagation. Based on our understanding of the mechanisms of new centromere formation, we discuss the perspectives of developing more stably segregating human artificial chromosomes to facilitate gene delivery in therapeutics and research.

Plasticity in Triticeae centromere DNA sequences: a wheat × tall wheatgrass (decaploid) model

Centromeres mediate chromosome attachment to microtubules and maintain the integrity of chromosomes for proper segregation of the sister chromatids during cell division. Advances in the assembly of Triticeae genome sequences combined with the capacity to recover hybrid species derived from very distantly related species provides potential experimental systems for linking retrotransposon amplification and repositioning of centromeres via non-mendelian inheritance in partial amphiploid breeds. The decaploid tall wheatgrass (Thinopyrum ponticum) is one of the most successfully used perennial species in wheat breeding for generating translocation lines with valuable agronomic traits. We found that wheat centromere retrotransposons CRW and Quinta widely occur within the tall wheatgrass genome. In addition, one of the genome donors to Th. ponticum, Pseudoroegneria stipifolia (StSt), has been shown to have Abigail and a satellite repeat, CentSt. We also found two other centromeric retrotransposons, Abia and CL135 in Th. ponticum by ChIP-seq. Examination of partial amphiploid lines that were generated in the 1970s demonstrated extensive modification in centromere sequences using CentSt, Abigail and Abia as probes. We also detected that St-genome chromosomes were more enriched with Abigail and CentSt, whereas E-genome chromosomes were enriched with CRW and Quinta in tall wheatgrass and its closer relatives. It can be concluded that bursts of transposition of retrotransposons and repositioning of centromeres via non-mendelian segregation are common in partial amphiploids derived from interspecific hybrids. Practically speaking, our study reveals that the existence of homologous centromere functional sequences in both a donor and its receptor can substantially contribute to the successful transfer of alien genes into crop species. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/sra/SRR9089557; https://www.ncbi.nlm.nih.gov/sra/SRR9089558; https://www.ncbi.nlm.nih.gov/sra/SRR9089559; https://www.ncbi.nlm.nih.gov/sra/SRR9089560; https://www.ncbi.nlm.nih.gov/sra/SRR9089561; https://www.ncbi.nlm.nih.gov/sra/SRR9089562; https://www.ncbi.nlm.nih.gov/sra/SRR9089563; https://www.ncbi.nlm.nih.gov/sra/SRR9089564; https://www.ncbi.nlm.nih.gov/nuccore/MK999394; https://www.ncbi.nlm.nih.gov/nuccore/MK999395; https://www.ncbi.nlm.nih.gov/nuccore/MK999396.

Dynamic turnover of centromeres drives karyotype evolution in Drosophila

Centromeres are the basic unit for chromosome inheritance, but their evolutionary dynamics is poorly understood. We generate high-quality reference genomes for multiple Drosophila obscura group species to reconstruct karyotype evolution. All chromosomes in this lineage were ancestrally telocentric and the creation of metacentric chromosomes in some species was driven by de novo seeding of new centromeres at ancestrally gene-rich regions, independently of chromosomal rearrangements. The emergence of centromeres resulted in a drastic size increase due to repeat accumulation, and dozens of genes previously located in euchromatin are now embedded in pericentromeric heterochromatin. Metacentric chromosomes secondarily became telocentric in the pseudoobscura subgroup through centromere repositioning and a pericentric inversion. The former (peri)centric sequences left behind shrunk dramatically in size after their inactivation, yet contain remnants of their evolutionary past, including increased repeat-content and heterochromatic environment. Centromere movements are accompanied by rapid turnover of the major satellite DNA detected in (peri)centromeric regions.

The Genomic Landscape of Centromeres in Cancers

Centromere genomics remain poorly characterized in cancer, due to technologic limitations in sequencing and bioinformatics methodologies that make high-resolution delineation of centromeric loci difficult to achieve. We here leverage a highly specific and targeted rapid PCR methodology to quantitatively assess the genomic landscape of centromeres in cancer cell lines and primary tissue. PCR-based profiling of centromeres revealed widespread heterogeneity of centromeric and pericentromeric sequences in cancer cells and tissues as compared to healthy counterparts. Quantitative reductions in centromeric core and pericentromeric markers (α-satellite units and HERV-K copies) were observed in neoplastic samples as compared to healthy counterparts. Subsequent phylogenetic analysis of a pericentromeric endogenous retrovirus amplified by PCR revealed possible gene conversion events occurring at numerous pericentromeric loci in the setting of malignancy. Our findings collectively represent a more comprehensive evaluation of centromere genetics in the setting of malignancy, providing valuable insight into the evolution and reshuffling of centromeric sequences in cancer development and progression.

Cellular Dynamics and Genomic Identity of Centromeres in Cereal Blast Fungus

Magnaporthe oryzae is an important fungal pathogen that causes a loss of 10% to 30% of the annual rice crop due to the devastating blast disease. In most organisms, kinetochores are clustered together or arranged at the metaphase plate to facilitate synchronized anaphase separation of sister chromatids in mitosis. In this study, we showed that the initially clustered kinetochores separate and position randomly prior to anaphase in M. oryzae. Centromeres in M. oryzae occupy large genomic regions and form on AT-rich DNA without any common sequence motifs. Overall, this study identified atypical kinetochore dynamics and mapped functional centromeres in M. oryzae to define the roles of centromeric and pericentric boundaries in kinetochore assembly on epigenetically specified centromere loci. This study should pave the way for further understanding of the contribution of heterochromatin in genome stability and virulence of the blast fungus and its related species of high economic importance.

Precise kinetochore-microtubule interactions ensure faithful chromosome segregation in eukaryotes. Centromeres, identified as scaffolding sites for kinetochore assembly, are among the most rapidly evolving chromosomal loci in terms of the DNA sequence and length and organization of intrinsic elements. Neither the centromere structure nor the kinetochore dynamics is well studied in plant-pathogenic fungi. Here, we sought to understand the process of chromosome segregation in the rice blast fungus Magnaporthe oryzae. High-resolution imaging of green fluorescent protein (GFP)-tagged inner kinetochore proteins CenpA and CenpC revealed unusual albeit transient declustering of centromeres just before anaphase separation of chromosomes in M. oryzae. Strikingly, the declustered centromeres positioned randomly at the spindle midzone without an apparent metaphase plate per se. Using CenpA chromatin immunoprecipitation followed by deep sequencing, all seven centromeres in M. oryzae were found to be regional, spanning 57-kb to 109-kb transcriptionally poor regions. Highly AT-rich and heavily methylated DNA sequences were the only common defining features of all the centromeres in rice blast. Lack of centromere-specific DNA sequence motifs or repetitive elements suggests an epigenetic specification of centromere function in M. oryzae. PacBio genome assemblies and synteny analyses facilitated comparison of the centromeric/pericentromeric regions in distinct isolates of rice blast and wheat blast and in Magnaporthiopsis poae. Overall, this study revealed unusual centromere dynamics and precisely identified the centromere loci in the top model fungal pathogens that belong to Magnaporthales and cause severe losses in the global production of food crops and turf grasses.

Human Artificial Chromosomes that Bypass Centromeric DNA

Recent breakthroughs with synthetic budding yeast chromosomes expedite the creation of synthetic mammalian chromosomes and genomes. Mammals, unlike budding yeast, depend on the histone H3 variant, CENP-A, to epigenetically specify the location of the centromere-the locus essential for chromosome segregation. Prior human artificial chromosomes (HACs) required large arrays of centromeric α-satellite repeats harboring binding sites for the DNA sequence-specific binding protein, CENP-B. We report the development of a type of HAC that functions independently of these constraints. Formed by an initial CENP-A nucleosome seeding strategy, a construct lacking repetitive centromeric DNA formed several self-sufficient HACs that showed no uptake of genomic DNA. In contrast to traditional α-satellite HAC formation, the non-repetitive construct can form functional HACs without CENP-B or initial CENP-A nucleosome seeding, revealing distinct paths to centromere formation for different DNA sequence types. Our developments streamline the construction and characterization of HACs to facilitate mammalian synthetic genome efforts.

Islands of retroelements are major components of Drosophila centromeres

Centromeres are essential chromosomal regions that mediate kinetochore assembly and spindle attachments during cell division. Despite their functional conservation, centromeres are among the most rapidly evolving genomic regions and can shape karyotype evolution and speciation across taxa. Although significant progress has been made in identifying centromere-associated proteins, the highly repetitive centromeres of metazoans have been refractory to DNA sequencing and assembly, leaving large gaps in our understanding of their functional organization and evolution. Here, we identify the sequence composition and organization of the centromeres of Drosophila melanogaster by combining long-read sequencing, chromatin immunoprecipitation for the centromeric histone CENP-A, and high-resolution chromatin fiber imaging. Contrary to previous models that heralded satellite repeats as the major functional components, we demonstrate that functional centromeres form on islands of complex DNA sequences enriched in retroelements that are flanked by large arrays of satellite repeats. Each centromere displays distinct size and arrangement of its DNA elements but is similar in composition overall. We discover that a specific retroelement, G2/Jockey-3, is the most highly enriched sequence in CENP-A chromatin and is the only element shared among all centromeres. G2/Jockey-3 is also associated with CENP-A in the sister species D. simulans, revealing an unexpected conservation despite the reported turnover of centromeric satellite DNA. Our work reveals the DNA sequence identity of the active centromeres of a premier model organism and implicates retroelements as conserved features of centromeric DNA.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Comparison of Oryza sativa and Oryza brachyantha Genomes Reveals Selection-Driven Gene Escape from the Centromeric Regions

Centromeres are dynamic chromosomal regions, and the genetic and epigenetic environment of the centromere is often regarded as oppressive to protein-coding genes. Here, we used comparative genomic and phylogenomic approaches to study the evolution of centromeres and centromere-linked genes in the genus Oryza We report a 12.4-Mb high-quality BAC-based pericentromeric assembly for Oryza brachyantha, which diverged from cultivated rice (Oryza sativa) ∼15 million years ago. The synteny analyses reveal seven medium (>50 kb) pericentric inversions in O. sativa and 10 in O. brachyantha Of these inversions, three resulted in centromere movement (Chr1, Chr7, and Chr9). Additionally, we identified a potential centromere-repositioning event, in which the ancestral centromere on chromosome 12 in O. brachyantha jumped ∼400 kb away, possibly mediated by a duplicated transposition event (>28 kb). More strikingly, we observed an excess of syntenic gene loss at and near the centromeric regions (P < 2.2 × 10^-16). Most (33/47) of the missing genes moved to other genomic regions; therefore such excess could be explained by the selective loss of the copy in or near centromeric regions after gene duplication. The pattern of gene loss immediately adjacent to centromeric regions suggests centromere chromatin dynamics (e.g., spreading or microrepositioning) may drive such gene loss.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male–female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Sex-specific meiotic crossover (CO) landscapes have been identified in multiple species. Here, the authors show that male and female meioses in maize have similar CO landscapes, and differences between COs in the two sexes only exists in their location relative to transcription start sites and some chromatin marks.

Meiotic Studies on Combinations of Chromosomes With Different Sized Centromeres in Maize

Multiple centromere misdivision derivatives of a translocation between the supernumerary B chromosome and the short arm of chromosome 9 (TB-9Sb) permit investigation of how centromeres of different sizes behave in meiosis in opposition or in competition with each other. In the first analysis, heterozygotes were produced between the normal TB-9Sb and derivatives of it that resulted from centromere misdivision that reduced the amounts of centromeric DNA. These heterozygotes could test whether these drastic differences would result in meiotic drive of the larger chromosome in female meiosis. Cytological determinations of the segregation of large and small centromeres among thousands of progeny of four combinations were made. The recovery of the larger centromere was at a few percent higher frequency in two of four combinations. However, examination of phosphorylated histone H2A-Thr133, a characteristic of active centromeres, showed a lack of correlation with the size of the centromeric DNA, suggesting an expansion of the basal protein features of the kinetochore in two of the three cases despite the reduction in the size of the underlying DNA. In the second analysis, plants containing different sizes of the B chromosome centromere were crossed to plants with TB-9Sb with a foldback duplication of 9S (TB-9Sb-Dp9). In the progeny, plants containing large and small versions of the B chromosome centromere were selected by FISH. A meiotic "tug of war" occurred in hybrid combinations by recombination between the normal 9S and the foldback duplication in those cases in which pairing occurred. Such pairing and recombination produce anaphase I bridges but in some cases the large and small centromeres progressed to the same pole. In one combination, new dicentric chromosomes were found in the progeny. Collectively, the results indicate that the size of the underlying DNA of a centromere does not dramatically affect its segregation properties or its ability to progress to the poles in meiosis potentially because the biochemical features of centromeres adjust to the cellular conditions.