Genetic definition and sequence analysis of Arabidopsis centromeres

Structure and evolution of metapolycentromeres

Metapolycentromeres consist of multiple sequential domains of centromeric chromatin associated with a centromere-specific variant of histone H3 (CENP-A), functioning collectively as a single centromere. To date, they have been revealed in nine flowering plant, five insect and six vertebrate species. In this paper, we focus on their structure and possible mechanisms of emergence and evolution. The metapolycentromeres may vary in the number of centromeric domains and in their genetic content and epigenetic modifications. However, these variations do not seem to affect their function. The emergence of metapolycentromeres has been attributed to multiple Robertsonian translocations and segmental duplications. Conditions of genomic instability, such as interspecific hybridization and malignant neoplasms, are suggested as triggers for the de novo emergence of metapolycentromeres. Addressing the "centromere paradox" - the rapid evolution of centromeric DNA and proteins despite their conserved cellular function - we explore the centromere drive hypothesis as a plausible explanation for the dynamic evolution of centromeres in general, and in particular the emergence of metapolycentromeres and holocentromeres. Apparently, metapolycentromeres are more common across different species than it was believed until recently. Indeed, a systematic review of the available cytogenetic publications allowed us to identify 27 candidate species with metapolycentromeres. Тhe list of the already established and newly revealed candidate species thus spans 27 species of flowering plants and eight species of gymnosperm plants, five species of insects, and seven species of vertebrates. This indicates an erratic phylogenetic distribution of the species with metapolycentromeres and may suggest an independent emergence of the metapolycentromeres in the course of evolution. However, the current catalog of species with identified and likely metapolycentromeres remains too short to draw reliable conclusions about their evolution, particularly in the absence of knowledge about related species without metapolycentromeres for comparative analysis. More studies are necessary to shed light on the mechanisms of metapolycentromere formation and evolution.

Retrotransposon addiction promotes centromere function via epigenetically activated small RNAs

Retrotransposons have invaded eukaryotic centromeres in cycles of repeat expansion and purging, but the function of centromeric retrotransposons has remained unclear. In Arabidopsis, centromeric ATHILA retrotransposons give rise to epigenetically activated short interfering RNAs in mutants in DECREASE IN DNA METHYLATION1 (DDM1). Here we show that mutants that lose both DDM1 and RNA-dependent RNA polymerase have pleiotropic developmental defects and mis-segregate chromosome 5 during mitosis. Fertility and segregation defects are epigenetically inherited with centromere 5, and can be rescued by directing artificial small RNAs to ATHILA5 retrotransposons that interrupt tandem satellite repeats. Epigenetically activated short interfering RNAs promote pericentromeric condensation, chromosome cohesion and chromosome segregation in mitosis. We propose that insertion of ATHILA silences centromeric transcription, while simultaneously making centromere function dependent on retrotransposon small RNAs in the absence of DDM1. Parallels are made with the fission yeast Schizosaccharomyces pombe, where chromosome cohesion depends on RNA interference, and with humans, where chromosome segregation depends on both RNA interference and HELLSDDM1.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

Centromere-proximal suppression of meiotic crossovers in Drosophila is robust to changes in centromere number, repetitive DNA content, and centromere-clustering

Accurate segregation of homologous chromosomes during meiosis depends on both the presence and the regulated placement of crossovers (COs). The centromere effect, or CO exclusion in pericentromeric regions of the chromosome, is a meiotic CO patterning phenomenon that helps prevent nondisjunction, thereby protecting against chromosomal disorders and other meiotic defects. Despite being identified nearly a century ago, the mechanisms behind this fundamental cellular process remain unknown, with most studies of the Drosophila centromere effect focusing on local influences of the centromere and pericentric heterochromatin. In this study, we sought to investigate whether dosage changes in centromere number and repetitive DNA content affect the strength of the centromere effect, using phenotypic recombination mapping. Additionally, we studied the effects of repetitive DNA function on centromere effect strength using satellite DNA-binding protein mutants displaying defective centromere-clustering in meiotic nuclei. Despite what previous studies suggest, our results show that the Drosophila centromere effect is robust to changes in centromere number, repetitive DNA content, as well as repetitive DNA function. Our study suggests that the centromere effect is unlikely to be spatially controlled, providing novel insight into the mechanisms behind the Drosophila centromere effect.

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

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

MSH2 stimulates interfering and inhibits non-interfering crossovers in response to genetic polymorphism

Meiotic crossovers can be formed through the interfering pathway, in which one crossover prevents another from forming nearby, or by an independent non-interfering pathway. In Arabidopsis, local sequence polymorphism between homologs can stimulate interfering crossovers in a MSH2-dependent manner. To understand how MSH2 regulates crossovers formed by the two pathways, we combined Arabidopsis mutants that elevate non-interfering crossovers with msh2 mutants. We demonstrate that MSH2 blocks non-interfering crossovers at polymorphic loci, which is the opposite effect to interfering crossovers. We also observe MSH2-independent crossover inhibition at highly polymorphic sites. We measure recombination along the chromosome arms in lines differing in patterns of heterozygosity and observe a MSH2-dependent crossover increase at the boundaries between heterozygous and homozygous regions. Here, we show that MSH2 is a master regulator of meiotic DSB repair in Arabidopsis, with antagonistic effects on interfering and non-interfering crossovers, which shapes the crossover landscape in relation to interhomolog polymorphism.

Retrotransposon addiction promotes centromere function via epigenetically activated small RNAs

Retrotransposons have invaded eukaryotic centromeres in cycles of repeat expansion and purging, but the function of centromeric retrotransposons, if any, has remained unclear. In Arabidopsis, centromeric ATHILA retrotransposons give rise to epigenetically activated short interfering RNAs (easiRNAs) in mutants in DECREASE IN DNA METHYLATION1 (DDM1), which promote histone H3 lysine-9 di-methylation (H3K9me2). Here, we show that mutants which lose both DDM1 and RNA dependent RNA polymerase (RdRP) have pleiotropic developmental defects and mis-segregation of chromosome 5 during mitosis. Fertility defects are epigenetically inherited with the centromeric region of chromosome 5, and can be rescued by directing artificial small RNAs to a single family of ATHILA5 retrotransposons specifically embedded within this centromeric region. easiRNAs and H3K9me2 promote pericentromeric condensation, chromosome cohesion and proper chromosome segregation in mitosis. Insertion of ATHILA silences transcription, while simultaneously making centromere function dependent on retrotransposon small RNAs, promoting the selfish survival and spread of centromeric retrotransposons. Parallels are made with the fission yeast S. pombe, where chromosome segregation depends on RNAi, and with humans, where chromosome segregation depends on both RNAi and HELLSDDM1.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

The regulation and potential functions of intronic satellite DNA

Satellite DNAs are arrays of tandem repeats found in the eukaryotic genome. They are mainly found in pericentromeric heterochromatin and have been believed to be mostly inert, leading satellite DNAs to be erroneously regarded as junk. Recent studies have started to elucidate the function of satellite DNA, yet little is known about the peculiar case where satellite DNA is found within the introns of protein coding genes, resulting in incredibly large introns, a phenomenon termed intron gigantism. Studies in Drosophila demonstrated that satellite DNA-containing introns are transcribed with the gene and require specialized mechanisms to overcome the burdens imposed by the extremely long stretches of repetitive DNA. Whether intron gigantism confers any benefit or serves any functional purpose for cells and/or organisms remains elusive. Here we review our current understanding of intron gigantism: where it is found, the challenges it imposes, how it is regulated and what purpose it may serve.

De novo centromere formation in pericentromeric region of rice chromosome 8

Neocentromeres develop when kinetochores assemble de novo at DNA loci that are not previously associated with CenH3 nucleosomes, and can rescue rearranged chromosomes that have lost a functional centromere. The molecular mechanisms associated with neocentromere formation in plants have been elusive. Here, we developed a Xian (indica) rice line with poor growth performance in the field due to approximately 272 kb deletion that spans centromeric DNA sequences, including the centromeric satellite repeat CentO, in the centromere of chromosome 8 (Cen8). The CENH3-binding domains were expanded downstream of the original CentO position in Cen8, which revealed a de novo centromere formation in rice. The neocentromere formation avoids chromosomal regions containing functional genes. Meanwhile, canonical histone H3 was replaced by CENH3 in the regions with low CENH3 levels, and the CenH3 nucleosomes in these regions became more periodic. In addition, we identified active genes in the deleted centromeric region, which are essential for chloroplast growth and development. In summary, our results provide valuable insights into neocentromere formation and show that functional genes exist in the centromeric regions of plant chromosomes.

Variation in the Number and Position of rDNA Loci Contributes to the Diversification and Speciation in Nigella (Ranunculaceae)

Nigella is a small genus belonging to the Ranunculaceae family which is presumably originated and distributed in Aegean and the adjacent Western-Irano-Turanian region. Comparative repeat analysis of N. sativa, N. damascena and N. bucharica was performed using low-pass Illumina genomic reads followed by karyotyping and FISH mapping of seven Nigella species using the in silico identified repeats and ribosomal DNA (rDNA) probes. High- and moderate-copy repeat sequences occupy 57.52, 59.01, and 64.73% of N. sativa, N. damascena and N. bucharica genomes, respectively. Roughly, half of the genomes are retrotransposons (class I transposons), while DNA transposons (class II transposons) contributed to only about 2% of the genomes. The analyzed Nigella species possess large genomes of about 7.4 to 12.4 Gbp/1C. Only two satellite repeats in N. sativa, one in N. damascena and four in N. bucharica were identified, which were mostly (peri)centromeric and represented about 1% of each genome. A high variation in number and position of 45S rDNA loci were found among Nigella species. Interestingly, in N. hispanica, each chromosome revealed at least one 45S rDNA site and one of them occurs in hemizygous condition. Based on the chromosome numbers, genome size and (peri)centromeric satellites, three karyotype groups were observed: Two with 2n = 2x = 12 and a karyotype formula of 10m + 2t (including N. sativa, N. arvensis, N. hispanica as the first group and N. damascena and N. orientalis as the second group) and a more distant group with 2n = 2x = 14 and a karyotype formula of 8m + 2st + 4t (including N. integrifolia and N. bucharica). These karyotype groups agreed with the phylogenetic analysis using ITS and rbcL sequences. We conclude that variation in (peri)centromeric sequences, number and localization of rDNA sites as well as chromosome number (dysploidy) are involved in the diversification of the genus Nigella.

Identification of a biomass unaffected pale green mutant gene in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Chlorophyll (Chl) is an essential component of the photosynthetic apparatus and pigments in plant greening. Leaf color is an important agronomic and commercial trait of Chinese cabbage. In this study, we identified a pale green mutant pgm created by ethyl methane sulfonate (EMS) mutagenesis in Chinese cabbage. Compared with wild-type (FT), pgm had a lower Chl content with a higher Chl a/b ratio, imperfect chloroplast structure, and lower non-photochemical quenching. However, its net photosynthetic rate and biomass showed no significant differences. Genetic analysis revealed that the pale green phenotype of pgm was controlled by a recessive nuclear gene, designated as Brpgm. We applied BSR-Seq, linkage analysis, and whole-genome resequencing to map Brpgm and predicted that the target gene was BraA10g007770.3C (BrCAO), which encodes chlorophyllide a oxygenase (CAO). Brcao sequencing results showed that the last nucleotide of its first intron changed from G to A, causing the deletion of the first nucleotide in its second CDS and termination of the protein translation. The expression of BrCAO in pgm was upregulated, and the enzyme activity of CAO in pgm was significantly decreased. These results provide an approach to explore the function of BrCAO and create a pale green variation in Chinese cabbage.

High-throughput estimation of allele frequencies using combined pooled-population sequencing and haplotype-based data processing

BACKGROUND

In addition to heterogeneity and artificial selection, natural selection is one of the forces used to combat climate change and improve agrobiodiversity in evolutionary plant breeding. Accurate identification of the specific genomic effects of natural selection will likely accelerate transfer between populations. Thus, insights into changes in allele frequency, adequate population size, gene flow and drift are essential. However, observing such effects often involves a trade-off between costs and resolution when a large sample of genotypes for many loci is analysed. Pool genotyping approaches achieve high resolution and precision in estimating allele frequency when sequence coverage is high. Nevertheless, high-coverage pool sequencing of large genomes is expensive.

RESULTS

Three pool samples (n = 300, 300, 288) from a barley backcross population were generated to assess the population's allele frequency. The tested population (BC2F21) has undergone 18 generations of natural adaption to conventional farming practice. The accuracies of estimated pool-based allele frequencies and genome coverage yields were compared using three next-generation sequencing genotyping methods. To achieve accurate allele frequency estimates with low sequence coverage, we employed a haplotyping approach. Low coverage allele frequencies of closely located single polymorphisms were aggregated into a single haplotype allele frequency, yielding 2-to-271-times higher depth and increased precision. When we combined different haplotyping tactics, we found that gene and chip marker-based haplotype analyses performed equivalently or better compared with simple contig haplotype windows. Comparing multiple pool samples and referencing against an individual sequencing approach revealed that whole-genome pool re-sequencing (WGS) achieved the highest correlation with individual genotyping (≥ 0.97). In contrast, transcriptome-based genotyping (MACE) and genotyping by sequencing (GBS) pool replicates were significantly associated with higher error rates and lower correlations, but are still valuable to detect large allele frequency variations.

CONCLUSIONS

The proposed strategy identified the allele frequency of populations with high accuracy at low cost. This is particularly relevant to evolutionary plant breeding of crops with very large genomes, such as barley. Whole-genome low coverage re-sequencing at 0.03 × coverage per genotype accurately estimated the allele frequency when a loci-based haplotyping approach was applied. The implementation of annotated haplotypes capitalises on the biological background and statistical robustness.

Characterization of Repetitive DNA in Saccharum officinarum and Saccharum spontaneum by Genome Sequencing and Cytological Assays

In most plant species, DNA repeated elements such as satellites and retrotransposons are composing the majority of their genomes. Saccharum officinarum (2n = 8x = 80) and S. spontaneum (2n = 40-128) are the two fundamental donors of modern sugarcane cultivars. These two species are polyploids with large genome sizes and are enriched in repetitive elements. In this work, we adopted a de novo strategy to isolate highly repetitive and abundant sequences in S. officinarum LA Purple and S. spontaneum SES208. The findings obtained from alignment to the genome assemblies revealed that the vast majority of the repeats (97.9% in LA Purple and 96.5% in SES208) were dispersed in the respective genomes. Fluorescence in situ hybridization assays were performed on 27 representative repeats to investigate their distributions and abundances. The results showed that the copies of some highly repeated sequences, including rDNA and centromeric or telomeric repeats, were underestimated in current genome assemblies. The analysis of the raw read mapping strategy showed more copy numbers for all studied repeats, suggesting that copy number underestimation is common for highly repeated sequences in current genome assemblies of LA Purple and SES208. In addition, the data showed that the centromeric retrotransposons in all SES208 centromeres were absent in certain S. spontaneum clones with different ploidies. This rapid turnover of centromeric DNA in sugarcane provides new clues regarding the pattern of centromeric retrotransposon formation and accumulation.

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.

Differentially Amplified Repetitive Sequences Among Aegilops tauschii Subspecies and Genotypes

Genomic repetitive sequences commonly show species-specific sequence type, abundance, and distribution patterns, however, their intraspecific characteristics have been poorly described. We quantified the genomic repetitive sequences and performed single nucleotide polymorphism (SNP) analysis between 29 Ae. tauschii genotypes and subspecies using publicly available raw genomic Illumina sequence reads and used fluorescence in situ hybridization (FISH) to experimentally analyze some repeats. The majority of the identified repetitive sequences had similar contents and proportions between anathera, meyeri, and strangulata subspecies. However, two Ty3/gypsy retrotransposons (CL62 and CL87) showed significantly higher abundances, and CL1, CL119, CL213, CL217 tandem repeats, and CL142 retrotransposon (Ty1/copia type) showed significantly lower abundances in subspecies strangulata compared with the subspecies anathera and meyeri. One tandem repeat and 45S ribosomal DNA (45S rDNA) abundances showed a high variation between genotypes but their abundances were not subspecies specific. Phylogenetic analysis using the repeat abundances of the aforementioned clusters placed the strangulata subsp. in a distinct clade but could not discriminate anathera and meyeri. A near complete differentiation of anathera and strangulata subspecies was observed using SNP analysis; however, var. meyeri showed higher genetic diversity. FISH using major tandem repeats couldn't detect differences between subspecies, although (GAA)₁₀ signal patterns generated two different karyotype groups. Taken together, the different classes of repetitive DNA sequences have differentially accumulated between strangulata and the other two subspecies of Ae. tauschii that is generally in agreement with spike morphology, implying that factors affecting repeatome evolution are variable even among highly closely related lineages.

The formation and evolution of centromeric satellite repeats in Saccharum species

Centromeres in eukaryotes are composed of tandem DNAs and retrotransposons. However, centromeric repeats exhibit considerable diversity, even among closely related species, and their origin and evolution are largely unknown. We conducted a genome-wide characterization of the centromeric sequences in sugarcane (Saccharum officinarum). Four centromeric tandem repeat sequences, So1, So103, So137 and So119, were isolated. So1 has a monomeric length of 137 bp, typical of a centromeric satellite, and has evolved four variants. However, these So1 variants had distinct centromere distributions and some were unique to an individual centromere. The distributions of the So1 variants were unexpectedly consistent among the Saccharum species that had different basic chromosome numbers or ploidy levels, thus suggesting evolutionary stability for approximately 7 million years in sugarcane. So103, So137 and So119 had unusually longer monomeric lengths that ranged from 327 to 1371 bp and lacked translational phasing on the CENH3 nucleosomes. Moreover, So103, So137 and So119 seemed to be highly similar to retrotransposons, which suggests that they originated from these mobile elements. Notably, all three repeats were flanked by direct repeats, and formed extrachromosomal circular DNAs (eccDNAs). The presence of circular molecules for these retrotransposon-derived centromeric satellites suggests an eccDNA-mediated centromeric satellite formation pathway in sugarcane.

Role of H1 and DNA methylation in selective regulation of transposable elements during heat stress

Heat-stressed Arabidopsis plants release heterochromatin-associated transposable element (TE) silencing, yet it is not accompanied by major reductions of epigenetic repressive modifications. In this study, we explored the functional role of histone H1 in repressing heterochromatic TEs in response to heat stress. We generated and analyzed RNA and bisulfite-sequencing data of wild-type and h1 mutant seedlings before and after heat stress. Loss of H1 caused activation of pericentromeric Gypsy elements upon heat treatment, despite these elements remaining highly methylated. By contrast, nonpericentromeric Copia elements became activated concomitantly with loss of DNA methylation. The same Copia elements became activated in heat-treated chromomethylase 2 (cmt2) mutants, indicating that H1 represses Copia elements through maintaining DNA methylation under heat. We discovered that H1 is required for TE repression in response to heat stress, but its functional role differs depending on TE location. Strikingly, H1-deficient plants treated with the DNA methyltransferase inhibitor zebularine were highly tolerant to heat stress, suggesting that both H1 and DNA methylation redundantly suppress the plant response to heat stress.

A dCas9-Based System Identifies a Central Role for Ctf19 in Kinetochore-Derived Suppression of Meiotic Recombination

In meiosis, crossover (CO) formation between homologous chromosomes is essential for faithful segregation. However, misplaced meiotic recombination can have catastrophic consequences on genome stability. Within pericentromeres, COs are associated with meiotic chromosome missegregation. In organisms ranging from yeast to humans, pericentromeric COs are repressed. We previously identified a role for the kinetochore-associated Ctf19 complex (Ctf19c) in pericentromeric CO suppression. Here, we develop a dCas9/CRISPR-based system that allows ectopic targeting of Ctf19c-subunits. Using this approach, we query sufficiency in meiotic CO suppression, and identify Ctf19 as a mediator of kinetochore-associated CO control. The effect of Ctf19 is encoded in its NH2-terminal tail, and depends on residues important for the recruitment of the Scc2-Scc4 cohesin regulator. This work provides insight into kinetochore-derived control of meiotic recombination. We establish an experimental platform to investigate and manipulate meiotic CO control. This platform can easily be adapted in order to investigate other aspects of chromosome biology.

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.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha

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

Transcriptome sequencing reveals hotspot mutation regions and dwarfing mechanisms in wheat mutants induced by γ-ray irradiation and EMS

Induced mutation is an important approach for creating novel plant germplasms. The introduction of dwarf or semi-dwarf genes into wheat has led to great advancements in yield improvement. In this study, four elite dwarf wheat mutants, named dm1-dm4, induced from γ-ray irradiation or ethyl methanesulfonate (EMS) mutagenesis, were used to identify transcriptome variations and dwarfing mechanisms. The results showed that the hotspot regions of mutations distributed on the chromosomes were consistent among the four mutant lines and these regions were mainly located around the 50, 360 and 400 Mb positions of chromosome 1A and the distal regions of chromosomes 2A and 2BL. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that 'protein processing in endoplasmic reticulum' was the most common significantly enriched pathway based on the differentially expressed genes (DEGs) between wildtype (WT) and the mutants. Notably, 18 out of 20 genes involved in this process encode heat shock proteins (HSPs). The results implied that HSPs might participate in wheat dwarfism response and function in the dwarfism process through protein folding and/or degradation. Moreover, seven genes in dm4 involved in modulating auxin levels were down-regulated and dm4 was more sensitive to auxin treatment compared with WT, indicating the important roles of auxin in regulation of dwarf phenotype in dm4. This study not only identified transcriptome sequence variation induced by physical and chemical mutagenesis but also revealed potential dwarfing mechanisms in the wheat mutant lines.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

Heterogeneous transposable elements as silencers, enhancers and targets of meiotic recombination

During meiosis, DNA double-strand breaks are initiated by the topoisomerase-like enzyme SPO11 and are repaired by inter-sister chromatid and inter-homologue DNA repair pathways. Genome-wide maps of initiating DNA double-strand breaks and inter-homologue repair events are now available for a number of mammalian, fungal and plant species. In mammals, PRDM9 specifies the location of meiotic recombination initiation via recognition of specific DNA sequence motifs by its C2H2 zinc finger array. In fungi and plants, meiotic recombination appears to be initiated less discriminately in accessible chromatin, including at gene promoters. Generally, meiotic crossover is suppressed in highly repetitive genomic regions that are made up of transposable elements (TEs), to prevent deleterious non-allelic homologous recombination events. However, recent and older studies have revealed intriguing relationships between meiotic recombination initiation and repair, and transposable elements. For instance, gene conversion events have been detected in maize centromeric retroelements, mouse MULE-MuDR DNA transposons undergo substantial meiotic recombination initiation, Arabidopsis Helitron TEs are among the hottest of recombination initiation hotspots, and human TE sequences can modify the crossover rate at adjacent PRDM9 motifs in cis. Here, we summarize the relationship between meiotic recombination and TEs, discuss recent insights from highly divergent eukaryotes and highlight outstanding questions in the field.

A New Portrait of Constitutive Heterochromatin: Lessons from Drosophila melanogaster

Constitutive heterochromatin represents a significant portion of eukaryotic genomes, but its functions still need to be elucidated. Even in the most updated genetics and molecular biology textbooks, constitutive heterochromatin is portrayed mainly as the 'silent' component of eukaryotic genomes. However, there may be more complexity to the relationship between heterochromatin and gene expression. In the fruit fly Drosophila melanogaster, a model for heterochromatin studies, about one-third of the genome is heterochromatic and is concentrated in the centric, pericentric, and telomeric regions of the chromosomes. Recent findings indicate that hundreds of D. melanogaster genes can 'live and work' properly within constitutive heterochromatin. The genomic size of these genes is generally larger than that of euchromatic genes and together they account for a significant fraction of the entire constitutive heterochromatin. Thus, this peculiar genome component in spite its ability to induce silencing, has in fact the means for being quite dynamic. A major scope of this review is to revisit the 'dogma of silent heterochromatin'.

Wavelet-Based Genomic Signal Processing for Centromere Identification and Hypothesis Generation

Various ‘omics data types have been generated for Populus trichocarpa, each providing a layer of information which can be represented as a density signal across a chromosome. We make use of genome sequence data, variants data across a population as well as methylation data across 10 different tissues, combined with wavelet-based signal processing to perform a comprehensive analysis of the signature of the centromere in these different data signals, and successfully identify putative centromeric regions in P. trichocarpa from these signals. Furthermore, using SNP (single nucleotide polymorphism) correlations across a natural population of P. trichocarpa, we find evidence for the co-evolution of the centromeric histone CENH3 with the sequence of the newly identified centromeric regions, and identify a new CENH3 candidate in P. trichocarpa.

Centromere Transcription: Means and Motive

The chromosome biology field at large has benefited from studies of the cell cycle components, protein cascades and genomic landscape that are required for centromere identity, assembly and stable transgenerational inheritance. Research over the past 20 years has challenged the classical descriptions of a centromere as a stable, unmutable, and transcriptionally silent chromosome component. Instead, based on studies from a broad range of eukaryotic species, including yeast, fungi, plants, and animals, the centromere has been redefined as one of the more dynamic areas of the eukaryotic genome, requiring coordination of protein complex assembly, chromatin assembly, and transcriptional activity in a cell cycle specific manner. What has emerged from more recent studies is the realization that the transcription of specific types of nucleic acids is a key process in defining centromere integrity and function. To illustrate the transcriptional landscape of centromeres across eukaryotes, we focus this review on how transcripts interact with centromere proteins, when in the cell cycle centromeric transcription occurs, and what types of sequences are being transcribed. Utilizing data from broadly different organisms, a picture emerges that places centromeric transcription as an integral component of centromere function.

Meiotic recombination within plant centromeres

Meiosis is a conserved eukaryotic cell division that increases genetic diversity in sexual populations. During meiosis homologous chromosomes pair and undergo recombination that can result in reciprocal genetic exchange, termed crossover. The frequency of crossover is highly variable along chromosomes, with hot spots and cold spots. For example, the centromeres that contain the kinetochore, which attach chromosomes to the microtubular spindle, are crossover cold spots. Plant centromeres typically consist of large tandemly repeated arrays of satellite sequences and retrotransposons, a subset of which assemble CENH3-variant nucleosomes, which bind to kinetochore proteins. Although crossovers are suppressed in centromeres, there is abundant evidence for gene conversion and homologous recombination between repeats, which plays a role in satellite array change. We review the evidence for recombination within plant centromeres and the implications for satellite sequence evolution. We speculate on the genetic and epigenetic features of centromeres that may influence meiotic recombination in these regions. We also highlight unresolved questions relating to centromere function and sequence change and how the advent of new technologies promises to provide insights.

Non-B-Form DNA Is Enriched at Centromeres

Animal and plant centromeres are embedded in repetitive "satellite" DNA, but are thought to be epigenetically specified. To define genetic characteristics of centromeres, we surveyed satellite DNA from diverse eukaryotes and identified variation in <10-bp dyad symmetries predicted to adopt non-B-form conformations. Organisms lacking centromeric dyad symmetries had binding sites for sequence-specific DNA-binding proteins with DNA-bending activity. For example, human and mouse centromeres are depleted for dyad symmetries, but are enriched for non-B-form DNA and are associated with binding sites for the conserved DNA-binding protein CENP-B, which is required for artificial centromere function but is paradoxically nonessential. We also detected dyad symmetries and predicted non-B-form DNA structures at neocentromeres, which form at ectopic loci. We propose that centromeres form at non-B-form DNA because of dyad symmetries or are strengthened by sequence-specific DNA binding proteins. This may resolve the CENP-B paradox and provide a general basis for centromere specification.

Kinetochores, cohesin, and DNA breaks: Controlling meiotic recombination within pericentromeres

In meiosis, DNA break formation and repair are essential for the formation of crossovers between homologous chromosomes. Without crossover formation, faithful meiotic chromosome segregation and sexual reproduction cannot occur. Crossover formation is initiated by the programmed, meiosis-specific introduction of numerous DNA double-strand breaks, after which specific repair pathways promote recombination between homologous chromosomes. Despite its crucial nature, meiotic recombination is fraud with danger: When positioned or repaired inappropriately, DNA breaks can have catastrophic consequences on genome stability of the resulting gametes. As such, DNA break formation and repair needs to be carefully controlled. Within centromeres and surrounding regions (i.e., pericentromeres), meiotic crossover recombination is repressed in organisms ranging from yeast to humans, and a failure to do so is implicated in chromosome missegregation and developmental aneuploidy. (Peri)centromere sequence identity and organization diverge considerably across eukaryotes, yet suppression of meiotic DNA break formation and repair appear universal. Here, we discuss emerging work that has used budding and fission yeast systems to study the mechanisms underlying pericentromeric suppression of DNA break formation and repair. We particularly highlight a role for the kinetochore, a universally conserved, centromere-associated structure essential for chromosome segregation, in suppressing (peri)centromeric DNA break formation and repair. We discuss the current understanding of kinetochore-associated and chromosomal factors involved in this regulation and suggest future avenues of research.

Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing

A male sterile mutant was created by ⁶⁰Co γ-rays of microspores isolated from Chinese cabbage DH line 'FT'. A candidate gene for the male sterile trait was identified as Bra010198. Male sterility is used for hybrid seed production in Chinese cabbage. In this study, we derived a male sterile mutant (ftms) from Chinese cabbage DH line 'FT' by irradiating microspores with ⁶⁰Co γ-rays and realized the rapid trait transformation from male fertility to sterility for creating valuable breeding materials. Genetic analysis indicated that the male sterile trait is controlled by a single recessive nuclear gene, ftms. Microspore development in mutant ftms was aborted at the tetrad stage and associated with severely retarded degeneration and vacuolation of tapetum. Using BSR-seq analysis, the candidate region for ftms was mapped on chromosome A05. A large F₂ population was created, and the region was narrowed to approximately 1.7-Mb between markers Indel20 and Indel14 via linkage analysis. The recombination frequency was extremely suppressed because the region was located on the chromosome A05 centromere. Whole-genome resequencing of mutant ftms and wild-type 'FT' aligned only one nonsynonymous SNP to Bra010198; this gene is a homolog of Arabidopsis KNS4/UPEX1, which encodes a putative β-(1,3)-galactosyltransferase that controls pollen exine development. Comparative sequencing verified the SNP position on the fifth exon of Bra010198 in mutant ftms. Further genotyping revealed that the male sterile phenotype was fully co-segregated with this SNP. Quantitative real-time PCR indicated that Bra0101918 specifically expressed in stamen. The data presented herein suggested that Bra010198 is a strong candidate gene for ftms. Hence, we developed a male sterile line for potential application in breeding and expanded the knowledge about the molecular mechanism underlying male sterility in Chinese cabbage.

Where to Cross Over? Defining Crossover Sites in Plants

It is believed that recombination in meiosis serves to reshuffle genetic material from both parents to increase genetic variation in the progeny. At the same time, the number of crossovers is usually kept at a very low level. As a consequence, many organisms need to make the best possible use from the one or two crossovers that occur per chromosome in meiosis. From this perspective, the decision of where to allocate rare crossover events becomes an important issue, especially in self-pollinating plant species, which experience limited variation due to inbreeding. However, the freedom in crossover allocation is significantly limited by other, genetic and non-genetic factors, including chromatin structure. Here we summarize recent progress in our understanding of those processes with a special emphasis on plant genomes. First, we focus on factors which influence the distribution of recombination initiation sites and discuss their effects at both, the single hotspot level and at the chromosome scale. We also briefly explain the aspects of hotspot evolution and their regulation. Next, we analyze how recombination initiation sites translate into the development of crossovers and their location. Moreover, we provide an overview of the sequence polymorphism impact on crossover formation and chromosomal distribution.

The Genomic Architecture of a Rapid Island Radiation: Recombination Rate Variation, Chromosome Structure, and Genome Assembly of the Hawaiian Cricket Laupala

Phenotypic evolution and speciation depend on recombination in many ways. Within populations, recombination can promote adaptation by bringing together favorable mutations and decoupling beneficial and deleterious alleles. As populations diverge, crossing over can give rise to maladapted recombinants and impede or reverse diversification. Suppressed recombination due to genomic rearrangements, modifier alleles, and intrinsic chromosomal properties may offer a shield against maladaptive gene flow eroding coadapted gene complexes. Both theoretical and empirical results support this relationship. However, little is known about this relationship in the context of behavioral isolation, where coevolving signals and preferences are the major hybridization barrier. Here we examine the genomic architecture of recently diverged, sexually isolated Hawaiian swordtail crickets (Laupala). We assemble a de novo genome and generate three dense linkage maps from interspecies crosses. In line with expectations based on the species' recent divergence and successful interbreeding in the laboratory, the linkage maps are highly collinear and show no evidence for large-scale chromosomal rearrangements. Next, the maps were used to anchor the assembly to pseudomolecules and estimate recombination rates across the genome to test the hypothesis that loci involved in behavioral isolation (song and preference divergence) are in regions of low interspecific recombination. Contrary to our expectations, the genomic region where a male song and female preference QTL colocalize is not associated with particularly low recombination rates. This study provides important novel genomic resources for an emerging evolutionary genetics model system and suggests that trait-preference coevolution is not necessarily facilitated by locally suppressed recombination.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Epigenetic activation of meiotic recombination near Arabidopsis thaliana centromeres via loss of H3K9me2 and non-CG DNA methylation

Eukaryotic centromeres contain the kinetochore, which connects chromosomes to the spindle allowing segregation. During meiosis, centromeres are suppressed for inter-homolog crossover, as recombination in these regions can cause chromosome missegregation and aneuploidy. Plant centromeres are surrounded by transposon-dense pericentromeric heterochromatin that is epigenetically silenced by histone 3 lysine 9 dimethylation (H3K9me2), and DNA methylation in CG and non-CG sequence contexts. However, the role of these chromatin modifications in control of meiotic recombination in the pericentromeres is not fully understood. Here, we show that disruption of Arabidopsis thaliana H3K9me2 and non-CG DNA methylation pathways, for example, via mutation of the H3K9 methyltransferase genes KYP/SUVH4 SUVH5 SUVH6, or the CHG DNA methyltransferase gene CMT3, increases meiotic recombination in proximity to the centromeres. Using immunocytological detection of MLH1 foci and genotyping by sequencing of recombinant plants, we observe that H3K9me2 and non-CG DNA methylation pathway mutants show increased pericentromeric crossovers. Increased pericentromeric recombination in H3K9me2/non-CG mutants occurs in hybrid and inbred backgrounds and likely involves contributions from both the interfering and noninterfering crossover repair pathways. We also show that meiotic DNA double-strand breaks (DSBs) increase in H3K9me2/non-CG mutants within the pericentromeres, via purification and sequencing of SPO11-1-oligonucleotides. Therefore, H3K9me2 and non-CG DNA methylation exert a repressive effect on both meiotic DSB and crossover formation in plant pericentromeric heterochromatin. Our results may account for selection of enhancer trap Dissociation (Ds) transposons into the CMT3 gene by recombination with proximal transposon launch-pads.

Massive crossover elevation via combination of HEI10 and recq4a recq4b during Arabidopsis meiosis

During meiosis, homologous chromosomes undergo reciprocal crossovers, which generate genetic diversity and underpin classical crop improvement. Meiotic recombination initiates from DNA double-strand breaks (DSBs), which are processed into single-stranded DNA that can invade a homologous chromosome. The resulting joint molecules can ultimately be resolved as crossovers. In Arabidopsis, competing pathways balance the repair of ∼100-200 meiotic DSBs into ∼10 crossovers per meiosis, with the excess DSBs repaired as noncrossovers. To bias DSB repair toward crossovers, we simultaneously increased dosage of the procrossover E3 ligase gene HEI10 and introduced mutations in the anticrossovers helicase genes RECQ4A and RECQ4B As HEI10 and recq4a recq4b increase interfering and noninterfering crossover pathways, respectively, they combine additively to yield a massive meiotic recombination increase. Interestingly, we also show that increased HEI10 dosage increases crossover coincidence, which indicates an effect on interference. We also show that patterns of interhomolog polymorphism and heterochromatin drive recombination increases distally towards the subtelomeres in both HEI10 and recq4a recq4b backgrounds, while the centromeres remain crossover suppressed. These results provide a genetic framework for engineering meiotic recombination landscapes in plant genomes.

H3K23me1 is an evolutionarily conserved histone modification associated with CG DNA methylation in Arabidopsis

Amino-terminal tails of histones are targets for diverse post-translational modifications whose combinatorial action may constitute a code that will be read and interpreted by cellular proteins to define particular transcriptional states. Here, we describe monomethylation of histone H3 lysine 23 (H3K23me1) as a histone modification not previously described in plants. H3K23me1 is an evolutionarily conserved mark in diverse species of flowering plants. Chromatin immunoprecipitation followed by high-throughput sequencing in Arabidopsis thaliana showed that H3K23me1 was highly enriched in pericentromeric regions and depleted from chromosome arms. In transposable elements it co-localized with CG, CHG and CHH DNA methylation as well as with the heterochromatic histone mark H3K9me2. Transposable elements are often rich in H3K23me1 but different families vary in their enrichment: LTR-Gypsy elements are most enriched and RC/Helitron elements are least enriched. The histone methyltransferase KRYPTONITE and normal DNA methylation were required for normal levels of H3K23me1 on transposable elements. Immunostaining experiments confirmed the pericentromeric localization and also showed mild enrichment in less condensed regions. Accordingly, gene bodies of protein-coding genes had intermediate H3K23me1 levels, which coexisted with CG DNA methylation. Enrichment of H3K23me1 along gene bodies did not correlate with transcription levels. Together, this work establishes H3K23me1 as a so far undescribed component of the plant histone code.

Genetic and epigenetic variation of transposable elements in Arabidopsis

Transposable elements are mobile genetic elements that are prevalent in plant genomes and are silenced by epigenetic modification. Different epigenetic modification pathways play distinct roles in the control of transposable element transcription, replication and recombination. The Arabidopsis genome contains families of all of the major transposable element classes, which are differentially enriched in particular genomic regions. Whole genome sequencing and DNA methylation profiling of hundreds of natural Arabidopsis accessions has revealed that transposable elements exhibit significant intraspecific genetic and epigenetic variation, and that genetic variation often underlies epigenetic variation. Together, epigenetic modification and the forces of selection define the scope within which transposable elements can contribute to, and control, genome evolution.

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.

Natural variation and dosage of the HEI10 meiotic E3 ligase control Arabidopsis crossover recombination

During meiosis, homologous chromosomes undergo crossover recombination, which creates genetic diversity and balances homolog segregation. Despite these critical functions, crossover frequency varies extensively within and between species. Although natural crossover recombination modifier loci have been detected in plants, causal genes have remained elusive. Using natural Arabidopsis thaliana accessions, we identified two major recombination quantitative trait loci (rQTLs) that explain 56.9% of crossover variation in Col×Ler F2 populations. We mapped rQTL1 to semidominant polymorphisms in HEI10, which encodes a conserved ubiquitin E3 ligase that regulates crossovers. Null hei10 mutants are haploinsufficient, and, using genome-wide mapping and immunocytology, we show that transformation of additional HEI10 copies is sufficient to more than double euchromatic crossovers. However, heterochromatic centromeres remained recombination-suppressed. The strongest HEI10-mediated crossover increases occur in subtelomeric euchromatin, which is reminiscent of sex differences in Arabidopsis recombination. Our work reveals that HEI10 naturally limits Arabidopsis crossovers and has the potential to influence the response to selection.

Centromeres Drive a Hard Bargain

Centromeres are essential chromosomal structures that mediate the accurate distribution of genetic material during meiotic and mitotic cell divisions. In most organisms, centromeres are epigenetically specified and propagated by nucleosomes containing the centromere-specific H3 variant, centromere protein A (CENP-A). Although centromeres perform a critical and conserved function, CENP-A and the underlying centromeric DNA are rapidly evolving. This paradox has been explained by the centromere drive hypothesis, which proposes that CENP-A is undergoing an evolutionary tug-of-war with selfish centromeric DNA. Here, we review our current understanding of CENP-A evolution in relation to centromere drive and discuss classical and recent advances, including new evidence implicating CENP-A chaperones in this conflict.

Quantification and Sequencing of Crossover Recombinant Molecules from Arabidopsis Pollen DNA

During meiosis, homologous chromosomes undergo recombination, which can result in formation of reciprocal crossover molecules. Crossover frequency is highly variable across the genome, typically occurring in narrow hotspots, which has a significant effect on patterns of genetic diversity. Here we describe methods to measure crossover frequency in plants at the hotspot scale (bp-kb), using allele-specific PCR amplification from genomic DNA extracted from the pollen of F1 heterozygous plants. We describe (1) titration methods that allow amplification, quantification and sequencing of single crossover molecules, (2) quantitative PCR methods to more rapidly measure crossover frequency, and (3) application of high-throughput sequencing for study of crossover distributions within hotspots. We provide detailed descriptions of key steps including pollen DNA extraction, prior identification of hotspot locations, allele-specific oligonucleotide design, and sequence analysis approaches. Together, these methods allow the rate and recombination topology of plant hotspots to be robustly measured and compared between varied genetic backgrounds and environmental conditions.

Visualization of specific repetitive genomic sequences with fluorescent TALEs in Arabidopsis thaliana

Live imaging of the dynamics of nuclear organization provides the opportunity to uncover the mechanisms responsible for four-dimensional genome architecture. Here, we describe the use of fluorescent protein (FP) fusions of transcription activator-like effectors (TALEs) to visualize endogenous genomic sequences in Arabidopsis thaliana. The ability to engineer sequence-specific TALEs permits the investigation of precise genomic sequences. We could detect TALE-FP signals associated with centromeric, telomeric, and rDNA repeats and the signal distribution was consistent with that observed by fluorescent in situ hybridization. TALE-FPs are advantageous because they permit the observation of intact tissues. We used our TALE-FP method to investigate the nuclei of several multicellular plant tissues including roots, hypocotyls, leaves, and flowers. Because TALE-FPs permit live-cell imaging, we successfully observed the temporal dynamics of centromeres and telomeres in plant organs. Fusing TALEs to multimeric FPs enhanced the signal intensity when observing telomeres. We found that the mobility of telomeres was different in sub-nuclear regions. Transgenic plants stably expressing TALE-FPs will provide new insights into chromatin organization and dynamics in multicellular organisms.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.

Hybrid Sterility in Rice (Oryza sativa L.) Involves the Tetratricopeptide Repeat Domain Containing Protein

Intersubspecific hybrid sterility is a common form of reproductive isolation in rice (Oryza sativa L.), which significantly hampers the utilization of heterosis between indica and japonica varieties. Here, we elucidated the mechanism of S7, which specially causes Aus-japonica/indica hybrid female sterility, through cytological and genetic analysis, map-based cloning, and transformation experiments. Abnormal positioning of polar nuclei and smaller embryo sac were observed in F1 compared with male and female parents. Female gametes carrying S7(cp) and S7(i) were aborted in S7(ai)/S7(cp) and S7(ai)/S7(i), respectively, whereas they were normal in both N22 and Dular possessing a neutral allele, S7(n) S7 was fine mapped to a 139-kb region in the centromere region on chromosome 7, where the recombination was remarkably suppressed due to aggregation of retrotransposons. Among 16 putative open reading frames (ORFs) localized in the mapping region, ORF3 encoding a tetratricopeptide repeat domain containing protein was highly expressed in the pistil. Transformation experiments demonstrated that ORF3 is the candidate gene: downregulated expression of ORF3 restored spikelet fertility and eliminated absolutely preferential transmission of S7(ai) in heterozygote S7(ai)/S7(cp); sterility occurred in the transformants Cpslo17-S7(ai) Our results may provide implications for overcoming hybrid embryo sac sterility in intersubspecific hybrid rice and utilization of hybrid heterosis for cultivated rice improvement.

A Gossypium BAC clone contains key repeat components distinguishing sub-genome of allotetraploidy cottons

Background

Dissecting genome organization is indispensable for further functional and applied studies. As genome sequences data shown, cotton genomes contain more than 60 % repetitive sequences, so study on repetitive sequences composition, structure, and distribution is the key step to dissect cotton genome.

Results

In this study, a bacterial artificial chromosome (BAC) clone enriched in repetitive sequences, was discovered initiatively by fluorescence in situ hybridization (FISH). FISHing with allotetraploidy cotton as target DNA, dispersed signals on most regions of all A sub-genome chromosomes, and only middle regions of all D sub-genome chromosomes were detected. Further FISHing with other cotton species bearing A or D genome as target DNA, specific signals were viewed. After BAC sequencing and bioinformational analysis, 129 repeat elements, size about 57,172 bp were found, accounting for more than 62 % of the BAC sequence (91,238 bp). Among them, a type of long terminal repeat-retrotransposon (LTR-RT), LTR/Gypsy was the key element causing the specific FISH results. Using the fragments of BAC matching with the identified Gypsy-like LTR as probes, the BAC-57I23-like FISH signals were reappeared. Running BLASTN, the fragments had good match with all chromosomes of G. arboreum (A₂) genome and A sub-genome of G. hirsutum (AD₁), and had relatively inferior match with all chromosomes of D sub-genome of AD₁, but had little match with the chromosomes of G. raimondii (D₅) genome, which was consistent with the FISH results.

Conclusion

A repeats-enriched cytogenetic marker to identify A and D sub-genomes of Gossypium was discovered by FISH. Combined sequences analysis with FISH verification, the assembly quality of repetitive sequences in the allotetraploidy cotton draft genome was assessed, and better chromosome belonging was verified. We also found the genomic distribution of the identified Gypsy-LTR-RT was similar to the distribution of heterochromatin. The expansion of this type of Gypsy-LTR-RT in heterochromatic regions may be one of the major reasons for the size gap between A and D genome. The findings showed here will help to understand the composition, structure, and evolution of cotton genome, and contribute to the further perfection of the draft genomes of cotton.

Centromere Locations in Brassica A and C Genomes Revealed Through Half-Tetrad Analysis

Locating centromeres on genome sequences can be challenging. The high density of repetitive elements in these regions makes sequence assembly problematic, especially when using short-read sequencing technologies. It can also be difficult to distinguish between active and recently extinct centromeres through sequence analysis. An effective solution is to identify genetically active centromeres (functional in meiosis) by half-tetrad analysis. This genetic approach involves detecting heterozygosity along chromosomes in segregating populations derived from gametes (half-tetrads). Unreduced gametes produced by first division restitution mechanisms comprise complete sets of nonsister chromatids. Along these chromatids, heterozygosity is maximal at the centromeres, and homologous recombination events result in homozygosity toward the telomeres. We genotyped populations of half-tetrad-derived individuals (from Brassica interspecific hybrids) using a high-density array of physically anchored SNP markers (Illumina Brassica 60K Infinium array). Mapping the distribution of heterozygosity in these half-tetrad individuals allowed the genetic mapping of all 19 centromeres of the Brassica A and C genomes to the reference Brassica napus genome. Gene and transposable element density across the B. napus genome were also assessed and corresponded well to previously reported genetic map positions. Known centromere-specific sequences were located in the reference genome, but mostly matched unanchored sequences, suggesting that the core centromeric regions may not yet be assembled into the pseudochromosomes of the reference genome. The increasing availability of genetic markers physically anchored to reference genomes greatly simplifies the genetic and physical mapping of centromeres using half-tetrad analysis. We discuss possible applications of this approach, including in species where half-tetrads are currently difficult to isolate.