Maize Centromere Structure and Evolution: Sequence Analysis of Centromeres 2 and 5 Reveals Dynamic Loci Shaped Primarily by Retrotransposons

We describe a comprehensive and general approach for mapping centromeres and present a detailed characterization of two maize centromeres. Centromeres are difficult to map and analyze because they consist primarily of repetitive DNA sequences, which in maize are the tandem satellite repeat CentC and interspersed centromeric retrotransposons of maize (CRM). Centromeres are defined epigenetically by the centromeric histone H3 variant, CENH3. Using novel markers derived from centromere repeats, we have mapped all ten centromeres onto the physical and genetic maps of maize. We were able to completely traverse centromeres 2 and 5, confirm physical maps by fluorescence in situ hybridization (FISH), and delineate their functional regions by chromatin immunoprecipitation (ChIP) with anti-CENH3 antibody followed by pyrosequencing. These two centromeres differ substantially in size, apparent CENH3 density, and arrangement of centromeric repeats; and they are larger than the rice centromeres characterized to date. Furthermore, centromere 5 consists of two distinct CENH3 domains that are separated by several megabases. Succession of centromere repeat classes is evidenced by the fact that elements belonging to the recently active recombinant subgroups of CRM1 colonize the present day centromeres, while elements of the ancestral subgroups are also found in the flanking regions. Using abundant CRM and non-CRM retrotransposons that inserted in and near these two centromeres to create a historical record of centromere location, we show that maize centromeres are fluid genomic regions whose borders are heavily influenced by the interplay of retrotransposons and epigenetic marks. Furthermore, we propose that CRMs may be involved in removal of centromeric DNA (specifically CentC), invasion of centromeres by non-CRM retrotransposons, and local repositioning of the CENH3.     

The size and sequence organization of the centromeric region of arabidopsis thaliana chromosome 5.

We have determined the size of the centromeric region of Arabidopsis thaliana chromosome 5, which corresponds to the genetically defined centromere by Copenhaver et al. (Science, 286, 2468-2474, 1999) on the basis of restriction analysis. As a large clone gap was present in the previously constructed contig map of the centromeric region, the restriction map of this region was constructed using Asc I, Not I, Apa I and Pme I and genomic DNA from a hypomethylated strain. The size of the centromeric region finally estimated by combination with the sequence data of cloned regions at both sides was 4.35 megabases (Mb). This value is over 2 Mb longer than those estimated in our previous work and also by Copenhaver et al. Combing this centromeric region with the physical map previously constructed, the entire length of chromosome 5 becomes 31 Mb. Although the internal moiety of the centromeric region has not been sequenced yet because of extremely high repetition, the result of sequence analysis from both sides toward the inside strongly suggests that the centromeric region is composed of the central 2.9-Mb domain consisting of mainly 180-bp repeats and Athila retrotransposons and flanking regions containing various types of transposons. On the basis of these observations, a structural model for the centromeric region is discussed.     

Centromeric DNA from Saccharomyces uvarum is functional in Saccharomyces cerevisiae

Previous comparisons of centromeric DNA sequences in laboratory strains of Saccharomyces cerevisiae have revealed conserved sequences within 120 base pairs (bp) which appear to be essential for centromere function. We wanted to find out whether centromeric DNA in Saccharomyces strains with different degrees of DNA sequence divergence carry the same conserved sequences or not. Bam HI DNA fragments from two S. cerevisiae strains and one Saccharomyces uvarum strain were cloned into a centromere selection vector and tested for centromere function in a S. cerevisiae laboratory strain. Fragments having centromere function were obtained at approximately equal frequencies from all three strains. Two of the S. uvarum centromeric DNAs and two of the S. cerevisiae centromeric DNAs were sequenced and shown to carry in a 120 bp region sequences essentially like those of centromeric DNA in S. cerevisiae laboratory strains. DNA hybridization to separated chromosomal DNAs revealed that the two newly determined S. cerevisiae centromeric DNA sequences belong to chromosomes V and XIII, respectively. On leave from: Department of Cell and Tumor Biology, Roswell Park Memorial Institute, Buffalo, NY 14263, USA; On leave from: The Biological Laboratories, University of Leiden, The Netherlands     

Analysis of DNA restriction fragments greater than 5.7 Mb in size from the centromeric region of human chromosomes

Pulsed electrophoresis was used to study the organization of the human centromeric region. Genomic DNA was digested with rare-cutting enzymes. DNA fragments from 0.2 to greater than 5.7 Mb were separated by electrophoresis and hybridized with alphoid and simple DNA repeats. Rare-cutting enzymes (Mlu I, Nar I, Not I, Nru I, Sal I, Sfi I, Sst II) demonstrated fewer restriction sites at centromeric regions than elsewhere in the genome. The enzyme Not I had the fewest restriction sites at centromeric regions. As much as 70% of these sequences from the centromeric region are present in Not I DNA fragments greater than 5.7 and estimated to be as large as 10 Mb in size. Other repetitive sequences such as short interspersed repeated segments (SINEs), long interspersed repeated segments (LINEs), ribosomal DNA, and minisatellite DNA that are not enriched at the centromeric region, are not enriched in Not I fragments of greater than 5.7 Mb in size.     

Physical Mapping of Meiotic Crossover Events in a 200-Kb Region of Neurospora Crassa Linkage Group I

We propose a general restriction fragment length polymorphism-based strategy to analyze the distribution of meiotic crossover events throughout specific genetic intervals. We have isolated 64 recombinant chromosomes carrying independent meiotic crossover events in the genetic interval eth-1-un-2 on linkage group I of Neurospora crassa. Thirty-eight crossover events were physically mapped with reference to a 200-kb region cloned by chromosome walking, using N. crassa λ and cosmid libraries. Crossovers were homogeneously distributed at intervals of 5.0 +/- 2.3 kb along the entire cloned interval. The ratio of physical to genetic distance appears to be higher in the region than in the overall N. crassa genome, suggesting that recombinational activity is less in large chromosomes than in small ones. The present work provides a method for defining the centromeric-telomeric orientation of single cloned DNA fragments. Their physical distance can also be estimated with respect to linked loci, provided that crossover events are distributed homogeneously in the interval. This strategy overcomes typical difficulties in defining the position and direction of chromosome walking steps on conventional linkage maps.     

Polymorphic Chromosomal Specificity of Centromere Satellite Families in Arabidopsis halleri ssp. gemmifera

The chromosomal localizations of repetitive DNA clusters (ribosomal DNA and centromere satellites) were analyzed by fluorescent in situ hybridization in five strains of Arabidopsis halleri ssp. gemmifera. All five A. gemmifera strains have three chromosome pairs with 45S (5.8S-16S-26S) rDNA loci, and one pair with both 5S and 45S rDNA loci. These localizations are different from that of A. thaliana. Very unusually, there are three families of centromeric satellite DNAs (pAa, pAge1, and pAge2), and they showed polymorphism among the five strains studied. Overall, we found four different centromere satellite compositions. A plant from Fumuro was heterozygous for the chromosome specificities of centromere satellite families, possibly due to a reciprocal translocation involving centromere regions. Changes of centromeric satellite repeats appear to be rapid and frequent events in the history of A. gemmifera, and seem to occur by exchanging clusters as units.     

A BAC library of <Emphasis Type="Italic">Beta vulgaris</Emphasis> L. for the targeted isolation of centromeric DNA and molecular cytogenetics of <Emphasis Type="Italic">Beta</Emphasis> species

We constructed a sugar beet (Beta vulgaris) bacterial artificial chromosome (BAC) library of the monosomic addition line PAT2. This chromosomal mutant carries a single additional chromosome fragment (minichromosome) derived from the wild beet Beta patellaris. Restriction analysis of the mutant line by pulsed-field gel electrophoresis was used to determine HindIII as a suitable enzyme for partial digestion of genomic DNA to generate large-insert fragments which were cloned into the vector pCC1. The library consists of 36,096 clones with an average insert size of 120 kb, and 2.2% of the clones contain mitochondrial or chloroplast DNA. Based on a haploid genome size of 758 Mbp, the library represents 5.7 genome equivalents providing the probability of 99.67% that any sequence of the PAT2 genome can be found in the library. Hybridization to high-density filters was used to isolate 89 BACs containing arrays of the centromere-associated satellite repeats pTS5 and pTS4.1. Using the identified BAC clones in fluorescent in situ hybridization experiments with PAT2 and Beta patellaris chromosome spreads their wild beet origin and centromeric localization was demonstrated. Multi-colour FISH with differently labelled satellite repeats pTS5 and pTS4.1 was used to investigate the large-scale organization of the centromere of the PAT2 minichromosome in detail. FISH studies showed that the centromeric satellite pTS5 is flanked on both sides by pTS4.1 arrays and the arms of the minichromosome are terminated by the Arabidopsis-type telomeric sequences. FISH with a BAC, selected from high-density filters after hybridization with an RFLP marker of the genetic linkage group I, demonstrated that it is feasible to correlate genetic linkage groups with chromosomes. Therefore, the PAT2 BAC library provides a useful tool for the characterization of Beta centromeres and a valuable resource for sugar beet genome analysis.     

Structure and dynamics of retrotransposons at wheat centromeres and pericentromeres

Little is known of the dynamics of centromeric DNA in polyploid plants. We report the sequences of two centromere-associated bacterial artificial chromosome clones from a Triticum boeoticum library. Both autonomous and non-autonomous wheat centromeric retrotransposons (CRWs) were identified, both being closely associated with the centromeres of wheat. Fiber-fluorescence in situ hybridization and chromatin immunoprecipitation analysis showed that wheat centromeric retrotransposons (CRWs) represent a dominant component of the wheat centromere and are associated with centromere function. CRW copy number showed variation among different genomes: the D genome chromosomes contained fewer copies than either the A or B genome chromosomes. The frequency of lengthy continuous CRW arrays was higher than that in either rice or maize. The dynamics of CRWs and other retrotransposons at centromeric and pericentromeric regions during diploid speciation and polyploidization of wheat and its related species are discussed. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Zhao Liu and Wei Yue made an equal contribution to this work.     

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres – an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half‐tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination – a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high‐density maps in species with strong recombination interference and will enrich many existing and future mapping resources.     

Cloning of human centromeres by transformation-associated recombination in yeast and generation of functional human artificial chromosomes

Human centromeres remain poorly characterized regions of the human genome despite their importance for the maintenance of chromosomes. In part this is due to the difficulty of cloning of highly repetitive DNA fragments and distinguishing chromosome-specific clones in a genomic library. In this work we report the highly selective isolation of human centromeric DNA using transformation-associated recombination (TAR) cloning. A TAR vector with alphoid DNA monomers as targeting sequences was used to isolate large centromeric regions of human chromosomes 2, 5, 8, 11, 15, 19, 21 and 22 from human cells as well as monochromosomal hybrid cells. The alphoid DNA array was also isolated from the 12 Mb human mini-chromosome ΔYq74 that contained the minimum amount of alphoid DNA required for proper chromosome segregation. Preliminary results of the structural analyses of different centromeres are reported in this paper. The ability of the cloned human centromeric regions to support human artificial chromosome (HAC) formation was assessed by transfection into human HT1080 cells. Centromeric clones from ΔYq74 did not support the formation of HACs, indicating that the requirements for the existence of a functional centromere on an endogenous chromosome and those for forming a de novo centromere may be distinct. A construct with an alphoid DNA array from chromosome 22 with no detectable CENP-B motifs formed mitotically stable HACs in the absence of drug selection without detectable acquisition of host DNAs. In summary, our results demonstrated that TAR cloning is a useful tool for investigating human centromere organization and the structural requirements for formation of HAC vectors that might have a potential for therapeutic applications.     

Yeast centromere DNA is in a unique and highly ordered structure in chromosomes and small circular minichromosomes

We have examined the chromatin structure of the centromere regions of chromosomes III and XI in yeast by using cloned functional centromere DNAs (CEN3 and CEN11) as labeled probes. When chromatin from isolated nuclei is digested with micrococcal nuclease and the resulting DNA fragments separated electrophoretically and blotted to nitrocellulose filters, the centromeric nucleosomal sub-units are resolved into significantly more distinct ladders than are those from the bulk of the chromatin. A discrete protected region of 220–250 bp of CEN sequence flanked by highly nuclease-sensitive sites was revealed by mapping the exact nuclease cleavage sites within the centromeric chromatin. On both sides of this protected region, highly phased and specific nuclease cutting sites exist at nucleosomal intervals (160 bp) for a total length of 12–15 nucleosomal subunits. The central protected region in the chromatin of both centromeres spans the 130 bp segment that exhibits the highest degree of sequence homology (71%) between functional CEN3 and CEN11 DNAs. This unique chromatin structure is maintained on CEN sequences introduced into yeast on autonomously replicating plasmids, but is not propagated through foreign DNA sequences flanking the inserted yeast DNA.     

Using microsatellites to understand the physical distribution of recombination on soybean chromosomes

Soybean is a major crop that is an important source of oil and proteins. A number of genetic linkage maps have been developed in soybean. Specifically, hundreds of simple sequence repeat (SSR) markers have been developed and mapped. Recent sequencing of the soybean genome resulted in the generation of vast amounts of genetic information. The objectives of this investigation were to use SSR markers in developing a connection between genetic and physical maps and to determine the physical distribution of recombination on soybean chromosomes. A total of 2,188 SSRs were used for sequence-based physical localization on soybean chromosomes. Linkage information was used from different maps to create an integrated genetic map. Comparison of the integrated genetic linkage maps and sequence based physical maps revealed that the distal 25% of each chromosome was the most marker-dense, containing an average of 47.4% of the SSR markers and 50.2% of the genes. The proximal 25% of each chromosome contained only 7.4% of the markers and 6.7% of the genes. At the whole genome level, the marker density and gene density showed a high correlation (R²) of 0.64 and 0.83, respectively with the physical distance from the centromere. Recombination followed a similar pattern with comparisons indicating that recombination is high in telomeric regions, though the correlation between crossover frequency and distance from the centromeres is low (R² = 0.21). Most of the centromeric regions were low in recombination. The crossover frequency for the entire soybean genome was 7.2%, with extremes much higher and lower than average. The number of recombination hotspots varied from 1 to 12 per chromosome. A high correlation of 0.83 between the distribution of SSR markers and genes suggested close association of SSRs with genes. The knowledge of distribution of recombination on chromosomes may be applied in characterizing and targeting genes.     

Centromeric DNA characterization in the model grass Brachypodium distachyon provides insights on the evolution of the genus

Brachypodium distachyon is a well‐established model monocot plant, and its small and compact genome has been used as an accurate reference for the much larger and often polyploid genomes of cereals such as Avena sativa (oats), Hordeum vulgare (barley) and Triticum aestivum (wheat). Centromeres are indispensable functional units of chromosomes and they play a core role in genome polyploidization events during evolution. As the Brachypodium genus contains about 20 species that differ significantly in terms of their basic chromosome numbers, genome size, ploidy levels and life strategies, studying their centromeres may provide important insight into the structure and evolution of the genome in this interesting and important genus. In this study, we isolated the centromeric DNA of the B. distachyon reference line Bd21 and characterized its composition via the chromatin immunoprecipitation of the nucleosomes that contain the centromere‐specific histone CENH3. We revealed that the centromeres of Bd21 have the features of typical multicellular eukaryotic centromeres. Strikingly, these centromeres contain relatively few centromeric satellite DNAs; in particular, the centromere of chromosome 5 (Bd5) consists of only ~40 kb. Moreover, the centromeric retrotransposons in B. distachyon (CRBds) are evolutionarily young. These transposable elements are located both within and adjacent to the CENH3 binding domains, and have similar compositions. Moreover, based on the presence of CRBds in the centromeres, the species in this study can be grouped into two distinct lineages. This may provide new evidence regarding the phylogenetic relationships within the Brachypodium genus.     

Strong nucleosomes of mouse genome including recovered centromeric sequences

Recently discovered strong nucleosomes (SNs) characterized by visibly periodical DNA sequences have been found to concentrate in centromeres of Arabidopsis thaliana and in transient meiotic centromeres of Caenorhabditis elegans. To find out whether such affiliation of SNs to centromeres is a more general phenomenon, we studied SNs of the Mus musculus. The publicly available genome sequences of mouse, as well as of practically all other eukaryotes do not include the centromere regions which are difficult to assemble because of a large amount of repeat sequences in the centromeres and pericentromeric regions. We recovered those missing sequences using the data from MNase-seq experiments in mouse embryonic stem cells, where the sequence of DNA inside nucleosomes, including missing regions, was determined by 100-bp paired-end sequencing. Those nucleosome sequences, which are not matching to the published genome sequence, would largely belong to the centromeres. By evaluating SN densities in centromeres and in non-centromeric regions, we conclude that mouse SNs concentrate in the centromeres of telocentric mouse chromosomes, with ~3.9 times excess compared to their density in the rest of the genome. The remaining non-centromeric SNs are harbored mainly by introns and intergenic regions, by retro-transposons, in particular. The centromeric involvement of the SNs opens new horizons for the chromosome and centromere structure studies.     

Diversity of centromeric repeats in two closely related wild rice species, Oryza officinalis and Oryza rhizomatis

Oryza officinalis (CC, 2n=24) and Oryza rhizomatis (CC, 2n=24) belong to the Oryza genus, which contains more than 20 identified wild rice species. Although much has been known about the molecular composition and organization of centromeres in Oryza sativa, relatively little is known of its wild relatives. In the present study, we isolated and characterized a 126-bp centromeric satellite (CentO-C) from three bacterial artificial chromosomes of O. officinalis. In addition to CentO-C, low abundance of CentO satellites is also present in O. officinalis. In order to determine the chromosomal locations and distributions of CentO-C (126-bp), CentO (155 bp) and TrsC (366 bp) satellite within O. officinalis, fluorescence in situ hybridization examination was done on pachytene or metaphase I chromosomes. We found that only ten centromeres (excluding centromere 7 and 2) contain CentO-C arrays in O. officinalis, while centromere 7 comprises CentO satellites, and centromere 2 is devoid of any detectable satellites. For TrsC satellites, it was detected at multiple subtelomeric regions in O. officinalis, however, in O. rhizomatis, TrsC sequences were detected both in the four centromeric regions (CEN 3, 4, 10, 11) and the multiple subtelomeric regions. Therefore, these data reveal the evolutionary diversification pattern of centromere DNA within/or between close related species, and could provide an insight into the dynamic evolutionary processes of rice centromere.     

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.     

Position-effect variegation and intercalary heterochromatin: a comparative study

The behaviour of IH (intercalary heterochromatin) regions of Drosophila melanogaster polytene chromosomes was compared with that of euchromatin condensed as a result of position-effect variegation. Normally replicating regions, when subject to such an effect, were found to become among the last regions in the genome to replicate. It is shown that the factors which enhance position effect (low temperature, the removal of the Y chromosome, genetic enhancers of position effect) increase the weak point frequency in the IH, i.e. enhance DNA underreplication in these regions. We suggest that the similarity in the properties of IH, CH (centromeric heterochromatin) and the dense blocks induced by position effect is due to strong genetic inactivation and supercondensation caused by specific proteins in early development. The primary DNA structure is not likely to play a key role in this process.     

Robertsonian Fusion and Centromere Repositioning Contributed to the Formation of Satellite-free Centromeres During the Evolution of Zebras

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated with highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and Equus asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras Equus burchelli (15 of 22) and Equus grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At 11 chromosomes, whose primary constrictions seemed to be associated with satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres; therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.