Dissection of complicate genetic architecture and breeding perspective of cottonseed traits by genome-wide association study.

BACKGROUND: Cottonseed is one of the most important raw materials for plant protein, oil and alternative biofuel for diesel engines. Understanding the complex genetic basis of cottonseed traits is requisite for achieving efficient genetic improvement of the traits. However, it is not yet clear about their genetic architecture in genomic level. GWAS has been an effective way to explore genetic basis of quantitative traits in human and many crops. This study aims to dissect genetic mechanism seven cottonseed traits by a GWAS for genetic improvement. RESULTS: A genome-wide association study (GWAS) based on a full gene model with gene effects as fixed and gene-environment interaction as random, was conducted for protein, oil and 5 fatty acids using 316 accessions and ~ 390 K SNPs. Totally, 124 significant quantitative trait SNPs (QTSs), consisting of 16, 21, 87 for protein, oil and fatty acids (palmitic, linoleic, oleic, myristic, stearic), respectively, were identified and the broad-sense heritability was estimated from 71.62 to 93.43%; no QTS-environment interaction was detected for the protein, the palmitic and the oleic contents; the protein content was predominantly controlled by epistatic effects accounting for 65.18% of the total variation, but the oil content and the fatty acids except the palmitic were mainly determined by gene main effects and no epistasis was detected for the myristic and the stearic. Prediction of superior pure line and hybrid revealed the potential of the QTSs in the improvement of cottonseed traits, and the hybrid could achieve higher or lower genetic values compared with pure lines. CONCLUSIONS: This study revealed complex genetic architecture of seven cottonseed traits at whole genome-wide by mixed linear model approach; the identified genetic variants and estimated genetic component effects of gene, gene-gene and gene-environment interaction provide cotton geneticist or breeders new knowledge on the genetic mechanism of the traits and the potential molecular breeding design strategy.

Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements.

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F2 individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

Updated foxtail millet genome assembly and gene mapping of nine key agronomic traits by resequencing a RIL population.

Foxtail millet (Setaria italica) provides food and fodder in semi-arid regions and infertile land. Resequencing of 184 foxtail millet recombinant inbred lines (RILs) was carried out to aid essential research on foxtail millet improvement. A total 483 414 single nucleotide polymorphisms were determined. Bin maps were constructed based on the RILs' recombination data. Based on the high-density bin map, we updated Zhanggu reference with 416 Mb after adding 16 Mb unanchored scaffolds and Yugu reference with some assembly error correction and 3158 gaps filled. Quantitative trait loci (QTL) mapping of nine agronomic traits was done based on this RIL population, five of which were controlled by a single gene. Meanwhile, two QTLs were found for plant height, and a candidate gene showed 89% identity to the known rice gibberellin-synthesis gene sd1. Three QTLs were found for the trait of heading date. The whole genome resequencing and QTL mapping provided important tools for foxtail millet research and breeding. Resequencing of the RILs could also provide an effective way for high-quality genome assembly and gene identification.

Updated sesame genome assembly and fine mapping of plant height and seed coat color QTLs using a new high-density genetic map.

BACKGROUND: Sesame is an important high-quality oil seed crop. The sesame genome was de novo sequenced and assembled in 2014 (version 1.0); however, the number of anchored pseudomolecules was higher than the chromosome number (2n = 2x = 26) due to the lack of a high-density genetic map with 13 linkage groups. RESULTS: We resequenced a permanent population consisting of 430 recombinant inbred lines and constructed a genetic map to improve the sesame genome assembly. We successfully anchored 327 scaffolds onto 13 pseudomolecules. The new genome assembly (version 2.0) included 97.5 % of the scaffolds greater than 150 kb in size present in assembly version 1.0 and increased the total pseudomolecule length from 233.7 to 258.4 Mb with 94.3 % of the genome assembled and 97.2 % of the predicted gene models anchored. Based on the new genome assembly, a bin map including 1,522 bins spanning 1090.99 cM was generated and used to identified 41 quantitative trait loci (QTLs) for sesame plant height and 9 for seed coat color. The plant height-related QTLs explained 3-24 % the phenotypic variation (mean value, 8 %), and 29 of them were detected in at least two field trials. Two major loci (qPH-8.2 and qPH-3.3) that contributed 23 and 18 % of the plant height were located in 350 and 928-kb spaces on Chr8 and Chr3, respectively. qPH-3.3, is predicted to be responsible for the semi-dwarf sesame plant phenotype and contains 102 candidate genes. This is the first report of a sesame semi-dwarf locus and provides an interesting opportunity for a plant architecture study of the sesame. For the sesame seed coat color, the QTLs of the color spaces L*, a*, and b* were detected with contribution rates of 3-46 %. qSCb-4.1 contributed approximately 39 % of the b* value and was located on Chr4 in a 199.9-kb space. A list of 32 candidate genes for the locus, including a predicted black seed coat-related gene, was determined by screening the newly anchored genome. CONCLUSIONS: This study offers a high-density genetic map and an improved assembly of the sesame genome. The number of linkage groups and pseudomolecules in this assembly equals the number of sesame chromosomes for the first time. The map and updated genome assembly are expected to serve as a platform for future comparative genomics and genetic studies.

Genetic analysis for rice grain quality traits in the YVB stable variant line using RAD-seq.

The future of rice breeding will likely be built on the basis of the further utilization of heterosis between elite cultivars and genetic resources from distant subspecies of rice. Previous studies have proved that exogenous genomic DNA transformation methods can be used to transfer genetic information from distant relatives (donor) into cultivated rice (recipient). However, the mechanism underlying this form of genetic transfer is poorly characterized, and the genes that cause the phenotypic changes in these variants are typically difficult to identify. This study examined YVB, a stable variant line with greatly improved grain quality traits that was derived from an indica variety (V20B) by transferring genomic DNA of O.minuta through the "spike-stalk injection method (SIM)". We used restriction-site associated DNA sequencing technology (RAD-seq) to evaluate a population of BC1F5 backcross lines (YVB × V20B); the RAD-seq data were used to construct a genetic linkage map with high-density SNPs for use in association analysis exploring genotype-phenotype relationships at the whole-genome level. A total of 17 quantitative trait loci (QTLs) for rice quality traits were mapped to chromosomes 3, 5, 6, 8, and 9. 8 major QTLs controlling different phenotypic variations were mapped to the same region of chromosome 5. This region contained the GS5 gene for grain weight and the qSW5/GW5 gene for grain width. This study provides new resources and insights into the molecular mechanisms of grain trait phenotypic variation and the transmission of genetic information via the introduction of genomic DNA to a distantly related crop relative species.

A precise chloroplast genome of Nelumbo nucifera (Nelumbonaceae) evaluated with Sanger, Illumina MiSeq, and PacBio RS II sequencing platforms: insight into the plastid evolution of basal eudicots.

BACKGROUND: The chloroplast genome is important for plant development and plant evolution. Nelumbo nucifera is one member of relict plants surviving from the late Cretaceous. Recently, a new sequencing platform PacBio RS II, known as 'SMRT (Single Molecule, Real-Time) sequencing', has been developed. Using the SMRT sequencing to investigate the chloroplast genome of N. nucifera will help to elucidate the plastid evolution of basal eudicots. RESULTS: The sizes of the de novo assembled complete chloroplast genome of N. nucifera were 163,307 bp, 163,747 bp and 163,600 bp with average depths of coverage of 7×, 712× and 105× sequenced by Sanger, Illumina MiSeq and PacBio RS II, respectively. The precise chloroplast genome of N. nucifera was obtained from PacBio RS II data proofread by Illumina MiSeq reads, with a quadripartite structure containing a large single copy region (91,846 bp) and a small single copy region (19,626 bp) separated by two inverted repeat regions (26,064 bp). The genome contains 113 different genes, including four distinct rRNAs, 30 distinct tRNAs and 79 distinct peptide-coding genes. A phylogenetic analysis of 133 taxa from 56 orders indicated that Nelumbo with an age of 177 million years is a sister clade to Platanus, which belongs to the basal eudicots. Basal eudicots began to emerge during the early Jurassic with estimated divergence times at 197 million years using MCMCTree. IR expansions/contractions within the basal eudicots seem to have occurred independently. CONCLUSIONS: Because of long reads and lack of bias in coverage of AT-rich regions, PacBio RS II showed a great promise for highly accurate 'finished' genomes, especially for a de novo assembly of genomes. N. nucifera is one member of basal eudicots, however, evolutionary analyses of IR structural variations of N. nucifera and other basal eudicots suggested that IR expansions/contractions occurred independently in these basal eudicots or were caused by independent insertions and deletions. The precise chloroplast genome of N. nucifera will present new information for structural variation of chloroplast genomes and provide new insight into the evolution of basal eudicots at the primary sequence and structural level.

Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing.

We explored genetic variation by sequencing a selection of 84 tomato accessions and related wild species representative of the Lycopersicon, Arcanum, Eriopersicon and Neolycopersicon groups, which has yielded a huge amount of precious data on sequence diversity in the tomato clade. Three new reference genomes were reconstructed to support our comparative genome analyses. Comparative sequence alignment revealed group-, species- and accession-specific polymorphisms, explaining characteristic fruit traits and growth habits in the various cultivars. Using gene models from the annotated Heinz 1706 reference genome, we observed differences in the ratio between non-synonymous and synonymous SNPs (dN/dS) in fruit diversification and plant growth genes compared to a random set of genes, indicating positive selection and differences in selection pressure between crop accessions and wild species. In wild species, the number of single-nucleotide polymorphisms (SNPs) exceeds 10 million, i.e. 20-fold higher than found in most of the crop accessions, indicating dramatic genetic erosion of crop and heirloom tomatoes. In addition, the highest levels of heterozygosity were found for allogamous self-incompatible wild species, while facultative and autogamous self-compatible species display a lower heterozygosity level. Using whole-genome SNP information for maximum-likelihood analysis, we achieved complete tree resolution, whereas maximum-likelihood trees based on SNPs from ten fruit and growth genes show incomplete resolution for the crop accessions, partly due to the effect of heterozygous SNPs. Finally, results suggest that phylogenetic relationships are correlated with habitat, indicating the occurrence of geographical races within these groups, which is of practical importance for Solanum genome evolution studies.

The Tarenaya hassleriana genome provides insight into reproductive trait and genome evolution of crucifers.

The Brassicaceae, including Arabidopsis thaliana and Brassica crops, is unmatched among plants in its wealth of genomic and functional molecular data and has long served as a model for understanding gene, genome, and trait evolution. However, genome information from a phylogenetic outgroup that is essential for inferring directionality of evolutionary change has been lacking. We therefore sequenced the genome of the spider flower (Tarenaya hassleriana) from the Brassicaceae sister family, the Cleomaceae. By comparative analysis of the two lineages, we show that genome evolution following ancient polyploidy and gene duplication events affect reproductively important traits. We found an ancient genome triplication in Tarenaya (Th-α) that is independent of the Brassicaceae-specific duplication (At-α) and nested Brassica (Br-α) triplication. To showcase the potential of sister lineage genome analysis, we investigated the state of floral developmental genes and show Brassica retains twice as many floral MADS (for minichromosome maintenance1, AGAMOUS, DEFICIENS and serum response factor) genes as Tarenaya that likely contribute to morphological diversity in Brassica. We also performed synteny analysis of gene families that confer self-incompatibility in Brassicaceae and found that the critical serine receptor kinase receptor gene is derived from a lineage-specific tandem duplication. The T. hassleriana genome will facilitate future research toward elucidating the evolutionary history of Brassicaceae genomes.

The sacred lotus genome provides insights into the evolution of flowering plants.

Sacred lotus (Nelumbo nucifera) is an ornamental plant that is also used for food and medicine. This basal eudicot species is especially important from an evolutionary perspective, as it occupies a critical phylogenetic position in flowering plants. Here we report the draft genome of a wild strain of sacred lotus. The assembled genome is 792 Mb, which is approximately 85-90% of genome size estimates. We annotated 392 Mb of repeat sequences and 36,385 protein-coding genes within the genome. Using these sequence data, we constructed a phylogenetic tree and confirmed the basal location of sacred lotus within eudicots. Importantly, we found evidence for a relatively recent whole-genome duplication event; any indication of the ancient paleo-hexaploid event was, however, absent. Genomic analysis revealed evidence of positive selection within 28 embryo-defective genes and one annexin gene that may be related to the long-term viability of sacred lotus seed. We also identified a significant expansion of starch synthase genes, which probably elevated starch levels within the rhizome of sacred lotus. Sequencing this strain of sacred lotus thus provided important insights into the evolution of flowering plant and revealed genetic mechanisms that influence seed dormancy and starch synthesis.

Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation.

About 8,000 years ago in the Fertile Crescent, a spontaneous hybridization of the wild diploid grass Aegilops tauschii (2n = 14; DD) with the cultivated tetraploid wheat Triticum turgidum (2n = 4x = 28; AABB) resulted in hexaploid wheat (T. aestivum; 2n = 6x = 42; AABBDD). Wheat has since become a primary staple crop worldwide as a result of its enhanced adaptability to a wide range of climates and improved grain quality for the production of baker's flour. Here we describe sequencing the Ae. tauschii genome and obtaining a roughly 90-fold depth of short reads from libraries with various insert sizes, to gain a better understanding of this genetically complex plant. The assembled scaffolds represented 83.4% of the genome, of which 65.9% comprised transposable elements. We generated comprehensive RNA-Seq data and used it to identify 43,150 protein-coding genes, of which 30,697 (71.1%) were uniquely anchored to chromosomes with an integrated high-density genetic map. Whole-genome analysis revealed gene family expansion in Ae. tauschii of agronomically relevant gene families that were associated with disease resistance, abiotic stress tolerance and grain quality. This draft genome sequence provides insight into the environmental adaptation of bread wheat and can aid in defining the large and complicated genomes of wheat species.

Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential.

Foxtail millet (Setaria italica), a member of the Poaceae grass family, is an important food and fodder crop in arid regions and has potential for use as a C(4) biofuel. It is a model system for other biofuel grasses, including switchgrass and pearl millet. We produced a draft genome (∼423 Mb) anchored onto nine chromosomes and annotated 38,801 genes. Key chromosome reshuffling events were detected through collinearity identification between foxtail millet, rice and sorghum including two reshuffling events fusing rice chromosomes 7 and 9, 3 and 10 to foxtail millet chromosomes 2 and 9, respectively, that occurred after the divergence of foxtail millet and rice, and a single reshuffling event fusing rice chromosome 5 and 12 to foxtail millet chromosome 3 that occurred after the divergence of millet and sorghum. Rearrangements in the C(4) photosynthesis pathway were also identified.

Development of novel EST-SSRs from sacred lotus (Nelumbo nucifera Gaertn) and their utilization for the genetic diversity analysis of N. nucifera.

Expressed sequence tags (ESTs) provide a valuable resource for the development of simple sequence repeat (SSR) or microsatellite markers. This study identified SSRs within ESTs from Nelumbo nucifera (lotus or sacred lotus), developed markers from them, and assessed the potential of those markers for diversity analysis. Within 2207 ESTs from N. nucifera downloaded from GenBank, 1483 unigenes (303 contigs and 1180 singletons) were identified. After eliminating for redundancy, 125 SSR-containing ESTs were derived, and 71 unique SSRs were detected with an average density of one SSR per 13.04 kb. Dinucleotide repeats were the dominant motif in N. nucifera, whereas the sequences AG/TC/GA/CT, AAG/TTC/GAT/AGA, and AAAGCC were the most frequent of di-, tri-, and hexanucleotide motifs, respectively. The AG/TC (40.85%) and AAG (5.63%) motifs were predominant for the di- and trinucleotide repeats, respectively. Sixty-two SSR-containing ESTs were suitable for primer design. From these sequences, 23 EST-SSR markers were developed and were applied to 39 cultivated varieties of N. nucifera, 10 accessions of wild N. nucifera, and 1 accession of Nelumbo lutea (American lotus). Genetic diversity and genetic relationships were examined by constructing unweighted pair-group method with arithmetic average dendrograms and principal coordinates analysis plots based on SSR polymorphisms. Results indicated genetic differentiation between cultivated and wild lotus and between seed lotus cultivars and rhizome lotus cultivars. These EST-SSR markers will be useful for further studies of the evolution and diversity of Nelumbo.

Polymorphic microsatellite markers in Euryale ferox Salisb. (Nymphaeaceae).

Eleven polymorphic microsatellite markers were isolated and identified in the aquatic plant Euryale ferox Salisb. (Nymphaeaceae). This species, which belongs to basal Magnoliophyta, reproduces sexually. All of these 11 microsatellite markers yielded 25 alleles in a survey of a wild population of 34 individuals. Two or three alleles per locus were detected, with expected heterozygosity ranging from 0.056 to 0.634 and observed heterozygosity from 0.000 to 0.088. These simple sequence repeat markers will be useful for evaluating the genetic structure of the E. ferox population in the future.

Sixteen polymorphic microsatellite markers from Zizania latifolia Turcz. (Poaceae).

Sixteen polymorphic microsatellite markers were isolated and identified in Zizania latifolia Turcz. (Poaceae), a perennial aquatic plant widespread in Eastern Asia. The microsatellite-enriched library was constructed using the fast isolation by AFLP of sequences containing repeats method. These markers revealed two to 14 alleles, with an average of 5.6 alleles per locus. The observed and expected heterozygosities varied from 0.071 to 0.690 and from 0.174 to 0.812, respectively. These markers will be useful for studying of gene flow and evaluating the genetic diversity of the Zizania latifolia population.