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 holocentromeres influence genome architecture and karyotype evolution.

The centromere represents a single region in most eukaryotic chromosomes. However, several plant and animal lineages assemble holocentromeres along the entire chromosome length. Here, we compare genome organization and evolution as a function of centromere type by assembling chromosome-scale holocentric genomes with repeat-based holocentromeres from three beak-sedge (Rhynchospora pubera, R. breviuscula, and R. tenuis) and their closest monocentric relative, Juncus effusus. We demonstrate that transition to holocentricity affected 3D genome architecture by redefining genomic compartments, while distributing centromere function to thousands of repeat-based centromere units genome-wide. We uncover a complex genome organization in R. pubera that hides its unexpected octoploidy and describe a marked reduction in chromosome number for R. tenuis, which has only two chromosomes. We show that chromosome fusions, facilitated by repeat-based holocentromeres, promoted karyotype evolution and diploidization. Our study thus sheds light on several important aspects of genome architecture and evolution influenced by centromere organization.

Nuclear variations and tapetum polyploidy related to pollen grain development in Passiflora L. (Passifloraceae).

Tapetal cells comprise an anther tissue fundamental to pollen grain development. They are associated with endoreduplication events, which culminate in polyploid and multinucleated cells, high metabolic activity, and different organelle arrangements to support all the development of the pollen grains. Passiflora species present a secretory tapetum, with diversity in the number and size of nuclei. Tapetal cells undergo numerous changes in a short period of development when compared to the plant's life span. To improve our knowledge of tapetum development, tests assessing ploidy levels, anatomy, cytochemistry, transmission electron microscopy, flow cytometry, as well as conventional and molecular cytogenetics were used in Passiflora actinia and P. elegans. The current data show striking differences in nuclear organisation during tapetal cell development, including mono to quadrinucleate cells, and ploidy levels from 2n to 32n. One of the most peculiar features was the atypical behaviour of the endoplasmic reticulum (ER), which accumulated in the cell border, similar to a 'cER', as well as large dictyosomes. This endomembrane configuration may be related to the tapetum nutritional network and secretion of compounds at the end of meiosis. Another atypical feature of the ER was the formation of an invagination to establish 'binucleated' polyploid cells. This membrane projection appears when the nuclei form two lobes, as well as when it organises a nucleoplasmic reticulum. These data demonstrate that there are important ultrastructural changes in tapetal cells, including organelle arrangements, ploidy levels, and nuclear activity, common to P. actinia and P. elegans, but different from the plant model A. thaliana.

Holocentric Karyotype Evolution in Rhynchospora Is Marked by Intense Numerical, Structural, and Genome Size Changes.

Cyperaceae is a family of Monocotyledons comprised of species with holocentric chromosomes that are associated with intense dysploidy and polyploidy events. Within this family the genus Rhynchospora has recently become the focus of several studies that characterize the organization of the holocentric karyotype and genome structures. To broaden our understanding of genome evolution in this genus, representatives of Rhynchospora were studied to contrast chromosome features, C-CMA/DAPI band distribution and genome sizes. Here, we carried out a comparative analysis for 35 taxa of Rhynchospora, and generated new genome size estimates for 20 taxa. The DNA 2C-values varied up to 22-fold, from 2C = 0.51 pg to 11.32 pg, and chromosome numbers ranged from 2n = 4 to 61. At least 37% of our sampling exhibited 2n different from the basic number x = 5, and chromosome rearrangements were also observed. A large variation in C-CMA/DAPI band accumulation and distribution was observed as well. We show that genome variation in Rhynchospora is much larger than previously reported. Phylogenetic analysis showed that most taxa were grouped in clades corresponding to previously described taxonomic sections. Basic chromosome numbers are the same within every section, however, changes appeared in all the clades. Ancestral chromosome number reconstruction revealed n = 5 as the most likely ancestral complements, but n = 10 appears as a new possibility. Chromosome evolution models point to polyploidy as the major driver of chromosome evolution in Rhynchospora, followed by dysploidy. A negative correlation between chromosome size and diploid number open the discussion for holokinetic drive-based genome evolution. This study explores relationships between karyotype differentiation and genome size variation in Rhynchospora, and contrasts it against the phylogeny of this holocentric group.

Capsidiol-related genes are highly expressed in response to Colletotrichum scovillei during Capsicum annuum fruit development stages.

Capsicum annuum is one of the most important horticultural crops worldwide. Anthracnose disease (Colletotrichum spp.) is a major constraint for chili production, causing substantial losses. Capsidiol is a sesquiterpene phytoalexin present in pepper fruits that can enhance plant resistance. The genetic mechanisms involved in capisidiol biosynthesis are still poorly understood. In this study, a 3' RNA sequencing approach was used to develop the transcriptional profile dataset of C. annuum genes in unripe (UF) and ripe fruits (RF) in response to C. scovillei infection. Results showed 4,845 upregulated and 4,720 downregulated genes in UF, and 2,560 upregulated and 1,762 downregulated genes in RF under fungus inoculation. Four capsidiol-related genes were selected for RT-qPCR analysis, two 5-epi-aristolochene synthase (CA12g05030, CA02g09520) and two 5-epi-aristolochene-1,3-dihydroxylase genes (CA12g05070, CA01g05990). CA12g05030 and CA01g05990 genes showed an early response to fungus infection in RF (24 h post-inoculation-HPI), being 68-fold and 53-fold more expressed at 96 HPI, respectively. In UF, all genes showed a late response, especially CA12g05030, which was 700-fold more expressed at 96 HPI compared to control plants. We are proving here the first high-throughput expression dataset of pepper fruits in response to anthracnose disease in order to contribute for future pepper breeding programs.

Are unusual ultrastructural features occurring in the pollen endomembrane system of Cyperaceae and other angiosperms?

Cyperaceae representatives present peculiar microsporogenesis and microgametogenesis, which raises the question of how regular the sedge pollen is. In order to answer this question, preanthesis pollen grains of Cyperaceae individuals were analyzed under different tools such as light and transmission electron microscopy, which included cytochemistry and immunogold procedures. The results showed that maturing pollen in Cyperaceae presents some unusual endomembrane behaviors. Endoplasmic reticulum and dictyosomes are concerned in classic secretion pathways in vegetative cells, and possibly the late breakdown of degenerative microspores. However, cortical and concentric endoplasmic reticulum are also present and are possibly related to other functions aside the biosynthetic pathway. Unconventional secretion of large membrane-bound bodies containing cell wall precursors was also observed and confirmed by immunogold. However, since these bodies most likely receive material from dictyosomes, as observed in silver nitrate reaction, the "unconventional" status of this secretion is debatable. Reports of the literature show that these unusual endomembrane formations are not exclusive of the sedge pollen, but little attention have been given to them so far. This could represent an opportunity to re-examine our understanding on the endomembrane system in pollen cells in general.

Analysis of retrotransposon abundance, diversity and distribution in holocentric Eleocharis (Cyperaceae) genomes.

Background and Aims: Long terminal repeat-retrotransposons (LTR-RTs) comprise a large portion of plant genomes, with massive repeat blocks distributed across the chromosomes. Eleocharis species have holocentric chromosomes, and show a positive correlation between chromosome numbers and the amount of nuclear DNA. To evaluate the role of LTR-RTs in karyotype diversity in members of Eleocharis (subgenus Eleocharis), the occurrence and location of different members of the Copia and Gypsy superfamilies were compared, covering interspecific variations in ploidy levels (considering chromosome numbers), DNA C-values and chromosomal arrangements. Methods: The DNA C-value was estimated by flow cytometry. Genomes of Eleocharis elegans and E. geniculata were partially sequenced using Illumina MiSeq assemblies, which were a source for searching for conserved proteins of LTR-RTs. POL domains were used for recognition, comparing families and for probe production, considering different families of Copia and Gypsy superfamilies. Probes were obtained by PCR and used in fluorescence in situ hybridization (FISH) against chromosomes of seven Eleocharis species. Key Results: A positive correlation between ploidy levels and the amount of nuclear DNA was observed, but with significant variations between samples with the same ploidy levels, associated with repetitive DNA fractions. LTR-RTs were abundant in E. elegans and E. geniculata genomes, with a predominance of Copia Sirevirus and Gypsy Athila/Tat clades. FISH using LTR-RT probes exhibited scattered and clustered signals, but with differences in the chromosomal locations of Copia and Gypsy. The diversity in LTR-RT locations suggests that there is no typical chromosomal distribution pattern for retrotransposons in holocentric chromosomes, except the CRM family with signals distributed along chromatids. Conclusions: These data indicate independent fates for each LTR-RT family, including accumulation between and within chromosomes and genomes. Differential activity and small changes in LTR-RTs suggest a secondary role in nuclear DNA variation, when compared with ploidy changes.

Structure and Distribution of Centromeric Retrotransposons at Diploid and Allotetraploid Coffea Centromeric and Pericentromeric Regions.

Centromeric regions of plants are generally composed of large array of satellites from a specific lineage of Gypsy LTR-retrotransposons, called Centromeric Retrotransposons. Repeated sequences interact with a specific H3 histone, playing a crucial function on kinetochore formation. To study the structure and composition of centromeric regions in the genus Coffea, we annotated and classified Centromeric Retrotransposons sequences from the allotetraploid C. arabica genome and its two diploid ancestors: Coffea canephora and C. eugenioides. Ten distinct CRC (Centromeric Retrotransposons in Coffea) families were found. The sequence mapping and FISH experiments of CRC Reverse Transcriptase domains in C. canephora, C. eugenioides, and C. arabica clearly indicate a strong and specific targeting mainly onto proximal chromosome regions, which can be associated also with heterochromatin. PacBio genome sequence analyses of putative centromeric regions on C. arabica and C. canephora chromosomes showed an exceptional density of one family of CRC elements, and the complete absence of satellite arrays, contrasting with usual structure of plant centromeres. Altogether, our data suggest a specific centromere organization in Coffea, contrasting with other plant genomes.

Centromeric and non-centromeric satellite DNA organisation differs in holocentric Rhynchospora species.

Satellite DNA repeats (or satDNA) are fast-evolving sequences usually associated with condensed heterochromatin. To test whether the chromosomal organisation of centromeric and non-centromeric satDNA differs in species with holocentric chromosomes, we identified and characterised the major satDNA families in the holocentric Cyperaceae species Rhynchospora ciliata (2n = 10), R. globosa (2n = 50) and R. tenuis (2n = 2x = 4 and 2n = 4x = 8). While conserved centromeric repeats (present in R. ciliata and R. tenuis) revealed linear signals at both chromatids, non-centromeric, species-specific satDNAs formed distinct clusters along the chromosomes. Colocalisation of both repeat types resulted in a ladder-like hybridisation pattern at mitotic chromosomes. In interphase, the centromeric satDNA was dispersed while non-centromeric satDNA clustered and partly colocalised to chromocentres. Despite the banding-like hybridisation patterns of the clustered satDNA, the identification of chromosome pairs was impaired due to the irregular hybridisation patterns of the homologues in R. tenuis and R. ciliata. These differences are probably caused by restricted or impaired meiotic recombination as reported for R. tenuis, or alternatively by complex chromosome rearrangements or unequal condensation of homologous metaphase chromosomes. Thus, holocentricity influences the chromosomal organisation leading to differences in the distribution patterns and condensation dynamics of centromeric and non-centromeric satDNA.

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.

Interstitial telomeric sites and Robertsonian translocations in species of Ipheion and Nothoscordum (Amaryllidaceae).

The genera Nothoscordum and Ipheion (Allioideae, Amaryllidaceae) are cytologically characterized by a dysploid series with variable numbers of metacentric and acrocentric chromosomes typical of karyotypes rearranged by Robertsonian translocations (RT). Since they have large chromosomes, low diploid numbers, and possess two telomeric motifs [the vertebrate-type (TTAGGG) n and the Arabidopsis-type (TTTAGGG) n ] they are suitable for investigating the occurrence and possible role of interstitial telomeric sites (ITS) associated with RT. We analyzed the distributions of telomeric sites in 12 species of Nothoscordum and Ipheion and found that both telomeric probes colocalized in all chromosome termini. Cloning and sequencing PCR products obtained using both telomeric primers simultaneously revealed long stretches of (TTAGGG) n and (TTTAGGG) n sequences together with degenerated telomeric sequences. Most acrocentric chromosomes have a 45S rDNA site at the terminal region of the short arms adjacent to the most distal telomeric sites. Telomeric signals were found at all chromosome ends, but ITS were also detected in a few proximal and subterminal regions in some Nothoscordum species. Although RT are common in this group of plants, our findings suggest that proximal positioning of telomeric motifs are not necessarily related to that kind of rearrangement. Rather, transposition of telomeric sequences followed by amplification, could better explain the presence of (TTAGGG) n and (TTTAGGG) n repeats at those sites. Furthermore, a few small interstitial sites found in some Nothoscordum species indicate that dispersion of these sequences was not restricted to the proximal region.

Holocentromeres in Rhynchospora are associated with genome-wide centromere-specific repeat arrays interspersed among euchromatin.

Holocentric chromosomes lack a primary constriction, in contrast to monocentrics. They form kinetochores distributed along almost the entire poleward surface of the chromatids, to which spindle fibers attach. No centromere-specific DNA sequence has been found for any holocentric organism studied so far. It was proposed that centromeric repeats, typical for many monocentric species, could not occur in holocentrics, most likely because of differences in the centromere organization. Here we show that the holokinetic centromeres of the Cyperaceae Rhynchospora pubera are highly enriched by a centromeric histone H3 variant-interacting centromere-specific satellite family designated "Tyba" and by centromeric retrotransposons (i.e., CRRh) occurring as genome-wide interspersed arrays. Centromeric arrays vary in length from 3 to 16 kb and are intermingled with gene-coding sequences and transposable elements. We show that holocentromeres of metaphase chromosomes are composed of multiple centromeric units rather than possessing a diffuse organization, thus favoring the polycentric model. A cell-cycle-dependent shuffling of multiple centromeric units results in the formation of functional (poly)centromeres during mitosis. The genome-wide distribution of centromeric repeat arrays interspersing the euchromatin provides a previously unidentified type of centromeric chromatin organization among eukaryotes. Thus, different types of holocentromeres exist in different species, namely with and without centromeric repetitive sequences.

Transcriptional activity, chromosomal distribution and expression effects of transposable elements in Coffea genomes.

Plant genomes are massively invaded by transposable elements (TEs), many of which are located near host genes and can thus impact gene expression. In flowering plants, TE expression can be activated (de-repressed) under certain stressful conditions, both biotic and abiotic, as well as by genome stress caused by hybridization. In this study, we examined the effects of these stress agents on TE expression in two diploid species of coffee, Coffea canephora and C. eugenioides, and their allotetraploid hybrid C. arabica. We also explored the relationship of TE repression mechanisms to host gene regulation via the effects of exonized TE sequences. Similar to what has been seen for other plants, overall TE expression levels are low in Coffea plant cultivars, consistent with the existence of effective TE repression mechanisms. TE expression patterns are highly dynamic across the species and conditions assayed here are unrelated to their classification at the level of TE class or family. In contrast to previous results, cell culture conditions per se do not lead to the de-repression of TE expression in C. arabica. Results obtained here indicate that differing plant drought stress levels relate strongly to TE repression mechanisms. TEs tend to be expressed at significantly higher levels in non-irrigated samples for the drought tolerant cultivars but in drought sensitive cultivars the opposite pattern was shown with irrigated samples showing significantly higher TE expression. Thus, TE genome repression mechanisms may be finely tuned to the ideal growth and/or regulatory conditions of the specific plant cultivars in which they are active. Analysis of TE expression levels in cell culture conditions underscored the importance of nonsense-mediated mRNA decay (NMD) pathways in the repression of Coffea TEs. These same NMD mechanisms can also regulate plant host gene expression via the repression of genes that bear exonized TE sequences.

Developmental and cytogenetic analyses of pollen sterility in Valeriana scandens L.

Valeriana scandens presents perfect and pistillate flowers, the latter with sterile anthers. The species is composed of two varieties with different ploidy; V. scandens var. scandens (2n = 28) and V. scandens var. candolleana (2n = 56), both of which occur in RS, Brazil. Crosses between these varieties may give rise to hybrids with pollen sterility. In this study, we analyzed the microsporogenesis and microgametogenesis of sterile and fertile anthers, and also investigate whether pollen sterility is caused by an irregular meiotic process. Developmental analysis using light microscopy and scanning electron microscopy showed that sterile anthers develop similarly to fertile anthers until the end of meiosis. After this stage, sterile tetrads do not separate as a consequence of exine fusion between adjacent microspores, which is similar to sterile pollen of Brassica ms-cdl1 mutants. In addition, vacuolated immature pollen grains degenerate after separation. The cytogenetic analysis of the microspore mother cell (MMC) showed that the diploid population of V. scandens var. scandens (2n = 28) has pollen sterility that is not caused by a cytogenetic disturbance. The MMCs analyzed from prophase I to tetrad stage showed a regular meiotic process, indicating the phenotype of V. scandens sterile pollen is a postmeiotic process formed by fusion of exine between opposite microspores.

Molecular and cytogenetic characterization of an AT-rich satellite DNA family in Urvillea chacoensis Hunz. (Paullinieae, Sapindaceae).

Urvillea chacoensis is a climber with 2n = 22 and some terminal AT-rich heterochromatin blocks that differentiate it from other species of the genus. The AT-rich highly repeated satellite DNA was isolated from U. chacoensis by the digestion of total nuclear DNA with HindIII and XbaI and cloned in Escherichia coli. Satellite DNA structure and chromosomal distribution were investigated. DNA sequencing revealed that the repeat length of satDNA ranges between 721 and 728 bp, the percentage of AT-base pairs was about 72-73% and the studied clones showed an identity of 92.5-95.9%. Although this monomer has a tetranucleosomal size, direct imperfect repetitions of ~180 bp subdividing it in four nucleosomal subregions were observed. The results obtained with FISH indicate that this monomer usually appears distributed in the terminal regions of most chromosomes and is associated to heterochromatin blocks observed after DAPI staining. These observations are discussed in relation to the satellite DNA evolution and compared with other features observed in several plant groups.

Chromosomal organization and phylogenetic relationships in Hypochaeris species (Asteraceae) from Brazil

The association of cytogenetic and molecular techniques has contributed to the analysis of chromosome organization and phylogeny in plants. The fluorochrome GC-specific CMA3, fluorescent in situ hybridization (FISH) and RAPD (Random Amplified Polymorphic DNA) markers were used to investigate chromosome structure and genetic relationships in Hypochaeris (Asteraceae). Seven species native to South America, and two species introduced from Europe (H. glabra and Hypochaeris sp) were studied. FISH with rDNA probes identified one or two loci of 18S-5.8S-25S rDNA in the South American Hypochaeris species and one locus in the European species. Only one 5S rDNA locus was seen in all species studied. Blocks of GC-rich heterochromatin (CMA-positive bands) associated to 18S-5.8S-25SrDNA loci were detected in all species investigated. Co-location of 5S rDNA and CMA bands was also observed, except for three South American species and Hypochaeris sp. In two South American species, additional CMA bands not related to rDNA were observed on the long arm of chromosome 2, near to the centromere. Hypochaeris glabra exhibited additional CMA-positive signals distributed at pericentromeric regions, on the short arms of all chromosomes. A total of 122 RAPD markers were used to determine the genetic relationships among species. The level of polymorphism was very high, revealing two genetic groups comprising the South American and the European species, thus supporting a previous hypothesis of monophyly of the South American Hypochaeris species. The coefficients of genetic similarity between European and South American species were 0.35, on average. Polymorphism was also high within the two groups. The genetic associations observed with RAPD markers were consistent with chromosome characteristics. Species carrying similar distribution of 45S rDNA loci and CMA-positive signals were included in the same group revealed by RAPDs. Cytogenetic and molecular data support the view that not only chromosome rearrangements, but also changes in DNA sequence took place during the diversification of the South American Hypochaeris species.

Localization of 45S rDNA and telomeric sites on holocentric chromosomes of Rhynchospora tenuis Link (Cyperaceae)

Rhynchospora tenuis Link (Cyperaceae) is a weed widely distributed in Brazil that presents a small number of holocentric chromosomes (2n = 4) with some autopolyploid populations (2n = 8). The haploid number n = 2 is considered as a derivative of the base number x = 5. 45S rDNA probes and telomeric DNA were hybridized in both chromosome races of R. tenuis, looking for indications of chromosome fusions. The results showed that hybridization sites of the telomeric probe were restricted to end chromosome regions whereas rDNA sites were terminally located. The chromosome race with n = 4 exhibited a doubled number of sites, with similar size and location to the hybridized sequences, confirming its autopolyploid origin. Furthermore, the terminal location of the single 45S rDNA site in the 2n = 4 race suggested that disploid reduction in Rhynchospora, from n = 5 to n = 2, was followed by elimination or reorganization events, keeping the terminal distribution of these sites, as in an others species of the genus

Characterization of diploid, tetraploid and hexaploid Helianthus species by chromosome banding and FISH with 45S rDNA probe.

Comparative karyotype analyses of five diploid, two tetraploid, and three hexaploid species of Helianthus were performed using Feulgen staining, Giemsa C and CMA3 (C-CMA) staining, and FISH with 45S rDNA probe. The karyotypes are composed by a basic number of x = 17 with a predominance of meta- and submetacentric chromosome types. A polyploid series is associated with the basic number. Giemsa C- and C-CMA banding revealed terminal or interstitial heterochromatin according to the species, suggesting the existence of a mechanism that may be acting in the dispersion of heterochromatic segments in Helianthus. The nucleolar organizer regions were located at terminal chromosome positions by FISH with 45S rDNA probe. Diploid species presented four, six, and eight rDNA sites, tetraploid species showed eight sites and hexaploid species presented 12 rDNA sites. Karyomorphological differences include variation in number, size and chromosome morphology, suggesting that rearrangements involving small heterochromatic and rDNA segments played a major role in karyotype evolution.