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.

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.

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.

Genetic and Sequence Analysis of Markers Tightly Linked to the Soybean mosaic virus Resistance Gene, Rsv3.

Soybean mosaic virus (SMV) is a major viral pathogen, affecting soybean [Glycine max (L.) Merr.] production worldwide. The Rsv3 gene of soybean confers resistance to three of the most virulent strains (G5-G7) of SMV. The objectives of this study were to map Rsv3 and develop polymerase chain reaction (PCR) based markers for marker-assisted selection (MAS) purposes. Disease-response data were collected from two F(2) mapping populations, L29 (Rsv3) x Lee68 (rsv3) and Tousan 140 (Rsv3) x Lee68 (rsv3). Bulk segregant analysis based on amplified fragment length polymorphism (AFLP) markers demonstrated that the Rsv3 locus maps to the soybean molecular linkage group (MLG) B2 between restriction fragment length polymorphism (RFLP) markers A519 and Mng247. These two tightly linked RFLP markers were converted to PCR-based markers to expedite MAS. Sequence analysis of the Mng247 genomic region revealed similarity to the consensus sequence of a leucine-rich repeat (LRR) characteristic of the extracellular LRR class of disease resistance genes. Results from this study will be useful in pyramiding viral resistance genes and in cloning the Rsv3 gene.

Processing and localization of bovine beta-casein expressed in transgenic soybean seeds under control of a soybean lectin expression cassette.

We have examined the processing and subcellular localization of a chimeric gene consisting of the bovine milk protein, beta-casein, under the control of a soybean seed lectin promoter and its 32 amino acid signal sequence in the seeds of transgenic soybean plants. The beta-casein expressed in developing soybean seeds is a doublet with apparent molecular weight slightly smaller than the bovine beta-casein and expression of the protein was highest in immature cotyledons. The casein proteins were purified from the immature soybean seeds by immunoaffinity chromatography and were analyzed by two-dimensional gel electrophoresis, blotting, and amino terminal sequencing. The N-terminal sequences of both of the doublet soybean casein polypeptides were identical to the N-terminal sequence of the bovine beta-casein indicating that the 32 amino acid lectin signal sequence was cleaved precisely from the chimeric protein in developing soybean seeds. Analysis of the purified soybean beta-casein polypeptides by mass spectrometry (MALDI-MS) showed that they are not phosphorylated. Absence of added phosphate groups is the cause of the size difference between the soybean beta-casein and native bovine beta-casein protein. Immunolocalization experiments showed that the casein protein was found in the protein storage vacuoles (PSV) in developing and mature soybean seeds. The precise removal of the 32 amino acid lectin amino terminal sequence from the chimeric lectin-casein fusion suggests that the lectin expression cassette can be used for production of pharmaceutical or other recombinant proteins of added value in the developing soybean seed.

Copia-like retrotransposons in rice: sequence heterogeneity, species distribution and chromosomal locations.

Degenerated oligonucleotide primers were used to amplify, clone, and analyze sequence heterogeneity and chromosomal distribution of 23 PCR fragments corresponding to the reverse transcriptase domain of copia-like retrotransposons in rice. Of the 23 fragments 22 could be aligned by their deduced amino acid sequences and were divided into 6 groups according to the phylogenetic and Southern blot analyses. Amino acid sequence differences among the 22 aligned fragments ranged from 1 to 64%. Southern blot analysis of 10 rice accessions including indica, japonica and common wild rice, using these 23 fragments as probes, showed that copia-like retrotransposons were present in moderate to high copy numbers in all the rice genome although the exact copy number cannot be determined. The major difference revealed by southern analysis is a differentiation between the four indica varieties as one group and the four japonica varieties and the two wild rice accessions as another group. Polymorphisms were also detected among the indica and japonica varieties by major bands and repeatable minor bands. Five hybridization bands were mapped to chromosomes 3, 4, 8, and 9, respectively. All the five bands were inherited in a dominant Mendelian fashion and were not allelic with each other, indicating that the same element did not reside on the same location in different rice accessions. No transcript of the copia-like reverse transcriptase was detected on northern blot. The results suggest that the sequence heterogeneity and distributional variability of retrotransposons may be one of contributory factors causing genetic diversity in rice.

Amplified fragment length polymorphism (AFLP) in soybean: species diversity, inheritance, and near-isogenic line analysis.

Amplified fragment length polymorphism (AFLP) analysis is a PCR-based technique capable of detecting more than 50 independent loci in a single PCR reaction. The objectives of the present study were to: (1) assess the extent of AFLP variation in cultivated (Gycine max L. Merr.) and wild soybean (G. soja Siebold & Zucc.), (2) determine genetic relationships among soybean accessions using AFLP data, and (3) evaluate the usefulness of AFLPs as genetic markers. Fifteen AFLP primer pairs detected a total of 759 AFLP fragments in a sample of 23 accessions of wild and cultivated soybean, with an average of 51 fragments produced per primer pair per accession. Two-hundred and seventy four fragments (36% of the total observed) were polymorphic, among which 127 (17%) were polymorphic in G. max and 237 (31%) were polymorphic in G. soja. F2 segregation analysis of six AFLP fragments indicated that they segregate as stable Mendelian loci. The number of polymorphic loci detected per AFLP primer pair in a sample of 23 accessions ranged from 9 to 27. The AFLP phenotypic diversity values were greater in wild than in cultivated soybean. Cluster and principal component analyses using AFLP data clearly separated G. max and G. soja accessions. Within the G. max group, adapted soybean cultivars were tightly clustered, illustrating the relatively low genetic diversity present in cultivated soybean. AFLP analysis of four soybean near-isogenic lines (NILs) identified three AFLP markers putatively linked to a virus resistance gene from two sources. The capacity of AFLP analysis to detect thousands of independent genetic loci with minimal cost and time requirements makes them an ideal marker for a wide array of genetic investigations.

Molecular-marker analysis of seed-weight: genomic locations, gene action, and evidence for orthologous evolution among three legume species.

The objectives of this study were to use molecular markers to: (1) identify quantitative trait loci (QTL) controlling seed-weight in soybean, (2) characterize the genetic basis of seed-weight expression, and (3) determine whether soybean shares orthologous seed-weight genes with cowpea and/or mung bean. An F2 population was developed between a large-seeded Glycine max breeding line and a small-seeded G. soja plant introduction. DNA samples from 150 F2 individuals were analyzed with 91 polymorphic genetic markers, including RFLPs, RAPDs and SSRs. Seed-weight was analyzed by randomly sampling 100 seeds from each of 150 greenhouse-grown F2 individuals, and their 150 F2∶3 lines, from a replicated field trial. Markers associated with seed-weight were identified using the computer program MapMaker-QTL and a one-way analysis of variance. Three and five markers were significantly associated with seed-weight variation (P<0.01) in the F2 and F2∶3 generations, respectively. Tests for digenic epistasis revealed three significant interactions in both generations. In a combined analysis, these markers and interactions explained 50 and 60% of the phenotypic variation for seed-weight in the F2 and F2∶3 generations, respectively. Comparison of our results in soybean (Glycine) with those previously reported in cowpea and mung bean (Vigna) indicated that soybean and cowpea share an orthologous seed-weight gene. In both species, a genomic region significantly associated with seed-weight spanned the same RFLP markers in the same linkage order. A significant digenic interaction involving this genomic region was conserved in all three species. These results suggest that the exploitation of "comparative QTL mapping" is an invaluable tool for quantitative geneticists working with poorly characterized plant systems.

Analysis of the barley and rice genomes by comparative RFLP linkage mapping.

Comparative genetic mapping of rice and barley, both major crop species with extensive genetic resources, offers the possibility of uniting two well-established and characterized genetic systems. In the present study, we screened 229 molecular markers and utilized 110 polymorphic orthologous loci to construct comparative maps of the rice and barley genomes. While extensive chromosomal rearrangements, including inversions and intrachromosomal translocations, differentiate the rice and barley genomes, several syntenous chromosomes are evident. Indeed, several chromosomes and chromosome arms appear to share nearly identical gene content and gene order. Seventeen regions of conserved organization were detected, spanning 287 cM (24%) and 321 cM (31%) of the rice and barley genomes, respectively. The results also indicate that most (72%) of the single-copy sequences in barley are also single copy in rice, suggesting that the large barley genome arose by unequal crossing over and amplification of repetitive DNA sequences and not by the duplication of single-copy sequences. Combining these results with those previously reported for comparative analyses of rice and wheat identified nine putatively syntenous chromosomes among barley, wheat and rice. The high degree of gene-order conservation as detected by comparative mapping has astonishing implications for interpreting genetic information among species and for elucidating chromosome evolution and speciation.

Microsatellite and amplified sequence length polymorphisms in cultivated and wild soybean.

The objectives of this study were to (i) assess the extent of genetic variation in soybean microsatellites (simple sequence repeats or SSRs), (ii) assay for amplified sequence length polymorphisms (ASLPs), and (iii) evaluate the usefulness of SSRs and ASLPs as genetic markers. Five microsatellites detected a total of 79 variants (alleles) in a sample of 94 accessions of wild (Glycine soja) and cultivated soybean (G. max). F2 segregation analysis of four of the five microsatellites identified these variants (alleles) with four loci located in independent linkage groups. The number of alleles per microsatellite locus ranged from 5 to 21; to our knowledge these are the largest numbers of alleles for single Mendelian loci reported in soybean. Allelic diversity for the SSR loci was greater in wild than in cultivated soybean. Overall, 43 more SSR alleles were detected in wild than in cultivated soybean. These results indicate that SSRs are the marker of choice, especially for species with low levels of variation as detected by other types of markers. Two alleles were detected at each of the three ASLP loci examined. A total of six ASLP alleles were observed in cultivated soybean and five were observed in wild soybean; all alleles detected in wild soybean were present in cultivated soybean. Allelic diversity values for the ASLP loci were near previous estimates for restriction fragment length polymorphisms and therefore ASLPs may be useful as genetic markers in site-directed mapping.