Genome size variation in the genus Avena

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.

Differential expression and global analysis of miR156/SQUAMOSA promoter binding-like proteins (SPL) module in oat

SQUAMOSA promoter binding-like proteins (SPLs) are important transcription factors that influence growth phase transition and reproduction in plants. SPLs are targeted by miR156 but the SPL/miR156 module is completely unknown in oat. We identified 28 oat SPL genes (AsSPLs) distributed across all 21 oat chromosomes except for 4C and 6D. The oat- SPL gene family represented six of eight SPL phylogenetic groups, with no AsSPLs in groups 3 and 7. A novel oat miR156 (AsmiR156) family with 21 precursors divided into 7 groups was characterized. A total of 16 AsSPLs were found to be targeted by AsmiR156. Intriguingly, AsSPL3s showed high transcript abundance during early inflorescence (GS-54), as compared to the lower abundance of AsmiR156, indicating their role in reproductive development. Unravelling the SPL/miR156 regulatory hub and alterations in expression patterns of AsSPLs could provide an essential toolbox for genetic improvement in the cultivated oat.

The comparison of polymorphism among Avena species revealed by retrotransposon-based DNA markers and soluble carbohydrates in seeds

Here, we compared the polymorphism among 13 Avena species revealed by the iPBS markers and soluble carbohydrate profiles in seeds. The application of seven iPBS markers generated 83 bands, out of which 20.5% were polymorphic. No species-specific bands were scored. Shannon's information index (I) and expected heterozygosity (He) revealed low genetic diversity, with the highest values observed for A. nuda (I = 0.099; He = 0.068). UPGMA clustering of studied Avena accessions and PCoA results showed that the polyploidy level is the main grouping criterion. High-resolution gas chromatography revealed that the studied Avena accessions share the same composition of soluble carbohydrates, but significant differences in the content of total (5.30-22.38 mg g-1 of dry weight) and particular sugars among studied samples were observed. Sucrose appeared as the most abundant sugar (mean 61.52% of total soluble carbohydrates), followed by raffinose family oligosaccharides (31.23%), myo-inositol and its galactosides (6.16%), and monosaccharides (1.09%). The pattern of interspecific variation in soluble carbohydrates, showed by PCA, was convergent to that revealed by iPBS markers. Thus, both methods appeared as a source of valuable data useful in the characterization of Avena resources or in the discussion on the evolution of this genus.

Genome Size Variation across a Cypriot Fabeae Tribe Germplasm Collection

DNA content is an important trait linked to the evolutionary routes of taxa and often connected to speciation. In the present study, we studied C-values variation across the Cypriot Fabeae gene pool. Several hundred plants (Vicia spp., Lens spp., Pisum spp.) were sampled across Cyprus. Accurate estimates were established by flow cytometry and propidium iodine staining for 155 discrete populations/accessions. A ten-fold variation was detected across lineages with 1C DNA content varying from 1.584 pg for V. cretica (ARI02420) to 13.983 pg for V. faba (ARI00187). In general, flow cytometry was precise for the characterization of species, even though there were instances of genome overlapping across taxa. Most analyses in the current work refer to species that have not been characterized before by flow cytometry (or any other DNA content estimation method). Still, a correlation to C-values previously reported in Kew Plant DNA C-values database was attempted. A high degree of correlation except for V. dalmatica was established. The evaluation of genome size trait in relation with the Fabeae phylogeny, revealed that Pisum and Lens genera were rather homogenous, but an astonishing fluctuation was shown for Vicia spp. Moreover, it was established that genome up- or down-scaling was not directly linked to speciation drivers. The genomic size measurements presented here could deliver extra quality control for the identification and characterization of taxa in germplasm collections, particularly in cases where species share morphological characters.

Population structure analysis and genome-wide association study of a hexaploid oat landrace and cultivar collection

INTRODUCTION

Oat (Avena sativa L.) is an important cereal crop grown worldwide for grain and forage, owing to its high adaptability to diverse environments. However, the genetic and genomics research of oat is lagging behind that of other staple cereal crops.

METHODS

In this study, a collection of 288 oat lines originating worldwide was evaluated using 2,213 single nucleotide polymorphism (SNP) markers obtained from an oat iSelect 6K-beadchip array to study its genetic diversity, population structure, and linkage disequilibrium (LD) as well as the genotype-phenotype association for hullessness and lemma color.

RESULTS

The average gene diversity and polymorphic information content (PIC) were 0.324 and 0.262, respectively. The first three principal components (PCs) accounted for 30.33% of the genetic variation, indicating that the population structure of this panel of oat lines was stronger than that reported in most previous studies. In addition, accessions could be classified into two subpopulations using a Bayesian clustering approach, and the clustering pattern of accessions was closely associated with their region of origin. Additionally, evaluation of LD decay using 2,143 mapped markers revealed that the intrachromosomal whole-genome LD decayed rapidly to a critical r² value of 0.156 for marker pairs separated by a genetic distance of 1.41 cM. Genome-wide association study (GWAS) detected six significant associations with the hullessness trait. Four of these six markers were located on the Mrg21 linkage group between 194.0 and 205.7 cM, while the other two significant markers mapped to Mrg05 and Mrg09. Three significant SNPs, showing strong association with lemma color, were located on linkage groups Mrg17, Mrg18, and Mrg20.

DISCUSSION

Our results discerned relevant patterns of genetic diversity, population structure, and LD among members of a worldwide collection of oat landraces and cultivars proposed to be 'typical' of the Qinghai-Tibetan Plateau. These results have important implications for further studies on association mapping and practical breeding in high-altitude oat.

Changes in physio-biochemical parameters and expression of metallothioneins in Avena sativa L. in response to drought

Drought is one of the major threats to food security. Among several mechanisms involved in plant stress tolerance, one protein family-the plant metallothioneins (MTs)-shows great promise for enhancing drought resistance. Plant metallothioneins in oat (Avena sativa L.) have not yet been deeply analysed, and the literature lacks a comprehensive study of the whole family of plant MTs in response to drought. In this study, we showed that the number and nature of cis-elements linked with stress response in promoters of AsMTs1-3 differed depending on the MT type. Drought stress in oat plants caused an increase in the expression of AsMT2 and AsMT3 and a decrease in the expression of AsMT1 compared to well-watered plants. Moreover, the low values of relative water content, water use efficiency, net photosynthesis (P_N), transpiration (E), stomatal conductance (g_s), chlorophyll a, and carotenoid were accompanied by high levels of electrolyte leakage, internal CO₂ concentration (C_i) and abscisic acid content, and high activity of antioxidants enzymes in plants under drought stress. The present study puts forward the idea that AsMTs are crucial for oat response to drought stress not only by regulating antioxidant activity but also by changing the plant water regime and photosynthesis. Our results support the hypothesis that structural differences among types of plant MTs reflect their diversified physiological roles.

Combining genome size and pollen morphology data to study species relationships in the genus Daucus (Apiaceae)

BACKGROUND

The genus Daucus (Apiaceae) comprises about 40 wild species and the cultivated carrot, a crop of great economic and nutritional importance. The rich genetic diversity of wild Daucus species makes them a valuable gene pool for carrot improvement breeding programs. Therefore, it is essential to have good knowledge of the genome structure and relationships among wild Daucus species. To broaden such knowledge, in this research, the nuclear DNA content for 14 Daucus accessions and four closely related species was estimated by flow cytometry and their pollen morphology was analyzed by light and scanning electron microscopy (SEM).

RESULTS

The flow cytometric analysis showed a 3.2-fold variation in the mean 2C values among Daucus taxa, ranging from 0.999 (D. carota subsp. sativus) to 3.228 pg (D. littoralis). Among the outgroup species, the mean 2C values were 1.775-2.882 pg. The pollen grains of Daucus were tricolporate, mainly prolate or perprolate (rarely) in shape, and mainly medium or small (rarely) in size (21.19-40.38 µm), whereas the outgroup species had tricolporate, perprolate-shaped, and medium-sized (26.01-49.86 µm) pollen grains. In the studied taxa, SEM analysis revealed that exine ornamentation was striate, rugulate, perforate, or the ornamentation pattern was mixed. At the time of shedding, all pollen grains were three-celled, as evidenced by DAPI staining. We also found high positive correlations between the length of the polar axis (P) and the length of the equatorial diameter (E) of pollen grains, as well as between P and P/E. However, when comparing cytogenetic information with palynological data, no significant correlations were observed.

CONCLUSIONS

This study complements the information on the nuclear DNA content in Daucus and provides comprehensive knowledge of the pollen morphology of its taxa. These findings may be important in elucidating the taxonomic relationships among Daucus species and can help in the correct identification of gene bank accessions. In a broader view, they could also be meaningful for the interpretation of evolutionary trends in the genus.

Reference genome assemblies reveal the origin and evolution of allohexaploid oat

Common oat (Avena sativa) is an important cereal crop serving as a valuable source of forage and human food. Although reference genomes of many important crops have been generated, such work in oat has lagged behind, primarily owing to its large, repeat-rich polyploid genome. Here, using Oxford Nanopore ultralong sequencing and Hi-C technologies, we have generated a reference-quality genome assembly of hulless common oat, comprising 21 pseudomolecules with a total length of 10.76 Gb and contig N50 of 75.27 Mb. We also produced genome assemblies for diploid and tetraploid Avena ancestors, which enabled the identification of oat subgenomes and provided insights into oat chromosomal evolution. The origin of hexaploid oat is inferred from whole-genome sequencing, chloroplast genomes and transcriptome assemblies of different Avena species. These findings and the high-quality reference genomes presented here will facilitate the full use of crop genetic resources to accelerate oat improvement.

Analysis of oat seed transcriptome with regards to proteins involved in celiac disease

Oat (Avena sativa L.) is considered to be a healthy food. In contrast to other grain crops, oat is high in protein, lipids, dietary fiber, antioxidants, and uniquely in avenanthramides. The question of whether it can also be consumed by people suffering from celiac disease is still unresolved. The main aim of this study was to extract and sequence genes for potentially harmful avenins, globulins, and α-amylase/trypsin inhibitors in six oat varieties and to establish their variability using PacBio sequencing technology of enriched libraries. The results were compared with sequences of the genes already present in databases. In total, 21 avenin, 75 globulin, and 25 α-amylase/trypsin inhibitor genes were identified and mapped in the hexaploid oat chromosomes. In all of the three gene families, only marginal sequence differences were found between the oat varieties within the individual genes. Avenin epitopes were found in all four types of avenin genes occurring in all oat varieties tested within this study. However, the number of avenin genes was nearly four times lower than of globulin genes and, on the protein level, formed only 10% of storage proteins. Therefore, the question of whether oat is safe to celiac disease people is a question of boundary values.

Toward the development of Ac/Ds transposon-mediated gene tagging system for functional genomics in oat (Avena sativa L.)

Cultivated oat (Avena sativa L.) is an important cereal grown worldwide due to its multifunctional uses for animal feed and human food. Oat has lagged behind other cereals in the genetic and genomic studies attributed to its large and complex genomes. Transposon-based genome characterization has been utilized successfully for identifying and determining gene function in large genome cereals. To develop gene tagging and gene-editing resources for oat, maize Activator (Ac) and Dissociation (Ds) transposons were introduced into the oat genome using the biolistic delivery system. A total of 2035 oat calli were bombarded and twenty-four independent, stable transgenic events were obtained. Transformation frequencies were up to 19.0%, and 1.9% for bialaphos and hygromycin selection, respectively. Re-mobilization of the non-autonomous Ds element, by introducing Ac transposase source, led to a transposition frequency up to 16.8%. The properties of ten unique flanking sequences have been characterized to reveal the Ds-tagged sites in the oat genome. Genes at Ds insertion sites showed homology to gibberellin 20-oxidase 3, (1,3;1,4)-beta-D-glucan synthase, and aspartate kinase. This Ac/Ds transposon-based gene tagging system could facilitate and expedite functional genomic studies in oat.

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.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

Discovery and Chromosomal Location a Highly Effective Oat Crown Rust Resistance Gene Pc50-5

Crown rust, caused by Puccinia coronata f. sp. avenae, is one of the most destructive fungal diseases of oat worldwide. Growing disease-resistant oat cultivars is the preferred method of preventing the spread of rust and potential epidemics. The object of the study was Pc50-5, a race-specific seedling crown rust resistant gene, highly effective at all growth stages, selected from the differential line Pc50 (Avena sterilis L. CW 486-1 × Pendek). A comparison of crown rust reaction as well as an allelism test showed the distinctiveness of Pc50-5, whereas the proportions of phenotypes in segregating populations derived from a cross with two crown rust-susceptible Polish oat cultivars, Kasztan × Pc50-5 and Bingo × Pc50-5, confirmed monogenic inheritance of the gene, indicating its usefulness in oat breeding programs. Effective gene introgression depends on reliable gene identification in the early stages of plant development; thus, the aim of the study was to develop molecular markers that are tightly linked to Pc50-5. Segregating populations of Kasztan × Pc50-5 were genotyped using DArTseq technology based on next-generation Illumina short-read sequencing. Markers associated with Pc50-5 were located on chromosome 6A of the current version of the oat reference genome (Avena sativa OT3098 v2, PepsiCo) in the region between 434,234,214 and 440,149,046 bp and subsequently converted to PCR-based SCAR (sequence-characterized amplified region) markers. Furthermore, 5426978_SCAR and 24031809_SCAR co-segregated with the Pc50-5 resistance allele and were mapped to the partial linkage group at 0.6 and 4.0 cM, respectively. The co-dominant 58163643_SCAR marker was the best diagnostic and it was located closest to Pc50-5 at 0.1 cM. The newly discovered, very strong monogenic crown rust resistance may be useful for oat improvement. DArTseq sequences converted into specific PCR markers will be a valuable tool for marker-assisted selection in breeding programs.

Subtelomeric assembly of a multi-gene pathway for antimicrobial defense compounds in cereals

Non-random gene organization in eukaryotes plays a significant role in genome evolution. Here, we investigate the origin of a biosynthetic gene cluster for production of defence compounds in oat-the avenacin cluster. We elucidate the structure and organisation of this 12-gene cluster, characterise the last two missing pathway steps, and reconstitute the entire pathway in tobacco by transient expression. We show that the cluster has formed de novo since the divergence of oats in a subtelomeric region of the genome that lacks homology with other grasses, and that gene order is approximately colinear with the biosynthetic pathway. We speculate that the positioning of the late pathway genes furthest away from the telomere may mitigate against a 'self-poisoning' scenario in which toxic intermediates accumulate as a result of telomeric gene deletions. Our investigations reveal a striking example of adaptive evolution underpinned by remarkable genome plasticity.

Sympatric genome size variation and hybridization of four oak species as determined by flow cytometry genome size variation and hybridization

The Quercus species serve as a powerful model for studying introgression in relation to species boundaries and adaptive processes. Coexistence of distant relatives, or lack of coexistence of closely relative oak species, introgression may play a role. In the current study, four closely related oak species were found in Zijinshan, China. We generated a comprehensive genome size (GS) database for 120 individuals of four species using flow cytometry-based approaches. We examined GS variability within and among the species and hybridization events among the four species. The mean GSs of Q. acutissima, Q. variabilis, Q. fabri, and Q. serrata var. brevipetiolata were estimated to be 1.87, 1.92, 1.97, and 1.97 pg, respectively. The intraspecific and interspecific variations of GS observed among the four oak species indicated adaptation to the environment. Hybridization occurred both within and between the sections. A hybrid offspring was produced from Q. fabri and Q. variabilis, which belonged to different sections. The GS evolutionary pattern for hybrid species was expansion. Hybridization between the sections may be affected by habitat disturbance. This study increases our understanding of the evolution of GS in Quercus and will help establish guidelines for the ecological protection of oak trees.

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.

Genome-Wide Association Study Reveals the Genetic Architecture of Seed Vigor in Oats

Seed vigor is crucial for crop early establishment in the field and is particularly important for forage crop production. Oat (Avena sativa L.) is a nutritious food crop and also a valuable forage crop. However, little is known about the genetics of seed vigor in oats. To investigate seed vigor-related traits and their genetic architecture in oats, we developed an easy-to-implement image-based phenotyping pipeline and applied it to 650 elite oat lines from the Collaborative Oat Research Enterprise (CORE). Root number, root surface area, and shoot length were measured in two replicates. Variables such as growth rate were derived. Using a genome-wide association (GWA) approach, we identified 34 and 16 unique loci associated with root traits and shoot traits, respectively, which corresponded to 41 and 16 unique SNPs at a false discovery rate < 0.1. Nine root-associated loci were organized into four sets of homeologous regions, while nine shoot-associated loci were organized into three sets of homeologous regions. The context sequences of five trait-associated markers matched to the sequences of rice, Brachypodium and maize (E-value < 10-10), including three markers matched to known gene models with potential involvement in seed vigor. These were a glucuronosyltransferase, a mitochondrial carrier protein domain containing protein, and an iron-sulfur cluster protein. This study presents the first GWA study on oat seed vigor and data of this study can provide guidelines and foundation for further investigations.

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.

Defensins of Grasses: A Systematic Review

The grass family (Poaceae) is one of the largest families of flowering plants, growing in all climatic zones of all continents, which includes species of exceptional economic importance. The high adaptability of grasses to adverse environmental factors implies the existence of efficient resistance mechanisms that involve the production of antimicrobial peptides (AMPs). Of plant AMPs, defensins represent one of the largest and best-studied families. Although wheat and barley seed γ-thionins were the first defensins isolated from plants, the functional characterization of grass defensins is still in its infancy. In this review, we summarize the current knowledge of the characterized defensins from cultivated and selected wild-growing grasses. For each species, isolation of defensins or production by heterologous expression, peptide structure, biological activity, and structure-function relationship are described, along with the gene expression data. We also provide our results on in silico mining of defensin-like sequences in the genomes of all described grass species and discuss their potential functions. The data presented will form the basis for elucidation of the mode of action of grass defensins and high adaptability of grasses to environmental stress and will provide novel potent molecules for practical use in medicine and agriculture.

A Chromosome-Scale Assembly of the Garden Orach (Atriplex hortensis L.) Genome Using Oxford Nanopore Sequencing

Atriplex hortensis (2n = 2x = 18, 1C genome size ∼1.1 gigabases), also known as garden orach and mountain-spinach, is a highly nutritious, broadleaf annual of the Amaranthaceae-Chenopodiaceae alliance (Chenopodiaceae sensu stricto, subfam. Chenopodioideae) that has spread in cultivation from its native primary domestication area in Eurasia to other temperate and subtropical regions worldwide. Atriplex L. is a highly complex but, as understood now, a monophyletic group of mainly halophytic and/or xerophytic plants, of which A. hortensis has been a vegetable of minor importance in some areas of Eurasia (from Central Asia to the Mediterranean) at least since antiquity. Nonetheless, it is a crop with tremendous nutritional potential due primarily to its exceptional leaf and seed protein quantities (approaching 30%) and quality (high levels of lysine). Although there is some literature describing the taxonomy and production of A. hortensis, there is a general lack of genetic and genomic data that would otherwise help elucidate the genetic variation, phylogenetic positioning, and future potential of the species. Here, we report the assembly of the first high-quality, chromosome-scale reference genome for A. hortensis cv. "Golden." Long-read data from Oxford Nanopore's MinION DNA sequencer was assembled with the program Canu and polished with Illumina short reads. Contigs were scaffolded to chromosome scale using chromatin-proximity maps (Hi-C) yielding a final assembly containing 1,325 scaffolds with a N50 of 98.9 Mb - with 94.7% of the assembly represented in the nine largest, chromosome-scale scaffolds. Sixty-six percent of the genome was classified as highly repetitive DNA, with the most common repetitive elements being Gypsy-(32%) and Copia-like (11%) long-terminal repeats. The annotation was completed using MAKER which identified 37,083 gene models and 2,555 tRNA genes. Completeness of the genome, assessed using the Benchmarking Universal Single Copy Orthologs (BUSCO) metric, identified 97.5% of the conserved orthologs as complete, with only 2.2% being duplicated, reflecting the diploid nature of A. hortensis. A resequencing panel of 21 wild, unimproved and cultivated A. hortensis accessions revealed three distinct populations with little variation within subpopulations. These resources provide vital information to better understand A. hortensis and facilitate future study.

Identification of molecular markers for the Pc39 gene conferring resistance to crown rust in oat

KEY MESSAGE

Six new PCR-based markers for the Pc39 crown rust resistance gene in Avena sativa L. were developed. Pc39 was mapped to Mrg11 of the oat consensus map using BLASTn analysis. The aim of this study was the identification of molecular markers for the Pc39 gene in cultivated oat (Avena sativa L.). Pc39 is a major race-specific crown rust resistance gene originally found in an Israeli accession of the wild hexaploid Avena sterilis. The effectiveness of this gene in Europe has decreased in recent years, but is still relatively high and breeding programs would benefit from the availability of molecular markers to aid in its mapping and deployment. The complexity of the oat genome poses a significant obstacle to genetic research. No oat rust resistance genes have yet been cloned, and even the number of relevant molecular markers is very limited. Here, genotyping of a segregating population derived from a cross 'Celer' (Pc39)/STH9210 (susceptible) was conducted using RAPD- and SRAP-PCR-based methods, as well as microarray-based DArT™ and next-generation sequencing DArTseq™ techniques. Markers associated with Pc39 were placed on the hexaploid oat consensus linkage group Mrg11 at 3.7-6.7 cM. Six new PCR-based markers were developed to allow identification of the resistant Pc39 allele. These tightly linked markers will be useful in marker-assisted selection, with the closest, SCAR_3456624, being within 0.37 cM of Pc39. The newly developed markers could find applications in the fine mapping or positional cloning of this gene. Moreover, easy-to-use PCR-based markers linked to Pc39 could facilitate the utilization of this gene in oat breeding programs, especially as a component of crown rust resistance gene pyramids.

Implementing within-cross genomic prediction to reduce oat breeding costs

A barrier to the adoption of genomic prediction in small breeding programs is the initial cost of genotyping material. Although decreasing, marker costs are usually higher than field trial costs. In this study we demonstrate the utility of stratifying a narrow-base biparental oat population genotyped with a modest number of markers to employ genomic prediction at early and later generations. We also show that early generation genotyping data can reduce the number of lines for later phenotyping based on selections of siblings to progress. Using sets of small families selected at an early generation could enable the use of genomic prediction for adaptation to multiple target environments at an early stage in the breeding program. In addition, we demonstrate that mixed marker data can be effectively integrated to combine cheap dominant marker data (including legacy data) with more expensive but higher density codominant marker data in order to make within generation and between lineage predictions based on genotypic information. Taken together, our results indicate that small programs can test and initiate genomic predictions using sets of stratified, narrow-base populations and incorporating low density legacy genotyping data. This can then be scaled to include higher density markers and a broadened population base.

Genomic insights from the first chromosome-scale assemblies of oat (Avena spp.) diploid species

BACKGROUND

Cultivated hexaploid oat (Common oat; Avena sativa) has held a significant place within the global crop community for centuries; although its cultivation has decreased over the past century, its nutritional benefits have garnered increased interest for human consumption. We report the development of fully annotated, chromosome-scale assemblies for the extant progenitor species of the As- and Cp-subgenomes, Avena atlantica and Avena eriantha respectively. The diploid Avena species serve as important genetic resources for improving common oat's adaptive and food quality characteristics.

RESULTS

The A. atlantica and A. eriantha genome assemblies span 3.69 and 3.78 Gb with an N50 of 513 and 535 Mb, respectively. Annotation of the genomes, using sequenced transcriptomes, identified ~ 50,000 gene models in each species-including 2965 resistance gene analogs across both species. Analysis of these assemblies classified much of each genome as repetitive sequence (~ 83%), including species-specific, centromeric-specific, and telomeric-specific repeats. LTR retrotransposons make up most of the classified elements. Genome-wide syntenic comparisons with other members of the Pooideae revealed orthologous relationships, while comparisons with genetic maps from common oat clarified subgenome origins for each of the 21 hexaploid linkage groups. The utility of the diploid genomes was demonstrated by identifying putative candidate genes for flowering time (HD3A) and crown rust resistance (Pc91). We also investigate the phylogenetic relationships among other A- and C-genome Avena species.

CONCLUSIONS

The genomes we report here are the first chromosome-scale assemblies for the tribe Poeae, subtribe Aveninae. Our analyses provide important insight into the evolution and complexity of common hexaploid oat, including subgenome origin, homoeologous relationships, and major intra- and intergenomic rearrangements. They also provide the annotation framework needed to accelerate gene discovery and plant breeding.

Comparative linkage mapping of diploid, tetraploid, and hexaploid Avena species suggests extensive chromosome rearrangement in ancestral diploids

The genus Avena (oats) contains diploid, tetraploid and hexaploid species that evolved through hybridization and polyploidization. Four genome types (named A through D) are generally recognized. We used GBS markers to construct linkage maps of A genome diploid (Avena strigosa x A. wiestii, 2n = 14), and AB genome tetraploid (A. barbata 2n = 28) oats. These maps greatly improve coverage from older marker systems. Seven linkage groups in the tetraploid showed much stronger homology and synteny with the A genome diploids than did the other seven, implying an allopolyploid hybrid origin of A. barbata from distinct A and B genome diploid ancestors. Inferred homeologies within A. barbata revealed that the A and B genomes are differentiated by several translocations between chromosomes within each subgenome. However, no translocation exchanges were observed between A and B genomes. Comparison to a consensus map of ACD hexaploid A. sativa (2n = 42) revealed that the A and D genomes of A. sativa show parallel rearrangements when compared to the A genomes of the diploids and tetraploids. While intergenomic translocations are well known in polyploid Avena, our results are most parsimoniously explained if translocations also occurred in the A, B and D genome diploid ancestors of polyploid Avena.

Polyploid evolution: The ultimate way to grasp the nettle

Polyploidy is one of the major forces of plant evolution and widespread mixed-ploidy species offer an opportunity to evaluate its significance. We therefore selected the cosmopolitan species Urtica dioica (stinging nettle), examined its cytogeography and pattern of absolute genome size, and assessed correlations with bioclimatic and ecogeographic data (latitude, longitude, elevation). We evaluated variation in ploidy level using an extensive dataset of 7012 samples from 1317 populations covering most of the species' distribution area. The widespread tetraploid cytotype (87%) was strongly prevalent over diploids (13%). A subsequent analysis of absolute genome size proved a uniform Cx-value of core U. dioica (except for U. d. subsp. cypria) whereas other closely related species, namely U. bianorii, U. kioviensis and U. simensis, differed significantly. We detected a positive correlation between relative genome size and longitude and latitude in the complete dataset of European populations and a positive correlation between relative genome size and longitude in a reduced dataset of diploid accessions (the complete dataset of diploids excluding U. d. subsp. kurdistanica). In addition, our data indicate an affinity of most diploids to natural and near-natural habitats and that the tetraploid cytotype and a small part of diploids (population from the Po river basin in northern Italy) tend to inhabit synanthropic sites. To sum up, the pattern of ploidy variation revealed by our study is in many aspects unique to the stinging nettle, being most likely first of all driven by the greater ecological plasticity and invasiveness of the tetraploid cytotype.

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

BACKGROUND

Repetitive DNA motifs - not coding genetic information and repeated millions to hundreds of times - make up the majority of many genomes. Here, we identify the nature, abundance and organization of all the repetitive DNA families in oats (Avena sativa, 2n = 6x = 42, AACCDD), a recognized health-food, and its wild relatives.

RESULTS

Whole-genome sequencing followed by k-mer and RepeatExplorer graph-based clustering analyses enabled assessment of repetitive DNA composition in common oat and its wild relatives' genomes. Fluorescence in situ hybridization (FISH)-based karyotypes are developed to understand chromosome and repetitive sequence evolution of common oat. We show that some 200 repeated DNA motifs make up 70% of the Avena genome, with less than 20 families making up 20% of the total. Retroelements represent the major component, with Ty3/Gypsy elements representing more than 40% of all the DNA, nearly three times more abundant than Ty1/Copia elements. DNA transposons are about 5% of the total, while tandemly repeated, satellite DNA sequences fit into 55 families and represent about 2% of the genome. The Avena species are monophyletic, but both bioinformatic comparisons of repeats in the different genomes, and in situ hybridization to metaphase chromosomes from the hexaploid species, shows that some repeat families are specific to individual genomes, or the A and D genomes together. Notably, there are terminal regions of many chromosomes showing different repeat families from the rest of the chromosome, suggesting presence of translocations between the genomes.

CONCLUSIONS

The relatively small number of repeat families shows there are evolutionary constraints on their nature and amplification, with mechanisms leading to homogenization, while repeat characterization is useful in providing genome markers and to assist with future assemblies of this large genome (c. 4100 Mb in the diploid). The frequency of inter-genomic translocations suggests optimum strategies to exploit genetic variation from diploid oats for improvement of the hexaploid may differ from those used widely in bread wheat.

Genome Size Unaffected by Variation in Morphological Traits, Temperature, and Precipitation in Turnip

Genome size (GS) was proposed as proxy for gross phenotypic and environmental changes in plants. GS organismal complexity is an enigma in evolutionary biology. While studies pertaining to intraspecific GS variation are abundant, literatures reporting the adaptive significance of GS are largelymissing. During food shortage, Brassica rapa var. rapa (turnip) is used as food and fodder for sustaining the livelihood of residents in the Qinghai Tibetan Plateau (QTP), which is also known as “the roof of the world”. Thus, climatic extremities make this region a natural environment to test adaptive significance of GS variation in turnip landraces. Therefore, from the QTP and its adjacent regions (the Hengduanshan and the Himalayas), we investigated adaptive evolution of GS in turnip landraces. Tuber diameter of turnip landraces was found to be significantly correlated with most of the environmental factors. GS was also shown not to be associated with morphological traits, temperature, and precipitation. Moreover, principal component analyses based on the whole dataset trisected the landraces into three distinct populations based on landrace usage—Hengduanshan, QTP, and the Himalayas. Nonetheless, our cumulative dataset showed evidence of adaptation of turnip landrace to different environments throughnonassociated genomic and phenomic plasticity.

Phylogenetic relationships in the genus Avena based on the nuclear Pgk1 gene

The phylogenetic relationships among 76 Avena taxa, representing 14 diploids, eight tetraploids, and four hexaploids were investigated by using the nuclear plastid 3-phosphoglycerate kinase gene (Pgk1). A significant deletion (131 bp) was detected in all the C genome homoeologues which reconfirmed a major structural divergence between the A and C genomes. Phylogenetic analysis indicated the Cp genome is more closely related to the polyploid species than is the Cv genome. Two haplotypes of Pgk1 gene were obtained from most of the AB genome tetraploids. Both types of the barbata group showed a close relationship with the As genome diploid species, supporting the hypothesis that both the A and B genomes are derived from an As genome. Two haplotypes were also detected in A. agadiriana, which showed close relationships with the As genome diploid and the Ac genome diploid, respectively, emphasizing the important role of the Ac genome in the evolution of A. agadiriana. Three homoeologues of the Pgk1 gene were detected in five hexaploid accessions. The homoeologues that might represent the D genome were tightly clustered with the tetraploids A. maroccana and A. murphyi, but did not show a close relationship with any extant diploid species.

Haplotype-based genotyping-by-sequencing in oat genome research

In a de novo genotyping-by-sequencing (GBS) analysis of short, 64-base tag-level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag-level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635-line diversity panel were used to infer chromosome-level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high-resolution genome analysis and genomic selection in oats. A combined genome-wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component-based genome-wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS-derived markers facilitate genome analysis and genomic selection in oat.

Genetic mapping and a new PCR-based marker linked to a dwarfing gene in oat (Avena sativa L.)

Short straw is a desired trait in cultivated hexaploid oat (Avena sativa L.) for some production environments. Marker-assisted selection, a key tool for achieving this objective, is limited by a lack of mapping data and available markers. Here, bulked-segregant analysis was used to identify PCR-based markers associated with a dwarfing gene. Genetic analysis identified a monogenic dominant inheritance of one dwarfing gene from WAOAT2132, temporarily designated DwWA. A simple sequence repeat (SSR) marker (AME117) that was already available and a new codominant PCR-based marker (bi17) developed by homologous cloning in the present study were both associated with the dwarfing gene. The two markers were located 21 and 1.2 cM from DwWA, respectively. The bi17 marker was mapped to neighboring SNP markers on chromosome 18D of the oat consensus map. Since Dw6 was previously mapped on chromosome 18, and since our new marker bi17 is also diagnostic for NILs generated for Dw6, there is strong evidence that the dwarfing gene identified in WAOAT2132 is Dw6. The newly developed markers could find applications in the identification of this gene in oat germplasm and in the fine mapping or positional cloning of the gene.

Genomic relationships among sixteen species of Avena based on (ACT)6 trinucleotide repeat FISH

Knowledge of the locations of repeat elements could be very important in the assembly of genome sequences and their assignment to physical chromosomes. Genomic and species relationships among 16 species were investigated using fluorescence in situ hybridization (FISH) with the Am1 and (ACT)6 probes. The Am1 oligonucleotide probe was particularly enriched in the C genomes, whereas the (ACT)6 trinucleotide repeat probe showed a diverse distribution of hybridization patterns in the A, AB, C, AC, and ACD genomes but might not be present in the B and D genomes. The hybridization pattern of Avena sativa was very similar to that of A. insularis, indicating that this species most likely originated from A. insularis as a tetraploid ancestor. Although the two FISH probes failed to identify relationships of more species, this proof-of-concept approach opens the way to the use of FISH probes in assigning other signature elements from genomic sequence to physical chromosomes.

Unraveling the evolutionary dynamics of ancient and recent polyploidization events in Avena (Poaceae)

Understanding the diversification of polyploid crops in the circum-Mediterranean region is a challenging issue in evolutionary biology. Sequence data of three nuclear genes and three plastid DNA fragments from 109 accessions of Avena L. (Poaceae) and the outgroups were used for maximum likelihood and Bayesian analyses. The evolution of cultivated oat (Avena sativa L.) and its close relatives was inferred to have involved ancient allotetraploidy and subsequent recent allohexaploidy events. The crown ages of two infrageneric lineages (Avena sect. Ventricosa Baum ex Romero-Zarco and Avena sect. Avena) were estimated to be in the early to middle Miocene, and the A. sativa lineages were dated to the late Miocene to Pliocene. These periods coincided with the mild seasonal climatic contrasts and the Mediterranean climate established in the Mediterranean Basin. Our results suggest that polyploidy, lineage divergence, and complex reticulate evolution have occurred in Avena, exemplifying the long-term persistence of tetraploids and the multiple origins of hexaploids related to paleoclimatic oscillations during the Miocene-Pliocene interval in the circum-Mediterranean region. This newly-resolved infrageneric phylogenetic framework represents a major step forward in understanding the origin of the cultivated oat.

De novo transcriptome assembly, development of EST-SSR markers and population genetic analyses for the desert biomass willow, Salix psammophila

Salix psammophila, a sandy shrub known as desert willow, is regarded as a potential biomass feedstock and plays an important role in maintaining local ecosystems. However, a lack of genomic data and efficient molecular markers limit the study of its population evolution and genetic breeding. In this study, chromosome counts, flow cytometry and SSR analyses indicated that S. psammophila is tetraploid. A total of 6,346 EST-SSRs were detected based on 71,458 de novo assembled unigenes from transcriptome data. Twenty-seven EST-SSR markers were developed to evaluate the genetic diversity and population structure of S. psammophila from eight natural populations in Northern China. High levels of genetic diversity (mean 10.63 alleles per locus; mean HE 0.689) were dectected in S. psammophila. The weak population structure and little genetic differentiation (pairwise FST = 0.006-0.016) were found among Population 1-Population 7 (Pop1-Pop7; Inner Mongolia and Shaanxi), but Pop8 (Ningxia) was clearly separated from Pop1-Pop7 and moderate differentiation (pairwise FST = 0.045-0.055) was detected between them, which may be influenced by local habitat conditions. Molecular variance analyses indicated that most of the genetic variation (94.27%) existed within populations. These results provide valuable genetic informations for natural resource conservation and breeding programme optimisation of S. psammophila.

High-density marker profiling confirms ancestral genomes of Avena species and identifies D-genome chromosomes of hexaploid oat

KEY MESSAGE

Genome analysis of 27 oat species identifies ancestral groups, delineates the D genome, and identifies ancestral origin of 21 mapped chromosomes in hexaploid oat. We investigated genomic relationships among 27 species of the genus Avena using high-density genetic markers revealed by genotyping-by-sequencing (GBS). Two methods of GBS analysis were used: one based on tag-level haplotypes that were previously mapped in cultivated hexaploid oat (A. sativa), and one intended to sample and enumerate tag-level haplotypes originating from all species under investigation. Qualitatively, both methods gave similar predictions regarding the clustering of species and shared ancestral genomes. Furthermore, results were consistent with previous phylogenies of the genus obtained with conventional approaches, supporting the robustness of whole genome GBS analysis. Evidence is presented to justify the final and definitive classification of the tetraploids A. insularis, A. maroccana (=A. magna), and A. murphyi as containing D-plus-C genomes, and not A-plus-C genomes, as is most often specified in past literature. Through electronic painting of the 21 chromosome representations in the hexaploid oat consensus map, we show how the relative frequency of matches between mapped hexaploid-derived haplotypes and AC (DC)-genome tetraploids vs. A- and C-genome diploids can accurately reveal the genome origin of all hexaploid chromosomes, including the approximate positions of inter-genome translocations. Evidence is provided that supports the continued classification of a diverged B genome in AB tetraploids, and it is confirmed that no extant A-genome diploids, including A. canariensis, are similar enough to the D genome of tetraploid and hexaploid oat to warrant consideration as a D-genome diploid.

A Consensus Map in Cultivated Hexaploid Oat Reveals Conserved Grass Synteny with Substantial Subgenome Rearrangement

Hexaploid oat ( L., 2 = 6 = 42) is a member of the Poaceae family and has a large genome (∼12.5 Gb) containing 21 chromosome pairs from three ancestral genomes. Physical rearrangements among parental genomes have hindered the development of linkage maps in this species. The objective of this work was to develop a single high-density consensus linkage map that is representative of the majority of commonly grown oat varieties. Data from a cDNA-derived single-nucleotide polymorphism (SNP) array and genotyping-by-sequencing (GBS) were collected from the progeny of 12 biparental recombinant inbred line populations derived from 19 parents representing oat germplasm cultivated primarily in North America. Linkage groups from all mapping populations were compared to identify 21 clusters of conserved collinearity. Linkage groups within each cluster were then merged into 21 consensus chromosomes, generating a framework consensus map of 7202 markers spanning 2843 cM. An additional 9678 markers were placed on this map with a lower degree of certainty. Assignment to physical chromosomes with high confidence was made for nine chromosomes. Comparison of homeologous regions among oat chromosomes and matches to orthologous regions of rice ( L.) reveal that the hexaploid oat genome has been highly rearranged relative to its ancestral diploid genomes as a result of frequent translocations among chromosomes. Heterogeneous chromosome rearrangements among populations were also evident, probably accounting for the failure of some linkage groups to match the consensus. This work contributes to a further understanding of the organization and evolution of hexaploid grass genomes.