Single-molecule sequencing and Hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution.

BACKGROUND: Amaranth (Amaranthus hypochondriacus) was a food staple among the ancient civilizations of Central and South America that has recently received increased attention due to the high nutritional value of the seeds, with the potential to help alleviate malnutrition and food security concerns, particularly in arid and semiarid regions of the developing world. Here, we present a reference-quality assembly of the amaranth genome which will assist the agronomic development of the species. RESULTS: Utilizing single-molecule, real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C) to close assembly gaps and scaffold contigs, respectively, we improved our previously reported Illumina-based assembly to produce a chromosome-scale assembly with a scaffold N50 of 24.4 Mb. The 16 largest scaffolds contain 98% of the assembly and likely represent the haploid chromosomes (n = 16). To demonstrate the accuracy and utility of this approach, we produced physical and genetic maps and identified candidate genes for the betalain pigmentation pathway. The chromosome-scale assembly facilitated a genome-wide syntenic comparison of amaranth with other Amaranthaceae species, revealing chromosome loss and fusion events in amaranth that explain the reduction from the ancestral haploid chromosome number (n = 18) for a tetraploid member of the Amaranthaceae. CONCLUSIONS: The assembly method reported here minimizes cost by relying primarily on short-read technology and is one of the first reported uses of in vivo Hi-C for assembly of a plant genome. Our analyses implicate chromosome loss and fusion as major evolutionary events in the 2n = 32 amaranths and clearly establish the homoeologous relationship among most of the subgenome chromosomes, which will facilitate future investigations of intragenomic changes that occurred post polyploidization.

The Amaranth Genome: Genome, Transcriptome, and Physical Map Assembly.

Amaranth ( L.) is an emerging pseudocereal native to the New World that has garnered increased attention in recent years because of its nutritional quality, in particular its seed protein and more specifically its high levels of the essential amino acid lysine. It belongs to the Amaranthaceae family, is an ancient paleopolyploid that shows disomic inheritance (2 = 32), and has an estimated genome size of 466 Mb. Here we present a high-quality draft genome sequence of the grain amaranth. The genome assembly consisted of 377 Mb in 3518 scaffolds with an N of 371 kb. Repetitive element analysis predicted that 48% of the genome is comprised of repeat sequences, of which -like elements were the most commonly classified retrotransposon. A de novo transcriptome consisting of 66,370 contigs was assembled from eight different amaranth tissue and abiotic stress libraries. Annotation of the genome identified 23,059 protein-coding genes. Seven grain amaranths (, , and ) and their putative progenitor () were resequenced. A single nucleotide polymorphism (SNP) phylogeny supported the classification of as the progenitor species of the grain amaranths. Lastly, we generated a de novo physical map for using the BioNano Genomics' Genome Mapping platform. The physical map spanned 340 Mb and a hybrid assembly using the BioNano physical maps nearly doubled the N of the assembly to 697 kb. Moreover, we analyzed synteny between amaranth and sugar beet ( L.) and estimated, using analysis, the age of the most recent polyploidization event in amaranth.

New Diversity Arrays Technology (DArT) markers for tetraploid oat (Avena magna Murphy et Terrell) provide the first complete oat linkage map and markers linked to domestication genes from hexaploid A. sativa L.

Nutritional benefits of cultivated oat (Avena sativa L., 2n = 6x = 42, AACCDD) are well recognized; however, seed protein levels are modest and resources for genetic improvement are scarce. The wild tetraploid, A. magna Murphy et Terrell (syn A. maroccana Gdgr., 2n = 4x = 28, CCDD), which contains approximately 31% seed protein, was hybridized with cultivated oat to produce a domesticated A. magna. Wild and cultivated accessions were crossed to generate a recombinant inbred line (RIL) population. Although these materials could be used to develop domesticated, high-protein oat, mapping and quantitative trait loci introgression is hindered by a near absence of genetic markers. Objectives of this study were to develop high-throughput, A. magna-specific markers; generate a genetic linkage map based on the A. magna RIL population; and map genes controlling oat domestication. A Diversity Arrays Technology (DArT) array derived from 10 A. magna genotypes was used to generate 2,688 genome-specific probes. These, with 12,672 additional oat clones, produced 2,349 polymorphic markers, including 498 (21.2%) from A. magna arrays and 1,851 (78.8%) from other Avena libraries. Linkage analysis included 974 DArT markers, 26 microsatellites, 13 SNPs, and 4 phenotypic markers, and resulted in a 14-linkage-group map. Marker-to-marker correlation coefficient analysis allowed classification of shared markers as unique or redundant, and putative linkage-group-to-genome anchoring. Results of this study provide for the first time a collection of high-throughput tetraploid oat markers and a comprehensive map of the genome, providing insights to the genome ancestry of oat and affording a resource for study of oat domestication, gene transfer, and comparative genomics.

Chromosomal localization of two novel repetitive sequences isolated from the Chenopodium quinoa Willd. genome.

The chromosomal organization of two novel repetitive DNA sequences isolated from the Chenopodium quinoa Willd. genome was analyzed across the genomes of selected Chenopodium species. Fluorescence in situ hybridization (FISH) analysis with the repetitive DNA clone 18-24J in the closely related allotetraploids C. quinoa and Chenopodium berlandieri Moq. (2n = 4x = 36) evidenced hybridization signals that were mainly present on 18 chromosomes; however, in the allohexaploid Chenopodium album L. (2n = 6x = 54), cross-hybridization was observed on all of the chromosomes. In situ hybridization with rRNA gene probes indicated that during the evolution of polyploidy, the chenopods lost some of their rDNA loci. Reprobing with rDNA indicated that in the subgenome labeled with 18-24J, one 35S rRNA locus and at least half of the 5S rDNA loci were present. A second analyzed sequence, 12-13P, localized exclusively in pericentromeric regions of each chromosome of C. quinoa and related species. The intensity of the FISH signals differed considerably among chromosomes. The pattern observed on C. quinoa chromosomes after FISH with 12-13P was very similar to GISH results, suggesting that the 12-13P sequence constitutes a major part of the repetitive DNA of C. quinoa.

A new chromosome nomenclature system for oat (Avena sativa L. and A. byzantina C. Koch) based on FISH analysis of monosomic lines.

Fluorescent in situ hybridization (FISH) with multiple probes was used to analyze mitotic and meiotic chromosome spreads of Avena sativa cv 'Sun II' monosomic lines, and of A. byzantina cv 'Kanota' monosomic lines from spontaneous haploids. The probes used were A. strigosa pAs120a (a repetitive sequence abundant in A-genome chromatin), A. murphyi pAm1 (a repetitive sequence abundant in C-genome chromatin), A. strigosa pITS (internal transcribed spacer of rDNA) and the wheat rDNA probes pTa71 (nucleolus organizer region or NOR) and pTa794 (5S). Simultaneous and sequential FISH employing pairs of these probes allowed the identification and genome assignation of all chromosomes. FISH mapping using mitotic and meiotic metaphases facilitated the genomic and chromosomal identification of the monosome in each line. Of the 17 'Sun II' lines analyzed, 13 distinct monosomic lines were found, corresponding to four monosomes of the A-genome, five of the C-genome and four of the D-genome. In addition, 12 distinct monosomic lines were detected among the 20 'Kanota' lines examined, corresponding to six monosomes of the A-genome, three of the C-genome and three of the D-genome. The results show that 19 chromosomes out of 21 of the complement are represented by monosomes between the two genetic backgrounds. The identity of the remaining chromosomes can be deduced either from one intergenomic translocation detected on both 'Sun II' and 'Kanota' lines, or from the single reciprocal, intergenomic translocation detected among the 'Sun II' lines. These results permit a new system to be proposed for numbering the 21 chromosome pairs of the hexaploid oat complement. Accordingly, the A-genome contains chromosomes 8A, 11A, 13A, 15A, 16A, 17A and 19A; the C-genome contains chromosomes 1C, 2C, 3C, 4C, 5C, 6C and 7C; and the D-genome consists of chromosomes 9D, 10D, 12D, 14D, 18D, 20D and 21D. Moreover, the FISH patterns of 16 chromosomes in 'Sun II' and 15 in 'Kanota' suggest that these chromosomes could be involved in intergenomic translocations. By comparing the identities of individually translocated chromosomes in the two hexaploid species with those of other hexaploids, we detected different types of intergenomic translocations.

Characterization of Salt Overly Sensitive 1 (SOS1) gene homoeologs in quinoa (Chenopodium quinoa Willd.).

Salt tolerance is an agronomically important trait that affects plant species around the globe. The Salt Overly Sensitive 1 (SOS1) gene encodes a plasma membrane Na+/H+ antiporter that plays an important role in germination and growth of plants in saline environments. Quinoa (Chenopodium quinoa Willd.) is a halophytic, allotetraploid grain crop of the family Amaranthaceae with impressive nutritional content and an increasing worldwide market. Many quinoa varieties have considerable salt tolerance, and research suggests quinoa may utilize novel mechanisms to confer salt tolerance. Here we report the cloning and characterization of two homoeologous SOS1 loci (cqSOS1A and cqSOS1B) from C. quinoa, including full-length cDNA sequences, genomic sequences, relative expression levels, fluorescent in situ hybridization (FISH) analysis, and a phylogenetic analysis of SOS1 genes from 13 plant taxa. The cqSOS1A and cqSOS1B genes each span 23 exons spread over 3477 bp and 3486 bp of coding sequence, respectively. These sequences share a high level of similarity with SOS1 homologs of other species and contain two conserved domains, a Nhap cation-antiporter domain and a cyclic-nucleotide binding domain. Genomic sequence analysis of two BAC clones (98 357 bp and 132 770 bp) containing the homoeologous SOS1 genes suggests possible conservation of synteny across the C. quinoa sub-genomes. This report represents the first molecular characterization of salt-tolerance genes in a halophytic species in the Amaranthaceae as well as the first comparative analysis of coding and non-coding DNA sequences of the two homoeologous genomes of C. quinoa.

Simple sequence repeat marker development and genetic mapping in quinoa (Chenopodium quinoa Willd.).

Quinoa is a regionally important grain crop in the Andean region of South America. Recently quinoa has gained international attention for its high nutritional value and tolerances of extreme abiotic stresses. DNA markers and linkage maps are important tools for germplasm conservation and crop improvement programmes. Here we report the development of 216 new polymorphic SSR (simple sequence repeats) markers from libraries enriched for GA, CAA and AAT repeats, as well as 6 SSR markers developed from bacterial artificial chromosome-end sequences (BES-SSRs). Heterozygosity (H) values of the SSR markers ranges from 0.12 to 0.90, with an average value of 0.57. A linkage map was constructed for a newly developed recombinant inbred lines (RIL) population using these SSR markers. Additional markers, including amplified fragment length polymorphisms (AFLPs), two 11S seed storage protein loci, and the nucleolar organizing region (NOR), were also placed on the linkage map. The linkage map presented here is the first SSR-based map in quinoa and contains 275 markers, including 200 SSR. The map consists of 38 linkage groups (LGs) covering 913 cM. Segregation distortion was observed in the mapping population for several marker loci, indicating possible chromosomal regions associated with selection or gametophytic lethality. As this map is based primarily on simple and easily-transferable SSR markers, it will be particularly valuable for research in laboratories in Andean regions of South America.

Molecular and cytological characterization of ribosomal RNA genes in Chenopodium quinoa and Chenopodium berlandieri.

The nucleolus organizer region (NOR) and 5S ribosomal RNA (rRNA) genes are valuable as chromosome landmarks and in evolutionary studies. The NOR intergenic spacers (IGS) and 5S rRNA nontranscribed spacers (NTS) were PCR-amplified and sequenced from 5 cultivars of the Andean grain crop quinoa (Chenopodium quinoa Willd., 2n = 4x = 36) and a related wild ancestor (C. berlandieri Moq. subsp. zschackei (Murr) A. Zobel, 2n = 4x = 36). Length heterogeneity observed in the IGS resulted from copy number difference in subrepeat elements, small re arrangements, and species-specific indels, though the general sequence composition of the 2 species was highly similar. Fifteen of the 41 sequence polymorphisms identified among the C. quinoa lines were synapomorphic and clearly differentiated the highland and lowland ecotypes. Analysis of the NTS sequences revealed 2 basic NTS sequence classes that likely originated from the 2 allopolyploid subgenomes of C. quinoa. Fluorescence in situ hybridization (FISH) analysis showed that C. quinoa possesses an interstitial and a terminal pair of 5S rRNA loci and only 1 pair of NOR, suggesting a reduction in the number of rRNA loci during the evolution of this species. C. berlandieri exhibited variation in both NOR and 5S rRNA loci without changes in ploidy.

Construction of a quinoa (Chenopodium quinoa Willd.) BAC library and its use in identifying genes encoding seed storage proteins.

Quinoa (Chenopodium quinoa Willd.) is adapted to the harsh environments of the Andean Altiplano region. Its seeds have a well-balanced amino acid composition and exceptionally high protein content with respect to human nutrition. Quinoa grain is a staple in the diet of some of the most impoverished people in the world. The plant is an allotetraploid displaying disomic inheritance (2n=4x=36) with a di-haploid genome of 967 Mbp (megabase pair), or 2C=2.01 pg. We constructed two quinoa BAC libraries using BamHI (26,880 clones) and EcoRI (48,000 clones) restriction endonucleases. Cloned inserts in the BamHI library average 113 kb (kilobase) with approximately 2% of the clones lacking inserts, whereas cloned inserts in the EcoRI library average 130 kb and approximately 1% lack inserts. Three plastid genes used as probes of high-density arrayed blots of 73,728 BACs identified approximately 2.8% of the clones as containing plastid DNA inserts. We estimate that the combined quinoa libraries represent at least 9.0 di-haploid nuclear genome equivalents. An average of 12.2 positive clones per probe were identified with 13 quinoa single-copy ESTs as probes of the high-density arrayed blots, suggesting that the estimate of 9.0x coverage of the genome is conservative. Utility of the BAC libraries for gene identification was demonstrated by probing the library with a partial sequence of the 11S globulin seed storage protein gene and identifying multiple positive clones. The presence of the 11S globulin gene in four of the clones was verified by direct comparison with quinoa genomic DNA on a Southern blot. Besides serving as a useful tool for gene identification, the quinoa BAC libraries will be an important resource for physical mapping of the quinoa genome.

An algorithm for analyzing linkages affected by heterozygous translocations: QuadMap.

Heterozygous chromosome rearrangements such as reciprocal translocations are most accurately displayed as two-dimensional linkage maps. Standard linkage mapping software packages, such as MapMaker, generate only one-dimensional maps and so reciprocal translocations appear as clusters of markers, even though they originate from two nonhomologous chromosomes. To more accurately map these regions, researchers have developed statistical methods that use the variance in map distance to distinguish among the four segments (two translocation, two interstitial) of the translocation. In this study, we describe modifications to one of these protocols, that proposed by Livingstone et al. (2000). We also introduce QuadMap, a new software application for dissecting heterozygous translocation-affected linkage maps.

A genetic linkage map of quinoa ( Chenopodium quinoa) based on AFLP, RAPD, and SSR markers.

Quinoa ( Chenopodium quinoa Willd.) is an important seed crop for human consumption in the Andean region of South America. It is the primary staple in areas too arid or saline for the major cereal crops. The objective of this project was to build the first genetic linkage map of quinoa. Selection of the mapping population was based on a preliminary genetic similarity analysis of four potential mapping parents. Breeding lines 'Ku-2' and '0654', a Chilean lowland type and a Peruvian Altiplano type, respectively, showed a low similarity coefficient of 0.31 and were selected to form an F(2) mapping population. The genetic map is based on 80 F(2) individuals from this population and consists of 230 amplified length polymorphism (AFLP), 19 simple-sequence repeat (SSR), and six randomly amplified polymorphic DNA markers. The map spans 1,020 cM and contains 35 linkage groups with an average marker density of 4.0 cM per marker. Clustering of AFLP markers was not observed. Additionally, we report the primer sequences and map locations for 19 SSR markers that will be valuable tools for future quinoa genome analysis. This map provides a key starting point for genetic dissection of agronomically important characteristics of quinoa, including seed saponin content, grain yield, maturity, and resistance to disease, frost, and drought. Current efforts are geared towards the generation of more than 200 mapped SSR markers and the development of several recombinant-inbred mapping populations.

C-banding variation in the Moroccan oat species Avena agadiriana (2n=4x=28).

The C-banding technique was used to describe the chromosomes of a relatively recently-discovered Moroccan oat species, Avena agadiriana (2n=4x=28). A substantial amount of polymorphism for arm ratios and C-banding patterns was observed among five accessions of this species. However a common set of ten putatively homologous chromosomes was identifiable among the five accessions. The chromosomes of A. Agadiriana do not closely match those of any of the previously described diploid or tetraploid oat species in terms of their arm ratios and C-banding patterns. However, their overall C-banded appearance generally resembles the A/B/D groups of chromosomes of Avena species, rather than the more hetrochromatic C genomes. Implications of these findings in terms of chromosome evolution in the genus Avena are discussed.

Molecular genetic identification of Avena chromosomes related to the group 1 chromosomes of the Triticeae.

A collection of 19 wheat (Triticum aestivum) probes, detecting sequences in the seven homoeologous groups of chromosomes, were hybridized to DNA from the 'Kanota' series of oat monosomic lines (Avena byzantina) to investigate their use for identifying groups of homoeologous oat chromosomes. Three probes from homoeologous group 1 of wheat, psr161, psr162, and psr121, mapped among the set of oat chromosomes 1C, 14, and 17. One homoeologous group 6 probe, psr167, mapped to oat chromosomes 1C and 17. Two oat probes that had previously been shown to map to oat chromosomes 1C, 14, and 17 were then hybridized to DNA from the 'Chinese Spring' wheat ditelosomics. They localized to homoeologous group 1 wheat chromosomes, one to the short arm and one to the long arm. These results reveal that in hexaploid oat there is a group of three chromosomes that correspond at least in part to homoeologous group 1 of wheat. The remaining wheat probes identifying other wheat homoeologous sets did not detect a complete series of homoeologous chromosomes in oat. This was presumably due to the incomplete status of the 'Kanota' monosomic series, chromosomal rearrangement in Avena, weak hybridization signals owing to low probe-target sequence homology, and (or) detection of only two hybridization bands by the wheat probe.

Genomic in situ hybridization differentiates between A/D- and C-genome chromatin and detects intergenomic translocations in polyploid oat species (genus Avena).

The genomic in situ hybridization (GISH) technique was used to discriminate between chromosomes of the C genome and those of the A and A/D genomes in allopolyploid oat species (genus Avena). Total biotinylated DNA from A. strigosa (2n = 2x = 14, AsAs genome) was mixed with sheared, unlabelled total DNA from A. eriantha (2n = 2x = 14, CpCp) at a ratio of 1:200 (labelled to unlabelled). The resulting hybridization pattern consisted of 28 mostly labelled and 14 mostly unlabelled chromosomes in the hexaploids. Attempts to discriminate between chromosomes of the A and D genomes in A. sativa (2n = 6x = 42, AACCDD) were unsuccessful using GISH. At least eight intergenomic translocation segments were detected in A. sativa 'Ogle', several of which were not observed in A. byzantina 'Kanota' (2n = 6x = 42, AACCDD) or in A. sterilis CW 439-2 (2n = 6x = 42, AACCDD). At least five intergenomic translocation segments were observed in A. maroccana CI 8330 'Magna' (2n = 4x = 28, AACC). In both 'Ogle' and 'Magna', positions of most of these translocations matched with C-banding patterns.

Chromosomal localization and polymorphisms of ribosomal DNA in oat (Avena spp.).

The 17S/5.8S/26S ribosomal DNA (rDNA) sequences were mapped to the three satellited (SAT) chromosomes in the common hexaploid cultivated oat Avena sativa (2n = 6x = 42, AACCDD genomes). In situ hybridization and Southern hybridization of maize and (or) wheat rDNA probes to DNA from nullisomics derived from the cultivar 'Sun II' allowed the placement of rDNA sequences to the physical chromosomes. A restriction map was produced for the rDNA sequences of 'Sun II' using a maize probe from the transcribed region of the 17S/26S rDNA repeat. The set of rDNA repeats on SAT 2 of 'Sun II' possesses a 10.5-kb EcoRI fragment not found in the rDNA repeats of SAT 1 and SAT 8. This 10.5-kb fragment results from the absence of an EcoRI site in the intergenic spacer (IGS) of SAT 2 repeats. Extensive polymorphisms were demonstrated for three hexaploid Avena species, namely, the Mediterranean-type cultivated oat A. byzantina and the wild species A. sterilis and A. fatua. However, geographically diverse A. sativa cultivars displayed little rDNA variation. In contrast with all of the A. sativa cultivars examined, the A. sterilis accessions generally lacked the 10.5-kb EcoRI fragment. The results support the hypothesis that A. sativa accessions descend from a limited ancestral cultivated population. The rDNA polymorphisms are attributed to differences in lengths and restriction sites of the IGS.

Identification of homoeologous chromosomes in hexaploid oat (A. byzantina cv Kanota) using monosomics and RFLP analysis.

The use of RFLP markers, together with a partial set of monosomics available in Avena byzantina cv Kanota, has enabled us to identify putative homoeologous chromosome sets in hexaploid Avena species (2n = 6x = 42, AACCDD). We first identified probes producing distinct three-band patterns on Southern blots that possibly reflect orthologous loci of the three genomes present in the hexaploid. Using monosomic analysis, 51 different restriction fragments that hybridized to 26 probes were localized to 12 different chromosomes for which monosomic stocks were available. These DNA restriction fragments were localized to specific monosomics using image analysis to quantify band intensity relative to other bands in the same lane. From these data, we have tentatively identified two complete homoeologous sets of three chromosomes each and two partial sets of two of the three chromosomes. The results indicate that RFLP dosage analysis is useful in the characterization of homoeologous chromosomes in hexaploid oat where nullisomics for many of the chromosomes are not available.Mention of a trademark or proprietary product does not constitute a guarantee or warranty by the USDA-ARS or the University of Minnesota and does not imply approval over other products that also may be suitable.

C-banded karyotypes and polymorphisms in hexaploid oat accessions (Avena spp.) using Wright's stain.

A chromosome C-banding protocol using Wright's stain was employed to compare chromosomes in cultivars and wild accessions of several hexaploid oat taxa (Avena spp.). This technique permits the identification of each of the 21 somatic hexaploid oat chromosomes. Digital images of C-banded cells were captured on computer and used to construct karyotypes of several oat accessions. Polymorphisms for C-bands among oat cultivars and wild accessions are described. These banding polymorphisms can be used to trace introgression of chromosomes from wild sources and to provide physical markers on the genetic map for oat. Although C-banding permits the identification of likely C-genome chromosomes based on comparisons with C-banding patterns in diploid and tetraploid Avena species, the A and D genomes cannot be readily differentiated based on their banding patterns.

Characterization of the hexaploid oat Avena byzantina cv. Kanota monosomic series using C-banding and RFLPs.

The establishment of a C-banded karyotype of hexaploid oat (Avena spp., 2n = 6x = 42) has facilitated the cytological characterization of a monosomic series in 'Kanota', an A. byzantina (C. Koch) cultivar. The 'Kanota' series of monosomics analyzed in this study consists of only 12 of the 21 different chromosome-deficient lines possible plus potential translocated segments of two or three additional chromosomes. These findings were confirmed by RFLP mapping data from studies in which oat probes were assigned to syntenic groups using the 'Kanota' set of monosomic lines. Among the remaining nine monosomic lines analyzed, eight are missing chromosomes represented in the set of 12 unique lines and one line, monosomic K13, is missing a chromosome from the unique set of 12 that possesses a cytologically detectable translocation. This same translocation, involving chromosomes 7C and 14, is found in 5 of the 21 'Kanota' monosomics. The incompleteness of the set of 'Kanota' monosomics might be due to (i) difficulty in identifying individual oat chromosomes without C-banding, (ii) plant genotypic and phenotypic variability in the original source population of the 'Kanota' monosomics, and (or) (iii) a high frequency of monosomic shifts in progency of the original 'Kanota' monosomic lines.