The repetitive DNA landscape in Avena (Poaceae): chromosome and genome evolution defined by major repeat classes in whole-genome sequence reads

Genome Studies in Amaranthus cruentus L. and A. hypochondriacus L. Based on Repeatomic and Cytogenetic Data

Amaranthus cruentus L. and Amaranthus hypochondriacus L. are valuable and promising food crops for multi-purpose use that are distributed worldwide in temperate, subtropical, and tropical zones. However, their karyotypes and genomic relationships still remain insufficiently studied. For the first time, a comparative repeatome analysis of A. cruentus and A. hypochondriacus was performed based on the available NGS data; bioinformatic analyses using RepeatExplorer/TAREAN pipelines; and chromosome FISH mapping of 45S rDNA, 5S rDNA, and the most abundant satellite DNAs. In the repeatomes of these species, interspecific variations in the amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were detected. In the repeatomes of both species, shared satDNAs with high sequence similarity were identified. The chromosome distribution patterns of four effective molecular markers, 45S rDNA, 5S rDNA, AmC4, and AmC9, allowed us to identify all chromosome pairs in the species karyotypes, construct unique karyograms of A. cruentus and A. hypochondriacus, and confirm the close relationship between their genomes. These results are important for comparative karyotypic studies within the genus Amaranthus. Our findings demonstrated that cytogenomic analyses might provide important data on genomic relationships within Amaranthus and increase knowledge on genome organization in these valuable crops.

Species of the Sections Hedysarum and Multicaulia of the Genus Hedysarum (Fabaceae): Taxonomy, Distribution, Chromosomes, Genomes, and Phylogeny

The genus Hedysarum L. (Fabaceae) includes about 200 species of annual and perennial herbs distributed in Asia, Europe, North Africa, and North America. Many species of this genus are valuable medicinal, melliferous, and forage resources. In this review, we consider the taxonomic history of the genus Hedysarum, the chromosomal organization of the species from the sections Hedysarum and Multicaulia, as well as phylogenetic relationships between these sections. According to morphological, genetic, and phylogenetic data, the genus Hedysarum is divided into three main sections: Hedysarum (= syn. Gamotion), Multicaulia, and Stracheya. In species of this genus, two basic chromosome numbers, x = 7 (section Hedysarum) and x = 8 (sections Multicaulia and Stracheya), were determined. The systematic positions of some species within the sections are still uncertain due to their morphological similarities. The patterns of distribution of molecular chromosomal markers (45S rDNA, 5S rDNA, and different satellite DNAs) in karyotypes of various Hedysarum species made it possible to determine their ploidy status and also specify genomic relationships within the sections Hedysarum and Multicaulia. Recent molecular phylogenetic studies clarified significantly the taxonomy and evolutionary development of the genus Hedysarum.

Genome Variability in Artificial Allopolyploid Hybrids of Avena sativa L. and Avena macrostachya Balansa ex Coss. et Durieu Based on Marker Sequences of Satellite DNA and the ITS1-5.8S rDNA Region

Artificial hybrids between cultivated Avena species and wild Avena macrostachya that possess genes for resistance to biotic and abiotic stresses can be important for oat breeding. For the first time, a comprehensive study of genomes of artificial fertile hybrids Avena sativa × Avena macrostachya and their parental species was carried out based on the chromosome FISH mapping of satellite DNA sequences (satDNAs) and also analysis of intragenomic polymorphism in the 18S-ITS1-5.8S rDNA region, using NGS data. Chromosome distribution patterns of marker satDNAs allowed us to identify all chromosomes in the studied karyotypes, determine their subgenomic affiliation, and detect several chromosome rearrangements. Based on the obtained cytogenomic data, we revealed differences between two A. macrostachya subgenomes and demonstrated that only one of them was inherited in the studied octoploid hybrids. Ribotype analyses showed that the second major ribotype of A. macrostachya was species-specific and was not represented in rDNA pools of the octoploids, which could be related to the allopolyploid origin of this species. Our results indicate that the use of marker satDNAs in cytogenomic studies can provide important data on genomic relationships within Avena allopolyploid species and hybrids, and also expand the potential for interspecific crosses for breeding.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

Satellite DNAs-From Localized to Highly Dispersed Genome Components

According to the established classical view, satellite DNAs are defined as abundant non-coding DNA sequences repeated in tandem that build long arrays located in heterochromatin. Advances in sequencing methodologies and development of specialized bioinformatics tools enabled defining a collection of all repetitive DNAs and satellite DNAs in a genome, the repeatome and the satellitome, respectively, as well as their reliable annotation on sequenced genomes. Supported by various non-model species included in recent studies, the patterns of satellite DNAs and satellitomes as a whole showed much more diversity and complexity than initially thought. Differences are not only in number and abundance of satellite DNAs but also in their distribution across the genome, array length, interspersion patterns, association with transposable elements, localization in heterochromatin and/or in euchromatin. In this review, we compare characteristic organizational features of satellite DNAs and satellitomes across different animal and plant species in order to summarize organizational forms and evolutionary processes that may lead to satellitomes' diversity and revisit some basic notions regarding repetitive DNA landscapes in genomes.

Complex polyploid and hybrid species in an apomictic and sexual tropical forage grass group: genomic composition and evolution in Urochloa (Brachiaria) species

BACKGROUND AND AIMS

Diploid and polyploid Urochloa (including Brachiaria, Panicum and Megathyrsus species) C4 tropical forage grasses originating from Africa are important for food security and the environment, often being planted in marginal lands worldwide. We aimed to characterize the nature of their genomes, the repetitive DNA and the genome composition of polyploids, leading to a model of the evolutionary pathways within the group including many apomictic species.

METHODS

Some 362 forage grass accessions from international germplasm collections were studied, and ploidy was determined using an optimized flow cytometry method. Whole-genome survey sequencing and molecular cytogenetic analysis were used to identify chromosomes and genomes in Urochloa accessions belonging to the 'brizantha' and 'humidicola' agamic complexes and U. maxima.

KEY RESULTS

Genome structures are complex and variable, with multiple ploidies and genome compositions within the species, and no clear geographical patterns. Sequence analysis of nine diploid and polyploid accessions enabled identification of abundant genome-specific repetitive DNA motifs. In situ hybridization with a combination of repetitive DNA and genomic DNA probes identified evolutionary divergence and allowed us to discriminate the different genomes present in polyploids.

CONCLUSIONS

We suggest a new coherent nomenclature for the genomes present. We develop a model of evolution at the whole-genome level in diploid and polyploid accessions showing processes of grass evolution. We support the retention of narrow species concepts for Urochloa brizantha, U. decumbens and U. ruziziensis, and do not consider diploids and polyploids of single species as cytotypes. The results and model will be valuable in making rational choices of parents for new hybrids, assist in use of the germplasm for breeding and selection of Urochloa with improved sustainability and agronomic potential, and assist in measuring and conserving biodiversity in grasslands.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes

BACKGROUND AND AIMS

Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution.

METHODS

Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs.

KEY RESULTS

Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements.

CONCLUSIONS

The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

Polyploidy: its consequences and enabling role in plant diversification and evolution

BACKGROUND

Most, if not all, green plant (Virdiplantae) species including angiosperms and ferns are polyploids themselves or have ancient polyploid or whole genome duplication signatures in their genomes. Polyploids are not only restricted to our major crop species such as wheat, maize, potato and the brassicas, but also occur frequently in wild species and natural habitats. Polyploidy has thus been viewed as a major driver in evolution, and its influence on genome and chromosome evolution has been at the centre of many investigations. Mechanistic models of the newly structured genomes are being developed that incorporate aspects of sequence evolution or turnover (low-copy genes and regulatory sequences, as well as repetitive DNAs), modification of gene functions, the re-establishment of control of genes with multiple copies, and often meiotic chromosome pairing, recombination and restoration of fertility.

SCOPE

World-wide interest in how green plants have evolved under different conditions - whether in small, isolated populations, or globally - suggests that gaining further insight into the contribution of polyploidy to plant speciation and adaptation to environmental changes is greatly needed. Forward-looking research and modelling, based on cytogenetics, expression studies, and genomics or genome sequencing analyses, discussed in this Special Issue of the Annals of Botany, consider how new polyploids behave and the pathways available for genome evolution. They address fundamental questions about the advantages and disadvantages of polyploidy, the consequences for evolution and speciation, and applied questions regarding the spread of polyploids in the environment and challenges in breeding and exploitation of wild relatives through introgression or resynthesis of polyploids.

CONCLUSION

Chromosome number, genome size, repetitive DNA sequences, genes and regulatory sequences and their expression evolve following polyploidy - generating diversity and possible novel traits and enabling species diversification. There is the potential for ever more polyploids in natural, managed and disturbed environments under changing climates and new stresses.

Non-B-form DNA tends to form in centromeric regions and has undergone changes in polyploid oat subgenomes

Centromeres are the specialized regions of the chromosomes that direct faithful chromosome segregation during cell division. Despite their functional conservation, centromeres display features of rapidly evolving DNA and wide evolutionary diversity in size and organization. Previous work found that the noncanonical B-form DNA structures are abundant in the centromeres of several eukaryotic species with a possible implication for centromere specification. Thus far, systematic studies into the organization and function of non-B-form DNA in plants remain scarce. Here, we applied the oat system to investigate the role of non-B-form DNA in centromeres. We conducted chromatin immunoprecipitation sequencing using an antibody to the centromere-specific histone H3 variant (CENH3); this accurately positioned oat centromeres with different ploidy levels and identified a series of centromere-specific sequences including minisatellites and retrotransposons. To define genetic characteristics of oat centromeres, we surveyed the repeat sequences and found that dyad symmetries were abundant in oat centromeres and were predicted to form non-B-DNA structures in vivo. These structures including bent DNA, slipped DNA, Z-DNA, G-quadruplexes, and R-loops were prone to form within CENH3-binding regions. Dynamic conformational changes of predicted non-B-DNA occurred during the evolution from diploid to tetraploid to hexaploid oat. Furthermore, we applied the single-molecule technique of AFM and DNA:RNA immunoprecipitation with deep sequencing to validate R-loop enrichment in oat centromeres. Centromeric retrotransposons exhibited strong associations with R-loop formation. Taken together, our study elucidates the fundamental character of non-B-form DNA in the oat genome and reveals its potential role in centromeres.

Plant Cytogenetics: From Chromosomes to Cytogenomics

Chromosomes have been studied since the late nineteenth century in the disciplines of cytology and cytogenetics. Analyzing their numbers, features, and dynamics has been tightly linked to the technical development of preparation methods, microscopes, and chemicals to stain them, with latest continuing developments described in this volume. At the end of the twentieth and beginning of the twenty-first centuries, DNA technology, genome sequencing, and bioinformatics have revolutionized how we see, use, and analyze chromosomes. The advent of in situ hybridization has shaped our understanding of genome organization and behavior by linking molecular sequence information with the physical location along chromosomes and genomes. Microscopy is the best technique to accurately determine chromosome number. Many features of chromosomes in interphase nuclei or pairing and disjunction at meiosis, involving physical movement of chromosomes, can only be studied by microscopy. In situ hybridization is the method of choice to characterize the abundance and chromosomal distribution of repetitive sequences that make up the majority of most plant genomes. These most variable components of a genome are found to be species- and occasionally chromosome-specific and give information about evolution and phylogeny. Multicolor fluorescence hybridization and large pools of BAC or synthetic probes can paint chromosomes and we can follow them through evolution involving hybridization, polyploidization, and rearrangements, important at a time when structural variations in the genome are being increasingly recognized. This volume discusses many of the most recent developments in the field of plant cytogenetics and gives carefully compiled protocols and useful resources.

Oat chromosome and genome evolution defined by widespread terminal intergenomic translocations in polyploids

Structural chromosome rearrangements involving translocations, fusions and fissions lead to evolutionary variation between species and potentially reproductive isolation and variation in gene expression. While the wheats (Triticeae, Poaceae) and oats (Aveneae) all maintain a basic chromosome number of x=7, genomes of oats show frequent intergenomic translocations, in contrast to wheats where these translocations are relatively rare. We aimed to show genome structural diversity and genome relationships in tetraploid, hexaploid and octoploid Avena species and amphiploids, establishing patterns of intergenomic translocations across different oat taxa using fluorescence in situ hybridization (FISH) with four well-characterized repetitive DNA sequences: pAs120, AF226603, Ast-R171 and Ast-T116. In A. agadiriana (2n=4x=28), the selected probes hybridized to all chromosomes indicating that this species originated from one (autotetraploid) or closely related ancestors with the same genomes. Hexaploid amphiploids were confirmed as having the genomic composition AACCDD, while octoploid amphiploids showed three different genome compositions: AACCCCDD, AAAACCDD or AABBCCDD. The A, B, C, and D genomes of oats differ significantly in their involvement in non-centromeric, intercalary translocations. There was a predominance of distal intergenomic translocations from the C- into the D-genome chromosomes. Translocations from A- to C-, or D- to C-genome chromosomes were less frequent, proving that at least some of the translocations in oat polyploids are non-reciprocal. Rare translocations from A- to D-, D- to A- and C- to B-genome chromosomes were also visualized. The fundamental research has implications for exploiting genomic biodiversity in oat breeding through introgression from wild species potentially with contrasting chromosomal structures and hence deleterious segmental duplications or large deletions in amphiploid parental lines.

Aegilops crassa Boiss. repeatome characterized using low-coverage NGS as a source of new FISH markers: Application in phylogenetic studies of the Triticeae

Aegilops crassa Boiss. is polyploid grass species that grows in the eastern part of the Fertile Crescent, Afghanistan, and Middle Asia. It consists of tetraploid (4x) and hexaploid (6x) cytotypes (2n = 4x = 28, D¹D (Abdolmalaki et al., 2019) X^crX^cr and 2n = 6x = 42, D¹D (Abdolmalaki et al., 2019) X^crX^crD²D (Adams and Wendel, 2005), respectively) that are similar morphologically. Although many Aegilops species were used in wheat breeding, the genetic potential of Ae. crassa has not yet been exploited due to its uncertain origin and significant genome modifications. Tetraploid Ae. crassa is thought to be the oldest polyploid Aegilops species, the subgenomes of which still retain some features of its ancient diploid progenitors. The D¹ and D² subgenomes of Ae. crassa were contributed by Aegilops tauschii (2n = 2x = 14, DD), while the X^cr subgenome donor is still unknown. Owing to its ancient origin, Ae. crassa can serve as model for studying genome evolution. Despite this, Ae. crassa is poorly studied genetically and no genome sequences were available for this species. We performed low-coverage genome sequencing of 4x and 6x cytotypes of Ae. crassa, and four Ae. tauschii accessions belonging to different subspecies; diploid wheatgrass Thinopyrum bessarabicum (J^b genome), which is phylogenetically close to D (sub)genome species, was taken as an outgroup. Subsequent data analysis using the pipeline RepeatExplorer2 allowed us to characterize the repeatomes of these species and identify several satellite sequences. Some of these sequences are novel, while others are found to be homologous to already known satellite sequences of Triticeae species. The copy number of satellite repeats in genomes of different species and their subgenome (D¹ or X^cr) affinity in Ae. crassa were assessed by means of comparative bioinformatic analysis combined with quantitative PCR (qPCR). Fluorescence in situ hybridization (FISH) was performed to map newly identified satellite repeats on chromosomes of common wheat, Triticum aestivum, 4x and 6x Ae. crassa, Ae. tauschii, and Th. bessarabicum. The new FISH markers can be used in phylogenetic analyses of the Triticeae for chromosome identification and the assessment of their subgenome affinities and for evaluation of genome/chromosome constitution of wide hybrids or polyploid species.

The repetitive DNA sequence landscape and DNA methylation in chromosomes of an apomictic tropical forage grass, Cenchrus ciliaris

Cenchrus ciliaris is an apomictic, allotetraploid pasture grass widely distributed in the tropical and subtropical regions of Africa and Asia. In this study, we aimed to investigate the genomic organization and characterize some of the repetitive DNA sequences in this species. Due to the apomictic propagation, various aneuploid genotypes are found, and here, we analyzed a 2n = 4x + 3 = 39 accession. The physical mapping of Ty1-copia and Ty3-gypsy retroelements through fluorescence in situ hybridization with a global assessment of 5-methylcytosine DNA methylation through immunostaining revealed the genome-wide distribution pattern of retroelements and their association with DNA methylation. Approximately one-third of Ty1-copia sites overlapped or spanned centromeric DAPI-positive heterochromatin, while the centromeric regions and arms of some chromosomes were labeled with Ty3-gypsy. Most of the retroelement sites overlapped with 5-methylcytosine signals, except for some Ty3-gypsy on the arms of chromosomes, which did not overlap with anti-5-mC signals. Universal retrotransposon probes did not distinguish genomes of C. ciliaris showing signals in pericentromeric regions of all 39 chromosomes, unlike highly abundant repetitive DNA motifs found in survey genome sequences of C. ciliaris using graph-based clustering. The probes developed from RepeatExplorer clusters gave strong in situ hybridization signals, mostly in pericentromeric regions of about half of the chromosomes, and we suggested that they differentiate the two ancestral genomes in the allotetraploid C. ciliaris, likely having different repeat sequence variants amplified before the genomes came together in the tetraploid.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Integration of Genomic and Cytogenetic Data on Tandem DNAs for Analyzing the Genome Diversity Within the Genus Hedysarum L. (Fabaceae)

The section Multicaulia is the largest clade in the genus Hedysarum L. (Fabaceae). Representatives of the sect. Multicaulia are valuable plants used for medicinal and fodder purposes. The taxonomy and phylogeny of the sect. Multicaulia are still ambiguous. To clarify the species relationships within sect. Multicaulia, we, for the first time, explored repeatomes of H. grandiflorum Pall., H. zundukii Peschkova, and H. dahuricum Turcz. using next-generation sequencing technologies and a subsequent bioinformatic analysis by RepeatExplorer/TAREAN pipelines. The comparative repeatome analysis showed that mobile elements made up 20-24% (Class I) and about 2-2.5% (Class II) of their repetitive DNAs. The amount of ribosomal DNA varied from 1 to 2.6%, and the content of satellite DNA ranged from 2.7 to 5.1%. For each species, five high confident putative tandem DNA repeats and 5-10 low confident putative DNA repeats were identified. According to BLAST, these repeats demonstrated high sequence similarity within the studied species. FISH-based mapping of 35S rDNA, 5S rDNA, and satDNAs made it possible to detect new effective molecular chromosome markers for Hedysarum species and construct the species karyograms. Comparison of the patterns of satDNA localization on chromosomes of the studied species allowed us to assess genome diversity within the sect. Multicaulia. In all studied species, we revealed intra- and interspecific variabilities in patterns of the chromosomal distribution of molecular chromosome markers. In H. gmelinii Ledeb. and H. setigerum Turcz. ex Fisch. et Meyer, similar subgenomes were detected, which confirmed the polyploid status of their genomes. Our findings demonstrated a close genomic relationship among six studied species indicating their common origin and confirmed the taxonomic status of H. setigerum as a subspecies of H. gmelinii as well as the validity of combining the sect. Multicaulia and Subacaulia into one sect. Multicaulia.

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.

New Insights into the Genomic Structure of Avena L.: Comparison of the Divergence of A-Genome and One C-Genome Oat Species

We used next-generation sequencing analysis of the 3′-part of 18S rDNA, ITS1, and a 5′-part of the 5.8S rDNA region to understand genetic variation among seven diploid A-genome Avena species. We used 4−49 accessions per species that represented the As genome (A. atlantica, A. hirtula, and wiestii), Ac genome (A. canariensis), Ad genome (A. damascena), Al genome (A. longiglumis), and Ap genome (A. prostrata). We also took into our analysis one C-genome species, A. clauda, which previously was found to be related to A-genome species. The sequences of 169 accessions revealed 156 haplotypes of which seven haplotypes were shared by two to five species. We found 16 ribotypes that consisted of a unique sequence with a characteristic pattern of single nucleotide polymorphisms and deletions. The number of ribotypes per species varied from one in A. longiglumis to four in A. wiestii. Although most ribotypes were species-specific, we found two ribotypes shared by three species (one for A. damascena, A. hirtula, and A. wiestii, and the second for A. longiglumis, A. atlantica, and A. wiestii), and a third ribotype shared between A. atlantica and A. wiestii. A characteristic feature of the A. clauda ribotype, a diploid C-genome species, is that two different families of ribotypes have been found in this species. Some of these ribotypes are characteristic of Cc-genome species, whereas others are closely related to As-genome ribotypes. This means that A. clauda can be a hybrid between As- and C-genome oats.

Draft Sequencing Crested Wheatgrass Chromosomes Identified Evolutionary Structural Changes and Genes and Facilitated the Development of SSR Markers

Crested wheatgrass (Agropyron cristatum), a wild relative of wheat, is an attractive source of genes and alleles for their improvement. Its wider use is hampered by limited knowledge of its complex genome. In this work, individual chromosomes were purified by flow sorting, and DNA shotgun sequencing was performed. The annotation of chromosome-specific sequences characterized the DNA-repeat content and led to the identification of genic sequences. Among them, genic sequences homologous to genes conferring plant disease resistance and involved in plant tolerance to biotic and abiotic stress were identified. Genes belonging to the important groups for breeders involved in different functional categories were found. The analysis of the DNA-repeat content identified a new LTR element, Agrocen, which is enriched in centromeric regions. The colocalization of the element with the centromeric histone H3 variant CENH3 suggested its functional role in the grass centromere. Finally, 159 polymorphic simple-sequence-repeat (SSR) markers were identified, with 72 of them being chromosome- or chromosome-arm-specific, 16 mapping to more than one chromosome, and 71 mapping to all the Agropyron chromosomes. The markers were used to characterize orthologous relationships between A. cristatum and common wheat that will facilitate the introgression breeding of wheat using A. cristatum.

Cytogenetic evidence supports Avena insularis being closely related to hexaploid oats

Cytogenetic observations, phylogenetic studies and genome analysis using high-density genetic markers have suggested a tetraploid Avena species carrying the C and D genomes (formerly C and A) to be the donor of all hexaploid oats (AACCDD). However, controversy surrounds which of the three extant CCDD tetraploid species—A. insularis, A. magna and A. murphyi—is most closely related to hexaploid oats. The present work describes a comparative karyotype analysis of these three CCDD tetraploid species and two hexaploid species, A. sativa and A. byzantina. This involved the use of FISH with six simple sequence repeats (SSRs) with the motifs CT, AAC, AAG, ACG, ATC and ACT, two repeated ribosomal sequences, and C genome-specific repetitive DNA. The hybridization pattern of A. insularis with oligonucleotide (AC)₁₀ was also determined and compared with those previously published for A. sativa and A. byzantina. Significant differences in the 5S sites and SSR hybridization patterns of A. murphyi compared to the other CCDD species rule out its being directly involved in the origin of the hexaploids. In contrast, the repetitive and SSR hybridization patterns shown by the D genome chromosomes, and by most of the C genome chromosomes of A. magna and A. insularis, can be equated with the corresponding chromosomes of the hexaploids. Several chromosome hybridization signals seen for A. insularis, but not for A. magna, were shared with the hexaploid oats species, especially with A. byzantina. These diagnostic signals add weight to the idea that the extant A. insularis, or a direct ancestor of it, is the most closely related progenitor of hexaploid oats. The similarity of the chromosome hybridization patterns of the hexaploids and CCDD tetraploids was taken as being indicative of homology. A common chromosome nomenclature for CCDD species based on that of the hexaploid species is proposed.

The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Magnolia officinalis, a representative tall aromatic tree of the Magnoliaceae family, is a medicinal plant that is widely used in diverse industries from medicine to cosmetics. We report a chromosome-scale draft genome of M. officinalis, in which ∼99.66% of the sequences were anchored onto 19 chromosomes with the scaffold N50 of 76.62 Mb. We found that a high proportion of repetitive sequences was a common feature of three Magnoliaceae with known genomic data. Magnoliids were a sister clade to eudicots-monocots, which provided more support for understanding the phylogenetic position among angiosperms. An ancient duplication event occurred in the genome of M. officinalis and was shared with Lauraceae. Based on RNA-seq analysis, we identified several key enzyme-coding gene families associated with the biosynthesis of lignans in the genome. The construction of the M. officinalis genome sequence will serve as a reference for further studies of Magnolia, as well as other Magnoliaceae.

Cytogenetic events in the endosperm of amphiploid Avena magna × A. longiglumis

This study analysed cytogenetic events occurring in the syncytial endosperm of the Avena magna H. C. Murphy & Terrell × Avena longiglumis Durieu amphiploid, which is a product of two wild species having different genomes. Selection through the elimination of chromosomes and their fragments, including those translocated, decreased the level of ploidy in the endosperm below the expected 3n, leading to the modal number close to 2n. During intergenomic translocations, fragments of the heterochromatin-rich C-genome were transferred to the D and Al genomes. Terminal and non-reciprocal exchanges dominated, whereas other types of translocations, including microexchanges, were less common. Using two probes and by counterstaining with DAPI, the A. longiglumis and the rare exchanges between the D and Al genomes were detected by GISH. The large discontinuity in the probe labelling in the C chromosomes demonstrated inequality in the distribution of repetitive sequences along the chromosome and probable intragenomic rearrangements. In the nucleus, the spatial arrangement of genomes was non-random and showed a sectorial-concentric pattern, which can vary during the cell cycle, especially in the less stable tissue like the hybrid endosperm.

The first draft genome of feather grasses using SMRT sequencing and its implications in molecular studies of Stipa

The Eurasian plant Stipa capillata is the most widespread species within feather grasses. Many taxa of the genus are dominants in steppe plant communities and can be used for their classification and in studies related to climate change. Moreover, some species are of economic importance mainly as fodder plants and can be used for soil remediation processes. Although large-scale molecular data has begun to appear, there is still no complete or draft genome for any Stipa species. Thus, here we present a single-molecule long-read sequencing dataset generated using the Pacific Biosciences Sequel System. A draft genome of about 1004 Mb was obtained with a contig N50 length of 351 kb. Importantly, here we report 81,224 annotated protein-coding genes, present 77,614 perfect and 58 unique imperfect SSRs, reveal the putative allopolyploid nature of S. capillata, investigate the evolutionary history of the genus, demonstrate structural heteroplasmy of the chloroplast genome and announce for the first time the mitochondrial genome in Stipa. The assembled nuclear, mitochondrial and chloroplast genomes provide a significant source of genetic data for further works on phylogeny, hybridisation and population studies within Stipa and the grass family Poaceae.

Comparative Repeat Profiling of Two Closely Related Conifers (Larix decidua and Larix kaempferi) Reveals High Genome Similarity With Only Few Fast-Evolving Satellite DNAs

In eukaryotic genomes, cycles of repeat expansion and removal lead to large-scale genomic changes and propel organisms forward in evolution. However, in conifers, active repeat removal is thought to be limited, leading to expansions of their genomes, mostly exceeding 10 giga base pairs. As a result, conifer genomes are largely littered with fragmented and decayed repeats. Here, we aim to investigate how the repeat landscapes of two related conifers have diverged, given the conifers' accumulative genome evolution mode. For this, we applied low-coverage sequencing and read clustering to the genomes of European and Japanese larch, Larix decidua (Lamb.) Carrière and Larix kaempferi (Mill.), that arose from a common ancestor, but are now geographically isolated. We found that both Larix species harbored largely similar repeat landscapes, especially regarding the transposable element content. To pin down possible genomic changes, we focused on the repeat class with the fastest sequence turnover: satellite DNAs (satDNAs). Using comparative bioinformatics, Southern, and fluorescent in situ hybridization, we reveal the satDNAs' organizational patterns, their abundances, and chromosomal locations. Four out of the five identified satDNAs are widespread in the Larix genus, with two even present in the more distantly related Pseudotsuga and Abies genera. Unexpectedly, the EulaSat3 family was restricted to L. decidua and absent from L. kaempferi, indicating its evolutionarily young age. Taken together, our results exemplify how the accumulative genome evolution of conifers may limit the overall divergence of repeats after speciation, producing only few repeat-induced genomic novelties.

Position Validation of the Dwarfing Gene Dw6 in Oat (Avena sativa L.) and Its Correlated Effects on Agronomic Traits

An F6 : 8 recombinant inbred line (RIL) population derived from the cross between WAOAT2132 (Dw6) and Caracas along with the two parents were used to evaluate the genetic effects of Dw6 dwarfing gene on plant height and other agronomic traits in oat (Avena sativa L.) across three environments, and develop closely linked markers for marker-assisted selection (MAS) for Dw6. The two parents differed in all investigated agronomic traits except for the number of whorls. The RIL lines showed a bimodal distribution for plant height in all three tested environments, supporting the height of this population was controlled by a single gene. Dw6 significantly reduced plant height (37.66∼44.29%) and panicle length (13.99∼22.10%) but without compromising the coleoptile length which was often positively associated with the reduced stature caused by dwarfing genes. Dw6 has also strong negative effects on hundred kernel weight (14.00∼29.55%), and kernel length (4.21∼9.47%), whereas the effects of Dw6 on the kernel width were not uniform across three environments. By contrast, lines with Dw6 produced more productive tillers (10.11∼10.53%) than lines without Dw6. All these together suggested the potential yield penalty associated with Dw6 might be partially due to the decrease of kernel weight which is attributed largely to the reduction of kernel length. Eighty-one simple sequence repeat (SSR) primer pairs from chromosome 6D were tested, five of them were polymorphic in two parents and in two contrasting bulks, confirming the 6D location of Dw6. By using the five polymorphic markers, Dw6 was mapped to an interval of 1.0 cM flanked by markers SSR83 and SSR120. Caution should be applied in using this information since maker order conflicts were observed. The close linkages of these two markers to Dw6 were further validated in a range of oat lines. The newly developed markers will provide a solid basis for future efforts both in the identification of Dw6 in oat germplasm and in the determination of the nature of the gene through positional cloning.

A universal karyotypic system for hexaploid and diploid Avena species brings oat cytogenetics into the genomics era

BACKGROUND

The identification of chromosomes among Avena species have been studied by C-banding and in situ hybridization. However, the complicated results from several cytogenetic nomenclatures for identifying oat chromosomes are often contradictory. A universal karyotyping nomenclature system for precise chromosome identification and comparative evolutionary studies would be essential for genus Avena based on the recently released genome sequences of hexaploid and diploid Avena species.

RESULTS

Tandem repetitive sequences were predicted and physically located on chromosomal regions of the released Avena sativa OT3098 genome assembly v1. Eight new oligonucleotide (oligo) probes for sequential fluorescence in situ hybridization (FISH) were designed and then applied for chromosome karyotyping on mitotic metaphase spreads of A. brevis, A. nuda, A. wiestii, A. ventricosa, A. fatua, and A. sativa species. We established a high-resolution standard karyotype of A. sativa based on the distinct FISH signals of multiple oligo probes. FISH painting with bulked oligos, based on wheat-barley collinear regions, was used to validate the linkage group assignment for individual A. sativa chromosomes. We integrated our new Oligo-FISH based karyotype system with earlier karyotype nomenclatures through sequential C-banding and FISH methods, then subsequently determined the precise breakage points of some chromosome translocations in A. sativa.

CONCLUSIONS

This new universal chromosome identification system will be a powerful tool for describing the genetic diversity, chromosomal rearrangements and evolutionary relationships among Avena species by comparative cytogenetic and genomic approaches.

Genome-Wide Analysis of Transposable Elements and Satellite DNAs in Spinacia Species to Shed Light on Their Roles in Sex Chromosome Evolution

Sex chromosome evolution has mostly been studied in species with heteromorphic sex chromosomes. The Spinacia genus serves as an ideal model for investigating evolutionary mechanisms underlying the transition from homomorphic to heteromorphic sex chromosomes. Among evolutionary factors, repetitive sequences play multiple roles in sex chromosome evolution while their forces have not been fully explored in Spinacia species. Here, we identified major repetitive sequence classes in male and female genomes of Spinacia species and their ancestral relative sugar beet to elucidate the evolutionary processes of sex chromosome evolution using next-generation sequencing (NGS) data. Comparative analysis revealed that the repeat elements of Spinacia species are considerably higher than of sugar beet, especially the Ty3/Gypsy and Ty1/Copia retrotransposons. The long terminal repeat retroelements (LTR) Angela, Athila, and Ogre may be accounted for the higher proportion of repeats in the spinach genome. Comparison of the repeats proportion between female and male genomes of three Spinacia species indicated the different representation in Spinacia tetrandra samples but not in the S. oleracea or S. turkestanica samples. From these results, we speculated that emergence of repetitive DNA sequences may correlate the formation of sex chromosome and the transition from homomorphic sex chromosomes to heteromorphic sex chromosomes as heteromorphic sex chromosomes exclusively existed in Spinacia tetrandra. Three novel sugar beet-specific satellites were identified and confirmed by fluorescence in situ hybridization (FISH); six out of eight new spinach-specific satellites were mapped to the short arm of sex chromosomes. A total of 141 copies of SolSat01-171-s were found in the sex determination region (SDR). Thus, the accumulation of satellite DNA on the short arm of chromosome 1 may be involved in the sex chromosome evolution in Spinacia species. Our study provides a fundamental resource for understanding repeat sequences in Spinacia species and their roles in sex chromosome evolution.

New evidence confirming the CD genomic constitutions of the tetraploid Avena species in the section Pachycarpa Baum

The tetraploid Avena species in the section Pachycarpa Baum, including A. insularis, A. maroccana, and A. murphyi, are thought to be involved in the evolution of hexaploid oats; however, their genome designations are still being debated. Repetitive DNA sequences play an important role in genome structuring and evolution, so understanding the chromosomal organization and distribution of these sequences in Avena species could provide valuable information concerning genome evolution in this genus. In this study, the chromosomal organizations and distributions of six repetitive DNA sequences (including three SSR motifs (TTC, AAC, CAG), one 5S rRNA gene fragment, and two oat A and C genome specific repeats) were investigated using non-denaturing fluorescence in situ hybridization (ND-FISH) in the three tetraploid species mentioned above and in two hexaploid oat species. Preferential distribution of the SSRs in centromeric regions was seen in the A and D genomes, whereas few signals were detected in the C genomes. Some intergenomic translocations were observed in the tetraploids; such translocations were also detected between the C and D genomes in the hexaploids. These results provide robust evidence for the presence of the D genome in all three tetraploids, strongly suggesting that the genomic constitution of these species is DC and not AC, as had been thought previously.

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.

Comparative chloroplast genome analyses of Avena: insights into evolutionary dynamics and phylogeny

BACKGROUND

Oat (Avena sativa L.) is a recognized health-food, and the contributions of its different candidate A-genome progenitor species remain inconclusive. Here, we report chloroplast genome sequences of eleven Avena species, to examine the plastome evolutionary dynamics and analyze phylogenetic relationships between oat and its congeneric wild related species.

RESULTS

The chloroplast genomes of eleven Avena species (size range of 135,889-135,998 bp) share quadripartite structure, comprising of a large single copy (LSC; 80,014-80,132 bp), a small single copy (SSC; 12,575-12,679 bp) and a pair of inverted repeats (IRs; 21,603-21,614 bp). The plastomes contain 131 genes including 84 protein-coding genes, eight ribosomal RNAs and 39 transfer RNAs. The nucleotide sequence diversities (Pi values) range from 0.0036 (rps19) to 0.0093 (rpl32) for ten most polymorphic genes and from 0.0084 (psbH-petB) to 0.0240 (petG-trnW-CCA) for ten most polymorphic intergenic regions. Gene selective pressure analysis shows that all protein-coding genes have been under purifying selection. The adjacent position relationships between tandem repeats, insertions/deletions and single nucleotide polymorphisms support the evolutionary importance of tandem repeats in causing plastome mutations in Avena. Phylogenomic analyses, based on the complete plastome sequences and the LSC intermolecular recombination sequences, support the monophyly of Avena with two clades in the genus.

CONCLUSIONS

Diversification of Avena plastomes is explained by the presence of highly diverse genes and intergenic regions, LSC intermolecular recombination, and the co-occurrence of tandem repeat and indels or single nucleotide polymorphisms. The study demonstrates that the A-genome diploid-polyploid lineage maintains two subclades derived from different maternal ancestors, with A. longiglumis as the first diverging species in clade I. These genome resources will be helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Avena and family Poaceae levels, and are potentially useful to exploit plastome variation in making hybrids for plant breeding.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Evolutionary Dynamics of the Repetitive DNA in the Karyotypes of Pipa carvalhoi and Xenopus tropicalis (Anura, Pipidae)

The large amphibian genomes contain numerous repetitive DNA components that have played an important role in the karyotypic diversification of this vertebrate group. Hypotheses based on the presumable primitive karyotype (2n = 20) of the anurans of the family Pipidae suggest that they have evolved principally through intrachromosomal rearrangements. Pipa is the only South American pipid, while all the other genera are found in Africa. The divergence of the South American lineages from the African ones occurred at least 136 million years ago and is thought to have had a strong biogeographic component. Here, we tested the potential of the repetitive DNA to enable a better understanding of the differentiation of the karyotype among the family Pipidae and to expand our capacity to interpret the chromosomal evolution in this frog family. Our results indicate a long history of conservation in the chromosome bearing the H3 histone locus, corroborating inferences on the chromosomal homologies between the species in pairs 6, 8, and 9. The chromosomal distribution of the microsatellite motifs also provides useful markers for comparative genomics at the chromosome level between Pipa carvalhoi and Xenopus tropicalis, contributing new insights into the evolution of the karyotypes of these species. We detected similar patterns in the distribution and abundance of the microsatellite arrangements, which reflect the shared organization in the terminal/subterminal region of the chromosomes between these two species. By contrast, the microsatellite probes detected a differential arrangement of the repetitive DNA among the chromosomes of the two species, allowing longitudinal differentiation of pairs that are identical in size and morphology, such as pairs 1, 2, 4, and 5. We also found evidence of the distinctive composition of the repetitive motifs of the centromeric region between the species analyzed in the present study, with a clear enrichment of the (CA) and (GA) microsatellite motifs in P. carvalhoi. Finally, microsatellite enrichment in the pericentromeric region of chromosome pairs 6, 8, and 9 in the P. carvalhoi karyotype, together with interstitial telomeric sequences (ITS), validate the hypothesis that pericentromeric inversions occurred during the chromosomal evolution of P. carvalhoi and reinforce the role of the repetitive DNA in the remodeling of the karyotype architecture of the Pipidae.

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co-evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn-over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3-6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.

Repetitive DNA landscape in essential A and supernumerary B chromosomes of Festuca pratensis Huds

Here, we characterized the basic properties of repetitive sequences in essential A and supernumerary B chromosomes of Festuca pratensis Huds. This was performed by comparative analysis of low-pass Illumina sequence reads of B chromosome lacking (-B) and B chromosome containing (+B) individuals of F. pratensis. 61% of the nuclear genome is composed of repetitive sequences. 43.1% of the genome are transposons of which DNA transposons and retrotransposons made up 2.3% and 40.8%, respectively. LTR retrotransposons are the most abundant mobile elements and contribute to 40.7% of the genome and divided into Ty3-gypsy and Ty1-copia super families with 32.97% and 7.78% of the genome, respectively. Eighteen different satellite repeats were identified making up 3.9% of the genome. Five satellite repeats were used as cytological markers for chromosome identification and genome analysis in the genus Festuca. Four satellite repeats were identified on B chromosomes among which Fp-Sat48 and Fp-Sat253 were specific to the B chromosome of F. pratensis.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

De novo identification of satellite DNAs in the sequenced genomes of Drosophila virilis and D. americana using the RepeatExplorer and TAREAN pipelines

Satellite DNAs are among the most abundant repetitive DNAs found in eukaryote genomes, where they participate in a variety of biological roles, from being components of important chromosome structures to gene regulation. Experimental methodologies used before the genomic era were insufficient, too laborious and time-consuming to recover the collection of all satDNAs from a genome. Today, the availability of whole sequenced genomes combined with the development of specific bioinformatic tools are expected to foster the identification of virtually all the "satellitome" of a particular species. While whole genome assemblies are important to obtain a global view of genome organization, most of them are incomplete and lack repetitive regions. We applied short-read sequencing and similarity clustering in order to perform a de novo identification of the most abundant satellite families in two Drosophila species from the virilis group: Drosophila virilis and D. americana, using the Tandem Repeat Analyzer (TAREAN) and RepeatExplorer pipelines. These species were chosen because they have been used as models to understand satDNA biology since the early 70's. We combined the computational approach with data from the literature and chromosome mapping to obtain an overview of the major tandem repeat sequences of these species. The fact that all of the abundant tandem repeats (TRs) we detected were previously identified in the literature allowed us to evaluate the efficiency of TAREAN in correctly identifying true satDNAs. Our results indicate that raw sequencing reads can be efficiently used to detect satDNAs, but that abundant tandem repeats present in dispersed arrays or associated with transposable elements are frequent false positives. We demonstrate that TAREAN with its parent method RepeatExplorer may be used as resources to detect tandem repeats associated with transposable elements and also to reveal families of dispersed tandem repeats.