Arabidopsis thaliana centromere regions: genetic map positions and repetitive DNA structure

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we show that SRF could reconstruct known satellites in human and well-studied model organisms. We also find satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress in genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

Emerging mechanisms and roles of meiotic crossover repression at centromeres

Crossover events during recombination in meiosis are essential for generating genetic diversity as well as crucial to allow accurate chromosomal segregation between homologous chromosomes. Spatial control for the distribution of crossover events along the chromosomes is largely a tightly regulated process and involves many facets such as interference, repression as well as assurance, to make sure that not too many or too few crossovers are generated. Repression of crossover events at the centromeres is a highly conserved process across all species tested. Failure to inhibit such recombination events can result in chromosomal mis-segregation during meiosis resulting in aneuploid gametes that are responsible for infertility or developmental disorders such as Down's syndrome and other trisomies in humans. In the past few decades, studies to understand the molecular mechanisms behind this repression have shown the involvement of a multitude of factors ranging from the centromere-specific proteins such as the kinetochore to the flanking pericentric heterochromatin as well as DNA double-strand break repair pathways. In this chapter, we review the different mechanisms of pericentric repression mechanisms known till date as well as highlight the importance of understanding this regulation in the context of chromosomal segregation defects. We also discuss the clinical implications of dysregulation of this process, especially in human reproductive health and genetic diseases.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

The genetic and epigenetic landscape of the Arabidopsis centromeres

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support CENTROMERE SPECIFIC HISTONE H3 (CENH3) occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites that show the least amount of divergence and occur in higher-order repeats. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization. Centromeric crossover recombination is suppressed, yet low levels of meiotic DNA double-strand breaks occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving through cycles of satellite homogenization and retrotransposon-driven diversification.

Sequencing and Chromosome-Scale Assembly of Plant Genomes, Brassica rapa as a Use Case

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Not only sequencing techniques, but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome scale is now feasible at the laboratory scale. Current technologies, in particular long-range technologies, are numerous, and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range scaffolding techniques, optical mapping, Omni-C, and Pore-C sequencing libraries, commercialized by Bionano Genomics, Dovetail Genomics, and Oxford Nanopore Technologies, respectively, or a combination of the three, in order to evaluate the capability of each technology.

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.

Centromeres under Pressure: Evolutionary Innovation in Conflict with Conserved Function

Centromeres are essential genetic elements that enable spindle microtubule attachment for chromosome segregation during mitosis and meiosis. While this function is preserved across species, centromeres display an array of dynamic features, including: (1) rapidly evolving DNA; (2) wide evolutionary diversity in size, shape and organization; (3) evidence of mutational processes to generate homogenized repetitive arrays that characterize centromeres in several species; (4) tolerance to changes in position, as in the case of neocentromeres; and (5) intrinsic fragility derived by sequence composition and secondary DNA structures. Centromere drive underlies rapid centromere DNA evolution due to the "selfish" pursuit to bias meiotic transmission and promote the propagation of stronger centromeres. Yet, the origins of other dynamic features of centromeres remain unclear. Here, we review our current understanding of centromere evolution and plasticity. We also detail the mutagenic processes proposed to shape the divergent genetic nature of centromeres. Changes to centromeres are not simply evolutionary relics, but ongoing shifts that on one side promote centromere flexibility, but on the other can undermine centromere integrity and function with potential pathological implications such as genome instability.

CENP-A binding domains and recombination patterns in horse spermatocytes

Centromeres exert an inhibitory effect on meiotic recombination, but the possible contribution of satellite DNA to this "centromere effect" is under debate. In the horse, satellite DNA is present at all centromeres with the exception of the one from chromosome 11. This organization of centromeres allowed us to investigate the role of satellite DNA on recombination suppression in horse spermatocytes at the stage of pachytene. To this aim we analysed the distribution of the MLH1 protein, marker of recombination foci, relative to CENP-A, marker of centromeric function. We demonstrated that the satellite-less centromere of chromosome 11 causes crossover suppression, similarly to satellite-based centromeres. These results suggest that the centromere effect does not depend on satellite DNA. During this analysis, we observed a peculiar phenomenon: while, as expected, the centromere of the majority of meiotic bivalent chromosomes was labelled with a single immunofluorescence centromeric signal, double-spotted or extended signals were also detected. Their number varied from 0 to 7 in different cells. This observation can be explained by positional variation of the centromeric domain on the two homologs and/or misalignment of pericentromeric satellite DNA arrays during homolog pairing confirming the great plasticity of equine centromeres.

Genomic inversions and GOLGA core duplicons underlie disease instability at the 15q25 locus

Human chromosome 15q25 is involved in several disease-associated structural rearrangements, including microdeletions and chromosomal markers with inverted duplications. Using comparative fluorescence in situ hybridization, strand-sequencing, single-molecule, real-time sequencing and Bionano optical mapping analyses, we investigated the organization of the 15q25 region in human and nonhuman primates. We found that two independent inversions occurred in this region after the fission event that gave rise to phylogenetic chromosomes XIV and XV in humans and great apes. One of these inversions is still polymorphic in the human population today and may confer differential susceptibility to 15q25 microdeletions and inverted duplications. The inversion breakpoints map within segmental duplications containing core duplicons of the GOLGA gene family and correspond to the site of an ancestral centromere, which became inactivated about 25 million years ago. The inactivation of this centromere likely released segmental duplications from recombination repression typical of centromeric regions. We hypothesize that this increased the frequency of ectopic recombination creating a hotspot of hominid inversions where dispersed GOLGA core elements now predispose this region to recurrent genomic rearrangements associated with disease.

Evolving Centromeres and Kinetochores

The genetic material, contained on chromosomes, is often described as the "blueprint for life." During nuclear division, the chromosomes are pulled into each of the two daughter nuclei by the coordination of spindle microtubules, kinetochores, centromeres, and chromatin. These four functional units must link the chromosomes to the microtubules, signal to the cell when the attachment is made so that division can proceed, and withstand the force generated by pulling the chromosomes to either daughter cell. To perform each of these functions, kinetochores are large protein complexes, approximately 5MDa in size, and they contain at least 45 unique proteins. Many of the central components in the kinetochore are well conserved, yielding a common core of proteins forming consistent structures. However, many of the peripheral subcomplexes vary between different taxonomic groups, including changes in primary sequence and gain or loss of whole proteins. It is still unclear how significant these changes are, and answers to this question may provide insights into adaptation to specific lifestyles or progression of disease that involve chromosome instability.

Centromere location in Arabidopsis is unaltered by extreme divergence in CENH3 protein sequence

During cell division, spindle fibers attach to chromosomes at centromeres. The DNA sequence at regional centromeres is fast evolving with no conserved genetic signature for centromere identity. Instead CENH3, a centromere-specific histone H3 variant, is the epigenetic signature that specifies centromere location across both plant and animal kingdoms. Paradoxically, CENH3 is also adaptively evolving. An ongoing question is whether CENH3 evolution is driven by a functional relationship with the underlying DNA sequence. Here, we demonstrate that despite extensive protein sequence divergence, CENH3 histones from distant species assemble centromeres on the same underlying DNA sequence. We first characterized the organization and diversity of centromere repeats in wild-type Arabidopsis thaliana. We show that A. thaliana CENH3-containing nucleosomes exhibit a strong preference for a unique subset of centromeric repeats. These sequences are largely missing from the genome assemblies and represent the youngest and most homogeneous class of repeats. Next, we tested the evolutionary specificity of this interaction in a background in which the native A. thaliana CENH3 is replaced with CENH3s from distant species. Strikingly, we find that CENH3 from Lepidium oleraceum and Zea mays, although specifying epigenetically weaker centromeres that result in genome elimination upon outcrossing, show a binding pattern on A. thaliana centromere repeats that is indistinguishable from the native CENH3. Our results demonstrate positional stability of a highly diverged CENH3 on independently evolved repeats, suggesting that the sequence specificity of centromeres is determined by a mechanism independent of CENH3.

Database of Periodic DNA Regions in Major Genomes

Summary. We analyzed several prokaryotic and eukaryotic genomes looking for the periodicity sequences availability and employing a new mathematical method. The method envisaged using the random position weight matrices and dynamic programming. Insertions and deletions were allowed inside periodicities, thus adding a novelty to the results we obtained. A periodicity length, one of the key periodicity features, varied from 2 to 50 nt. Totally over 60,000 periodicity sequences were found in 15 genomes including some chromosomes of the H. sapiens (partial), C. elegans, D. melanogaster, and A. thaliana genomes.

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D-genome) of the allopolyploid (AD-genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin-immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon-related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A-genome and related diploid species (B-, F- and G-genomes), indicating that they colonized the centromeres of D-genome lineage after the divergence of the A- and D- ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A- and D-subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D- and AD-genome species, yet localized to just the NORs in A-, B-, F-, and G-genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.

A De Novo Genome Sequence Assembly of the Arabidopsis thaliana Accession Niederzenz-1 Displays Presence/Absence Variation and Strong Synteny

Arabidopsis thaliana is the most important model organism for fundamental plant biology. The genome diversity of different accessions of this species has been intensively studied, for example in the 1001 genome project which led to the identification of many small nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels). In addition, presence/absence variation (PAV), copy number variation (CNV) and mobile genetic elements contribute to genomic differences between A. thaliana accessions. To address larger genome rearrangements between the A. thaliana reference accession Columbia-0 (Col-0) and another accession of about average distance to Col-0, we created a de novo next generation sequencing (NGS)-based assembly from the accession Niederzenz-1 (Nd-1). The result was evaluated with respect to assembly strategy and synteny to Col-0. We provide a high quality genome sequence of the A. thaliana accession (Nd-1, LXSY01000000). The assembly displays an N50 of 0.590 Mbp and covers 99% of the Col-0 reference sequence. Scaffolds from the de novo assembly were positioned on the basis of sequence similarity to the reference. Errors in this automatic scaffold anchoring were manually corrected based on analyzing reciprocal best BLAST hits (RBHs) of genes. Comparison of the final Nd-1 assembly to the reference revealed duplications and deletions (PAV). We identified 826 insertions and 746 deletions in Nd-1. Randomly selected candidates of PAV were experimentally validated. Our Nd-1 de novo assembly allowed reliable identification of larger genic and intergenic variants, which was difficult or error-prone by short read mapping approaches alone. While overall sequence similarity as well as synteny is very high, we detected short and larger (affecting more than 100 bp) differences between Col-0 and Nd-1 based on bi-directional comparisons. The de novo assembly provided here and additional assemblies that will certainly be published in the future will allow to describe the pan-genome of A. thaliana.

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear-cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200,000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high-density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F(ST)), but not absolute (d(XY)), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to 'cold-spots' of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.

Centromeric chromatin and its dynamics in plants

Centromeres are chromatin structures that are required for proper separation of chromosomes during mitosis and meiosis. The centromere is composed of centromeric DNA, often enriched in satellite repeats, and kinetochore complex proteins. To date, over 100 kinetochore components have been identified in various eukaryotes. Kinetochore assembly begins with incorporation of centromeric histone H3 variant CENH3 into centromeric nucleosomes. Protein components of the kinetochore are either present at centromeres throughout the cell cycle or localize to centromeres transiently, prior to attachment of microtubules to each kinetochore in prometaphase of mitotic cells. This is the case for the spindle assembly checkpoint (SAC) proteins in animal cells. The SAC complex ensures equal separation of chromosomes between daughter nuclei by preventing anaphase onset before metaphase is complete, i.e. the sister kinetochores of all chromosomes are attached to spindle fibers from opposite poles. In this review, we focus on the organization of centromeric DNA and the kinetochore assembly in plants. We summarize recent advances regarding loading of CENH3 into the centromere, and the subcellular localization and protein-protein interactions of Arabidopsis thaliana proteins involved in kinetochore assembly and function. We describe the transcriptional activity of corresponding genes based on in silico analysis of their promoters and cell cycle-dependent expression. Additionally, barley homologs of all selected A. thaliana proteins have been identified in silico, and their sequences and domain structures are presented.

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

Background

SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.

Results

We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.

Conclusions

We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0678-1) contains supplementary material, which is available to authorized users.

Naturally occurring differences in CENH3 affect chromosome segregation in zygotic mitosis of hybrids

The point of attachment of spindle microtubules to metaphase chromosomes is known as the centromere. Plant and animal centromeres are epigenetically specified by a centromere-specific variant of Histone H3, CENH3 (a.k.a. CENP-A). Unlike canonical histones that are invariant, CENH3 proteins are accumulating substitutions at an accelerated rate. This diversification of CENH3 is a conundrum since its role as the key determinant of centromere identity remains a constant across species. Here, we ask whether naturally occurring divergence in CENH3 has functional consequences. We performed functional complementation assays on cenh3-1, a null mutation in Arabidopsis thaliana, using untagged CENH3s from increasingly distant relatives. Contrary to previous results using GFP-tagged CENH3, we find that the essential functions of CENH3 are conserved across a broad evolutionary landscape. CENH3 from a species as distant as the monocot Zea mays can functionally replace A. thaliana CENH3. Plants expressing variant CENH3s that are fertile when selfed show dramatic segregation errors when crossed to a wild-type individual. The progeny of this cross include hybrid diploids, aneuploids with novel genetic rearrangements and haploids that inherit only the genome of the wild-type parent. Importantly, it is always chromosomes from the plant expressing the divergent CENH3 that missegregate. Using chimeras, we show that it is divergence in the fast-evolving N-terminal tail of CENH3 that is causing segregation errors and genome elimination. Furthermore, we analyzed N-terminal tail sequences from plant CENH3s and discovered a modular pattern of sequence conservation. From this we hypothesize that while the essential functions of CENH3 are largely conserved, the N-terminal tail is evolving to adapt to lineage-specific centromeric constraints. Our results demonstrate that this lineage-specific evolution of CENH3 causes inviability and sterility of progeny in crosses, at the same time producing karyotypic variation. Thus, CENH3 evolution can contribute to postzygotic reproductive barriers.

Identification of centromeric regions on the linkage map of cotton using centromere-related repeats

Centromere usually contains high-copy-number retrotransposons and satellite repeats, which are difficult to map, clone and sequence. Currently, very little is known about the centromere in cotton. Here, we sequenced a bacterial artificial chromosome (BAC) mapping to the centromeric region and predicted four long-terminal-repeat (LTR) retrotransposons. They were located in the heterochromatic centromeric regions of all 52 pachytene chromosomes in Gossypium hirsutum. Fiber-FISH mapping revealed that these retrotransposons span an area of at least 1.8Mb in the centromeric region. Comparative analysis showed that these retrotransposons generated similar, strong fluorescent signals in the D progenitor Gossypium raimondii but not in the A progenitor Gossypium herbaceum, suggesting that the centromere sequence of tetraploid cotton might be derived from the D progenitor. Centromeric regions were anchored on 13 chromosomes of D-genome sequence. Characterization of these centromere-related repeats and regions will enhance cotton centromere mapping, sequencing and evolutionary studies.

Identification and diversity of functional centromere satellites in the wild rice species Oryza brachyantha

The centromere is a key chromosomal component for sister chromatid cohesion and is the site for kinetochore assembly and spindle fiber attachment, allowing each sister chromatid to faithfully segregate to each daughter cell during cell division. It is not clear what types of sequences act as functional centromeres and how centromere sequences are organized in Oryza brachyantha, an FF genome species. In this study, we found that the three classes of centromere-specific CentO-F satellites (CentO-F1, CentO-F2, and CentOF3) in O. brachyantha share no homology with the CentO satellites in Oryza sativa. The three classes of CentO-F satellites are all located within the chromosomal regions to which the spindle fibers attach and are characterized by megabase tandem arrays that are flanked by centromere-specific retrotransposons, CRR-F, in the O. brachyantha centromeres. Although these CentO-F satellites are quantitatively variable among 12 O. brachyantha centromeres, immunostaining with an antibody specific to CENH3 indicates that they are colocated with CENH3 in functional centromere regions. Our results demonstrate that the three classes of CentO-F satellites may be the major components of functional centromeres in O. brachyantha.

Sequence organization and evolutionary dynamics of Brachypodium-specific centromere retrotransposons

Brachypodium distachyon is a wild annual grass belonging to the Pooideae, more closely related to wheat, barley, and forage grasses than rice and maize. As an experimental model, the completed genome sequence of B. distachyon provides a unique opportunity to study centromere evolution during the speciation of grasses. Centromeric satellite sequences have been identified in B. distachyon, but little is known about centromeric retrotransposons in this species. In the present study, bacterial artificial chromosome (BAC)-fluorescence in situ hybridization was conducted in maize, rice, barley, wheat, and rye using B. distachyon (Bd) centromere-specific BAC clones. Eight Bd centromeric BAC clones gave no detectable fluorescence in situ hybridization (FISH) signals on the chromosomes of rice and maize, and three of them also did not yield any FISH signals in barley, wheat, and rye. In addition, four of five Triticeae centromeric BAC clones did not hybridize to the B. distachyon centromeres, implying certain unique features of Brachypodium centromeres. Analysis of Brachypodium centromeric BAC sequences identified a long terminal repeat (LTR)-centromere retrotransposon of B. distachyon (CRBd1). This element was found in high copy number accounting for 1.6 % of the B. distachyon genome, and is enriched in Brachypodium centromeric regions. CRBd1 accumulated in active centromeres, but was lost from inactive ones. The LTR of CRBd1 appears to be specific to B. distachyon centromeres. These results reveal different evolutionary events of this retrotransposon family across grass species.

Fine-scale heterogeneity in crossover rate in the garnet-scalloped region of the Drosophila melanogaster X chromosome

Homologous recombination affects myriad aspects of genome evolution, from standing levels of nucleotide diversity to the efficacy of natural selection. Rates of crossing over show marked variability at all scales surveyed, including species-, population-, and individual-level differences. Even within genomes, crossovers are nonrandomly distributed in a wide diversity of taxa. Although intra- and intergenomic heterogeneities in crossover distribution have been documented in Drosophila, the scale and degree of crossover rate heterogeneity remain unclear. In addition, the genetic features mediating this heterogeneity are unknown. Here we quantify fine-scale heterogeneity in crossover distribution in a 2.1-Mb region of the Drosophila melanogaster X chromosome by localizing crossover breakpoints in 2500 individuals, each containing a single crossover in this specific X chromosome region. We show 90-fold variation in rates of crossing over at a 5-kb scale, place this variation in the context of several aspects of genome evolution, and identify several genetic features associated with crossover rates. Our results shed new light on the scale and magnitude of crossover rate heterogeneity in D. melanogaster and highlight potential features mediating this heterogeneity.

An Arabidopsis tissue-specific RNAi method for studying genes essential to mitosis

A large fraction of the genes in plants can be considered essential in the sense that when absent the plant fails to develop past the first few cell divisions. The fact that angiosperms pass through a haploid gametophyte stage can make it challenging to propagate such mutants even in the heterozygous condition. Here we describe a tissue-specific RNAi method that allows us to visualize cell division phenotypes in petals, which are large dispensable organs. Portions of the APETALA (AP3) and PISTILLATA (PI) promoters confer early petal-specific expression. We show that when either promoter is used to drive the expression of a beta-glucuronidase (GUS) RNAi transgene in plants uniformly expressing GUS, GUS expression is knocked down specifically in petals. We further tested the system by targeting the essential kinetochore protein CENPC and two different components of the Spindle Assembly Checkpoint (MAD2 and BUBR1). Plant lines expressing petal-specific RNAi hairpins targeting these genes exhibited an array of petal phenotypes. Cytological analyses of the affected flower buds confirmed that CENPC knockdown causes cell cycle arrest but provided no evidence that either MAD2 or BUBR1 are required for mitosis (although both genes are required for petal growth by this assay). A key benefit of the petal-specific RNAi method is that the phenotypes are not expressed in the lineages leading to germ cells, and the phenotypes are faithfully transmitted for at least four generations despite their pronounced effects on growth.

Variation of presence/absence genes among Arabidopsis populations

BACKGROUND

Gene presence/absence (P/A) polymorphisms are commonly observed in plants and are important in individual adaptation and species differentiation. Detecting their abundance, distribution and variation among individuals would help to understand the role played by these polymorphisms in a given species. The recently sequenced 80 Arabidopsis genomes provide an opportunity to address these questions.

RESULTS

By systematically investigating these accessions, we identified 2,407 P/A genes (or 8.9%) absent in one or more genomes, averaging 444 absent genes per accession. 50.6% of P/A genes belonged to multi-copy gene families, or 31.0% to clustered genes. However, the highest proportion of P/A genes, outnumbered in singleton genes, was observed in the regions near centromeres. In addition, a significant correlation was observed between the P/A gene frequency among the 80 accessions and the diversity level at P/A loci. Furthermore, the proportion of P/A genes was different among functional gene categories. Finally, a P/A gene tree showed a diversified population structure in the worldwide Arabidopsis accessions.

CONCLUSIONS

An estimate of P/A genes and their frequency distribution in the worldwide Arabidopsis accessions was obtained. Our results suggest that there are diverse mechanisms to generate or maintain P/A genes, by which individuals and functionally different genes can selectively maintain P/A polymorphisms for a specific adaptation.

The cotton centromere contains a Ty3-gypsy-like LTR retroelement

The centromere is a repeat-rich structure essential for chromosome segregation; with the long-term aim of understanding centromere structure and function, we set out to identify cotton centromere sequences. To isolate centromere-associated sequences from cotton, (Gossypium hirsutum) we surveyed tandem and dispersed repetitive DNA in the genus. Centromere-associated elements in other plants include tandem repeats and, in some cases, centromere-specific retroelements. Examination of cotton genomic survey sequences for tandem repeats yielded sequences that did not localize to the centromere. However, among the repetitive sequences we also identified a gypsy-like LTR retrotransposon (Centromere Retroelement Gossypium, CRG) that localizes to the centromere region of all chromosomes in domestic upland cotton, Gossypium hirsutum, the major commercially grown cotton. The location of the functional centromere was confirmed by immunostaining with antiserum to the centromere-specific histone CENH3, which co-localizes with CRG hybridization on metaphase mitotic chromosomes. G. hirsutum is an allotetraploid composed of A and D genomes and CRG is also present in the centromere regions of other AD cotton species. Furthermore, FISH and genomic dot blot hybridization revealed that CRG is found in D-genome diploid cotton species, but not in A-genome diploid species, indicating that this retroelement may have invaded the A-genome centromeres during allopolyploid formation and amplified during evolutionary history. CRG is also found in other diploid Gossypium species, including B and E2 genome species, but not in the C, E1, F, and G genome species tested. Isolation of this centromere-specific retrotransposon from Gossypium provides a probe for further understanding of centromere structure, and a tool for future engineering of centromere mini-chromosomes in this important crop species.

Differing requirements for RAD51 and DMC1 in meiotic pairing of centromeres and chromosome arms in Arabidopsis thaliana

During meiosis homologous chromosomes pair, recombine, and synapse, thus ensuring accurate chromosome segregation and the halving of ploidy necessary for gametogenesis. The processes permitting a chromosome to pair only with its homologue are not fully understood, but successful pairing of homologous chromosomes is tightly linked to recombination. In Arabidopsis thaliana, meiotic prophase of rad51, xrcc3, and rad51C mutants appears normal up to the zygotene/pachytene stage, after which the genome fragments, leading to sterility. To better understand the relationship between recombination and chromosome pairing, we have analysed meiotic chromosome pairing in these and in dmc1 mutant lines. Our data show a differing requirement for these proteins in pairing of centromeric regions and chromosome arms. No homologous pairing of mid-arm or distal regions was observed in rad51, xrcc3, and rad51C mutants. However, homologous centromeres do pair in these mutants and we show that this does depend upon recombination, principally on DMC1. This centromere pairing extends well beyond the heterochromatic centromere region and, surprisingly, does not require XRCC3 and RAD51C. In addition to clarifying and bringing the roles of centromeres in meiotic synapsis to the fore, this analysis thus separates the roles in meiotic synapsis of DMC1 and RAD51 and the meiotic RAD51 paralogs, XRCC3 and RAD51C, with respect to different chromosome domains.

Meiosis is a specialised cell division that acts to halve the chromosome complement, or ploidy, in the production of gametes for sexual reproduction in eukaryotes. To ensure that each gamete has a full complement of the genetic material, homologous chromosomes must pair and then separate in a coordinated manner during meiosis, and this is mediated by recombination in the majority of studied eukaryotes. To better understand the relationship between recombination and meiotic homologue pairing, we have analysed meiotic chromosome pairing in plant mutants lacking key recombination proteins. This work provides new insights into the homologous chromosome pairing mechanisms occurring in meiotic prophase of Arabidopsis thaliana: heterochromatic centromeres and 5S rDNA regions pair early, and their pairing has different requirements for recombination proteins than does that of the chromosome arms. These data raise a number of questions concerning the specificities and roles of recombination at different chromosome and/or chromatin regions in the synapsis of homologous chromosomes at meiosis.

Dysfunction of chromatin assembly factor 1 induces shortening of telomeres and loss of 45S rDNA in Arabidopsis thaliana

Chromatin Assembly Factor 1 (CAF1) is a three-subunit H3/H4 histone chaperone responsible for replication-dependent nucleosome assembly. It is composed of CAC 1-3 in yeast; p155, p60, and p48 in humans; and FASCIATA1 (FAS1), FAS2, and MULTICOPY SUPPRESSOR OF IRA1 in Arabidopsis thaliana. We report that disruption of CAF1 function by fas mutations in Arabidopsis results in telomere shortening and loss of 45S rDNA, while other repetitive sequences (5S rDNA, centromeric 180-bp repeat, CACTA, and Athila) are unaffected. Substantial telomere shortening occurs immediately after the loss of functional CAF1 and slows down at telomeres shortened to median lengths around 1 to 1.5 kb. The 45S rDNA loss is progressive, leaving 10 to 15% of the original number of repeats in the 5th generation of mutants affecting CAF1, but the level of the 45S rRNA transcripts is not altered in these mutants. Increasing severity of the fas phenotype is accompanied by accumulation of anaphase bridges, reduced viability, and plant sterility. Our results show that appropriate replication-dependent chromatin assembly is specifically required for stable maintenance of telomeres and 45S rDNA.

Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4-1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

A molecular-cytogenetic method for locating genes to pericentromeric regions facilitates a genomewide comparison of synteny between the centromeric regions of wheat and rice

Centromeres, because of their repeat structure and lack of sequence conservation, are difficult to assemble and compare across organisms. It was recently discovered that rice centromeres often contain genes. This suggested a method for studying centromere homologies between wheat and rice chromosomes by mapping rice centromeric genes onto wheat aneuploid stocks. Three of the seven cDNA clones of centromeric genes from rice centromere 8 (Cen8), 6729.t09, 6729.t10, and 6730.t11 which lie in the Cen8 kinetochore region, and three wheat ESTs, BJ301191, BJ305475, and BJ280500, with similarity to sequences of rice centromeric genes, were mapped to the centromeric regions of the wheat group-7 (W7) chromosomes. A possible pericentric inversion in chromosome 7D was detected. Genomewide comparison of wheat ESTs that mapped to centromeric regions against rice genome sequences revealed high conservation and a one-to-one correspondence of centromeric regions between wheat and rice chromosome pairs W1-R5, W2-R7, W3-R1, W5-R12, W6-R2, and W7-R8. The W4 centromere may share homology with R3 only or with R3 + R11. Wheat ESTs that mapped to the pericentromeric region of the group-5 long arm anchored to the rice BACs located in the recently duplicated region at the distal ends of the short arms of rice chromosomes 11 and 12. A pericentric inversion specific to the rice lineage was detected. The depicted framework provides a working model for further studies on the structure and evolution of cereal chromosome centromeres.

Divergence in centromere structure distinguishes related genomes in Coix lacryma-jobi and its wild relative

Knowledge about the composition and structure of centromeres is critical for understanding how centromeres perform their functional roles. Here, we report the sequences of one centromere-associated bacterial artificial chromosome clone from a Coix lacryma-jobi library. Two Ty3/gypsy-class retrotransposons, centromeric retrotransposon of C. lacryma-jobi (CRC) and peri-centromeric retrotransposon of C. lacryma-jobi, and a (peri)centromere-specific tandem repeat with a unit length of 153 bp were identified. The CRC is highly homologous to centromere-specific retrotransposons reported in grass species. An 80-bp DNA region in the 153-bp satellite repeat was found to be conserved to centromeric satellite repeats from maize, rice, and pearl millet. Fluorescence in situ hybridization showed that the three repetitive sequences were located in (peri-)centromeric regions of both C. lacryma-jobi and Coix aquatica. However, the 153-bp satellite repeat was only detected on 20 out of the 30 chromosomes in C. aquatica. Immunostaining with an antibody against rice CENH3 indicates that the 153-bp satellite repeat and CRC might be both the major components for functional centromeres, but not all the 153-bp satellite repeats or CRC sequences are associated with CENH3. The evolution of centromeric repeats of C. lacryma-jobi during the polyploidization was discussed.

RNA polymerase V functions in Arabidopsis interphase heterochromatin organization independently of the 24-nt siRNA-directed DNA methylation pathway

In Arabidopsis, pericentromeric repeats, retroelements, and silenced rRNA genes are assembled into heterochromatin within nuclear structures known as chromocenters. The mechanisms governing higher-order heterochromatin organization are poorly understood but 24-nt small interfering RNAs (siRNAs) are known to play key roles in heterochromatin formation. Nuclear RNA polymerase IV (Pol IV), RNA-DEPENDENT RNA POLYMERASE 2 (RDR2), and DICER-LIKE 3 (DCL3) are required for biogenesis of 24-nt siRNAs that associate with ARGONAUTE 4 (AGO4). Nuclear RNA polymerase V (Pol V) collaborates with DRD1 (DEFICIENT IN RNA-DEPENDENT DNA METHYLATION 1) to generate transcripts at heterochromatic loci that are hypothesized to bind to siRNA-AGO4 complexes and subsequently recruit the de-novo DNA methylation and/or histone modifying machinery. Here, we report that decondensation of the major pericentromeric repeats and depletion of the heterochromatic mark histone H3 lysine 9 dimethylation at chromocenters occurs specifically in pol V and drd1 mutants. Disruption of pericentromeric repeats condensation is coincident with transcriptional reactivation of specific classes of pericentromeric 180-bp repeats. We further demonstrate that Pol V functions independently of Pol IV, RDR2, and DCL3-mediated siRNA production to affect interphase heterochromatin organization, possibly by involving RNAs that recruit structural or chromatin-modifying proteins.

Coding DNA repeated throughout intergenic regions of the Arabidopsis thaliana genome: evolutionary footprints of RNA silencing

Pyknons are non-random sequence patterns significantly repeated throughout non-coding genomic DNA that also appear at least once among coding genes. They are interesting because they portend an unforeseen connection between coding and non-coding DNA. Pyknons have only been discovered in the human genome, so it is unknown whether pyknons have wider biological relevance or are simply a phenomenon of the human genome. To address this, DNA sequence patterns from the Arabidopsis thaliana genome were detected using a probability-based method. 24 654 statistically significant sequence patterns, 16 to 24 nucleotides long, repeating 10 or more times in non-coding DNA also appeared in 46% of A. thaliana protein-coding genes. A. thaliana pyknons exhibit features similar to human pyknons, including being distinct sequence patterns, having multiple instances in genes and having remarkable similarity to small RNA sequences with roles in gene silencing. Chromosomal position mapping revealed that genomic pyknon density has concordance with siRNA and transposable element positioning density. Because the A. thaliana and human genomes have approximately the same number of genes but drastically different amounts of non-coding DNA, these data reveal that pyknons represent a biologically important link between coding and non-coding DNA. Because of the association of pyknons with siRNAs and localization to silenced regions of heterochromatin, we postulate that RNA-mediated gene silencing leads to the accumulation of gene sequences in non-coding DNA regions.

Repetitive sequence-derived markers tag centromeres and telomeres and provide insights into chromosome evolution in Brassica napus

Centromeres and telomeres are obvious markers on chromosomes but their location on genetic maps is difficult to determine, which hampers many basic and applied research programmes. In this study, we used the characteristic distribution of five Brassica repeated sequences to generate physically anchored molecular markers tentatively tagging Brassica centromeres (84 markers) and telomeres (31 markers). These markers were mapped to the existing oilseed rape genetic map. Clusters of centromere-related loci were observed on 14 linkage groups; in addition to previous reports, we could thus provide information about the most likely position of centromeres on 17 of the 19 B. napus linkage groups. The location of centromeres on linkage groups usually matches their position on chromosomes and coincides with sites of evolutionary breakage between chromosomes. Most telomere sequence-derived markers mapped interstitially or in the proximity of centromeres; this result echoes previous reports on many eukaryote genomes and may reflect different forms of chromosome evolution. Seven telomere sequence-derived markers were located at the outermost positions of seven linkage groups and therefore probably tagged telomeres.

Progress in understanding and sequencing the genome of Brassica rapa

Brassica rapa, which is closely related to Arabidopsis thaliana, is an important crop and a model plant for studying genome evolution via polyploidization. We report the current understanding of the genome structure of B. rapa and efforts for the whole-genome sequencing of the species. The tribe Brassicaceae, which comprises ca. 240 species, descended from a common hexaploid ancestor with a basic genome similar to that of Arabidopsis. Chromosome rearrangements, including fusions and/or fissions, resulted in the present-day "diploid" Brassica species with variation in chromosome number and phenotype. Triplicated genomic segments of B. rapa are collinear to those of A. thaliana with InDels. The genome triplication has led to an approximately 1.7-fold increase in the B. rapa gene number compared to that of A. thaliana. Repetitive DNA of B. rapa has also been extensively amplified and has diverged from that of A. thaliana. For its whole-genome sequencing, the Brassica rapa Genome Sequencing Project (BrGSP) consortium has developed suitable genomic resources and constructed genetic and physical maps. Ten chromosomes of B. rapa are being allocated to BrGSP consortium participants, and each chromosome will be sequenced by a BAC-by-BAC approach. Genome sequencing of B. rapa will offer a new perspective for plant biology and evolution in the context of polyploidization.

Epigenetic modification of centromeric chromatin: hypomethylation of DNA sequences in the CENH3-associated chromatin in Arabidopsis thaliana and maize

The centromere in eukaryotes is defined by the presence of a special histone H3 variant, CENH3. Centromeric chromatin consists of blocks of CENH3-containing nucleosomes interspersed with blocks of canonical H3-containing nucleosomes. However, it is not known how CENH3 is precisely deposited in the centromeres. It has been suggested that epigenetic modifications of the centromeric chromatin may play a role in centromere identity. The centromeres of Arabidopsis thaliana are composed of megabase-sized arrays of a 178-bp satellite repeat. Here, we report that the 178-bp repeats associated with the CENH3-containing chromatin (CEN chromatin) are hypomethylated compared with the same repeats located in the flanking pericentromeric regions. A similar hypomethylation of DNA in CEN chromatin was also revealed in maize (Zea mays). Hypomethylation of the DNA in CEN chromatin is correlated with a significantly reduced level of H3K9me2 in Arabidopsis. We demonstrate that the 178-bp repeats from CEN chromatin display a distinct distribution pattern of the CG and CNG sites, which may provide a foundation for the differential methylation of these repeats. Our results suggest that DNA methylation plays an important role in epigenetic demarcation of the CEN chromatin.

Large-scale polymorphism of heterochromatic repeats in the DNA of Arabidopsis thaliana

BACKGROUND

The composition of the individual eukaryote's genome and its variation within a species remain poorly defined. Even for a sequenced genome such as that of the model plant Arabidopsis thaliana accession Col-0, the large arrays of heterochromatic repeats are incompletely sequenced, with gaps of uncertain size persisting in them.

RESULTS

Using geographically separate populations of A. thaliana, we assayed variation in the heterochromatic repeat arrays using two independent methods and identified significant polymorphism among them, with variation by as much as a factor of two in the centromeric 180 bp repeat, in the 45S rDNA arrays and in the Athila retroelements. In the accession with highest genome size as measured by flow cytometry, Loh-0, we found more than a two-fold increase in 5S RNA gene copies relative to Col-0; results from fluorescence in situ hybridization with 5S probes were consistent with the existence of size polymorphism between Loh-0 and Col-0 at the 5S loci. Comparative genomic hybridization results of Loh-0 and Col-0 did not support contiguous variation in copy number of protein-coding genes on the scale needed to explain their observed genome size difference. We developed a computational data model to test whether the variation we measured in the repeat fractions could account for the different genome sizes determined with flow cytometry, and found that this proposed relationship could account for about 50% of the variance in genome size among the accessions.

CONCLUSION

Our analyses are consistent with substantial repeat number polymorphism for 5S and 45S ribosomal genes among accession of A. thaliana. Differences are also suggested for centromeric and pericentromeric repeats. Our analysis also points to the difficulties in measuring the repeated fraction of the genome and suggests that independent validation of genome size should be sought in addition to flow cytometric measurements.

Large-scale dissociation and sequential reassembly of pericentric heterochromatin in dedifferentiated Arabidopsis cells

Chromocenters in Arabidopsis thaliana are discrete nuclear domains of mainly pericentric heterochromatin. They are characterized by the presence of repetitive sequences, methylated DNA and dimethylated histone H3K9. Here we show that dedifferentiation of specialized mesophyll cells into undifferentiated protoplasts is accompanied by the disruption of chromocenter structures. The dramatic reduction of heterochromatin involves the decondensation of all major repeat regions, also including the centromeric 180 bp tandem repeats. Only the 45S rDNA repeat remained in a partly compact state in most cells. Remarkably, the epigenetic indicators for heterochromatin, DNA methylation and H3K9 dimethylation, did not change upon decondensation. Furthermore, the decondensation of pericentric heterochromatin did not result in transcriptional reactivation of silent genomic elements. The decondensation process was reversible upon prolonged culturing. Strikingly, recondensation of heterochromatin into chromocenters is a stepwise process. Compaction of the tandemly arranged 45S rDNA regions occurs first, followed by the centromeric 180 bp and the 5S rDNA repeats and finally the dispersed repeats, including transposons. The sequence of reassembly seems to be correlated to the size of the repeat domains. Our results indicate that different types of pericentromeric repeats form different types of heterochromatin, which subsequently merge to form a chromocenter.

Plant centromere organization: a dynamic structure with conserved functions

Although the structural features of centromeres from most multicellular eukaryotes remain to be characterized, recent analyses of the complete sequences of two centromeric regions of rice, together with data from Arabidopsis thaliana and maize, have illuminated the considerable size variation and sequence divergence of plant centromeres. Despite the severe suppression of meiotic chromosomal exchange in centromeric and pericentromeric regions of rice, the centromere core shows high rates of unequal homologous recombination in the absence of chromosomal exchange, resulting in frequent and extensive DNA rearrangement. Not only is the sequence of centromeric tandem and non-tandem repeats highly variable but also the copy number, spacing, order and orientation, providing ample natural variation as the basis for selection of superior centromere performance. This review article focuses on the structural and evolutionary dynamics of plant centromere organization and the potential molecular mechanisms responsible for the rapid changes of centromeric components.

Characterization of the centromere and peri-centromere retrotransposons in Brassica rapa and their distribution in related Brassica species

We report the identification and characterization of the major repeats in the centromeric and peri-centromeric heterochromatin of Brassica rapa. The analysis involved the characterization of 88 629 bacterial artificial chromosomes (BAC) end sequences and the complete sequences of two BAC clones. We identified centromere-specific retrotransposons of Brassica (CRB) and various peri-centromere-specific retrotransposons (PCRBr). Three copies of the CRB were identified in one BAC clone as nested insertions within a tandem array of 24 copies of a 176 bp centromeric repeat, CentBr. A complex mosaic structure consisting of nine PCRBr elements and large blocks of 238 bp degenerate tandem repeats (TR238) were found in or near a derivative of 5S-25S rDNA sequences. The chromosomal positions of selected repeats were determined using in situ hybridization. These revealed that CRB is a major component of all centromeres in three diploid Brassica species and their allotetraploid relatives. However, CentBr was not detected in the most distantly related of the diploid species analyzed, B. nigra. PCRBr and TR238 were found to be major components in the peri-centromeric heterochromatin blocks of four chromosomes of B. rapa. These repetitive elements were not identified in B. oleracea or B. nigra, indicating that they are A-genome-specific. GenBank accession numbers: KBrH001P13 (AC 166739); KBrH015B20 (AC 166740); end sequences of KBrH BAC library (CW 978640 - CW 988843); end sequences of KBrS BAC library (DU 826965 - DU 835595); end sequences of KBrB BAC library (DX 010661 - DX 083363).

Ecotype-specific and chromosome-specific expansion of variant centromeric satellites in Arabidopsis thaliana

Despite the conserved roles and conserved protein machineries of centromeres, their nucleotide sequences can be highly diverse even among related species. The diversity reflects rapid evolution, but the underlying mechanism is largely unknown. One approach to monitor rapid evolution is examination of intra-specific variation. Here we report variant centromeric satellites of Arabidopsis thaliana found through survey of 103 natural accessions (ecotypes). Among them, a cluster of variant centromeric satellites was detected in one ecotype, Cape Verde Islands (Cvi). Recombinant inbred mapping revealed that the variant satellites are distributed in centromeric region of the chromosome 5 (CEN5) of this ecotype. This apparently recent variant accumulation is associated with large deletion of a pericentromeric region and the expansion of satellite region. The variant satellite was bound to HTR12 (centromeric variant histone H3), although expansion of the satellite was not associated with comparable increase in the HTR12 binding. The results suggest that variant satellites with centromere function can rapidly accumulate in one centromere, supporting the model that the satellite repeats in the array are homogenized by occasional unequal crossing-over, which has a potential to generate an expansion of local sequence variants within a centromere cluster.

Genomic and genetic characterization of rice Cen3 reveals extensive transcription and evolutionary implications of a complex centromere

The centromere is the chromosomal site for assembly of the kinetochore where spindle fibers attach during cell division. In most multicellular eukaryotes, centromeres are composed of long tracts of satellite repeats that are recalcitrant to sequencing and fine-scale genetic mapping. Here, we report the genomic and genetic characterization of the complete centromere of rice (Oryza sativa) chromosome 3. Using a DNA fiber-fluorescence in situ hybridization approach, we demonstrated that the centromere of chromosome 3 (Cen3) contains approximately 441 kb of the centromeric satellite repeat CentO. Cen3 includes an approximately 1,881-kb domain associated with the centromeric histone CENH3. This CENH3-associated chromatin domain is embedded within a 3,113-kb region that lacks genetic recombination. Extensive transcription was detected within the CENH3 binding domain based on comprehensive annotation of protein-coding genes coupled with empirical measurements of mRNA levels using RT-PCR and massively parallel signature sequencing. Genes <10 kb from the CentO satellite array were expressed in several rice tissues and displayed histone modification patterns consistent with euchromatin, suggesting that rice centromeric chromatin accommodates normal gene expression. These results support the hypothesis that centromeres can evolve from gene-containing genomic regions.