Search

Quinoa is emerging as a key seed crop for global food security due to its ability to grow in marginal environments and its excellent nutritional properties. Because quinoa is partially allogamous, we have developed quinoa inbred lines necessary for molecular genetic analysis. Our comprehensive genomic analysis showed that the quinoa inbred lines fall into three genetic subpopulations: northern highland, southern highland, and lowland. Lowland and highland quinoa are the same species, but have very different genotypes and phenotypes. Lowland quinoa has relatively small grains and a darker grain color, and is widely tested and grown around the world. In contrast, the white, large-grained highland quinoa is grown in the Andean highlands, including the region where quinoa originated, and is exported worldwide as high-quality quinoa. Recently, we have shown that viral vectors can be used to regulate endogenous genes in quinoa, paving the way for functional genomics to reveal the diversity of quinoa. However, although a high-quality assembly has recently been reported for a lowland quinoa line, genomic resources of the quality required for functional genomics are not available for highland quinoa lines. Here we present high-quality chromosome-level genome assemblies for two highland inbred quinoa lines, J075 representing the northern highland line and J100 representing the southern highland line, using PacBio HiFi sequencing and dpMIG-seq. In addition, we demonstrate the importance of verifying and correcting reference-based scaffold assembly with other approaches such as linkage maps. The assembled genome sizes of J075 and J100 are 1.29 and 1.32 Gb, with contigs N50 of 66.3 and 12.6 Mb, and scaffold N50 of 71.2 and 70.6 Mb, respectively, comprising 18 pseudochromosomes. The repetitive sequences of J075 and J100 represent 72.6% and 71.5% of the genome, the majority of which are long terminal repeats, representing 44.0% and 42.7% of the genome, respectively. The de novo assembled genomes of J075 and J100 were predicted to contain 65,303 and 64,945 protein-coding genes, respectively. The high quality genomes of these highland quinoa lines will facilitate quinoa functional genomics research on quinoa and contribute to the identification of key genes involved in environmental adaptation and quinoa domestication.

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), ‘415-1’, by using γ-irradiated pollen. ‘415-1’ carries the S-genotype S4dS5S5, with “d” indicating a duplication of S5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between ‘515-20’ (S1S3) and ‘415-1’. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

Given the importance of prioritizing genome-based breeding of sweet potato to enable the promotion of food and nutritional security for future human societies, here, we aimed to dissect the genetic basis of storage root starch content (SC) when associated with a complex set of breeding traits including dry matter (DM) rate, storage root fresh weight (SRFW), and anthocyanin (AN) content in a mapping population containing purple-fleshed sweet potato. A polyploid genome-wide association study (GWAS) was extensively exploited using 90,222 single-nucleotide polymorphisms (SNPs) obtained from a bi-parental 204 F1 population between ‘Konaishin’ (having high SC but no AN) and ‘Akemurasaki’ (having high AN content but moderate SC). Through the comparison of polyploid GWAS on the whole set of the 204 F1, 93 high-AN-containing F1, and 111 low-AN-containing F1 populations, a total of two (consists of six SNPs), two (14 SNPs), four (eight SNPs), and nine (214 SNPs) significantly associated signals were identified for the variations of SC, DM, SRFW, and the relative AN content, respectively. Of them, a novel signal associated with SC, which was most consistent in 2019 and 2020 in both the 204 F1 and 111 low-AN-containing F1 populations, was identified in homologous group 15. The five SNP markers associated with homologous group 15 could affect SC improvement with a degree of positive effect (~4.33) and screen high-starch-containing lines with higher efficiency (~68%). In a database search of 62 genes involved in starch metabolism, five genes including enzyme genes granule-bound starch synthase I (IbGBSSI), α-amylase 1D, α-amylase 1E, and α-amylase 3, and one transporter gene ATP/ADP-transporter were located on homologous group 15. In an extensive qRT-PCR of these genes using the storage roots harvested at 2, 3, and 4 months after field transplantation in 2022, IbGBSSI, which encodes the starch synthase isozyme that catalyzes the biosynthesis of amylose molecule, was most consistently elevated during starch accumulation in sweet potato. These results would enhance our understanding of the underlying genetic basis of a complex set of breeding traits in the starchy roots of sweet potato, and the molecular information, particularly for SC, would be a potential platform for molecular marker development for this trait.

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Abstract Eustoma grandiflorumE. grandiflorum2E. grandiflorumC. canephoraO. pumilaE. grandiflorumEustoma

Nakano et al. sequenced the genome of a model species for hexaploid cultivated chrysanthemum, an important cut flower. While they investigated evolution, duplication events, and a novel retrotransposon in this study, they also demonstrated the translational implication of the genome sequence in breeding the economically important plant.

Abstract Background Understanding mechanisms of sugar accumulation and composition is essential to determining fruit quality and maintaining a desirable balance of sugars in plant storage organs. The major sugars in mature Rosaceae fruits are sucrose, fructose, glucose, and sorbitol. Among these, sucrose and fructose have high sweetness, whereas glucose and sorbitol have low sweetness. Japanese pear has extensive variation in individual sugar contents in mature fruit. Increasing total sugar content and that of individual high-sweetness sugars is a major target of breeding programs. The objective of this study was to identify quantitative trait loci (QTLs) associated with fruit traits including individual sugar accumulation, to infer the candidate genes underlying the QTLs, and to assess the potential of genomic selection for breeding pear fruit traits. Results We evaluated 10 fruit traits and conducted genome-wide association studies (GWAS) for 106 cultivars and 17 breeding populations (1112 F1 individuals) using 3484 tag single-nucleotide polymorphisms (SNPs). By implementing a mixed linear model and a Bayesian multiple-QTL model in GWAS, 56 SNPs associated with fruit traits were identified. In particular, a SNP located close to acid invertase gene PPAIV3 on chromosome 7 and a newly identified SNP on chromosome 11 had quite large effects on accumulation of sucrose and glucose, respectively. We used ‘Golden Delicious’ doubled haploid 13 (GDDH13), an apple reference genome, to infer the candidate genes for the identified SNPs. In the region flanking the SNP on chromosome 11, there is a tandem repeat of early responsive to dehydration (ERD6)-like sugar transporter genes that might play a role in the phenotypes observed. Conclusions SNPs associated with individual sugar accumulation were newly identified at several loci, and candidate genes underlying QTLs were inferred using advanced apple genome information. The candidate genes for the QTLs are conserved across Pyrinae genomes, which will be useful for further fruit quality studies in Rosaceae. The accuracies of genomic selection for sucrose, fructose, and glucose with genomic best linear unbiased prediction (GBLUP) were relatively high (0.67–0.75), suggesting that it would be possible to select individuals having high-sweetness fruit with high sucrose and fructose contents and low glucose content.

Identifying genes/loci for resistance to pine wilt disease (PWD) caused by the pine wood nematode (PWN) is beneficial for improving resistance breeding in Pinus thunbergii, but to date, genetic information using molecular markers has been limited. Here, we constructed a high-density linkage map using genotyping-by-sequencing (GBS) and conducted quantitative trait loci (QTL) analysis for PWD resistance for the self-pollinated progeny of “Namikata 73,” which is the most resistant variety among resistant varieties of P. thunbergii, following inoculation tests with PWN. An S1 mapping population consisting of the 116 progenies derived from self-pollination of the resistant variety, “Namikata 73” (resistance rank 5 to PWN), was inoculated with PWN isolate Ka-4 and evaluated for disease symptoms. To construct a high-density linkage map, we used single-nucleotide polymorphisms (SNPs) identified by GBS based on next-generation sequencing technology and some anchor DNA markers, expressed sequence tag (EST)-derived SNP markers and EST-derived simple sequence repeat (SSR) markers, and genomic SSR markers. The linkage map had 13 linkage groups (LGs) consisting of 2,365 markers including 2,243 GBS-SNP markers over a total map distance of 1968.4 centimorgans (cM). Results from QTL analysis using phenotype data and the linkage map indicated that PWD resistance is controlled by a single locus located on LG-3, as identified in a previous study. This locus showed overdominant genetic action in the present study. With the confirmation of PWD1 in two different mapping populations (present study and a previous study), the locus associated with this region is thought to be a good target for marker-assisted selection in P. thunbergii breeding programs in order to obtain high levels of resistance to PWD caused by PWN.

The southern root-knot nematode (SRKN; Meloidogyne incognita) is a typical parasitic nematode that affects sweetpotato [Ipomoea batatas (L.) Lam.], causing a significant decrease in crop yield and commercial value. In Japan, the SRKN is classified into 10 races: SP1–SP5, SP6-1, SP6-2, and SP7–SP9, with the dominant race differing according to the cultivation area. Soil insecticides have previously been used to reduce the soil density of SRKNs; however, this practice is both costly and labor intensive. Therefore, the development of SRKN-resistant sweetpotato lines and cultivars is necessary. However, due to the complexity of polyploid inheritance and the highly heterogeneous genomic composition of sweetpotato, genetic information and research for this species are significantly lacking compared to those for other major diploid crop species. In this study, we utilized the recently developed genome-wide association approach, which uses multiple-dose markers to assess autopolyploid species. We performed an association analysis to investigate resistance toward SRKN-SP2, which is the major race in areas with high sweetpotato production in Japan. The segregation ratio of resistant and susceptible lines in the F1 mapping population derived from the resistant “J-Red” and susceptible “Choshu” cultivars was fitted to 1: 3, suggesting that resistance to SP2 may be regulated by two loci present in the simplex. By aligning the double digest restriction-site associated DNA sequencing reads to the published Ipomoea trifida reference sequence, 46,982 single nucleotide polymorphisms (SNPs) were identified (sequencing depth > 200). The association study yielded its highest peak on chromosome 7 (Chr07) and second highest peak on chromosome 3 (Chr03), presenting as a single-dose in both loci. Selective DNA markers were developed to screen for resistant plants using the SNPs identified on Chr03 and Chr07. Our results showed that SRKN-SP2-resistant plants were selected with a probability of approximately 70% when combining the two selective DNA markers. This study serves as a model for the identification of genomic regions that control agricultural traits and the elucidation of their effects, and is expected to greatly advance marker-assisted breeding and association studies in polyploid crop species.

Next-generation sequencing (NGS)-based genotyping methods can generate numerous genetic markers in a single experiment and have contributed to plant genetic mapping. However, for high precision genetic analysis, the complicated genetic segregation mode in polyploid organisms requires high-coverage NGS data and elaborate analytical algorithms. In the present study, we propose a simple strategy for the genetic mapping of polyploids using low-coverage NGS data. The validity of the strategy was investigated using simulated data. Previous studies indicated that accurate allele dosage estimation from low-coverage NGS data (read depth < 40) is difficult. Therefore, we used allele dosage probabilities calculated from read counts in association analyses to detect loci associated with phenotypic variations. The allele dosage probabilities showed significant detection power, although higher allele dosage estimation accuracy resulted in higher detection power. On the contrary, differences in the segregation patterns between the marker and causal genes resulted in a drastic decrease in detection power even if the marker and casual genes were in complete linkage and the allele dosage estimation was accurate. These results indicated that the use of a larger number of markers is advantageous, even if the accuracy of allele dosage estimation is low. Finally, we applied the strategy for the genetic mapping of autohexaploid sweet potato (Ipomoea batatas) populations to detect loci associated with agronomic traits. Our strategy could constitute a cost-effective approach for preliminary experiments done performed to large-scale studies.

Abstract Objective Low DNA sequence polymorphism despite enormous phenotypic variations in peanut indicates the possible role of epigenetic variations. An attempt was made to analyze genome-wide DNA methylation pattern and its influence on gene expression across 11 diverse genotypes of peanut. Results The genotypes were subjected to bisulfite sequencing after 21 days of sowing (DAS). CHG regions showed the highest (30,537,376) DNA methylation followed by CpG (30,356,066) and CHH (15,993,361) across 11 genotypes. The B sub-genome exhibited higher DNA methylation sites (46,294,063) than the A sub-genome (30,415,166). Overall, the DNA methylation was more frequent in inter-genic regions than in the genic regions. The genes showing altered methylation and expression between the parent (TMV 2) and its EMS-derived mutant (TMV 2-NLM) were identified. Foliar disease resistant genotypes showed significant differential DNA methylation at 766 sites corresponding to 25 genes. Of them, two genes (Arahy.1XYC2X on chromosome 01 and Arahy.00Z2SH on chromosome 17) coding for senescence-associated protein showed differential expression with resistant genotypes recording higher fragments per kilobase of transcript per million mapped reads (FPKM) at their epialleles. Overall, the study indicated the variation in the DNA methylation pattern among the diverse genotypes of peanut and its influence of gene expression.

Genetic variation in phenological traits is the key in expanding production areas of crops. Southern highbush blueberry (SHB) is a blueberry cultivar group adapted to warmer climates and has been developed by multiple interspecific hybridizations between elite northern highbush blueberry (NHB) (Vaccinium corymbosum L.) and low-chill Vaccinium species native to the southern United States. In this study, we employed a collection of diverse SHB accessions and performed a genome-wide association study (GWAS) for five phenology-related traits [chilling requirement (CR), flowering date, ripening date, fruit development period, and continuous flowering] using polyploid GWAS models. Phenology-related traits showed higher heritability and larger correlation coefficients between year replications, which resulted in the detection of robust phenotype–genotype association peaks. Notably, a single association peak for the CR was detected on Chromosome 4. Comparison of genotypes at the GWAS peaks between NHB and SHB revealed the putative introgression of low-chill and late-flowering alleles into the highbush genetic pool. Our results provide basic insights into the diversity of phenological traits in blueberry and the genetic establishment of current highbush cultivar groups.

To gain insights into the genetic mechanisms underlying blooming and petal movement in flowering cherry (Cerasus × yedoensis), we performed time-course RNA-seq analysis of the floral buds and open-flowers of the most popular flowering cherry cultivar, ‘Somei-Yoshino.’ Independent biological duplicate samples of floral buds and open-flowers were collected from ‘Somei-Yoshino’ trees grown at three different locations in Japan. RNA-seq reads obtained from floral bud and open-flower samples collected in the current study (in 2019) and in a previous study (in 2017) were aligned against the genome sequence of ‘Somei-Yoshino’ to quantify gene transcript levels. Clustering analysis of RNA-seq reads revealed dynamic changes in the transcriptome, with genes in seven modules predominantly expressed at specific time points, ranging from 5 weeks before flowering to 2 weeks after flowering. Based on the identified gene modules and Gene Ontology (GO) terms enriched at different floral stages, we speculate that the genetic mechanisms underlying petal movement and flower opening in cherry involve the processes of development, cell wall organization, reproduction, and metabolism, which are executed by genes encoding transcription factors, phytohormones, transporters, and polysaccharide metabolic enzymes. Furthermore, we established a statistical model for cherry bloom forecasting, based on gene expression levels as RNA markers at different time points before flowering.

Abstract Wild plants are often tolerant to biotic and abiotic stresses in their natural environments, whereas domesticated plants such as crops frequently lack such resilience. This difference is thought to be due to the high levels of genome heterozygosity in wild plant populations and the low levels of heterozygosity in domesticated crop species. In this study, common vetch (Vicia sativa) was used as a model to examine this hypothesis. The common vetch genome (2n = 14) was estimated as 1.8 Gb in size. Genome sequencing produced a reference assembly that spanned 1.5 Gb, from which 31,146 genes were predicted. Using this sequence as a reference, 24,118 single nucleotide polymorphisms were discovered in 1243 plants from 12 natural common vetch populations in Japan. Common vetch genomes exhibited high heterozygosity at the population level, with lower levels of heterozygosity observed at specific genome regions. Such patterns of heterozygosity are thought to be essential for adaptation to different environments. The resources generated in this study will provide insights into de novo domestication of wild plants and agricultural enhancement.

Horsegram (Macrotyloma uniflorum [Lam.] Verdc.) is an underutilized warm-season diploid legume (2n = 20, 22). Because of its ability to grow under water-deficient and marginal soil conditions, horsegram is a preferred choice in the era of global climate change. In recognition of its potential as a crop species, we generated and analyzed a draft genome sequence for a horsegram variety, HPK-4. Ten chromosome-scale pseudomolecules were created by aligning Illumina scaffold sequences onto a linkage map. The total length of the ten pseudomolecules was 259.2 Mbp, covering 89% of the total length of the assembled sequences. A total of 36,105 genes were predicted on the assembled sequences. Diversity analysis of 89 horsegram accessions by dd-RAD-Seq identified 277 single nucleotide polymorphisms (SNPs), suggesting narrow genetic diversity among the horsegram accessions. This is the first attempt to generate a draft genome sequence of horsegram and will provide a reference for sequence-based analysis of horsegram germplasm.

With an objective of identifying the genomic regions for productivity and quality traits in peanut, a recombinant inbred line (RIL) population developed from an elite variety, TMV 2 and its ethyl methane sulfonate (EMS)-derived mutant was phenotyped over six seasons and genotyped with genotyping-by-sequencing (GBS), Arachis hypogaea transposable element (AhTE) and simple sequence repeats (SSR) markers. The genetic map with 700 markers spanning 2,438.1 cM was employed for quantitative trait loci (QTL) analysis which identified a total of 47 main-effect QTLs for the productivity and oil quality traits with the phenotypic variance explained (PVE) of 10–52% over the seasons. A common QTL region (46.7–50.1 cM) on Ah02 was identified for the multiple traits, such as a number of pods per plant (NPPP), pod weight per plant (PWPP), shelling percentage (SP), and test weight (TW). Similarly, a QTL (7.1–18.0 cM) on Ah16 was identified for both SP and protein content (PC). Epistatic QTL (epiQTL) analysis revealed intra- and inter-chromosomal interactions for the main-effect QTLs and other genomic regions governing these productivity traits. The markers identified by a single marker analysis (SMA) mapped to the QTL regions for most of the traits. Among the five potential candidate genes identified for PC, SP and oil quality, two genes (Arahy.7A57YA and Arahy.CH9B83) were affected by AhMITE1 transposition, and three genes (Arahy.J5SZ1I, Arahy.MZJT69, and Arahy.X7PJ8H) involved functional single nucleotide polymorphisms (SNPs). With major and consistent effects, the genomic regions, candidate genes, and the associated markers identified in this study would provide an opportunity for gene cloning and genomics-assisted breeding for increasing the productivity and enhancing the quality of peanut.

The common fig (Ficus carica L.) has a gynodioecious breeding system, and its sex phenotype is an important trait for breeding because only female plant fruits are edible. During breeding to select for female plants, we analyzed the FcRAN1 genotype, which is strongly associated with the sex phenotype. In 12 F1 populations derived from 13 cross combinations, the FcRAN1 genotype segregation ratio was 1:1, whereas the M119-226 × H238-107 hybridization resulted in an extremely male-biased segregation ratio (178:7 = male:female). This finding suggests that the segregation distortion was caused by some genetic factor(s). A whole-genome resequencing of breeding parents (paternal and maternal lines) identified 9,061 high-impact SNPs in the parents. A genome-wide linkage analysis exploring the gene(s) responsible for the distortion revealed 194 high-impact SNPs specific to Caprifig6085 (i.e., seed parent ancestor) and 215 high-impact SNPs specific to H238-107 (i.e., pollen parent) in 201 annotated genes. A comparison between the annotated genes and the genes required for normal embryo or gametophyte development and function identified several candidate genes possibly responsible for the segregation distortion. This is the first report describing segregation distortion in F. carica.

Abstract Background Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD. Results A total F1 mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties ‘Tanabe 54’ (resistant rank 2 to PWD) and ‘Tosashimizu 63’ (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of ‘Tanabe 54’ consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of ‘Tosashimizu 63’ consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the ‘Tanabe 54’ female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1. Conclusions The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.

Abstract Background Plasmodiophora brassicae is a soil-borne plant pathogen that causes clubroot disease, which results in crop yield loss in cultivated Brassica species. Here, we investigated whether a quantitative trait locus (QTL) in B. rapa might confer resistance to a Korean P. brassicae pathotype isolate, Seosan. We crossed resistant and susceptible parental lines and analyzed the segregation pattern in a F2 population of 348 lines. We identified and mapped a novel clubroot resistance QTL using the same mapping population that included susceptible Chinese cabbage and resistant turnip lines. Forty-five resistant and 45 susceptible F2 lines along with their parental lines were used for double digest restriction site-associated DNA sequencing (ddRAD-seq). High resolution melting (HRM)-based validation of SNP positions was conducted to confirm the novel locus. Results A 3:1 ratio was observed for resistant: susceptible genotypes, which is in accordance with Mendelian segregation. ddRAD-seq identified a new locus, CRs, on chromosome A08 that was different from the clubroot resistance (CR) locus, Crr1. HRM analysis validated SNP positions and constricted CRs region. Four out of seventeen single nucleotide polymorphisms (SNPs) positions were within a 0.8-Mb region that included three NBS-LRR candidate genes but not Crr1. Conclusion The newly identified CRs locus is a novel clubroot resistance locus, as the cultivar Akimeki bears the previously known Crr1 locus but remains susceptible to the Seosan isolate. These results could be exploited to develop molecular markers to detect Seosan-resistant genotypes and develop resistant Chinese cabbage cultivars.

Cultivated strawberry is the most widely consumed fruit crop in the world, and therefore, many breeding programs are underway to improve its agronomic traits such as fruit quality. Strawberry cultivars were vegetatively propagated through runners and carried a high risk of infection with viruses and insects. To solve this problem, the development of F1 hybrid seeds has been proposed as an alternative breeding strategy in strawberry. In this study, we conducted a potential assessment of genomic selection (GS) in strawberry F1 hybrid breeding. A total of 105 inbred lines were developed as candidate parents of strawberry F1 hybrids. In addition, 275 parental combinations were randomly selected from the 105 inbred lines and crossed to develop test F1 hybrids for GS model training. These populations were phenotyped for petiole length, leaf area, Brix, fruit hardness, and pericarp color. Whole-genome shotgun sequencing of the 105 inbred lines detected 20,811 single nucleotide polymorphism sites that were provided for subsequent GS analyses. In a GS model construction, inclusion of dominant effects showed a slight advantage in GS accuracy. In the across population prediction analysis, GS models using the inbred lines showed predictability for the test F1 hybrids and vice versa, except for Brix. Finally, the GS models were used for phenotype prediction of 5,460 possible F1 hybrids from 105 inbred lines to select F1 hybrids with high fruit hardness or high pericarp color. These F1 hybrids were developed and phenotyped to evaluate the efficacy of the GS. As expected, F1 hybrids that were predicted to have high fruit hardness or high pericarp color expressed higher observed phenotypic values than the F1 hybrids that were selected for other objectives. Through the analyses in this study, we demonstrated that GS can be applied for strawberry F1 hybrid breeding.

AbstractPea (Pisum sativum

Abstract Tomato is an attractive fruiting vegetable crop that can be used as an ornamental plant. Agronomical traits have been subjected to extensive genetic dissection to enhance vegetable breeding programs. By contrast, there are few genetic studies of ornamental traits for the development of ornamental tomato varieties. To investigate genetic loci linked to desired ornamental traits, we performed genetic analyses using an intraspecific mapping population that segregated for fruit color (yellow or red), fruit shape (round or pear), and plant height (high or compact). A genetic map was constructed with 965 single nucleotide polymorphisms (SNPs) and 33 simple sequence repeat markers. Subsequent linkage analysis using quantitative locus analysis and genome‐wide association study detected four genetic loci for the three selected traits, all of which were located near the reported genes. We performed KASP—kompetitive allele‐specific PCR—to develop SNP markers that were tightly linked to the four loci. Highly accurate genotyping data were obtained from the four SNPs across 187 F2 plants, which enabled us to select two lines with homozygous alleles for compact plant size and yellow pear‐shaped fruits. These newly developed SNP markers and genetic strategies could be used to accelerate breeding programs for ornamental tomato plants.

[This corrects the article DOI: 10.1371/journal.pgen.1008566.].

Most angiosperms bear hermaphroditic flowers, but a few species have evolved outcrossing strategies, such as dioecy, the presence of separate male and female individuals. We previously investigated the mechanisms underlying dioecy in diploid persimmon (D. lotus) and found that male flowers are specified by repression of the autosomal gene MeGI by its paralog, the Y-encoded pseudo-gene OGI. This mechanism is thought to be lineage-specific, but its evolutionary path remains unknown. Here, we developed a full draft of the diploid persimmon genome (D. lotus), which revealed a lineage-specific whole-genome duplication event and provided information on the architecture of the Y chromosome. We also identified three paralogs, MeGI, OGI and newly identified Sister of MeGI (SiMeGI). Evolutionary analysis suggested that MeGI underwent adaptive evolution after the whole-genome duplication event. Transformation of tobacco plants with MeGI and SiMeGI revealed that MeGI specifically acquired a new function as a repressor of male organ development, while SiMeGI presumably maintained the original function. Later, a segmental duplication event spawned MeGI's regulator OGI on the Y-chromosome, completing the path leading to dioecy, and probably initiating the formation of the Y-chromosome. These findings exemplify how duplication events can provide flexible genetic material available to help respond to varying environments and provide interesting parallels for our understanding of the mechanisms underlying the transition into dieocy in plants.

Phytophthora infestans (Mont.) de Bary, a hemibiotrophic oomycete, has caused severe epidemics of late blight in tomato and potato crops around the world since the Irish Potato Famine in the 1840s. Breeding of late blight resistant cultivars is one of the most effective strategies to overcome this disruptive disease. However, P. infestans is able to break down host resistance and acquire resistance to various fungicides, possibly because of the existence of high genetic variability among P. infestans isolates via sexual and asexual reproduction. Therefore, to manage this disease, it is important to understand the genetic divergence of P. infestans isolates. In this study, we analyzed the genomes of P. infestans isolates collected from Egypt and Japan using various molecular approaches including the mating type assay and genotyping simple sequence repeats, mitochondria DNA, and effector genes. We also analyzed genome-wide single nucleotide polymorphisms using double-digest restriction-site associated DNA sequencing and whole genome resequencing (WGRS). The isolates were classified adequately using high-resolution genome-wide approaches. Moreover, these analyses revealed new clusters of P. infestans isolates in the Egyptian population. Monitoring the genetic divergence of P. infestans isolates as well as breeding of resistant cultivars would facilitate the elimination of the late blight disease.

Abstract Objective In peanut, the DNA polymorphism is very low despite enormous phenotypic variations. This limits the use of genomics-assisted breeding to enhance peanut productivity. This study aimed to develop and validate new AhMITE1 and cleaved amplified polymorphic sequences (CAPS) markers. Results In total, 2957 new AhMITE1 markers were developed in addition to identifying 465 already reported markers from the whole genome re-sequencing data (WGRS) of 33 diverse genotypes of peanut. The B sub-genome (1620) showed more number of markers than the A sub-genome (1337). Distribution also varied among the chromosomes of both the sub-genomes. Further, 52.6% of the markers were from genic regions; where 31.0% were from intronic regions and 5.2% were from exonic regions. Of the 343 randomly selected markers, 82.2% showed amplification validation, with up to 35.5% polymorphism. From the SNPs on the A03, B01, B02 and B03 chromosomes, 11,730 snip-SNPs (potential CAPS sites) were identified, and 500 CAPS markers were developed from chromosome A03. Of these markers, 30.0% showed validation and high polymorphism. This study demonstrated the potential of the WGRS data to develop AhMITE1 and CAPS markers, which showed high level of validation and polymorphism. These marker resources will be useful for various genetic studies and mapping in peanut.

Abstract Background Cultivated strawberry (Fragaria × ananassa Duch.) has homoeologous chromosomes because of allo-octoploidy. For example, two homoeologous chromosomes that belong to different sub-genome of allopolyploids have similar base sequences. Thus, when conducting de novo assembly of DNA sequences, it is difficult to determine whether these sequences are derived from the same chromosome. To avoid the difficulties associated with homoeologous chromosomes and demonstrate the possibility of sequencing allopolyploids using single chromosomes, we conducted sequence analysis using microdissected single somatic chromosomes of cultivated strawberry. Results Three hundred and ten somatic chromosomes of the Japanese octoploid strawberry ‘Reiko’ were individually selected under a light microscope using a microdissection system. DNA from 288 of the dissected chromosomes was successfully amplified using a DNA amplification kit. Using next-generation sequencing, we decoded the base sequences of the amplified DNA segments, and on the basis of mapping, we identified DNA sequences from 144 samples that were best matched to the reference genomes of the octoploid strawberry, F. × ananassa, and the diploid strawberry, F. vesca. The 144 samples were classified into seven pseudo-molecules of F. vesca. The coverage rates of the DNA sequences from the single chromosome onto all pseudo-molecular sequences varied from 3 to 29.9%. Conclusion We demonstrated an efficient method for sequence analysis of allopolyploid plants using microdissected single chromosomes. On the basis of our results, we believe that whole-genome analysis of allopolyploid plants can be enhanced using methodology that employs microdissected single chromosomes.

Abstract Background The strawberry, Fragaria × ananassa, is an allo-octoploid (2n = 8x = 56) and outcrossing species. Although it is the most widely consumed berry crop in the world, its complex genome structure has hindered its genetic and genomic analysis, and thus discrimination of subgenome-specific loci among the homoeologous chromosomes is needed. In the present study, we identified candidate subgenome-specific single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) loci, and constructed a linkage map using an S1 mapping population of the cultivar ‘Reikou’ with an IStraw90 Axiom® SNP array and previously published SSR markers. Results The ‘Reikou’ linkage map consisted of 11,574 loci (11,002 SNPs and 572 SSR loci) spanning 2816.5 cM of 31 linkage groups. The 11,574 loci were located on 4738 unique positions (bin) on the linkage map. Of the mapped loci, 8999 (8588 SNPs and 411 SSR loci) showed a 1:2:1 segregation ratio of AA:AB:BB allele, which suggested the possibility of deriving loci from candidate subgenome-specific sequences. In addition, 2575 loci (2414 SNPs and 161 SSR loci) showed a 3:1 segregation of AB:BB allele, indicating they were derived from homoeologous genomic sequences. Comparative analysis of the homoeologous linkage groups revealed differences in genome structure among the subgenomes. Conclusions Our results suggest that candidate subgenome-specific loci are randomly located across the genomes, and that there are small- to large-scale structural variations among the subgenomes. The mapped SNPs and SSR loci on the linkage map are expected to be seed points for the construction of pseudomolecules in the octoploid strawberry.

Castanea sativa is an important multipurpose species in Europe for nut and timber production as well as for its role in the landscape and in the forest ecosystem. This species has low tolerance to chestnut gall wasp (Dryocosmus kuriphilus Yasumatsu), which is a pest that was accidentally introduced into Europe in early 2000 and devastated forest and orchard trees. Resistance to the gall wasp was found in the hybrid cultivar ‘Bouche de Bétizac’ (C. sativa × C. crenata) and studied by developing genetic linkage maps using a population derived from a cross between ‘Bouche de Bétizac’ and the susceptible cultivar ‘Madonna’ (C. sativa). The high-density genetic maps were constructed using double-digest restriction site-associated DNA-seq and simple sequence repeat markers. The map of ‘Bouche de Bétizac’ consisted of 1459 loci and spanned 809.6 cM; the map of ‘Madonna’ consisted of 1089 loci and spanned 753.3 cM. In both maps, 12 linkage groups were identified. A single major QTL was recognized on the ‘Bouche de Bétizac’ map, explaining up to 67–69% of the phenotypic variance of the resistance trait (Rdk1). The Rdk1 quantitative trait loci (QTL) region included 11 scaffolds and two candidate genes putatively involved in the resistance response were identified. This study will contribute to C. sativa breeding programs and to the study of Rdk1 genes.

The aim of this study was to identify candidate resistance genes for late leaf spot (LLS) and rust diseases in peanut (Arachis hypogaea L.). We used a double-digest restriction-site associated DNA sequencing (ddRAD-Seq) technique based on next-generation sequencing (NGS) for genotyping analysis across the recombinant inbred lines (RILs) derived from a cross between a susceptible line, TAG 24, and a resistant line, GPBD 4. A total of 171 SNPs from the ddRAD-Seq together with 282 markers published in the previous studies were mapped on a genetic map covering 1510.1 cM. Subsequent quantitative trait locus (QTL) analysis revealed major genetic loci for LLS and rust resistance on chromosomes A02 and A03, respectively. Heterogeneous inbred family-derived near isogenic lines and the pedigree of the resistant gene donor, A. cardenasii Krapov. & W.C. Greg., including the resistant derivatives of ICGV 86855 and VG 9514 as well as GPBD 4, were employed for whole-genome resequencing analysis. The results indicated the QTL candidates for LLS and rust resistance were located in 1.4- and 2.7-Mb genome regions on A02 and A03, respectively. In these regions, four and six resistance-related genes with deleterious mutations were selected as candidates for LLS and rust resistance, respectively. These delimited genomic regions may be beneficial in breeding programs aimed at improving disease resistance and enhancing peanut productivity.

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion's large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources.

Mitigating methane production by ruminants is a significant challenge to global livestock production. This research offers a new paradigm to reduce methane emissions from ruminants by breeding climate-clever clovers. We demonstrate wide genetic diversity for the trait methanogenic potential in Australia’s key pasture legume, subterranean clover (Trifolium subterraneum L.). In a bi-parental population the broadsense heritability in methanogenic potential was moderate (H2 = 0.4) and allelic variation in a region of Chr 8 accounted for 7.8% of phenotypic variation. In a genome-wide association study we identified four loci controlling methanogenic potential assessed by an in vitro fermentation system. Significantly, the discovery of a single nucleotide polymorphism (SNP) on Chr 5 in a defined haplotype block with an upstream putative candidate gene from a plant peroxidase-like superfamily (TSub_g18548) and a downstream lectin receptor protein kinase (TSub_g18549) provides valuable candidates for an assay for this complex trait. In this way haplotype variation can be tracked to breed pastures with reduced methanogenic potential. Of the quantitative trait loci candidates, the DNA-damage-repair/toleration DRT100-like protein (TSub_g26967), linked to avoid the severity of DNA damage induced by secondary metabolites, is considered central to enteric methane production, as are disease resistance (TSub_g26971, TSub_g26972, and TSub_g18549) and ribonuclease proteins (TSub_g26974, TSub_g26975). These proteins are good pointers to elucidate the genetic basis of in vitro microbial fermentability and enteric methanogenic potential in subterranean clover. The genes identified allow the design of a suite of markers for marker-assisted selection to reduce rumen methane emission in selected pasture legumes. We demonstrate the feasibility of a plant breeding approach without compromising animal productivity to mitigate enteric methane emissions, which is one of the most significant challenges to global livestock production.

A mapping population of recombinant inbred lines (RILs) derived from TMV 2 and its mutant, TMV 2-NLM was employed for mapping important taxonomic and productivity traits using genic and non-genic transposable element markers in peanut. Single nucleotide polymorphism and copy number variation using RAD-Sequencing data indicated very limited polymorphism between TMV 2 and TMV 2-NLM. But phenotypically they differed significantly for many taxonomic and productivity traits. Also, the RIL population showed significant variation for a few additional agronomic traits. A genetic linkage map of 1,205.66 cM was constructed using 91 genic and non-genic Arachis hypogaea transposable element (AhTE) markers. Using single marker analysis and QTL analysis, the markers with high phenotypic variance explained (PVE) were identified for branching pattern (32.3%), number of primary and secondary branches (19.9% and 28.4%, respectively), protein content (26.4%), days to 50% flowering (22.0%), content of oleic acid (15.1%), test weight (13.6%) and pod width (12.0%). Three genic markers (AhTE0357, AhTE0391, AhTE0025) with Arachis hypogaea miniature inverted-repeat transposable element (AhMITE1) activity in the genes Araip.TG1BL (B02 chromosome), Aradu.7N61X (A09 chromosome) and Aradu.7065G (A07 chromosome), respectively showed strong linkage with these taxonomic, productivity and quality traits. Since TMV 2 and TMV 2-NLM differed subtly at DNA level, the background noise in detecting the marker-trait associations was minimum; therefore, the markers identified in this study for the taxonomic and productivity traits may be significant and useful in peanut molecular breeding.

Tomato late blight caused by Phytophthora infestans (Mont.) de Bary, also known as the Irish famine pathogen, is one of the most destructive plant diseases. Wild relatives of tomato possess useful resistance genes against this disease, and could therefore be used in breeding to improve cultivated varieties. In the genome of a wild relative of tomato, Solanum habrochaites accession LA1777, we identified a new quantitative trait locus for resistance against blight caused by an aggressive Egyptian isolate of P. infestans. Using double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology, we determined 6,514 genome-wide SNP genotypes of an F2 population derived from an interspecific cross. Subsequent association analysis of genotypes and phenotypes of the mapping population revealed that a 6.8 Mb genome region on chromosome 6 was a candidate locus for disease resistance. Whole-genome resequencing analysis revealed that 298 genes in this region potentially had functional differences between the parental lines. Among of them, two genes with missense mutations, Solyc06g071810.1 and Solyc06g083640.3, were considered to be potential candidates for disease resistance. SNP and SSR markers linking to this region can be used in marker-assisted selection in future breeding programs for late blight disease, including introgression of new genetic loci from wild species. In addition, the approach developed in this study provides a model for identification of other genes for attractive agronomical traits.

Genome-wide genotyping data regarding breeding materials are essential resources for improving breeding efficiency, especially in plants with complex genomes with a high degree of polyploidy. Several current breeding efforts in cultivated peanut ( L.), which has a tetraploid genome, are devoted to developing high oleic acid cultivars. Genetic maps for such breeding programs have been developed using simple-sequence repeat (SSR) markers, the use of which requires time-consuming electrophoretic analyses. Next-generation sequencing (NGS) technology can overcome this technical hurdle. Initially, we attempted double-digest restriction site-associated DNA sequencing on peanut breeding materials used to develop high oleic acid cultivars. However, this method was not effective because few single nucleotide polymorphism (SNPs) were available because of low genetic diversity of the lines. The genome sequences of the probable diploid ancestors of cultivated peanut, Krapov. & W. C. Greg. and Krapov. & W. C. Greg., are available. Therefore, we next employed whole-genome resequencing analysis to obtain genome-wide SNP data. In this analysis, we observed large biases in the numbers and genomic positions of interspecific and intraspecific SNPs. For genome-wide genotyping, we selected a subset of SNPs covering the peanut genome as the targets of amplicon sequencing analysis. Using this technique, genome-wide genotypes of the breeding materials were easily and rapidly determined. The SNP information and analytic methods developed in this study should accelerate genetics, genomics, and breeding in peanut.

Miniature inverted-repeat transposable elements (MITEs) are ubiquitous, non-autonomous class II transposable elements. Here, we conducted genome-wide comparative analysis of 20 MITE families in B. rapa, B. oleracea, and Arabidopsis thaliana. A total of 5894 and 6026 MITE members belonging to the 20 families were found in the whole genome pseudo-chromosome sequences of B. rapa and B. oleracea, respectively. Meanwhile, only four of the 20 families, comprising 573 members, were identified in the Arabidopsis genome, indicating that most of the families were activated in the Brassica genus after divergence from Arabidopsis. Copy numbers varied from 4 to 1459 for each MITE family, and there was up to 6-fold variation between B. rapa and B. oleracea. In particular, analysis of intact members showed that whereas eleven families were present in similar copy numbers in B. rapa and B. oleracea, nine families showed copy number variation ranging from 2- to 16-fold. Four of those families (BraSto-3, BraTo-3, 4, 5) were more abundant in B. rapa, and the other five (BraSto-1, BraSto-4, BraTo-1, 7 and BraHAT-1) were more abundant in B. oleracea. Overall, 54% and 51% of the MITEs resided in or within 2 kb of a gene in the B. rapa and B. oleracea genomes, respectively. Notably, 92 MITEs were found within the CDS of annotated genes, suggesting that MITEs might play roles in diversification of genes in the recently triplicated Brassica genome. MITE insertion polymorphism (MIP) analysis of 289 MITE members showed that 52% and 23% were polymorphic at the inter- and intra-species levels, respectively, indicating that there has been recent MITE activity in the Brassica genome. These recently activated MITE families with abundant MIP will provide useful resources for molecular breeding and identification of novel functional genes arising from MITE insertion.

Vigna reflexo-pilosa, which includes a neglected crop, is the only one tetraploid species in genus Vigna. The ancestral species that make up this allotetraploid species have not conclusively been identified, although previous studies suggested that a donor genome of V. reflexo-pilosa is V. trinervia. In this study, 1,429 azuki bean EST-SSR markers were developed of which 38 EST-SSR primer pairs that amplified one product in diploid species and two discrete products in tetraploid species were selected to analyze 268 accessions from eight taxa of seven Asian Vigna species including V. reflexo-pilosa var. glabra, V. reflexo-pilosa var. reflexo-pilosa, V. exilis, V. hirtella, V. minima, V. radiata var. sublobata, V. tenuicaulis and V. trinervia to identify genome donor of V. reflexo-pilosa. Since both diploid and tetraploid species were analyzed and each SSR primer pair detected two loci in the tetraploid species, we separated genomes of the tetraploid species into two different diploid types, viz. A and B. In total, 445 alleles were detected by 38 EST-SSR markers. The highest gene diversity was observed in V. hirtella. By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes, we were able to identify the two genome donor parents of créole bean. Phylogenetic and principal coordinate analyses suggested that V. hirtella is a species complex and may be composed of at least three distinct taxa. Both analyses also clearly demonstrated that V. trinervia and one taxon of V. hirtella are the genome donors of V. reflexo-pilosa. Gene diversity indicates that the evolution rate of EST-SSRs on genome B of créole bean might be faster than that on genome A. Species relationship among the Vigna species in relation to genetic data, morphology and geographical distribution are presented.

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), '415-1', by using g-irradiated pollen. '415-1' carries the S-genotype S4dS5S5, with "d" indicating a duplication of S5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between '515-20' (S1S3) and '415-1'. Amplicon sequencing of S-RNase and SFBB18 clearly detected Shaplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRADseq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals. [ABSTRACT FROM AUTHOR]