Heterozygous Genome Research Articles

Genome survey sequencing of purple elephant grass (Pennisetum purpureum Schum ‘Zise’) and identification of its SSR markers

Elephant grass (Pennisetum purpureum) is a perennial grass in the Poaceae family with high tolerance and one of the best forage plants. Despite its economic importance, the inheritance information of P. purpureum has remained largely unknown. To obtain the whole reference genome, we first conducted a genome survey of P. purpureum. Next-generation sequencing (NGS) was used to perform the de novo whole genome sequencing. As a result, the estimated genome size of elephant grass was 2.01 Gb, with 71.36% repetitive elements. The heterozygosity was 1.02%, which indicates a highly heterozygous genome. The retroelements (9.36%) were the most repetitive elements, followed by DNA transposons (3.66%). In the meantime, 83,706 high-quality genomic simple sequence repeat (SSR) markers, in which the greatest SSR unit length was 3, were developed. Thirty pairs of SSR markers were randomly selected to verify the efficiency and all of them yielded clear amplification products, among which 28 pairs (93.3%) of the primers showed polymorphism. The genome data obtained in this research provided a large amount of gene resources for further investigating Pennisetum species.

Hybridity has a greater effect than paternal genome dosage on heterosis in sugar beet (Beta vulgaris)

BackgroundThe phenomenon of heterosis is critical to plant breeding and agricultural productivity. Heterosis occurs when F1 hybrid offspring display quantitative improvements in traits to levels that do not occur in the parents. Increasing the genome dosage (i.e. ploidy level) of F1 offspring can contribute to heterosis effects. Sugar beet (Beta vulgaris) provides a model for investigating the relative effects of genetic hybridity and genome dosage on heterosis. Sugar beet lines of different ploidy levels were crossed to generate diploid and triploid F1 offspring to investigate the effect of; (1) paternal genome dosage increase on F1 heterosis, and; (2) homozygous versus heterozygous tetraploid male parents on F1 triploid heterosis. A range of traits of agronomic and commercial importance were analyzed for the extent of heterosis effects observed in the F1 offspring.ResultsComparisons of parental lines to diploid (EA, EB) and triploid (EAA, EBB) F1 hybrids for total yield, root yield, and sugar yield indicated that there was no effect of paternal genome dosage increases on heterosis levels, indicating that hybridity is the main contributor to the heterosis levels observed. For all traits measured (apart from seed viability), F1 triploid hybrids derived from heterozygous tetraploid male parents displayed equivalent levels of heterosis as F1 triploid hybrids generated with homozygous tetraploid male parents, suggesting that heterosis gains in F1 triploids do not arise by simply increasing the extent of multi-locus heterozygosity in sugar beet F1 offspring.ConclusionsOverall, our study indicates that; (1) increasing the paternal genome dosage does not enhance heterosis in F1 hybrids, and; (2) increasing multi-locus heterozygosity using highly heterozygous paternal genomes to generate F1 triploid hybrids does not enhance heterosis. Our findings have implications for the design of future F1 hybrid improvement programs for sugar beet.

Assembly of the Boechera retrofracta Genome and Evolutionary Analysis of Apomixis-Associated Genes.

Closely related to the model plant Arabidopsis thaliana, the genus Boechera is known to contain both sexual and apomictic species or accessions. Boechera retrofracta is a diploid sexually reproducing species and is thought to be an ancestral parent species of apomictic species. Here we report the de novo assembly of the B. retrofracta genome using short Illumina and Roche reads from 1 paired-end and 3 mate pair libraries. The distribution of 23-mers from the paired end library has indicated a low level of heterozygosity and the presence of detectable duplications and triplications. The genome size was estimated to be equal 227 Mb. N50 of the assembled scaffolds was 2.3 Mb. Using a hybrid approach that combines homology-based and de novo methods 27,048 protein-coding genes were predicted. Also repeats, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were annotated. Finally, genes of B. retrofracta and 6 other Brassicaceae species were used for phylogenetic tree reconstruction. In addition, we explored the histidine exonuclease APOLLO locus, related to apomixis in Boechera, and proposed model of its evolution through the series of duplications. An assembled genome of B. retrofracta will help in the challenging assembly of the highly heterozygous genomes of hybrid apomictic species.

Draft genome sequence of ramie, Boehmeria nivea (L.) Gaudich.

Ramie, Boehmeria nivea (L.) Gaudich, family Urticaceae, is a plant native to eastern Asia, and one of the world's oldest fibre crops. It is also used as animal feed and for the phytoremediation of heavy metal-contaminated farmlands. Thus, the genome sequence of ramie was determined to explore the molecular basis of its fibre quality, protein content and phytoremediation. For further understanding ramie genome, different paired-end and mate-pair libraries were combined to generate 134.31Gb of raw DNA sequences using the Illumina whole-genome shotgun sequencing approach. The highly heterozygous B. nivea genome was assembled using the Platanus Genome Assembler, which is an effective tool for the assembly of highly heterozygous genome sequences. The final length of the draft genome of this species was approximately 341.9Mb (contig N50=22.62kb, scaffold N50=1,126.36kb). Based on ramie genome annotations, 30,237 protein-coding genes were predicted, and the repetitive element content was 46.3%. The completeness of the final assembly was evaluated by benchmarking universal single-copy orthologous genes (BUSCO); 90.5% of the 1,440 expected embryophytic genes were identified as complete, and 4.9% were identified as fragmented. Phylogenetic analysis based on single-copy gene families and one-to-one orthologous genes placed ramie with mulberry and cannabis, within the clade of urticalean rosids. Genome information of ramie will be a valuable resource for the conservation of endangered Boehmeria species and for future studies on the biogeography and characteristic evolution of members of Urticaceae.

Comparison of three assembly strategies for a heterozygous seedless grapevine genome assembly

BackgroundDe novo heterozygous assembly is an ongoing challenge requiring improved assembly approaches. In this study, three strategies were used to develop de novo Vitis vinifera ‘Sultanina’ genome assemblies for comparison with the inbred V. vinifera (PN40024 12X.v2) reference genome and a published Sultanina ALLPATHS-LG assembly (AP). The strategies were: 1) a default PLATANUS assembly (PLAT_d) for direct comparison with AP assembly, 2) an iterative merging strategy using METASSEMBLER to combine PLAT_d and AP assemblies (MERGE) and 3) PLATANUS parameter modifications plus GapCloser (PLAT*_GC).ResultsThe three new assemblies were greater in size than the AP assembly. PLAT*_GC had the greatest number of scaffolds aligning with a minimum of 95% identity and ≥1000 bp alignment length to V. vinifera (PN40024 12X.v2) reference genome. SNP analysis also identified additional high quality SNPs. A greater number of sequence reads mapped back with zero-mismatch to the PLAT_d, MERGE, and PLAT*_GC (>94%) than was found in the AP assembly (87%) indicating a greater fidelity to the original sequence data in the new assemblies than in AP assembly. A de novo gene prediction conducted using seedless RNA-seq data predicted > 30,000 coding sequences for the three new de novo assemblies, with the greatest number (30,544) in PLAT*_GC and only 26,515 for the AP assembly. Transcription factor analysis indicated good family coverage, but some genes found in the VCOST.v3 annotation were not identified in any of the de novo assemblies, particularly some from the MYB and ERF families.ConclusionsThe PLAT_d and PLAT*_GC had a greater number of synteny blocks with the V. vinifera (PN40024 12X.v2) reference genome than AP or MERGE. PLAT*_GC provided the most contiguous assembly with only 1.2% scaffold N, in contrast to AP (10.7% N), PLAT_d (6.6% N) and Merge (6.4% N). A PLAT*_GC pseudo-chromosome assembly with chromosome alignment to the reference genome V. vinifera, (PN40024 12X.v2) provides new information for use in seedless grape genetic mapping studies. An annotated de novo gene prediction for the PLAT*_GC assembly, aligned with VitisNet pathways provides new seedless grapevine specific transcriptomic resource that has excellent fidelity with the seedless short read sequence data.

Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F1 derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F1 sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

Scarless Genome Editing of Human Pluripotent Stem Cells via Transient Puromycin Selection.

SummaryGenome-edited human pluripotent stem cells (hPSCs) have broad applications in disease modeling, drug discovery, and regenerative medicine. We present and characterize a robust method for rapid, scarless introduction or correction of disease-associated variants in hPSCs using CRISPR/Cas9. Utilizing non-integrated plasmid vectors that express a puromycin N-acetyl-transferase (PAC) gene, whose expression and translation is linked to that of Cas9, we transiently select for cells based on their early levels of Cas9 protein. Under optimized conditions, co-delivery with single-stranded donor DNA enabled isolation of clonal cell populations containing both heterozygous and homozygous precise genome edits in as little as 2 weeks without requiring cell sorting or high-throughput sequencing. Edited cells isolated using this method did not contain any detectable off-target mutations and displayed expected functional phenotypes after directed differentiation. We apply the approach to a variety of genomic loci in five hPSC lines cultured using both feeder and feeder-free conditions.

Draft Sequencing of the Heterozygous Diploid Genome of Satsuma (Citrus unshiu Marc.) Using a Hybrid Assembly Approach.

Satsuma (Citrus unshiu Marc.) is one of the most abundantly produced mandarin varieties of citrus, known for its seedless fruit production and as a breeding parent of citrus. De novo assembly of the heterozygous diploid genome of Satsuma (“Miyagawa Wase”) was conducted by a hybrid assembly approach using short-read sequences, three mate-pair libraries, and a long-read sequence of PacBio by the PLATANUS assembler. The assembled sequence, with a total size of 359.7 Mb at the N50 length of 386,404 bp, consisted of 20,876 scaffolds. Pseudomolecules of Satsuma constructed by aligning the scaffolds to three genetic maps showed genome-wide synteny to the genomes of Clementine, pummelo, and sweet orange. Gene prediction by modeling with MAKER-P proposed 29,024 genes and 37,970 mRNA; additionally, gene prediction analysis found candidates for novel genes in several biosynthesis pathways for gibberellin and violaxanthin catabolism. BUSCO scores for the assembled scaffold and predicted transcripts, and another analysis by BAC end sequence mapping indicated the assembled genome consistency was close to those of the haploid Clementine, pummel, and sweet orange genomes. The number of repeat elements and long terminal repeat retrotransposon were comparable to those of the seven citrus genomes; this suggested no significant failure in the assembly at the repeat region. A resequencing application using the assembled sequence confirmed that both kunenbo-A and Satsuma are offsprings of Kishu, and Satsuma is a back-crossed offspring of Kishu. These results illustrated the performance of the hybrid assembly approach and its ability to construct an accurate heterozygous diploid genome.

Genome Dynamics of Hybrid Saccharomyces cerevisiae During Vegetative and Meiotic Divisions.

Mutation and recombination are the major sources of genetic diversity in all organisms. In the baker’s yeast, all mutation rate estimates are in homozygous background. We determined the extent of genetic change through mutation and loss of heterozygosity (LOH) in a heterozygous Saccharomyces cerevisiae genome during successive vegetative and meiotic divisions. We measured genome-wide LOH and base mutation rates during vegetative and meiotic divisions in a hybrid (S288c/YJM789) S. cerevisiae strain. The S288c/YJM789 hybrid showed nearly complete reduction in heterozygosity within 31 generations of meioses and improved spore viability. LOH in the meiotic lines was driven primarily by the mating of spores within the tetrad. The S288c/YJM789 hybrid lines propagated vegetatively for the same duration as the meiotic lines, showed variable LOH (from 2 to 3% and up to 35%). Two of the vegetative lines with extensive LOH showed frequent and large internal LOH tracts that suggest a high frequency of recombination repair. These results suggest significant LOH can occur in the S288c/YJM789 hybrid during vegetative propagation presumably due to return to growth events. The average base substitution rates for the vegetative lines (1.82 × 10−10 per base per division) and the meiotic lines (1.22 × 10−10 per base per division) are the first genome-wide mutation rate estimates for a hybrid yeast. This study therefore provides a novel context for the analysis of mutation rates (especially in the context of detecting LOH during vegetative divisions), compared to previous mutation accumulation studies in yeast that used homozygous backgrounds.

Influence of genetic background of engineered xylose-fermenting industrial Saccharomyces cerevisiae strains for ethanol production from lignocellulosic hydrolysates.

An industrial ethanol-producing Saccharomyces cerevisiae strain with genes of fungal oxido-reductive pathway needed for xylose fermentation integrated into its genome (YRH1415) was used to obtain haploids and diploid isogenic strains. The isogenic strains were more effective in metabolizing xylose than YRH1415 strain and able to co-ferment glucose and xylose in the presence of high concentrations of inhibitors resulting from the hydrolysis of lignocellulosic biomass (switchgrass). The rate of xylose consumption did not appear to be affected by the ploidy of strains or the presence of two copies of the xylose fermentation genes but by heterozygosity of alleles for xylose metabolism in YRH1415. Furthermore, inhibitor tolerance was influenced by the heterozygous genome of the industrial strain, which also showed a marked influenced on tolerance to increasing concentrations of toxic compounds, such as furfural. In this work, selection of haploid derivatives was found to be a useful strategy to develop efficient xylose-fermenting industrial yeast strains.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination

BackgroundRoot and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected “orphan” crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop.ResultsTo accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage.ConclusionsGuinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

HaploMerger2: rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly.

SummaryDe novo assembly is a difficult issue for heterozygous diploid genomes. The advent of high-throughput short-read and long-read sequencing technologies provides both new challenges and potential solutions to the issue. Here, we present HaploMerger2 (HM2), an automated pipeline for rebuilding both haploid sub-assemblies from the polymorphic diploid genome assembly. It is designed to work on pre-existing diploid assemblies, which are typically created by using de novo assemblers. HM2 can process any diploid assemblies, but it is especially suitable for diploid assemblies with high heterozygosity (≥3%), which can be difficult for other tools. This pipeline also implements flexible and sensitive assembly error detection, a hierarchical scaffolding procedure and a reliable gap-closing method for haploid sub-assemblies. Using HM2, we demonstrate that two haploid sub-assemblies reconstructed from a real, highly-polymorphic diploid assembly show greatly improved continuity.Availability and ImplementationSource code, executables and the testing dataset are freely available at https://github.com/mapleforest/HaploMerger2/releases/.Supplementary information Supplementary data are available at Bioinformatics online.

Understanding potato with the help of genomics

Potato (Solanum tuberosum L.) is an important staple crop with a highly heterozygous and complex genome. Despite its cultural and economic significance, potato improvement efforts have been held back by the relative lack of genetic resources available to producers and breeders. The publication of the potato reference genome and advances in high-throughput sequencing technologies have led to the development of a wide range of genomic and transcriptomic resources. An overview these new tools, from the updated versions of the potato reference genome and transcriptome, to more recent gene expression, regulatory motif, re-sequencing and SNP genotyping analyses, paints a picture of modern potato research and how it will change our understanding of potato as well as other tuber producing Solanaceae.

A Genetic Algorithm for Diploid Genome Reconstruction Using Paired-End Sequencing.

The genome of many species in the biosphere is a diploid consisting of paternal and maternal haplotypes. The differences between these two haplotypes range from single nucleotide polymorphisms (SNPs) to large-scale structural variations (SVs). Existing genome assemblers for next-generation sequencing platforms attempt to reconstruct one consensus sequence, which is a mosaic of two parental haplotypes. Reconstructing paternal and maternal haplotypes is an important task in linkage analysis and association studies. This study designs and implemented HapSVAssembler on the basis of Genetic Algorithm (GA) and paired-end sequencing. The proposed method builds a consensus sequence, identifies various types of heterozygous variants, and reconstructs the paternal and maternal haplotypes by solving an optimization problem with a GA algorithm. Experimental results indicate that the HapSVAssembler has high accuracy and contiguity under various sequencing coverage, error rates, and insert sizes. The program is tested on pilot sequencing of a highly heterozygous genome, and 12,781 heterozygous SNPs and 602 hemizygous SVs are identified. We observe that, although the number of SVs is much less than that of SNPs, the genomic regions occupied by SVs are much larger, implying the heterozygosity computed using SNPs or k-mer spectrum may be under-estimated.

Doubled Haploid 'CUDH2107' as a Reference for Bulb Onion (Allium cepa L.) Research: Development of a Transcriptome Catalogue and Identification of Transcripts Associated with Male Fertility.

Researchers working on model plants have derived great benefit from developing genomic and genetic resources using ‘reference’ genotypes. Onion has a large and highly heterozygous genome making the sharing of germplasm and analysis of sequencing data complicated. To simplify the discovery and analysis of genes underlying important onion traits, we are promoting the use of the homozygous double haploid line ‘CUDH2107’ by the onion research community. In the present investigation, we performed transcriptome sequencing on vegetative and reproductive tissues of CUDH2107 to develop a multi-organ reference transcriptome catalogue. A total of 396 million 100 base pair paired reads was assembled using the Trinity pipeline, resulting in 271,665 transcript contigs. This dataset was analysed for gene ontology and transcripts were classified on the basis of putative biological processes, molecular function and cellular localization. Significant differences were observed in transcript expression profiles between different tissues. To demonstrate the utility of our CUDH2107 transcriptome catalogue for understanding the genetic and molecular basis of various traits, we identified orthologues of rice genes involved in male fertility and flower development. These genes provide an excellent starting point for studying the molecular regulation, and the engineering of reproductive traits.

Phased diploid genome assembly with single-molecule real-time sequencing.

While genome assembly projects have been successful in many haploid and inbred species, the assembly of noninbred or rearranged heterozygous genomes remains a major challenge. To address this challenge, we introduce the open-source FALCON and FALCON-Unzip algorithms (https://github.com/PacificBiosciences/FALCON/) to assemble long-read sequencing data into highly accurate, contiguous, and correctly phased diploid genomes. We generate new reference sequences for heterozygous samples including an F1 hybrid of Arabidopsis thaliana, the widely cultivated Vitis vinifera cv. Cabernet Sauvignon, and the coral fungus Clavicorona pyxidata, samples that have challenged short-read assembly approaches. The FALCON-based assemblies are substantially more contiguous and complete than alternate short- or long-read approaches. The phased diploid assembly enabled the study of haplotype structure and heterozygosities between homologous chromosomes, including the identification of widespread heterozygous structural variation within coding sequences.

Anthocyanin pigmentation study in flower petals of Asiatic hybrid lily

- Lilium is very important bulbous flower, cosmopolitan and has high demand in International flower market. Asiatic hybrid lily cultivars often have dark red spots on the interior surface of their petals. Pigments accumulated in petal spots are anthocyanins. In Lilium, the genetic basis of flower colour and spot formation are little understood because of the heterozygous genome structure. In this study 10 intravarietal Lilium hybrids along with their parents were tested and evaluated. The aim of this study was to estimate and analyse the flower colour of Asiatic Lilium hybrids. Highest value of anthocyanin was recorded in the hybrid PKLH-12 (43.48 mg/100g) and lowest in hybrid PKLH-6 (3.670 mg/100g). Anthocyanin content in the basal part was slightly higher (12.02 mg/100g) than the upper part (12.00mg/100g). Lilium petal during day of anthesis exhibits the highest amount of anthocyanin (13.96 mg/100g) than the 2 days after anthesis(10.05 mg/100g). Based on the variation in anthocyanin content in the present genotypes it may be concluded that this variability can be used for further development of indigenous hybrids with altered colours, hues, and patterns with long post harvest life. This study is important to understand the genetic basis of the traits and to apply molecular markers linked to map loci for MAS.

The draft genome of MD-2 pineapple using hybrid error correction of long reads.

The introduction of the elite pineapple variety, MD-2, has caused a significant market shift in the pineapple industry. Better productivity, overall increased in fruit quality and taste, resilience to chilled storage and resistance to internal browning are among the key advantages of the MD-2 as compared with its previous predecessor, the Smooth Cayenne. Here, we present the genome sequence of the MD-2 pineapple (Ananas comosus (L.) Merr.) by using the hybrid sequencing technology from two highly reputable platforms, i.e. the PacBio long sequencing reads and the accurate Illumina short reads. Our draft genome achieved 99.6% genome coverage with 27,017 predicted protein-coding genes while 45.21% of the genome was identified as repetitive elements. Furthermore, differential expression of ripening RNASeq library of pineapple fruits revealed ethylene-related transcripts, believed to be involved in regulating the process of non-climacteric pineapple fruit ripening. The MD-2 pineapple draft genome serves as an example of how a complex heterozygous genome is amenable to whole genome sequencing by using a hybrid technology that is both economical and accurate. The genome will make genomic applications more feasible as a medium to understand complex biological processes specific to pineapple.

An Improved Genome Assembly of Azadirachta indica A. Juss.

Neem (Azadirachta indica A. Juss.), an evergreen tree of the Meliaceae family, is known for its medicinal, cosmetic, pesticidal and insecticidal properties. We had previously sequenced and published the draft genome of a neem plant, using mainly short read sequencing data. In this report, we present an improved genome assembly generated using additional short reads from Illumina and long reads from Pacific Biosciences SMRT sequencer. We assembled short reads and error-corrected long reads using Platanus, an assembler designed to perform well for heterozygous genomes. The updated genome assembly (v2.0) yielded 3- and 3.5-fold increase in N50 and N75, respectively; 2.6-fold decrease in the total number of scaffolds; 1.25-fold increase in the number of valid transcriptome alignments; 13.4-fold less misassembly and 1.85-fold increase in the percentage repeat, over the earlier assembly (v1.0). The current assembly also maps better to the genes known to be involved in the terpenoid biosynthesis pathway. Together, the data represent an improved assembly of the A. indica genome.

Natural Allelic Variations in Highly Polyploidy Saccharum Complex.

Sugarcane (Saccharum spp.) is an important sugar and biofuel crop with high polyploid and complex genomes. The Saccharum complex, comprised of Saccharum genus and a few related genera, are important genetic resources for sugarcane breeding. A large amount of natural variation exists within the Saccharum complex. Though understanding their allelic variation has been challenging, it is critical to dissect allelic structure and to identify the alleles controlling important traits in sugarcane. To characterize natural variations in Saccharum complex, a target enrichment sequencing approach was used to assay 12 representative germplasm accessions. In total, 55,946 highly efficient probes were designed based on the sorghum genome and sugarcane unigene set targeting a total of 6 Mb of the sugarcane genome. A pipeline specifically tailored for polyploid sequence variants and genotype calling was established. BWA-mem and sorghum genome approved to be an acceptable aligner and reference for sugarcane target enrichment sequence analysis, respectively. Genetic variations including 1,166,066 non-redundant SNPs, 150,421 InDels, 919 gene copy number variations, and 1,257 gene presence/absence variations were detected. SNPs from three different callers (Samtools, Freebayes, and GATK) were compared and the validation rates were nearly 90%. Based on the SNP loci of each accession and their ploidy levels, 999,258 single dosage SNPs were identified and most loci were estimated as largely homozygotes. An average of 34,397 haplotype blocks for each accession was inferred. The highest divergence time among the Saccharum spp. was estimated as 1.2 million years ago (MYA). Saccharum spp. diverged from Erianthus and Sorghum approximately 5 and 6 MYA, respectively. The target enrichment sequencing approach provided an effective way to discover and catalog natural allelic variation in highly polyploid or heterozygous genomes.