How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae)

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

Oligo-barcode illuminates holocentric karyotype evolution in Rhynchospora (Cyperaceae)

Holocentric karyotypes are assumed to rapidly evolve through chromosome fusions and fissions due to the diffuse nature of their centromeres. Here, we took advantage of the recent availability of a chromosome-scale reference genome for Rhynchospora breviuscula, a model species of this holocentric genus, and developed the first set of oligo-based barcode probes for a holocentric plant. These probes were applied to 13 additional species of the genus, aiming to investigate the evolutionary dynamics driving the karyotype evolution in Rhynchospora. The two sets of probes were composed of 27,392 (green) and 23,968 (magenta) oligonucleotides (45-nt long), and generated 15 distinct FISH signals as a unique barcode pattern for the identification of all five chromosome pairs of the R. breviuscula karyotype. Oligo-FISH comparative analyzes revealed different types of rearrangements, such as fusions, fissions, putative inversions and translocations, as well as genomic duplications among the analyzed species. Two rounds of whole genome duplication (WGD) were demonstrated in R. pubera, but both analyzed accessions differed in the complex chain of events that gave rise to its large, structurally diploidized karyotypes with 2n = 10 or 12. Considering the phylogenetic relationships and divergence time of the species, the specificity and synteny of the probes were maintained up to species with a divergence time of ~25 My. However, karyotype divergence in more distant species hindered chromosome mapping and the inference of specific events. This barcoding system is a powerful tool to study chromosomal variations and genomic evolution in holocentric chromosomes of Rhynchospora species.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Investigating the diversification of holocentromeric satellite DNA Tyba in Rhynchospora (Cyperaceae)

BACKGROUND AND AIMS

Satellite DNAs (satDNAs) are repetitive sequences composed by tandemly arranged, often highly homogenized units called monomers. Although satDNAs are usually fast evolving, some satDNA families can be conserved across species separated by several millions of years, probably because of their functional roles in the genomes. Tyba was the first centromere-specific satDNA described for a holocentric organism, until now being characterized for only eight species of the genus Rhynchospora Vahl. (Cyperaceae). Here, we characterized Tyba across a broad sampling of the genus, analysing and comparing its evolutionary patterns with other satDNAs.

METHODS

We characterized the structure and sequence evolution of satDNAs across a robust dadated phylogeny based on Hybrid Target-Capture Sequencing (hyb-seq) of 70 species. We mined the repetitive fraction for Tyba-like satellites to compare its features with other satDNAs and to construct a Tyba-based phylogeny for the genus.

KEY RESULTS

Our results show that Tyba is present in the majority of examined species of the genus, spanning four of the five major clades and maintaining intrafamily pairwise identity of 70.9% over 31 Myr. In comparison, other satellite families presented higher intrafamily pairwise identity but are phylogenetically restricted. Furthermore, Tyba sequences could be divided into 12 variants grouped into three different clade-specific subfamilies, showing evidence of traditional models of satDNA evolution, such as the concerted evolution and library models. Besides, a Tyba-based phylogeny showed high congruence with the hyb-seq topology. Our results show structural indications of a possible relationship of Tyba with nucleosomes, given its high curvature peaks over conserved regions and overall high bendability values compared with other non-centromeric satellites.

CONCLUSIONS

Overall, Tyba shows a remarkable sequence conservation and phylogenetic significance across the genus Rhynchospora, which suggests that functional roles might lead to long-term stability and conservation for satDNAs in the genome.

Cytogenetic characterization and karyotype evolution in six Macroptilium species (Leguminosae)

Macroptilium (Benth.) Urb. is a neotropical legume genus from the subtribe Phaseolinae. The investigated species present a stable chromosome number (2n = 22), but differ in their karyotype formulae, suggesting the presence of chromosome rearrangements. In this work, we comparatively analysed the karyotypes of six species (M. atropurpureum, M. bracteatum, M. erythroloma, M. gracile, M. lathyroides and M. martii) from the two main clades that form the genus. Heterochromatin distribution was investigated with CMA/DAPI staining and fluorescent in situ hybridization (FISH) was used to localise the 5S and 35S ribosomal DNA (rDNA) sites. Single copy BACs (Bacterial Artificial Chromosomes) previously mapped in the related genera Phaseolus L. and Vigna Savi were used to establish chromosome orthologies and to investigate possible rearrangements among species. CMA+/DAPI- bands were observed, mostly associated with rDNA sites. Additional weak, pericentromeric bands were observed on several chromosomes. Although karyotypes were similar, species could be differentiated mainly by the number and position of the 5S and 35S rDNA sites. BAC markers demonstrated conserved synteny of the main rDNA sites on orthologous chromosomes 6 and 10, as previously observed for Phaseolus and Vigna. The karyotypes of the six species could be differentiated, shedding light on its karyotype evolution.

Differential Repeat Accumulation in the Bimodal Karyotype of Agave L

The genus Agave presents a bimodal karyotype with x = 30 (5L, large, +25S, small chromosomes). Bimodality within this genus is generally attributed to allopolyploidy in the ancestral form of Agavoideae. However, alternative mechanisms, such as the preferential accumulation of repetitive elements at the macrochromosomes, could also be important. Aiming to understand the role of repetitive DNA within the bimodal karyotype of Agave, genomic DNA from the commercial hybrid 11648 (2n = 2x = 60, 6.31 Gbp) was sequenced at low coverage, and the repetitive fraction was characterized. In silico analysis showed that ~67.6% of the genome is mainly composed of different LTR retrotransposon lineages and one satellite DNA family (AgSAT171). The satellite DNA localized at the centromeric regions of all chromosomes; however, stronger signals were observed for 20 of the macro- and microchromosomes. All transposable elements showed a dispersed distribution, but not uniform across the length of the chromosomes. Different distribution patterns were observed for different TE lineages, with larger accumulation at the macrochromosomes. The data indicate the differential accumulation of LTR retrotransposon lineages at the macrochromosomes, probably contributing to the bimodality. Nevertheless, the differential accumulation of the satDNA in one group of macro- and microchromosomes possibly reflects the hybrid origin of this Agave accession.

Plastome sequencing of South American Podocarpus species reveals low rearrangement rates despite ancient gondwanan disjunctions

BACKGROUND

Historical reconstructions within Podocarpaceae have provided valuable information to disentangle biogeographic scenarios that begun 65 Mya. However, early molecular phylogenies of Podocarpaceae failed to agree on the intergeneric relationships within the family. The aims of this study were to test whether plastome organization is stable within the genus Podocarpus, to estimate the selective regimes affecting plastome protein-coding genes, and to strengthen our understanding of the phylogenetic relationships and biogeographic history.

METHODS AND RESULTS

We sequenced the plastomes of four South American species from Patagonia, southern Yungas, and Brazilian subtropical forests. We compared their plastomes to those published from Brazil, Africa, New Zealand, and Southeast Asia, along with representatives from other genera within Podocarpaceae as outgroups. The four newly sequenced plastomes ranged in size between 133,791 bp and 133,991 bp. Gene content and order among chloroplasts from South American, African and Asian Podocarpus were conserved and different from the plastome of P. totara, from New Zealand. Most genes showed substitution patterns consistent with a conservative selective regime. Phylogenies inferred from either complete sequences or protein coding regions were mostly congruent with previous studies, but showed earlier branching of P. salignus, P. totara and P. sellowii.

CONCLUSIONS

Highly similar and conserved plastomes of African, South American and Asian species suggest that P. totara plastome should be revised and compared to other species from Oceanic distribution. Furthermore, given such structural conservation, we suggest plastome sequencing is not useful to test whether genomic order can be climatically or geologically structured.

Karyotype asymmetry in Cuscuta L. subgenus Pachystigma reflects its repeat DNA composition

Cuscuta is a cytogenetically diverse genus, with karyotypes varying 18-fold in chromosome number and 127-fold in genome size. Each of its four subgenera also presents particular chromosomal features, such as bimodal karyotypes in Pachystigma. We used low coverage sequencing of the Cuscuta nitida genome (subgenus Pachystigma), as well as chromosome banding and molecular cytogenetics of three subgenus representatives, to understand the origin of bimodal karyotypes. All three species, C. nitida, C. africana (2n = 28) and C. angulata (2n = 30), showed heterochromatic bands mainly in the largest chromosome pairs. Eighteen satellite DNAs were identified in C. nitida genome, two showing similarity to mobile elements. The most abundant were present at the largest pairs, as well as the highly abundant ribosomal DNAs. The most abundant Ty1/Copia and Ty3/Gypsy elements were also highly enriched in the largest pairs, except for the Ty3/Gypsy CRM, which also labelled the pericentromeric regions of the smallest chromosomes. This accumulation of repetitive DNA in the larger pairs indicates that these sequences are largely responsible for the formation of bimodal karyotypes in the subgenus Pachystigma. The repetitive DNA fraction is directly linked to karyotype evolution in Cuscuta.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.

Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA

BACKGROUND AND AIMS

With the advance of high-throughput sequencing, reduced-representation methods such as target capture sequencing (TCS) emerged as cost-efficient ways of gathering genomic information, particularly from coding regions. As the off-target reads from such sequencing are expected to be similar to genome skimming (GS), we assessed the quality of repeat characterization in plant genomes using these data.

METHODS

Repeat composition obtained from TCS datasets of five Rhynchospora (Cyperaceae) species were compared with GS data from the same taxa. In addition, a FISH probe was designed based on the most abundant satellite found in the TCS dataset of Rhynchospora cephalotes. Finally, repeat-based phylogenies of the five Rhynchospora species were constructed based on the GS and TCS datasets and the topologies were compared with a gene-alignment-based phylogenetic tree.

KEY RESULTS

All the major repetitive DNA families were identified in TCS, including repeats that showed abundances as low as 0.01 % in the GS data. Rank correlations between GS and TCS repeat abundances were moderately high (r = 0.58-0.85), increasing after filtering out the targeted loci from the raw TCS reads (r = 0.66-0.92). Repeat data obtained by TCS were also reliable in developing a cytogenetic probe of a new variant of the holocentromeric satellite Tyba. Repeat-based phylogenies from TCS data were congruent with those obtained from GS data and the gene-alignment tree.

CONCLUSIONS

Our results show that off-target TCS reads can be recycled to identify repeats for cyto- and phylogenomic investigations. Given the growing availability of TCS reads, driven by global phylogenomic projects, our strategy represents a way to recycle genomic data and contribute to a better characterization of plant biodiversity.

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

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

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics

The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Large vs small genomes in Passiflora: the influence of the mobilome and the satellitome

While two lineages of retrotransposons were more abundant in larger Passiflora genomes, the satellitome was more diverse and abundant in the smallest genome analysed. Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within Passiflora genus, a tenfold variation in genome size, not attributed to polyploidy, is known. Here, we applied a combined in silico and cytological approach to study the organization and diversification of repetitive elements in three species of this genus representing its known range in genome size variation. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons as the most abundant elements in the three genomes, showing a considerable variation among species. Satellite DNAs (satDNAs) were less representative, but highly diverse between subgenera. Our results clearly confirm that the largest genome species (Passiflora quadrangularis) presents a higher accumulation of repetitive DNA sequences, specially Angela and Tekay elements, making up most of its genome. Passiflora cincinnata, with intermediate genome and from the same subgenus, showed similarity with P. quadrangularis regarding the families of repetitive DNA sequences, but in different proportions. On the other hand, Passiflora organensis, the smallest genome, from a different subgenus, presented greater diversity and the highest proportion of satDNA. Altogether, our data indicates that while large genomes evolved by an accumulation of retrotransposons, the smallest genome known for the genus has evolved by diversification of different repeat types, particularly satDNAs.

Multiple and independent rearrangements revealed by comparative cytogenetic mapping in the dysploid Leptostachyus group (Phaseolus L., Leguminosae)

Polyploidy and dysploidy have been reported as the main events in karyotype evolution of plants. In the genus Phaseolus L. (2n = 22), a small monophyletic group of three species, the Leptostachyus group, presents a dysploid karyotype with 2n = 20. It was shown in Phaseolus leptostachyus that the dysploidy was caused by a nested chromosome fusion (NCF) accompanied by several translocations, suggesting a high rate of karyotype evolution in the group. To verify if this karyotype restructuring was a single event or occurred progressively during the evolution of this group, we analysed P. macvaughii, sister to Phaseolus micranthus + P. leptostachyus. Twenty-four genomic clones of P. vulgaris previously mapped on P. leptostachyus, in addition to the 5S and 35S rDNA probes, were used for fluorescence in situ hybridization. Only a single rearrangement was common to the two species: the nested chromosome fusion (NCF) involving chromosomes 10 and 11. The translocation of chromosome 2 is not the same found in P. leptostachyus, and pericentric inversions in chromosomed 3 and 4 were exclusive of P. macvaughii. The other rearrangements observed in P. leptostachyus were not shared with this species, suggesting that they occurred after the separation of these lineages. The presence of private rearrangements indicates a progressive accumulation of karyotype changes in the Leptostachyus group instead of an instant genome-wide repatterning.

Analysis of the small chromosomal Prionium serratum (Cyperid) demonstrates the importance of reliable methods to differentiate between mono- and holocentricity

For a long time, the Cyperid clade (Thurniceae-Juncaceae-Cyperaceae) was considered a group of species possessing holocentromeres exclusively. The basal phylogenetic position of Prionium serratum (Thunb.) Drège (Thurniceae) within Cyperids makes this species an important specimen to understand the centromere evolution within this clade. In contrast to the expectation, the chromosomal distribution of the centromere-specific histone H3 (CENH3), alpha-tubulin and different centromere-associated post-translational histone modifications (H3S10ph, H3S28ph and H2AT120ph) demonstrate a monocentromeric organisation of P. serratum chromosomes. Analysis of the high-copy repeat composition resulted in the identification of two centromere-localised satellite repeats. Hence, monocentricity was the ancestral condition for the Juncaceae-Cyperaceae-Thurniaceae Cyperid clade, and holocentricity in this clade has independently arisen at least twice after differentiation of the three families, once in Juncaceae and the other one in Cyperaceae. In this context, methods suitable for the identification of holocentromeres are discussed.

How diverse is heterochromatin in the Caesalpinia group? Cytogenomic characterization of Erythrostemon hughesii Gagnon & G.P. Lewis (Leguminosae: Caesalpinioideae)

Cytogenomic characterization of Erythrostemon hughesii reveals a heterogeneity of repeats in its subtelomeric heterochromatin. Comparative analyses with other Caesalpinia group species reveal a significant reduction in the abundance of Ty3-gypsy/Chromovirus Tekay retrotransposons during its evolution. In numerically stable karyotypes, repetitive DNA variability is one of the main causes of genome and chromosome variation and evolution. Species from the Caesalpinia group (Leguminosae) are karyotypically characterized by 2n = 24, with small chromosomes and highly variable CMA⁺ heterochromatin banding patterns that correlate with environmental variables. Erythrostemon hughesii differs from other species of the group examined to date for having subtelomeric CMA⁺ bands; this contrasts with most species in the group which have proximal bands. Here we analyse the repeatome of E. hughesii using genome skimming and chromosomal mapping approaches to characterize the identity of the most abundant repetitive elements and their physical location. The repetitive fraction of E. hughesii comprises 28.73% of the genome. The most abundant elements were retrotransposons (RT) with long terminal repeats (LTR-RT; 9.76%) and satellite DNAs (7.83%). Within the LTR-RTs, the most abundant lineages were: Ty1/copia-Ale (1%), Ty3/gypsy CRM (0.88%) and Ty3/gypsy Athila (0.75%). Using fluorescent in situ hybridization four satellite DNAs and several LTR-RT elements were shown to be present in most subtelomeric CMA⁺ bands. These results highlight how the repeatome in E. hughesii, a species from Oaxaca state in Mexico, is clearly distinct from Northeast Brazilian species of the Caesalpinia group, mainly due to its high diversity of repeats in its subtelomeric heterochromatic bands and low amount of LTR-RT Ty3/gypsy-Tekay elements. Comparative sequence analysis of Tekay elements from different species is congruent with a clade-specific origin of this LTR-RT after the divergence of the Caesalpinia group. We hypothesize that repeat-rich heterochromatin may play a role in leading to faster genomic divergence between individuals, increasing speciation and diversification.

Integration of mandarin (Citrus reticulata) cytogenetic map with its genome sequence

Citrus is an extremely important genus in terms of world fruit production. Despite its economic importance and the small genome sizes of its species (2n = 18, 1C = 430 ± 68 Mbp), entire genomic assemblies have only recently become available for some of its representatives. Together with the previous CMA/DAPI banding and fluorescence in situ hybridization (FISH) in the group, these data are important for understanding the complex relationships between its species and for assisting breeding programs. To anchor genomic data with the cytogenetic map of mandarin (Citrus reticulata), the parental species of several economically important hybrids such as sweet orange and clementine, 18 BAC (bacterial artificial chromosome) clones were used. Eleven clementine BACs were positioned by BAC-FISH, doubling the number of chromosome markers so far available for BAC-FISH in citrus. Additionally, six previously mapped BACs were end-sequenced, allowing, together with one BAC previously sequenced, their assignment to scaffolds and the subsequent integration of chromosomes and the genome assembly. This study therefore established correlations between mandarin scaffolds and chromosomes, allowing further structural genomic and comparative study with the sweet orange genome, as well as insights into the chromosomal evolution of the group.

Super-Resolution Microscopy Reveals Diversity of Plant Centromere Architecture

Centromeres are essential for proper chromosome segregation to the daughter cells during mitosis and meiosis. Chromosomes of most eukaryotes studied so far have regional centromeres that form primary constrictions on metaphase chromosomes. These monocentric chromosomes vary from point centromeres to so-called "meta-polycentromeres", with multiple centromere domains in an extended primary constriction, as identified in Pisum and Lathyrus species. However, in various animal and plant lineages centromeres are distributed along almost the entire chromosome length. Therefore, they are called holocentromeres. In holocentric plants, centromere-specific proteins, at which spindle fibers usually attach, are arranged contiguously (line-like), in clusters along the chromosomes or in bands. Here, we summarize findings of ultrastructural investigations using immunolabeling with centromere-specific antibodies and super-resolution microscopy to demonstrate the structural diversity of plant centromeres. A classification of the different centromere types has been suggested based on the distribution of spindle attachment sites. Based on these findings we discuss the possible evolution and advantages of holocentricity, and potential strategies to segregate holocentric chromosomes correctly.

Evolutionary convergence or homology? Comparative cytogenomics of Caesalpinia group species (Leguminosae) reveals diversification in the pericentromeric heterochromatic composition

We demonstrated by cytogenomic analysis that the proximal heterochromatin of the Northeast Brazilian species of Caesalpinia group is enriched with phylogenetically conserved Ty3/Gypsy-Tekay RT, but diverge in the presence of Ty3/Gypsy-Athila RT and satDNA. The Caesalpinia Group includes 225 species and 27 monophyletic genera of which four occur in Northeastern Brazil: Erythrostemon (1 sp.), Cenostigma (7 spp.), Libidibia (1 sp.), and Paubrasilia (1 sp.). The last three genera are placed in different clades in the Caesalpinia Group phylogeny, and yet they are characterized by having a numerically stable karyotype 2n = 24 (16 M+8A) and GC-rich heterochromatic bands (chromomycin A₃ positive/CMA⁺ bands) in the proximal chromosome regions. To characterize the composition of their heterochromatin and test for the homology of these chromosomal regions, genomic DNA was extracted from Cenostigma microphyllum, Libidibia ferrea, and Paubrasilia echinata, and sequenced at low coverage using the Illumina platform. The genomic repetitive fractions were characterized using a Galaxy/RepeatExplorer-Elixir platform. The most abundant elements of each genome were chromosomally located by fluorescent in situ hybridization (FISH) and compared to the CMA⁺ heterochromatin distribution. The repetitive fraction of the genomes of C. microphyllum, L. ferrea, and P. echinata were estimated to be 41.70%, 38.44%, and 72.51%, respectively. Ty3/Gypsy retrotransposons (RT), specifically the Tekay lineage, were the most abundant repeats in each of the three genomes. FISH mapping revealed species-specific patterns for the Tekay elements in the proximal regions of the chromosomes, co-localized with CMA⁺ bands. Other species-specific patterns were observed, e.g., for the Ty3/Gypsy RT Athila elements which were found in all the proximal heterochromatin of L. ferrea or restricted to the acrocentric chromosomes of C. microphyllum. This Athila labeling co-localized with satellite DNAs (satDNAs). Although the Caesalpinia Group diverged around 55 Mya, our results suggest an ancestral colonization of Tekay RT in the proximal heterochromatin. Thus, the present-day composition of the pericentromeric heterochromatin in these Northeast Brazilian species is a combination of the maintenance of an ancestral Tekay distribution with a species-specific accumulation of other repeats.

Diversity of repetitive sequences within compact genomes of Phaseolus L. beans and allied genera Cajanus L. and Vigna Savi

Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within tribe Phaseoleae (Fabaceae), some genera, such as Phaseolus, Vigna, and Cajanus, show small genome and mostly stable chromosome number. Here, we applied a combined computational and cytological approach to study the organization and diversification of repetitive elements in some species of these genera. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons, especially Ogre and Chromovirus elements, making up most of genomes, other than P. acutifolius and Vigna species. Satellite DNAs (SatDNAs) were less representative, but highly diverse among species, showing a clear phylogenetic relationship. In situ localization revealed preferential location at pericentromeres and centromeres for both types of sequences, suggesting a heterogeneous composition, especially for centromeres. Few elements showed subterminal accumulation. Copy number variation among chromosomes within and among species was observed for all nine identified SatDNAs. Altogether, our data pointed two main elements (Ty3/Gypsy retrotransponsons and SatDNAs) to the diversification on the repetitive landscape in Phaseoleae, with a typical set of repeats in each species. The high turnover of these sequences, however, did not affect total genome size.

Together But Different: The Subgenomes of the Bimodal Eleutherine Karyotypes Are Differentially Organized

Bimodal karyotypes are characterized by the presence of two sets of chromosomes of contrasting size. Eleutherine bulbosa (2n = 12) presents a bimodal karyotype with a large chromosome pair, which has a pericentric inversion in permanent heterozygosity with suppressed recombination, and five pairs of three to four times smaller chromosomes. Aiming to understand whether high copy number sequence composition differs between both chromosome sets, we investigated the repetitive DNA fraction of E. bulbosa and compared it to the chromosomal organization of the related Eleutherine latifolia species, not containing the pericentric inversion. We also compared the repetitive sequence proportions between the heteromorphic large chromosomes of E. bulbosa and between E. bulbosa and E. latifolia to understand the influence of the chromosome inversion on the dynamics of repetitive sequences. The most abundant repetitive families of the genome showed a similar chromosomal distribution in both homologs of the large pair and in both species, apparently not influenced by the species-specific inversions. The repeat families Ebusat1 and Ebusat4 are localized interstitially only on the large chromosome pair, while Ebusat2 is located in the centromeric region of all chromosomes. The four most abundant retrotransposon lineages are accumulated in the large chromosome pair. Replication timing and distribution of epigenetic and transcriptional marks differ between large and small chromosomes. The differential distribution of retroelements appears to be related to the bimodal condition and is not influenced by the nonrecombining chromosome inversions in these species. Thus, the large and small chromosome subgenomes of the bimodal Eleutherine karyotype are differentially organized and probably evolved by repetitive sequences accumulation on the large chromosome set.

Identification of passion fruit (Passiflora edulis) chromosomes using BAC-FISH

Passiflora edulis, the yellow passion fruit, is the main crop from the Passiflora genus, which comprises 525 species with its diversity center in South America. Genetic maps and a BAC (bacterial artificial chromosome) genomic library are available, but the nine chromosome pairs of similar size and morphology (2n = 18) hamper chromosome identification, leading to different proposed karyotypes. Thus, the aim of this study was to establish chromosome-specific markers for the yellow passion fruit using single-copy and repetitive sequences as probes in fluorescent in situ hybridizations (FISH) to allow chromosome identification and future integration with whole genome data. Thirty-six BAC clones harboring genes and three retrotransposons (Ty1-copy, Ty3-gypsy, and LINE) were selected. Twelve BACs exhibited a dispersed pattern similar to that revealed by retroelements, and one exhibited subtelomeric distribution. Twelve clones showed unique signals in terminal or subterminal regions of the chromosomes, allowing their genes to be anchored to six chromosome pairs that can be identified with single-copy markers. The markers developed herein will provide an important tool for genomic and evolutionary studies in the Passiflora genus.

Comparative cytogenetic mapping of chromosomes bearing rDNA sites in Passiflora subgenus Passiflora (Passifloracea)

The genus Passiflora L. comprises around 500 species with different haploid numbers (n = 6, 9, 10 e 12) and variability in the number of ribosomal DNA sites (rDNA). The subgenus Passiflora Juss. is one of the most diverse and includes several species of economic importance, such as Passiflora edulis and Passiflora alata Curtis, as well as Passiflora watsoniana Mast., which is endemic to Northeastern Brazil. These three species present 2n = 18 and two (P. edulis, P. alata) or three (P. watsoniana) pairs of 35S rDNA sites. In P. edulis, these sites are located in chromosomes 7 and 8, identified by fluorescent in situ hybridization (FISH) with chromosome-specific bacterial artificial chromosomes (BACs). With the aim of investigating the conservation of synteny between these species, we used FISH with BACs 173B16 and 164M13, markers for chromosome 7 and 8, respectively, followed by FISH with rDNA. As previously observed in P. edulis, 35S rDNA sites were observed in opposite arms to the BAC markers of chromosome 7 and 8 of both, P. alata and P. watsoniana. In P. alata, the 35S rDNA was observed in the short arm of chromosome 7 and in the long arm of chromosome 8, as in P. edulis. In P. watsoniana, however, the 35S rDNA was located in the long arm of three chromosome pairs, two of them confirmed as chromosomes 7 and 8. The change between short and long arms in the chromosome 8 of P. watsoniana, as well as the presence of an extra rDNA site, is compatible with its phylogenetic distance and indicates karyotypic changes even among species with n = 9 in the Passiflora subgenus.Financial support: CAPES, CNPq and FACEPE

Deciphering the repetitive DNA landscape in Phaseolus beans and Allied genera (Cajanus and Vigna) by a comparative cytogenomic approach

Several bean species (Fabaceae) have been undergoing intensive genomic research, which includes the generation of massive amount of DNA sequence data. We performed a comparative investigation of the repetitive DNA composition in some bean species (Phaseolus vulgaris - Pvu, P. acutifolius - Pac, P. coccineus - Pco, Cajanus cajan - Cca, Vigna angularis - Van and V. nakashimae - Vna). Genomic sequences were obtained from the NCBI public database Cca/SRX031567, Pac/SRX1182015), Pco/SRX1182014, Pvu-andean accession/SRX1182008 and mesoamerican accession/SRX1182013, Vna/SRX758247 and Van/SRX7575509) and run in the Repeat Explorer pipeline using an all-to-all comparison analysis. Physical mapping of the most abundant repetitive elements was performed by FISH using mitotic chromosomes. Overall, the repetitive DNA fraction makes up around 30% of the genomes, except for both Vigna species (9.7% and 14.6%, for Van and Vna, respectively). LTR-retrotransposons represented the major class, with Ogre-Tat elements showing a large variation, followed by LTR-Chromovirus, satellite DNA repeats, LINEs and DNA transposons. The composition of satellite repeats was variable and four new families were identified, together with five previously described. All satellites were genus- or species-specific. Chromosomal mapping of LTR sequences (Chromovirus and Ogre) showed a differential accumulation on the pericentromeres of most chromosomes, particularly on Phaseolus. Among the satellite repeats, five were used for FISH (four newly identified). The three satellite families from Phaseolus showed different species-specific distribution. PacSAT363 was detected in 12 pericentromeres. PcoSAT36 and PcoSAT336 were accumulated at pericentromes, but PcoSAT336 were also located in some chromosome termini co-localized with 45 rDNA, suggesting that this satellite have probably originated from ribosomal sequences. Cajanus cajan and Vigna showed only one satellite repeat each. CcaSAT100 was located in 18 centromeres, while VanSAT174 corresponded to the 45 rDNA in Vna, but appeared also as an independent centromeric satellite repeat in Van. Our analysis indicates a very dynamic landscape on the evolution of repetitive DNA among beans. Retrotransposons appeared as the main component of genomes and showed higher accumulation in the pericentromeric heterochromatin of Phaseolus species. Satellite repeats were mostly centromeric or pericentromeric and frequently species-specific, suggesting a fast evolution of this type of sequence in the bean genomes.

Sucessive karyotipic changes revealed for comparative cytogenetic mapping in the Leptostachyus group, Phaseolus (Fabaceae)

Polyploidy and disploidy have been reported as the main events in karyotype evolution of plants. In the genus Phaseolus L., numerically stable (2n = 22), a small monophyletic group of three species from Mexico and Central America, the Leptostachyus group, presents the dysploid karyotype (2n = 20) and one larger chromosome pair. Many rearrangements, including translocations, inversions and a Nested Chromosome Fusion (NCF), were observed in P. leptostachyus Benth, revealing a high rate of karyotype evolution. The analysis of the first diverging species of this clade, P. macvaughii Delgado, revealed that the NCF event is a synapomorphy of the Leptostachyus clade, resulting in the descending dysploidy. Thus, this study aimed to understand the karyotype evolution of this group, verifying if P. macvaughii shares with P. leptostachyus the same rearrangements in the homeologs to P. vulgaris chromosomes Pv1, Pv2, Pv4 and Pv8. For this, fluorescence in situ hybridizations (FISH) were performed on P. macvaughii mitotic slides using nine sequences of P. vulgaris cloned into bacterial artificial chromosomes (BACs) and previously mapped on P. leptostachyus chromosomes. Sequences selected for Pv1 hybridized collinearly in P. macvaughii, contrasting to the observed in P. leptostachyus, which showed synteny break. Markers from Pv2 showed loss of synteny in both species, due to a translocation. Part of the long arm of chromosome 2 translocated to the short arm of chromosome 1, altering the arm ratio of these chromosomes in P. macvaughii. The chromosomes Pv4 and Pv8 kept their collinearity to P. vulgaris in P. macvaughii, different to the observed in P. leptostachyus. The presence of shared and non-shared rearrangements indicates a progressive accumulation of karyotypic changes in the Leptostachyus group.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.

Are holocentrics doomed to change? Limited chromosome number variation in Rhynchospora Vahl (Cyperaceae)

Karyotype evolution in species with non-localised centromeres (holocentric chromosomes) is usually very dynamic and associated with recurrent fission and fusion (also termed agmatoploidy/symploidy) events. In Rhynchospora (Cyperaceae), one of the most species-rich sedge genera, all analysed species have holocentric chromosomes and their numbers range from 2n = 4 to 2n = 84. Agmatoploidy/symploidy and polyploidy were suggested as the main processes in the reshuffling of Rhynchospora karyotypes, although testing different scenarios of chromosome number evolution in a phylogenetic framework has not been attempted until now. Here, we used maximum likelihood and model-based analyses, in combination with genome size estimation and ribosomal DNA distribution, to understand chromosome evolution in Rhynchospora. Overall, chromosome number variation showed a significant phylogenetic signal and the majority of the lineages maintained a karyotype of 2n = 10 (~48% of the species), the most likely candidate for the ancestral number of the genus. Higher and lower chromosome numbers were restricted to specific clades, whilst polyploidy and/or fusion/fission events were present in specific branches. Variation in genome size and ribosomal DNA site number showed no correlation with ploidy level or chromosome number. Although different mechanisms of karyotype evolution (polyploidy, fusion and fission) seem to be acting in distinct lineages, the degree of chromosome variation and the main mechanisms involved are comparable to those found in some monocentric genera and lower than expected for a holocentric genus.

Evolutionary dynamics of satellite DNA repeats from Phaseolus beans

Common bean (Phaseolus vulgaris) subtelomeres are highly enriched for khipu, the main satellite DNA identified so far in this genome. Here, we comparatively investigate khipu genomic organization in Phaseolus species from different clades. Additionally, we identified and characterized another satellite repeat, named jumper, associated to khipu. A mixture of P. vulgaris khipu clones hybridized in situ confirmed the presence of khipu-like sequences on subterminal chromosome regions in all Phaseolus species, with differences in the number and intensity of signals between species and when species-specific clones were used. Khipu is present as multimers of ∼500 bp and sequence analyses of cloned fragments revealed close relationship among khipu repeats. The new repeat, named jumper, is a 170-bp satellite sequence present in all Phaseolus species and inserted into the nontranscribed spacer (NTS) of the 5S rDNA in the P. vulgaris genome. Nevertheless, jumper was found as a high-copy repeat at subtelomeres and/or pericentromeres in the Phaseolus microcarpus lineage only. Our data argue for khipu as an important subtelomeric satellite DNA in the genus and for a complex satellite repeat composition of P. microcarpus subtelomeres, which also contain jumper. Furthermore, the differential amplification of these repeats in subtelomeres or pericentromeres reinforces the presence of a dynamic satellite DNA library in Phaseolus.

Developmental programmed cell death during asymmetric microsporogenesis in holocentric species of Rhynchospora (Cyperaceae)

Members of the Cyperaceae family exhibit an asymmetric microsporogenesis that results in the degeneration of three out of four meiotic products. Efforts have been made previously to describe the resulting structure, named the pseudomonad, but mechanisms concerning the establishment of cell domains, nuclear development, and programmed cell death are largely unknown. Using the Rhynchospora genus as a model, evidence for cell asymmetry, cytoplasmic isolation, and programmed cell death was obtained by a combination of electron microscopic, cytochemical, immunocytochemical, in situ hybridization, and flow cytometric methods. Degenerative cells were identified at the abaxial region, with the cytoskeleton marking their delimitation from the functional domain after meiosis. After attempting to initiate cell division with an unreplicated genome and abnormal spindle assembly, these cells exhibited a gradual process of cytoplasmic contraction associated with hypermethylation of cytosines and differential loss of DNA. These results indicate that the asymmetric tetrad establishes a functional cell, where one nucleus is preferentially selected to survive. Degenerative haploid cells are then eliminated in a multistep process associated with mitotic disorder, non-random elimination of repetitive DNA, vacuolar cell death, and DNA fragmentation.

Holocentromere identity: from the typical mitotic linear structure to the great plasticity of meiotic holocentromeres

The centromere is the chromosomal site of kinetochore assembly and is responsible for the correct chromosome segregation during mitosis and meiosis in eukaryotes. Contrary to monocentrics, holocentric chromosomes lack a primary constriction, what is attributed to a kinetochore activity along almost the entire chromosome length during mitosis. This extended centromere structure imposes a problem during meiosis, since sister holocentromeres are not co-oriented during first meiotic division. Thus, regardless of the relatively conserved somatic chromosome structure of holocentrics, during meiosis holocentric chromosomes show different adaptations to deal with this condition. Recent findings in holocentrics have brought back the discussion of the great centromere plasticity of eukaryotes, from the typical CENH3-based holocentromeres to CENH3-less holocentric organisms. Here, we summarize recent and former findings about centromere/kinetochore adaptations shown by holocentric organisms during mitosis and meiosis and discuss how these adaptations are related to the type of meiosis found.

Speeding up chromosome evolution in Phaseolus: multiple rearrangements associated with a one-step descending dysploidy

The genus Phaseolus L. has been subject of extensive cytogenetic studies due to its global economic importance. It is considered karyotypically stable, with most of its ca. 75 species having 2n = 22 chromosomes, and only three species (Phaseolus leptostachyus, Phaseolus macvaughii, and Phaseolus micranthus), which form the Leptostachyus clade, having 2n = 20. To test whether a simple chromosomal fusion was the cause of this descending dysploidy, mitotic chromosomes of P. leptostachyus (2n = 20) were comparatively mapped by fluorescent in situ hybridization (FISH) using bacterial artificial chromosomes (BACs) and ribosomal DNA (rDNA) probes. Our results corroborated the conservation of the 5S and 45S rDNA sites on ancestral chromosomes 10 and 6, respectively. The reduction from x = 11 to x = 10 was the result of the insertion of chromosome 10 into the centromeric region of chromosome 11, supporting a nested chromosome fusion (NCF) as the main cause of this dysploidy. Additionally, the terminal region of the long arm of chromosome 6 was translocated to this larger chromosome. Surprisingly, the NCF was accompanied by several additional translocations and inversions previously unknown for the genus, suggesting that the dysploidy may have been associated to a burst of genome reorganization in this otherwise stable, diploid plant genus.

Low genetic diversity and high differentiation among relict populations of the neotropical gymnosperm Podocarpus sellowii (Klotz.) in the Atlantic Forest

Podocarpus sellowii (Podocarpaceae) is one of only a few gymnosperms native to Brazil and the sole species of the genus found in the northeastern region of that country. It has a very restricted distribution in this region, with only three known populations in highland forests (called Brejos de Altitude), which apparently have been isolated from each other since the Pleistocene. Due to this long-term isolation and the fact that these populations have few adult individuals and suffer great anthropogenic pressure, low genetic variability is expected, compromising their long-term viability. The present work assessed the genetic variability and structure of northeastern populations of P. sellowii to investigate the role of Pleistocene glaciations on the genetic relationships between them and to propose strategies for their conservation by analyzing the SSR and ISSR markers of adult and juvenile individuals. Low genetic diversity was found with both markers, associated with a high differentiation of the Brejo de Baturité population in relation to the others-suggesting their isolation at different points in time, probably during the Pleistocene. Actions directed towards increasing the genetic diversity of these populations will be needed, such as planting seedlings with high genetic variability-but the high degrees of differentiation observed between the populations must be taken into account.

Chiasmatic and achiasmatic inverted meiosis of plants with holocentric chromosomes

Meiosis is a specialized cell division in sexually reproducing organisms before gamete formation. Following DNA replication, the canonical sequence in species with monocentric chromosomes is characterized by reductional segregation of homologous chromosomes during the first and equational segregation of sister chromatids during the second meiotic division. Species with holocentric chromosomes employ specific adaptations to ensure regular disjunction during meiosis. Here we present the analysis of two closely related plant species with holocentric chromosomes that display an inversion of the canonical meiotic sequence, with the equational division preceding the reductional. In-depth analysis of the meiotic divisions of Rhynchospora pubera and R. tenuis reveals that during meiosis I sister chromatids are bi-oriented, display amphitelic attachment to the spindle and are subsequently separated. During prophase II, chromatids are connected by thin chromatin threads that appear instrumental for the regular disjunction of homologous non-sister chromatids in meiosis II.

The absence of a defined centromere in organisms with holocentric chromosomes presents particular problems for the control of chromosome segregation during meiosis. Cabral et al. present evidence that two plant species overcome this challenge by inverting the conventional sequence of meiotic divisions.

Begin at the beginning: A BAC-end view of the passion fruit (Passiflora) genome

Background

The passion fruit (Passiflora edulis) is a tropical crop of economic importance both for juice production and consumption as fresh fruit. The juice is also used in concentrate blends that are consumed worldwide. However, very little is known about the genome of the species. Therefore, improving our understanding of passion fruit genomics is essential and to some degree a pre-requisite if its genetic resources are to be used more efficiently. In this study, we have constructed a large-insert BAC library and provided the first view on the structure and content of the passion fruit genome, using BAC-end sequence (BES) data as a major resource.

Results

The library consisted of 82,944 clones and its levels of organellar DNA were very low. The library represents six haploid genome equivalents, and the average insert size was 108 kb. To check its utility for gene isolation, successful macroarray screening experiments were carried out with probes complementary to eight Passiflora gene sequences available in public databases. BACs harbouring those genes were used in fluorescent in situ hybridizations and unique signals were detected for four BACs in three chromosomes (n = 9). Then, we explored 10,000 BES and we identified reads likely to contain repetitive mobile elements (19.6% of all BES), simple sequence repeats and putative proteins, and to estimate the GC content (~42%) of the reads. Around 9.6% of all BES were found to have high levels of similarity to plant genes and ontological terms were assigned to more than half of the sequences analysed (940). The vast majority of the top-hits made by our sequences were to Populus trichocarpa (24.8% of the total occurrences), Theobroma cacao (21.6%), Ricinus communis (14.3%), Vitis vinifera (6.5%) and Prunus persica (3.8%).

Conclusions

We generated the first large-insert library for a member of Passifloraceae. This BAC library provides a new resource for genetic and genomic studies, as well as it represents a valuable tool for future whole genome study. Remarkably, a number of BAC-end pair sequences could be mapped to intervals of the sequenced Arabidopsis thaliana, V. vinifera and P. trichocarpa chromosomes, and putative collinear microsyntenic regions were identified.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-816) contains supplementary material, which is available to authorized users.

The Subtelomeric khipu Satellite Repeat from Phaseolus vulgaris: Lessons Learned from the Genome Analysis of the Andean Genotype G19833

Subtelomeric regions in eukaryotic organisms are known for harboring species-specific tandemly repeated satellite sequences. However, studies on the molecular organization and evolution of subtelomeric repeats are scarce, especially in plants. Khipu is a satellite DNA of 528-bp repeat unit, specific of the Phaseolus genus, with a subtelomeric distribution in common bean, P. vulgaris. To investigate the genomic organization and the evolution of khipu, we performed genome-wide analysis on the complete genome sequence of the common bean genotype G19833. We identified 2,460 khipu units located at most distal ends of the sequenced regions. Khipu units are arranged in discrete blocks of 2-55 copies and are heterogeneously distributed among the different chromosome ends of G19833 (from 0 to 555 khipus units per chromosome arm). Phylogenetically related khipu units are spread between numerous chromosome ends, suggesting frequent exchanges between non-homologous subtelomeres. However, most subclades contain numerous khipu units from only one or few chromosome ends indicating that local duplication is also driving khipu expansion. Unexpectedly, we also identified 81 khipu units located at centromeres. All the centromeric khipu units belong to a single divergent clade also comprised of a few units from several subtelomeres, suggesting that a few sequence exchanges between centromeres and subtelomeres took place in the common bean genome. The divergence and low copy number of these centromeric units from the subtelomeric units could explain why they were not detected by FISH (Fluorescence in situ Hybridization) although it can not be excluded that these centromeric units may have resulted from errors in the pseudomolecule assembly. Altogether our data highlight extensive sequence exchanges in subtelomeres between non-homologous chromosomes in common bean and confirm that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.

Identification and characterization of functional centromeres of the common bean

In higher eukaryotes, centromeres are typically composed of megabase-sized arrays of satellite repeats that evolve rapidly and homogenize within a species' genome. Despite the importance of centromeres, our knowledge is limited to a few model species. We conducted a comprehensive analysis of common bean (Phaseolus vulgaris) centromeric satellite DNA using genomic data, fluorescence in situ hybridization (FISH), immunofluorescence and chromatin immunoprecipitation (ChIP). Two unrelated centromere-specific satellite repeats, CentPv1 and CentPv2, and the common bean centromere-specific histone H3 (PvCENH3) were identified. FISH showed that CentPv1 and CentPv2 are predominantly located at subsets of eight and three centromeres, respectively. Immunofluorescence- and ChIP-based assays demonstrated the functional significance of CentPv1 and CentPv2 at centromeres. Genomic analysis revealed several interesting features of CentPv1 and CentPv2: (i) CentPv1 is organized into an higher-order repeat structure, named Nazca, of 528 bp, whereas CentPv2 is composed of tandemly organized monomers; (ii) CentPv1 and CentPv2 have undergone chromosome-specific homogenization; and (iii) CentPv1 and CentPv2 are not likely to be commingled in the genome. These findings suggest that two distinct sets of centromere sequences have evolved independently within the common bean genome, and provide insight into centromere satellite evolution.

Karyotype stability in the genus Phaseolus evidenced by the comparative mapping of the wild species Phaseolus microcarpus

The genus Phaseolus L. (Fabaceae) is monophyletic and comprises approximately 75 species distributed into two principal clades. The five cultivated species, including the common bean (Phaseolus vulgaris), were placed in clade B. Clade A comprises only wild species, with more limited distribution. In the present work, bacterial artificial chromosomes (BACs) previously mapped in common bean (2n = 22) were used as probes in fluorescent in situ hybridization (FISH) in this comparative study of Phaseolus microcarpus (2n = 22), a species from clade A. We also analyzed the chromomycin A3 (CMA)/4',6-diamidino-2-phenylindole (DAPI) banding pattern and the localization of rDNA and telomeric DNA sites. The single 45S rDNA site from P. microcarpus was mapped to chromosome 6, showing conservation to the P. vulgaris homeolog. Of the two 5S rDNA sites identified in both species, only the site on chromosome 10 appeared conserved. In spite of the phylogenetic distance between the two species, all of the single-copy BACs demonstrated conservation of synteny. However, four collinearity breaks were observed, probably caused by para- and pericentric inversions. Some variation in the repetitive fraction of the genome was also observed. Thus, a broader analysis of the genus confirms that few, rare inversions seem to represent the main karyotype changes during the evolution of this genus.

High macro-collinearity between lima bean (Phaseolus lunatus L.) and the common bean (P. vulgaris L.) as revealed by comparative cytogenetic mapping

Common bean (P. vulgaris) and lima bean (P. lunatus) are the most important crop species from the genus Phaseolus. Both species have the same chromosome number (2n = 22) and previous cytogenetic mapping of BAC clones suggested conserved synteny. Nevertheless, karyotype differences were observed, suggesting structural rearrangements. In this study, comparative cytogenetic maps for chromosomes 3, 4 and 7 were built and the collinearity between the common bean and lima bean chromosomes was investigated. Thirty-two markers (30 BACs and 2 bacteriophages) from P. vulgaris were hybridized in situ on mitotic chromosomes from P. lunatus. Nine BACs revealed a repetitive DNA pattern with pericentromeric distribution and 23 markers showed unique signals. Nine of these markers were mapped on chromosome 3, eight on chromosome 4 and six on chromosome 7. The order and position of all analyzed BACs were similar between the two species, indicating a high level of macro-collinearity. Thus, although few inversions have probably altered centromere position in other chromosomes, the main karyotypic differences were associated with the repetitive DNA fraction.

Karyotype stability in the genus Phaseolus evidenced by the comparative mapping of the wild species Phaseolus microcarpus

The genus Phaseolus L. (Fabaceae) is monophyletic and comprises approximately 75 species distributed into two principal clades. The five cultivated species, including the common bean (Phaseolus vulgaris), were placed in clade B. Clade A comprises only wild species, with more limited distribution. In the present work, bacterial artificial chromosomes (BACs) previously mapped in common bean (2n = 22) were used as probes in fluorescent in situ hybridization (FISH) in this comparative study of Phaseolus microcarpus (2n = 22), a species from clade A. We also analyzed the chromomycin A3 (CMA)/4',6-diamidino-2-phenylindole (DAPI) banding pattern and the localization of rDNA and telomeric DNA sites. The single 45S rDNA site from P. microcarpus was mapped to chromosome 6, showing conservation to the P. vulgaris homeolog. Of the two 5S rDNA sites identified in both species, only the site on chromosome 10 appeared conserved. In spite of the phylogenetic distance between the two species, all of the single-copy BACs demonstrated conservation of synteny. However, four collinearity breaks were observed, probably caused by para- and pericentric inversions. Some variation in the repetitive fraction of the genome was also observed. Thus, a broader analysis of the genus confirms that few, rare inversions seem to represent the main karyotype changes during the evolution of this genus.

Contrasting evolution of a satellite DNA and its ancestral IGS rDNA in Phaseolus (Fabaceae)

CC4 is a satellite DNA from common bean (Phaseolus vulgaris L.) that is similar to its intergenic spacer (IGS) rDNA. CC4 was originally hypothesized to be an old, fast evolving satellite family that has invaded common bean rDNA. To test this hypothesis and contribute to the understanding of IGS-like satellite DNA evolution, we have investigated its distribution in the genus Phaseolus and related species. CC4 was cloned and used as probe for Southern blot and FISH experiments. CC4 was observed as an independent satellite in common bean, forming two to three major and a few minor pericentromeric clusters. In Phaseolus coccineus, CC4 was present in four major clusters, also not co-localized with the 45S rDNA sites. Remarkably, in the less related species of the genus, signals were detected co-localized with the 45S rDNA sites, but co-localization was not observed in the species where CC4 is present as an independent satellite. No signal was detected in species from related genera. Altogether, the data suggest that CC4 has originated from the IGS rDNA in the P. vulgaris-P. coccineus lineage and has evolved slower than the IGS rDNA from this lineage.

Comparative cytogenetic mapping between the lima bean (Phaseolus lunatus L.) and the common bean (P. vulgaris L.)

The common bean (Phaseolus vulgaris) and lima bean (P. lunatus) are among the most important legumes in terms of direct human consumption. The present work establishes a comparative cytogenetic map of P. lunatus, using previously mapped markers from P. vulgaris, in association with analyses of heterochromatin distribution using the fluorochromes chromomycin A3 (CMA) and 4',6-diamidino-2-phenylindole (DAPI) and localization of the 5S and 45S ribosomal DNA (rDNA) probes. Seven BACs selected from different common bean chromosomes demonstrated a repetitive pericentromeric pattern corresponding to the heterochromatic regions revealed by CMA/DAPI and could not be mapped. The subtelomeric repetitive pattern observed for BAC 63H6 in most of the chromosome ends of common bean was not detected in lima bean, indicating lack of conservation of this subtelomeric repeat. All chromosomes could be identified and 16 single-copy clones were mapped. These results showed a significant conservation of synteny between species, although change in centromere position suggested the occurrence of pericentric inversions on chromosomes 2, 9 and 10. The low number of structural rearrangements reflects the karyotypic stability of the genus.