Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits.

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Rapid and easy construction of a simplified amplicon sequencing (simplified AmpSeq) library for marker-assisted selection.

Marker-assisted selection (MAS) is fundamental for plant breeding programs, as it can identify desirable seedlings at a young stage and reduce the cost, time and space needed for plant maintenance, especially for perennial crops. To facilitate the process of genotyping, which is time consuming and laborious, we developed a simplified amplicon sequencing (simplified AmpSeq) library construction method for next-generation sequencing that can be applied to MAS in breeding programs. The method is based on one-step PCR with a mixture of two primer sets: the first consisting of tailed target primers, the second of primers that contain flow-cell binding sites, indexes and tail sequences complementary to those in the first set. To demonstrate the process of MAS using s implified AmpSeq, we created databases of genotypes for important traits by using cultivar collections including triploid cultivars and segregating seedlings of Japanese pear (Pyrus pyrifolia Nakai), Japanese chestnut (Castanea crenata Sieb. et Zucc.) and apple (Malus domestica Borkh.). Simplified AmpSeq has the advantages of high repeatability, ability to estimate allele number in polyploid species and semi-automatic evaluation using target allele frequencies. Because this method provides high flexibility for designing primer sets and targeting any variant, it will be useful for plant breeding programs.

Gene prediction for leaf margin phenotype and fruit flesh color in pineapple (Ananas comosus) using haplotype-resolved genome sequencing.

Pineapple (Ananas comosus (L.) Merr.) is one of the most economically important tropical fruit species. The major aim of the breeding programs in several countries, including Japan, is quality improvement, mainly for the fresh market. 'Yugafu', a Japanese cultivar with distinctive pipe-type leaf margin phenotype and white flesh color, is popular for fresh consumption. Therefore, genome sequencing of 'Yugafu' is expected to assist pineapple breeding. Here, we developed a haplotype-resolved assembly for the heterozygous genome of 'Yugafu' using long-read sequencing technology and obtained a pair of 25 pseudomolecule sequences inherited from the parental accessions 'Cream pineapple' and 'HI101'. The causative genes for leaf margin and fruit flesh color were identified. Fine mapping revealed a 162-kb region on CLG23 for the leaf margin phenotype. In this region, 20 kb of inverted repeat was specifically observed in the 'Cream pineapple' derived allele, and the WUSCHEL-related homeobox 3 (AcWOX3) gene was predicted as the key gene for leaf margin morphogenesis. Dominantly repressed AcWOX3 via RNAi was suggested to be the cause of the pipe-type leaf margin phenotype. Quantitative trait locus (QTL) analysis revealed that the terminal region of CLG08 contributed to white flesh and low carotenoid content. Carotenoid cleaved dioxygenase 4 (AcCCD4), a key gene for carotenoid degradation underlying this QTL, was predicted as the key gene for white flesh color through expression analysis. These findings could assist in modern pineapple breeding and facilitate marker-assisted selection for important traits.

Development of SSR Databases Available for Both NGS and Capillary Electrophoresis in Apple, Pear and Tea.

Developing new varieties in fruit and tea breeding programs is very costly and labor-intensive. Thus, establishing a variety discrimination system is important for protecting breeders' rights and producers' profits. Simple sequence repeat (SSR) databases that can be utilized for both next-generation sequencing (SSR-GBS) and polymerase chain reaction-capillary electrophoresis (PCR-CE) would be very useful in variety discrimination. In the present study, SSRs with tri-, tetra- and pentanucleotide repeats were examined in apple, pear and tea. Out of 37 SSRs that showed clear results in PCR-CE, 27 were suitable for SSR-GBS. Among the remaining markers, there was allele dropout for some markers that caused differences between the results of PCR-CE and SSR-GBS. For the selected 27 markers, the alleles detected by SSR-GBS were comparable to those detected by PCR-CE. Furthermore, we developed a computational pipeline for automated genotyping using SSR-GBS by setting a value "α" for each marker, a criterion whether a genotype is homozygous or heterozygous based on allele frequency. The set of 27 markers contains 10, 8 and 9 SSRs for apple, pear and tea, respectively, that are useful for both PCR-CE and SSR-GBS and suitable for automation. The databases help researchers discriminate varieties in various ways depending on sample size, markers and methods.

DNA marker for resistance to Puccinia horiana in chrysanthemum (Chrysanthemum morifolium Ramat.) "Southern Pegasus".

White rust caused by Puccinia horiana Henn. adversely affects chrysanthemum (Chrysanthemum morifolium Ramat.) production. The breeding of resistant varieties is effective in controlling the disease. Here we aimed to develop DNA markers for the strong resistance to P. horiana. We conducted a linkage analysis based on the genome-wide association study (GWAS) method. We employed a biparental population for the GWAS, wherein the single nucleotide polymorphism (SNP) allele frequency could be predicted. The population was derived from crosses between a strong resistant "Southern Pegasus" and a susceptible line. The GWAS used simplex and double-simplex SNP markers selected out of SNP candidates mined from ddRAD-Seq data of an F1 biparental population. These F1 individuals segregated in a 1:1 ratio of resistant to susceptible. Twenty-one simplex SNPs were significantly associated with P. horiana resistance in "Southern Pegasus" and generated one linkage group. These results show the presence of a single resistance gene in "Southern Pegasus". We identified the nearest SNP marker located 2.2 cM from P. horiana resistance locus and demonstrated this SNP marker-resistance link using an independent population. This is the first report of an effective DNA marker linked to a gene for P. horiana resistance in chrysanthemum.

Tracing founder haplotypes of Japanese apple varieties: application in genomic prediction and genome-wide association study.

Haplotypes provide useful information for genomics-based approaches, genomic prediction, and genome-wide association study. As a small number of superior founders have contributed largely to the breeding history of fruit trees, the information of founder haplotypes may be relevant for performing the genomics-based approaches in these plants. In this study, we proposed a method to estimate 14 haplotypes from 7 founders and automatically trace the haplotypes forward to apple parental (185 varieties) and breeding (659 F1 individuals from 16 full-sib families) populations based on 11,786 single-nucleotide polymorphisms, by combining multiple algorithms. Overall, 92% of the single-nucleotide polymorphisms information in the parental and breeding populations was characterized by the 14 founder haplotypes. The use of founder haplotype information improved the accuracy of genomic prediction in 7 traits and the resolution of genome-wide association study in 13 out of 27 fruit quality traits analyzed in this study. We also visualized the significant propagation of the founder haplotype with the largest genetic effect in genome-wide association study over the pedigree tree of the parental population. These results suggest that the information of founder haplotypes can be useful for not only genetic improvement of fruit quality traits in apples but also for understanding the selection history of founder haplotypes in the breeding program of Japanese apple varieties.

SSR markers developed using next-generation sequencing technology in pineapple, Ananas comosus (L.) Merr.

Simple sequence repeat (SSR) markers provide a reliable tool for the identification of accessions and the construction of genetic linkage maps because of their co-dominant inheritance. In the present study, we developed new SSR markers with next-generation sequencing using the Roche 454 GS FLX+ platform. Five hundred SSR primer sets were tested for the genetic identification of pineapple, including 100 each for the di-, tri-, tetra-, penta-, and hexa-nucleotide motif SSRs. In total, 160 SSR markers successfully amplified fragments and exhibited polymorphism among accessions. The SSR markers revealed the number of alleles per locus (ranging from 2 to 13), the expected heterozygosity (ranging from 0.041 to 0.823), and the observed heterozygosity (ranging from 0 to 0.875). A total of 117 SSR markers with tri- or greater nucleotide motifs were shown to be effective at facilitating accurate genotyping. Using the SSR markers, 25 accessions were distinguished genetically, with the exception of accessions 'MD-2' and 'Yonekura'. The developed SSR markers could facilitate the establishment of efficient and accurate genetic identification systems and the construction of genetic linkage maps in the future.

The use of a fertile doubled haploid apple line for QTL analysis of fruit traits.

Apple is an economically important crop, and various approaches to genetic analysis in breeding programs have been attempted, including the production of doubled haploid (DH) lines, which are genetically homozygous. In this study, we used a DH line for QTL analyses, for the first time in a fruit tree, expecting it to simplify the analysis of the inheritance of quantitative traits and thus to enhance QTL detection power. Using an F1 population from 'Prima' × 'Apple Chukanbohon 95P6' (DH), we constructed a genetic map of 'Prima', and identified 19 QTLs for 13 traits. These QTLs had comparatively high LOD scores and explained a large part of the variation of the phenotypes. In particular, acidity, juice browning, and skin splitting clearly segregated at a 1:1 ratio, consistent with the segregation of the alleles at the detected QTLs in linkage group 16; these traits appeared to be regulated by single genes, despite general consideration that they are quantitative traits. Using this simple genetic composition of the F1 population, we concluded that the skin splitting of apple fruit has recessive inheritance, and that the allele for splitting is tightly linked with those for high acidity and low juice browning in 'Prima'.

Genomic dissection of a 'Fuji' apple cultivar: re-sequencing, SNP marker development, definition of haplotypes, and QTL detection.

'Fuji' is one of the most popular and highly-produced apple cultivars worldwide, and has been frequently used in breeding programs. The development of genotypic markers for the preferable phenotypes of 'Fuji' is required. Here, we aimed to define the haplotypes of 'Fuji' and find associations between haplotypes and phenotypes of five traits (harvest day, fruit weight, acidity, degree of watercore, and flesh mealiness) by using 115 accessions related to 'Fuji'. Through the re-sequencing of 'Fuji' genome, total of 2,820,759 variants, including single nucleotide polymorphisms (SNPs) and insertions or deletions (indels) were detected between 'Fuji' and 'Golden Delicious' reference genome. We selected mapping-validated 1,014 SNPs, most of which were heterozygous in 'Fuji' and capable of distinguishing alleles inherited from the parents of 'Fuji' (i.e., 'Ralls Janet' and 'Delicious'). We used these SNPs to define the haplotypes of 'Fuji' and trace their inheritance in relatives, which were shown to have an average of 27% of 'Fuji' genome. Analysis of variance (ANOVA) based on 'Fuji' haplotypes identified one quantitative trait loci (QTL) each for harvest time, acidity, degree of watercore, and mealiness. A haplotype from 'Delicious' chr14 was considered to dominantly cause watercore, and one from 'Ralls Janet' chr1 was related to low-mealiness.

Fine mapping of the gene for susceptibility to black spot disease in Japanese pear (Pyrus pyrifolia Nakai).

Black spot disease, which is caused by the Japanese pear pathotype of the filamentous fungus Alternaria alternata (Fries) Keissler, is one of the most harmful diseases in Japanese pear cultivation. We mapped a gene for susceptibility to black spot disease in the Japanese pear (Pyrus pyrifolia Nakai) cultivar 'Kinchaku' (Aki gene) at the top of linkage group 11, similar to the positions of the susceptibility genes Ani in 'Osa Nijisseiki' and Ana in 'Nansui'. Using synteny-based marker enrichment, we developed novel apple SSR markers in the target region. We constructed a fine map of linkage group 11 of 'Kinchaku' and localized the Aki locus within a 1.5-cM genome region between SSR markers Mdo.chr11.28 and Mdo.chr11.34. Marker Mdo.chr11.30 co-segregated with Aki in all 621 F1 plantlets of a 'Housui' × 'Kinchaku' cross. The physical size of the Aki region, which includes three markers (Mdo.chr11.28, Mdo.chr11.30, and Mdo.chr11.34), was estimated to be 250 Kb in the 'Golden Delicious' apple genome and 107 Kb in the 'Dangshansuli' Chinese pear genome. Our results will help to identify the candidate gene for susceptibility to black spot disease in Japanese pear.

S-genotype identification based on allele-specific PCR in Japanese pear.

Gametophytic self-incompatibility in Japanese pear (Pyrus pyrifolia Nakai) is controlled by the single, multi-allelic S-locus. Information about the S-genotypes is important for breeding and the selection of pollen donors for fruit production. Rapid and reliable S-genotype identification system is necessary for efficient breeding of new cultivars in Japanese pear. We designed S allele-specific PCR primer pairs for ten previously reported S-RNase alleles (S (1)-S (9) and S (k)) as simple and reliable method. Specific nucleotide sequences were chosen to design the primers to amplify fragments of only the corresponding S alleles. The developed primer pairs were evaluated by using homozygous S-genotypes (S (1)/S (1)-S (9)/S (9) and S (4sm)/S (4sm)) and 14 major Japanese pear cultivars, and found that S allele-specific primer pairs can identify S-genotypes effectively. The S allele-specific primer pairs developed in this study will be useful for efficient S-genotyping and for marker-assisted selection in Japanese pear breeding programs.

Identification of QTLs for fruit quality traits in Japanese apples: QTLs for early ripening are tightly related to preharvest fruit drop.

Many important apple (Malus × domestica Borkh.) fruit quality traits are regulated by multiple genes, and more information about quantitative trait loci (QTLs) for these traits is required for marker-assisted selection. In this study, we constructed genetic linkage maps of the Japanese apple cultivars 'Orin' and 'Akane' using F1 seedlings derived from a cross between these cultivars. The 'Orin' map consisted of 251 loci covering 17 linkage groups (LGs; total length 1095.3 cM), and the 'Akane' map consisted of 291 loci covering 18 LGs (total length 1098.2 cM). We performed QTL analysis for 16 important traits, and found that four QTLs related to harvest time explained about 70% of genetic variation, and these will be useful for marker-assisted selection. The QTL for early harvest time in LG15 was located very close to the QTL for preharvest fruit drop. The QTL for skin color depth was located around the position of MYB1 in LG9, which suggested that alleles harbored by 'Akane' are regulating red color depth with different degrees of effect. We also analyzed soluble solids and sugar component contents, and found that a QTL for soluble solids content in LG16 could be explained by the amount of sorbitol and fructose.

Identification of QTLs controlling harvest time and fruit skin color in Japanese pear (Pyrus pyrifolia Nakai).

Using an F1 population from a cross between Japanese pear (Pyrus pyrifolia Nakai) cultivars 'Akiakari' and 'Taihaku', we performed quantitative trait locus (QTL) analysis of seven fruit traits (harvest time, fruit skin color, flesh firmness, fruit weight, acid content, total soluble solids content, and preharvest fruit drop). The constructed simple sequence repeat-based genetic linkage map of 'Akiakari' consisted of 208 loci and spanned 799 cM; that of 'Taihaku' consisted of 275 loci and spanned 1039 cM. Out of significant QTLs, two QTLs for harvest time, one for fruit skin color, and one for flesh firmness were stably detected in two successive years. The QTLs for harvest time were located at the bottom of linkage group (LG) Tai3 (nearest marker: BGA35) and at the top of LG Tai15 (nearest markers: PPACS2 and MEST050), in good accordance with results of genome-wide association study. The PPACS2 gene, a member of the ACC synthase gene family, may control harvest time, preharvest fruit drop, and fruit storage potential. One major QTL associated with fruit skin color was identified at the top of LG 8. QTLs identified in this study would be useful for marker-assisted selection in Japanese pear breeding programs.

Genotyping of strawberry (Fragaria x ananassa Duch.) cultivars by DNA markers: interlaboratory study.

Fourteen Japanese laboratories validated the reproducibility of genotyping by 25 cleavage amplified polymorphic sequence (CAPS) markers for discrimination of strawberry (Fragaria x ananassa Duch.) cultivars. Both the sensitivity and specificity rate of 12 markers were 100%, those of another 12 were >95%, and those of 1 were >90%. These results indicate that the method of genotyping by the CAPS markers was highly reproducible and could provide a useful basis for practical identification of strawberry cultivars. This is the first report of the statistical validation of crop genotyping by DNA markers.

Identification of QTLs for resistance to powdery mildew and SSR markers diagnostic for powdery mildew resistance genes in melon (Cucumis melo L.).

Powdery mildew caused by Podosphaera xanthii is an important foliar disease in melon. To find molecular markers for marker-assisted selection, we constructed a genetic linkage map of melon based on a population of 93 recombinant inbred lines derived from crosses between highly resistant AR 5 and susceptible 'Earl's Favourite (Harukei 3)'. The map spans 877 cM and consists of 167 markers, comprising 157 simple sequence repeats (SSRs), 7 sequence characterized amplified region/cleavage amplified polymorphic sequence markers and 3 phenotypic markers segregating into 20 linkage groups. Among them, 37 SSRs and 6 other markers were common to previous maps. Quantitative trait locus (QTL) analysis identified two loci for resistance to powdery mildew. The effects of these QTLs varied depending on strain and plant stage. The percentage of phenotypic variance explained for resistance to the pxA strain was similar between QTLs (R (2) = 22-28%). For resistance to pxB strain, the QTL on linkage group (LG) XII was responsible for much more of the variance (41-46%) than that on LG IIA (12-13%). The QTL on LG IIA was located between two SSR markers. Using an independent population, we demonstrated the effectiveness of these markers. This is the first report of universal and effective markers linked to a gene for powdery mildew resistance in melon.

Geranyl diphosphate:4-hydroxybenzoate geranyltransferase from Lithospermum erythrorhizon. Cloning and characterization of a ket enzyme in shikonin biosynthesis.

Two cDNAs encoding geranyl diphosphate:4-hy- droxybenzoate 3-geranyltransferase were isolated from Lithospermum erythrorhizon by nested PCR using the conserved amino acid sequences among polyprenyl- transferases for ubiquinone biosynthesis. They were functionally expressed in yeast COQ2 disruptant and showed a strict substrate specificity for geranyl diphosphate as the prenyl donor, in contrast to ubiquinone biosynthetic enzymes, suggesting that they are involved in the biosynthesis of shikonin, a naphthoquinone secondary metabolite. Regulation of their expression by various culture conditions coincided with that of geranyltransferase activity and the secondary metabolites biosynthesized via this enzyme. This is the first established plant prenyltransferase that transfers the prenyl chain to an aromatic substrate.