Characterization of some satellite DNA families in Deschampsia antarctica (Poaceae)

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.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

A Pipeline NanoTRF as a New Tool for De Novo Satellite DNA Identification in the Raw Nanopore Sequencing Reads of Plant Genomes

High-copy tandemly organized repeats (TRs), or satellite DNA, is an important but still enigmatic component of eukaryotic genomes. TRs comprise arrays of multi-copy and highly similar tandem repeats, which makes the elucidation of TRs a very challenging task. Oxford Nanopore sequencing data provide a valuable source of information on TR organization at the single molecule level. However, bioinformatics tools for de novo identification of TRs in raw Nanopore data have not been reported so far. We developed NanoTRF, a new python pipeline for TR repeat identification, characterization and consensus monomer sequence assembly. This new pipeline requires only a raw Nanopore read file from low-depth (<1×) genome sequencing. The program generates an informative html report and figures on TR genome abundance, monomer sequence and monomer length. In addition, NanoTRF performs annotation of transposable elements (TEs) sequences within or near satDNA arrays, and the information can be used to elucidate how TR−TE co-evolve in the genome. Moreover, we validated by FISH that the NanoTRF report is useful for the evaluation of TR chromosome organization—clustered or dispersed. Our findings showed that NanoTRF is a robust method for the de novo identification of satellite repeats in raw Nanopore data without prior read assembly. The obtained sequences can be used in many downstream analyses including genome assembly assistance and gap estimation, chromosome mapping and cytogenetic marker development.

Repeatome Analyses and Satellite DNA Chromosome Patterns in Deschampsia sukatschewii, D. cespitosa, and D. antarctica (Poaceae)

Subpolar and polar ecotypes of Deschampsia sukatschewii (Popl.) Roshev, D. cespitosa (L.) P. Beauv, and D. antarctica E. Desv. are well adapted to stressful environmental conditions, which make them useful model plants for genetic research and breeding. For the first time, the comparative repeatome analyses of subpolar and polar D. sukatschewii, D. cespitosa, and D. antarctica was performed using RepeatExplorer/TAREAN pipelines and FISH-based chromosomal mapping of the identified satellite DNA families (satDNAs). In the studied species, mobile genetic elements of class 1 made up the majority of their repetitive DNA; interspecific variations in the total amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were revealed; 12–18 high confident and 7–9 low confident putative satDNAs were identified. According to BLAST, most D. sukatschewii satDNAs demonstrated sequence similarity with satDNAs of D. antarctica and D. cespitosa indicating their common origin. Chromosomal mapping of 45S rDNA, 5S rDNA, and satDNAs of D. sukatschewii allowed us to construct the species karyograms and detect new molecular chromosome markers important for Deschampsia species. Our findings confirmed that genomes of D. sukatschewii and D. cespitosa were more closely related compared to D. antarctica according to repeatome composition and patterns of satDNA chromosomal distribution.

Genetic divergence of Deschampsia antarctica (Poaceae) population groups in the maritime Antarctic

We used inter-simple sequence repeats and inter-retrotransposon amplified polymorphism markers to assess genetic variation in Deschampsia antarctica populations in the context of its uneven distribution in the northern and central maritime Antarctic. Genetic diversity and population structure and differentiation were assessed in nine populations from geographically isolated population groups of D. antarctica, including the South Shetland Islands, Anvers Island and the Argentine Islands regions. In total, 265 amplified DNA fragments were scored, of which 220 (83.0%) were polymorphic. The total sample showed low genetic diversity (unbiased expected heterozygosity = 0.081 and Shannon diversity index = 0.115) and high population differentiation (molecular variance among populations = 0.659). We also found a trend toward a decrease in genetic diversity and an increase in population differentiation toward the southern edge of the species range. Principal coordinates analysis of polymerase chain reaction data and Bayesian population structure analysis showed three main clusters, which included plants originating from three spatially isolated population groups. The unweighted pair group method with arithmetic mean clustering of populations based on Nei’s genetic distances was mainly in agreement with this pattern. Testing of isolation by distance using the Mantel test demonstrated a significant correlation between genetic and ln-transformed geographical distance (r = 0.703). The data obtained indicate that the geographically isolated D. antarctica populations in the maritime Antarctic might form genetic clusters within the total range.

Chromosomal Differentiation of Deschampsia (Poaceae) Based on Four Satellite DNA Families

Diverse families of satellite DNA (satDNA) were detected in heterochromatin regions of Deschampsia. This kind of repetitive DNA consists of tandem repeat sequences forming big arrays in genomes, and can contribute to lineages differentiation. The differentiation between types of satDNA is related to their sequence identity, the size and number of monomers forming the array, and their chromosomal location. In this work, four families of satDNA (D2, D3, D12, D13), previously isolated by genomic analysis, were studied on chromosomal preparations of 12 species of Deschampsia (D. airiformis, D. antarctica, D. cespitosa, D. cordillerarum, D. elongata, D. kingii, D. laxa, D. mendocina, D. parvula, D. patula, D. venustula, and Deschampsia sp) and one of Deyeuxia (D. eminens). Despite the number of satDNA loci showing interspecific variation, the general distribution pattern of each satDNA family is maintained. The four satDNA families are AT-rich and associated with DAPI + heterochromatin regions. D2, D3, and D12 have mainly subterminal distribution, while D13 is distributed in intercalary regions. Such conservation of satDNA patterns suggests a not random distribution in genomes, where the variation between species is mainly associated with the array size and the loci number. The presence of satDNA in all species studied suggests a low genetic differentiation of sequences. On the other hand, the variation of the distribution pattern of satDNA has no clear association with phylogeny. This may be related to high differential amplification and contraction of sequences between lineages, as explained by the library model.

Cytogenomics of Deschampsia P. Beauv. (Poaceae) Species Based on Sequence Analyses and FISH Mapping of CON/COM Satellite DNA Families

The genus Deschampsia P. Beauv. (Poaceae) involves a group of widespread polymorphic species, and many of them are highly tolerant to stressful environmental conditions. Genome diversity and chromosomal phylogeny within the genus are still insufficiently studied. Satellite DNAs, including CON/COM families, are the main components of the plant repeatome, which contribute to chromosome organization. For the first time, using PCR-based (Polymerase Chain Reaction) techniques and sequential BLAST (Basic Local Alignment Search Tool) and MSA (Multiple Sequence Alignment) analyses, we identified and classified CON/COM repeats in genomes of eleven Deschampsia accessions and three accessions from related genera. High homology of CON/COM sequences were revealed in the studied species though differences in single-nucleotide alteration profiles detected in homologous CON/COM regions indicated that they tended to diverge independently. The performed chromosome mapping of 45S rDNA, 5S rDNA, and CON/COM repeats in six Deschampsia species demonstrated interspecific variability in localization of these cytogenetic markers and facilitated the identification of different chromosomal rearrangements. Based on the obtained data, the studied Deschampsia species were distinguished into karyological groups, and MSA-based schematic trees were built, which could clarify the relationships within the genus. Our findings can be useful for further genetic and phylogenetic studies.

Intense proliferation of rDNA sites and heterochromatic bands in two distantly related Cuscuta species (Convolvulaceae) with very large genomes and symmetric karyotypes

The genome size varies widely among angiosperms but only a few clades present huge variation at a low phylogenetic level. Among diploid species of the genus Cuscuta the genome size increased enormously in at least two independent lineages: in species of subgenus Monogynella and in at least one species (C. indecora) of the subgenus Grammica. Curiously, the independent events lead to similar karyotypes, with 2n = 30 mostly metacentric chromosomes. In this paper we compared the patterns of heterochromatic bands and rDNA sites of C. indecora and C. monogyna, aiming to evaluate the role of these repetitive fractions in these karyotypes. We found out that the large genomes of these species were incremented by a huge number of small heterochromatic CMA+ and DAPI+ bands and 5S and 35 rDNA sites, most of them clearly colocalized with CMA+ bands. Silver nitrate impregnation revealed that the maximum number of nucleoli per nucleus was low in both species, suggesting that some of these sites may be inactive. Noteworthy, the tandem repeats did not generate large bands or sites but rather dozens of small blocks dispersed throughout the chromosomes, apparently contributing to conserve the original karyotype symmetry.

Differential expression after UV-B radiation and characterization of chalcone synthase from the Patagonian hairgrass Deschampsia antarctica

Deschampsiaantarctica inhabits the maritime territory of Antarctica and South Patagonia. It grows under very harsh environmental conditions. The survival of this species in low freezing temperatures and under high levels of UV-B radiation may constitute some of the most remarkable adaptive plant responses and suggests that this plant possesses genes associated with cold and UV tolerance. Frequently, increased levels of flavonoids have been linked to highly UV-B irradiated plants. Studies examining the biosynthesis of flavonoids in D. antarctica may provide clues to its success in this extreme environment. In this study, we characterized the family of genes encoding chalcone synthase, a key enzyme of the flavonoid biosynthetic pathway. DaCHS was cloned, sequenced and characterized by using software tools. CHS contains two domains, the N-terminal domain ranges from amino acid 8 to 231 and the C-terminal domain ranges from amino acid 241 to 391. Sequence analysis of the three family members revealed a high degree of identity after comparison with other monocotyledons such as Oryza sativa L., Zea mays L. and Hordeum vulgare L. According to these results, DaCHS can be grouped together with H. vulgare CHS1 in the same branch. The phylogenetic tree was built using MEGA software and the neighbour join method with 1000 bootstrap replicates. A model of DaCHS was constructed by way of structural tools and key amino acid residues were identified at the active motif site.

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.

Characterization of repeated DNA sequences in genomes of blue-flowered flax

BACKGROUND

Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method.

RESULTS

It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes.

CONCLUSIONS

The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.

The Repetitive DNA Composition in the Natural Pesticide Producer Tanacetum cinerariifolium: Interindividual Variation of Subtelomeric Tandem Repeats

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevir.) Sch. Bip.), a plant species endemic to the east Adriatic coast, is used worldwide for production of the organic insecticide, pyrethrin. Most studies concerning Dalmatian pyrethrum have focused on its morphological and biochemical traits relevant for breeding. However, little is known about the chromosomal evolution and genome organization of this species. Our study aims are to identify, classify, and characterize repetitive DNA in the T. cinerariifolium genome using clustering analysis of a low coverage genomic dataset. Repetitive DNA represents about 71.63% of the genome. T. cinerariifolium exhibits linked 5S and 35S rDNA configuration (L-type). FISH reveals amplification of interstitial telomeric repeats (ITRs) in T. cinerariifolium. Of the three newly identified satellite DNA families, TcSAT1 and TcSAT2 are located subterminally on most of T. cinerariifolium chromosomes, while TcSAT3 family is located intercalary within the longer arm of two chromosome pairs. FISH reveals high levels of polymorphism of the TcSAT1 and TcSAT2 sites by comparative screening of 28 individuals. TcSAT2 is more variable than TcSAT1 regarding the number and position of FISH signals. Altogether, our data highlights the dynamic nature of DNA sequences associated with subtelomeres in T. cinerariifolium and suggests that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.