Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids.

Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost-efficient genotyping using next-generation sequencing. However, accurate genotype calling from next-generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome-wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype-calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum, and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype-calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost-efficiently implement GWAS and genomic prediction in polyploid crops.

Managing flowering time in Miscanthus and sugarcane to facilitate intra- and intergeneric crosses.

Miscanthus is a close relative of Saccharum and a potentially valuable genetic resource for improving sugarcane. Differences in flowering time within and between Miscanthus and Saccharum hinders intra- and interspecific hybridizations. A series of greenhouse experiments were conducted over three years to determine how to synchronize flowering time of Saccharum and Miscanthus genotypes. We found that day length was an important factor influencing when Miscanthus and Saccharum flowered. Sugarcane could be induced to flower in a central Illinois greenhouse using supplemental lighting to reduce the rate at which days shortened during the autumn and winter to 1 min d-1, which allowed us to synchronize the flowering of some sugarcane genotypes with Miscanthus genotypes primarily from low latitudes. In a complementary growth chamber experiment, we evaluated 33 Miscanthus genotypes, including 28 M. sinensis, 2 M. floridulus, and 3 M. ×giganteus collected from 20.9° S to 44.9° N for response to three day lengths (10 h, 12.5 h, and 15 h). High latitude-adapted M. sinensis flowered mainly under 15 h days, but unexpectedly, short days resulted in short, stocky plants that did not flower; in some cases, flag leaves developed under short days but heading did not occur. In contrast, for M. sinensis and M. floridulus from low latitudes, shorter day lengths typically resulted in earlier flowering, and for some low latitude genotypes, 15 h days resulted in no flowering. However, the highest ratio of reproductive shoots to total number of culms was typically observed for 12.5 h or 15 h days. Latitude of origin was significantly associated with culm length, and the shorter the days, the stronger the relationship. Nearly all entries achieved maximal culm length under the 15 h treatment, but the nearer to the equator an accession originated, the less of a difference in culm length between the short-day treatments and the 15 h day treatment. Under short days, short culms for high-latitude accessions was achieved by different physiological mechanisms for M. sinensis genetic groups from the mainland in comparison to those from Japan; for mainland accessions, the mechanism was reduced internode length, whereas for Japanese accessions the phyllochron under short days was greater than under long days. Thus, for M. sinensis, short days typically hastened floral induction, consistent with the expectations for a facultative short-day plant. However, for high latitude accessions of M. sinensis, days less than 12.5 h also signaled that plants should prepare for winter by producing many short culms with limited elongation and development; moreover, this response was also epistatic to flowering. Thus, to flower M. sinensis that originates from high latitudes synchronously with sugarcane, the former needs day lengths >12.5 h (perhaps as high as 15 h), whereas that the latter needs day lengths <12.5 h.

Training Population Optimization for Genomic Selection in Miscanthus.

Miscanthus is a perennial grass with potential for lignocellulosic ethanol production. To ensure its utility for this purpose, breeding efforts should focus on increasing genetic diversity of the nothospecies Miscanthus × giganteus (M×g) beyond the single clone used in many programs. Germplasm from the corresponding parental species M. sinensis (Msi) and M. sacchariflorus (Msa) could theoretically be used as training sets for genomic prediction of M×g clones with optimal genomic estimated breeding values for biofuel traits. To this end, we first showed that subpopulation structure makes a substantial contribution to the genomic selection (GS) prediction accuracies within a 538-member diversity panel of predominately Msi individuals and a 598-member diversity panels of Msa individuals. We then assessed the ability of these two diversity panels to train GS models that predict breeding values in an interspecific diploid 216-member M×g F2 panel. Low and negative prediction accuracies were observed when various subsets of the two diversity panels were used to train these GS models. To overcome the drawback of having only one interspecific M×g F2 panel available, we also evaluated prediction accuracies for traits simulated in 50 simulated interspecific M×g F2 panels derived from different sets of Msi and diploid Msa parents. The results revealed that genetic architectures with common causal mutations across Msi and Msa yielded the highest prediction accuracies. Ultimately, these results suggest that the ideal training set should contain the same causal mutations segregating within interspecific M×g populations, and thus efforts should be undertaken to ensure that individuals in the training and validation sets are as closely related as possible.

Contrasting allelic distribution of CO/Hd1 homologues in Miscanthus sinensis from the East Asian mainland and the Japanese archipelago.

The genus Miscanthus is a perennial C4 grass native to eastern Asia and is a promising candidate bioenergy crop for cool temperate areas. Flowering time is a crucial factor governing regional and seasonal adaptation; in addition, it is also a key target trait for extending the vegetative phase to improve biomass potential. Homologues of CONSTANS (CO)/Heading date 1(Hd1) were cloned from Miscanthus sinensis and named MsiHd1. Sequences of MsiHd1 homologues were compared among 24 wild M. sinensis accessions from Japan, 14 from China, and three from South Korea. Two to five MsiHd1 alleles in each accession were identified, suggesting that MsiHd1 consists of at least three loci in the Miscanthus genome. Verifying the open reading frame in MsiHd1, they were classified as putative functional alleles without mutations or non-functional alleles caused by indels. The Neighbor-Joining tree indicated that one of the multiple MsiHd1 loci is a pseudogene locus without any functional alleles. The pseudogene locus was named MsiHd1b, and the other loci were considered to be part of the MsiHd1a multi-locus family. Interestingly, in most Japanese accessions 50% or more of the MsiHd1a alleles were non-functional, whereas accessions from the East Asian mainland harboured only functional alleles. Five novel miniature inverted transposable elements (MITEs) (MsiMITE1-MsiMITE5) were observed in MsiHd1a/b. MsiMITE1, detected in exon 1 of MsiHd1a, was only observed in Japanese accessions and its revertant alleles derived from retransposition were predominantly in Chinese accessions. These differences in MsiHd1a show that the dependency on functional MsiHd1a alleles is different between accessions from the East Asian mainland and Japan.

Comparison of the enzymatic digestibility of physically and chemically pretreated selected line of diploid-Miscanthus sinensis Shiozuka and triploid-M.×giganteus.

The diploid Miscanthus sinensis "Shiozuka" which was selected as a high-biomass producing line, and the triploid M. × giganteus (M×G) were treated by ball milling (physical treatment) and alkaline hydrogen peroxide treatment (AHP; chemical treatment), and their structural sugar compositions and enzymatic digestibility were compared. The structural sugar content of Shiozuka was moderate and lower than that of M×G. The Klason lignin content of Shiozuka was also lower than that of M×G. However, Shiozuka was sensitive to ball milling and AHP treatment; ball milled and AHP-treated Shiozuka had higher enzymatic digestibility than ball milled and AHP-treated M×G. Shiozuka would be promising feedstock to obtain fermentable sugars with low energy consumption. Finally, enzymes for the hydrolysis of chemically treated Miscanthus were isolated from Trichoderma reesei ATCC 66589 and Penicillium pinophilum. The sugar yield could be increased by enzymatic hydrolysis of AHP-treated samples with NaOH and H2O2 and the isolated enzymes.