Published papers

Key message: A putative ketoacyl-ACP reductase (CaKR1) that was not previously known to be associated with pungency of Capsicum was identified from map-based cloning and functional characterization. Abstract: The pungency of chili pepper fruits is due to the presence of capsaicinoids, which are synthesized through the convergence of the phenylpropanoid and branched-chain fatty acid pathways. The extensive, global use of pungent and non-pungent peppers underlines the importance of understanding the genetic mechanism underlying capsaicinoid biosynthesis for breeding pepper cultivars. Although Capsicum is one of the earliest domesticated plant genera, the only reported genetic causes of its loss of pungency are mutations in acyltransferase (Pun1) and putative aminotransferase (pAMT). In this study, a single recessive gene responsible for the non-pungency of pepper No.3341 (C. chinense) was identified on chromosome 10 using an F 2 population derived from a cross between Habanero and No.3341. Five candidate genes were identified in the target region, within a distance of 220 kb. A candidate gene, a putative ketoacyl-ACP reductase (CaKR1), of No.3341 had an insertion of a 4.5-kb transposable element (TE) sequence in the first intron, resulting in the production of a truncated transcript missing the region coding the catalytic domain. Virus-induced gene silencing of CaKR1 in pungent peppers resulted in the decreased accumulation of capsaicinoids, a phenotype consistent with No.3341. Moreover, GC–MS analysis of 8-methyl-6-nonenoic acid, which is predicted to be synthesized during the elongation cycle of branched-chain fatty acid biosynthesis, revealed that its deficiency in No.3341. Genetic, genomic, transcriptional, silencing, and biochemical precursor analyses performed in combination provide a solid ground for the conclusion that CaKR1 is involved in capsaicinoid biosynthesis and that its disruption results in a loss of pungency.

An increasing number of field studies have shown that the phenotype of an individual plant depends not only on its genotype but also on those of neighboring plants; however, this fact is not taken into consideration in genome-wide association studies (GWAS). Based on the Ising model of ferromagnetism, we incorporated neighbor genotypic identity into a regression model, named “Neighbor GWAS”. Our simulations showed that the effective range of neighbor effects could be estimated using an observed phenotype when the proportion of phenotypic variation explained (PVE) by neighbor effects peaked. The spatial scale of the first nearest neighbors gave the maximum power to detect the causal variants responsible for neighbor effects, unless their effective range was too broad. However, if the effective range of the neighbor effects was broad and minor allele frequencies were low, there was collinearity between the self and neighbor effects. To suppress the false positive detection of neighbor effects, the fixed effect and variance components involved in the neighbor effects should be tested in comparison with a standard GWAS model. We applied neighbor GWAS to field herbivory data from 199 accessions of Arabidopsis thaliana and found that neighbor effects explained 8% more of the PVE of the observed damage than standard GWAS. The neighbor GWAS method provides a novel tool that could facilitate the analysis of complex traits in spatially structured environments and is available as an R package at CRAN (https://cran.rproject.org/package=rNeighborGWAS).

Phenotypes of sessile organisms, such as plants, rely not only on their own genotypes but also on those of neighboring individuals. Previously, we incorporated such neighbor effects into a single-marker regression using the Ising model of ferromagnetism. However, little is known regarding how neighbor effects should be incorporated in quantitative trait locus (QTL) mapping. In this study, we propose a new method for interval QTL mapping of neighbor effects, designated “neighbor QTL,” the algorithm of which includes: (1) obtaining conditional self-genotype probabilities with recombination fraction between flanking markers; (2) calculating conditional neighbor genotypic identity using the self-genotype probabilities; and (3) estimating additive and dominance deviations for neighbor effects. Our simulation using F2 and backcross lines showed that the power to detect neighbor effects increased as the effective range decreased. The neighbor QTL was applied to insect herbivory on Col × Kas recombinant inbred lines of Arabidopsis thaliana. Consistent with previous results, the pilot experiment detected a self-QTL effect on the herbivory at the GLABRA1 locus. Regarding neighbor QTL effects on herbivory, we observed a weak QTL on the top of chromosome 4, at which a weak self-bolting QTL was also identified. The neighbor QTL method is available as an R package (https://cran.r-project.org/package¼rNeighborQTL), providing a novel tool to investigate neighbor effects in QTL studies.

Three new taxa, Rhododendron sohayakiense Y.Watan. & T.Yukawa (Ericaceae), and its two varieties, var. kiusianum Y.Watan., T.Yukawa & T.Minamitani and var. koreanum Y.Watan. & T.Yukawa are described and illustrated from Japan and South Korea. They can be distinguished from each other and from the other members of the R. tschonoskii alliance, i.e. R. tschonoskii, R. tetramerum, R. trinerve and R. tsusiophyllum, through their combination of leaf size, leaf morphologies including lateral nerves on abaxial leaf surface, corolla morphologies including number of corolla lobes, style length and anther form. Phylogenetic inferences based on chloroplast DNA and genome-wide sequences revealed that each of the three new taxa is monophyletic and they further form a clade. Distributions of the three taxa are also clearly separated from each other and also from the other members of the R. tschonoskii alliance.

Niche conservatism is the tendency of lineages to retain the same niche as their ancestors. It constrains biological groups and prevents ecological divergence. However, theory predicts that niche conservatism can hinder gene flow, strengthen drift and increase local adaptation: does it mean that it also can facilitate speciation? Why does this happen? We aim to answer these questions. We examined the variation of chloroplast DNA, genome-wide single nucleotide polymorphisms, morphological traits and environmental variables across the Dioon merolae cycad populations. We tested geographical structure, scenarios of demographic history, and niche conservatism between population groups. Lineage divergence is associated with the presence of a geographical barrier consisting of unsuitable habitats for cycads. There is a clear genetic and morphological distinction between the geographical groups, suggesting allopatric divergence. However, even in contrasting available environmental conditions, groups retain their ancestral niche, supporting niche conservatism. Niche conservatism is a process that can promote speciation. In D. merolae, lineage divergence occurred because unsuitable habitats represented a barrier against gene flow, incurring populations to experience isolated demographic histories and disparate environmental conditions. This study explains why cycads, despite their ancient lineage origin and biological stasis, have been able to diversify into modern ecosystems worldwide.

In 2017, a leaf sample from a single chili pepper (Capsicum annuum) plant exhibiting yellowing was collected from Aceh province, Indonesia. Total RNA was extracted from this sample, and RNA-Seq analysis was conducted. Putative infecting viruses were detected by mapping the obtained reads to the full-length viral genome sequences available in the GenBank database (7457 sequences) and the de novo-assembled contigs. RNA-Seq analysis detected polerovirus, begomovirus, and amalgavirus sequences, and the polerovirus-like sequences showed strong similarity to those of previously reported pepper vein yellows viruses (PeVYVs). The complete viral genome sequence obtained by RT-PCR had a length of 6023 nt, had the typical genome organization of a polerovirus and showed a high degree of sequence similarity to PeVYV-2 from Israel. Moreover, the predicted amino acid sequence of the P0 protein of the Indonesian isolate was 85.1% to 88.8% identical to those of other PeVYVs. In accordance with the polerovirus species demarcation criteria, this isolate should be assigned to a new polerovirus species, and we propose the name “pepper vein yellows virus 9” (PeVYV-9) for this virus.

The economically and nutritionally important genus Citrus belongs to the subfamily Aurantioideae in the family Rutaceae. Here, we analyzed the phylogenetic relationships of the subfamily Aurantioideae based on RAD-Seq. The RAD-Seq data produced phylogenetic trees with high support values, clear discriminations based on branch length, and elucidations of early branching events. Our genetic classification corresponded well with the classical morphological classification system and supported the subdivision of Citreae, one of two tribes of the Aurantioideae, into three subtribes—Triphasiinae, Citrinae, and Balsamocitrinae. Additionally, it was largely consistent with the subdivision of Clauseneae, the other tribe of the Aurantioideae, into three subtribes—Micromelinae, Clauseninae, and Merrillinae; the exception was Murraya paniculata. With the exception of members of primitive citrus fruit trees, namely, Severinia buxifolia and Hesperethusa crenulata, lower-level morphological groupings under subtribes based on genetic and morphological classifications corresponded well. The phylogenetic relationship between Asian “true citrus fruit trees” (genera Citrus, Poncirus, and Fortunella) and Australian/New Guinean citrus fruit trees (genera Microcitrus, Eremocitrus, and Clymenia) was inconsistent between present classification based mainly on the nuclear genome and the previous classification based on the chloroplast genome. This inconsistency may be explained by chloroplast capture. Our findings provide a valuable insight into the genetic relationships of the subfamily Aurantioideae in the family Rutaceae.

Phenotypic polymorphism within a species is a notable phenomenon in evolutionary biology to understand the process of adaptive speciation and other historical events. The Saxifraga fortunei complex is a widespread herb found in East Asia. It includes several ecotypic taxa corresponding to their habitat environments. The distribution of the various ecotypes in a limited area of the Japanese Archipelago makes the species a suitable model to investigate the impact of population demographic history and natural selection on lineage diversification. Here, Sanger-based sequencing was used to estimate the divergence timeframe between populations of the Eurasian continent and Japan. Genome-wide SNPs obtained by ddRAD sequencing were used to investigate the phylogeographic origins of ecotypic taxa. The phylogenetic analyses revealed the divergence of the Japanese population from the continental population in the late Miocene. Two distinct regional clades of North and South Japan were identified; phenotypic diversification was evident only in the southern clade. The South Japan clades displayed a historical distribution expansion from north to south. The phenotypic variations appeared to have generated during the expansion. The ecotypic boundaries were incongruent with the genetic grouping. We propose that morphological and ecological specialization in Japanese populations was repeatedly generated by local natural selection.

• Background and Aims: The processes and mechanisms underlying lineage diversification are major topics in evolutionary biology. Eurasian goldenrod species of the Solidago virgaurea complex show remarkable morphological and ecological diversity in the Japanese Archipelago, with ecotypic taxa well adapted to specific environments (climate, edaphic conditions and disturbance regimes). The species complex is a suitable model to investigate the evolutionary processes of actively speciating plant groups, due to its ability to evolve in relation to environmental adaptation and its historical population dynamics. • Methods: Two chloroplast markers, 18 nuclear microsatellite markers and ddRAD-sequencing were used to infer population genetic demography of S. virgaurea complex with its related species/genera. • Key Results: Our analysis showed that populations in Japan form an evolutionary unit, which was genetically diverged from adjacent continental populations. The phylogenetic structure within the archipelago strongly corresponds to the geography, but interestingly there is no concordance between genetic structure and ecotypic boundaries; neighbouring populations of distinct ecotypes share a genetic background. • Conclusions: We propose that the traits specific to the ecotypic entities are maintained by natural selection or are very recently generated and have little effect on the genomes, making genome-wide genetic markers unsuitable for detecting ecotypic differentiation. Furthermore, some sporadically distributed taxa (found as rheophytes and alpine plants) were repeatedly generated from a more widespread taxon in geographically distant areas by means of selection. Overall, this study showed that the goldenrod complex has a high ability to evolve, enabling rapid ecological diversification over a recent timeframe.

Rhodeus ocellatus kurumeus is a small cyprinid fish endemic to the Japanese archipelago. This fish frequently hybridizes with an exotic subspecies Rhodeus ocellatus ocellatus, and non-introgressed populations of R. o. kurumeus have greatly decreased. Previous studies based on mitochondrial and nuclear DNA identified introgressed populations and inferred phylogenetic relationships, but these approaches may lead to underestimates of introgression or ambiguous results owing to the small number of genetic markers used. The present study applied double digest restriction site-associated DNA (ddRAD) sequencing to assignment tests and phylogenetic inferences. Data matrices can vary in size and content depending on the strategies used to process ddRAD sequences; therefore, we prepared 25 data matrices with various processing strategies. The assignment tests based on the 25 data matrices resulted in similar assignment patterns and provided evidence for introgression in some populations that had been identified as non-introgressed populations in previous studies. The maximum-likelihood (ML) phylogenetic trees differed in robustness and topology among the 25 matrices; interestingly, ML trees with low statistical support reflected geographical distributions of fish better than ML trees with high statistical support. These results indicate that ddRAD sequencing can detect introgression with greater sensitivity than conventional DNA markers, that ddRAD sequencing is useful for phylogenetic inference among closely related populations within a subspecies, and that statistically robust ML trees do not necessarily reflect the true phylogeny. Application of ddRAD data to conservation genetics and evolutionary history of this subspecies are also discussed.

Streptocarpus shows a great diversity of shoot and flower morphology, but lacks resources for genetic studies. RAD sequencing with NGS technologies provides large amounts of genotyping data quickly and cost effectively and is very useful for genetic map construction. In this study, three approaches for constructing RAD loci were combined to build the first genetic map for Streptocarpus (Gesneriaceae). A backcross population between Streptocarpus grandis and S. rexii was used as the mapping population. The parents are both diploid (2n = 32) and largely autogamous. Streptocarpus grandis is unifoliate, retaining the enlarged macrocotyledon as sole foliar organ which is monocarpic, and S. rexii is rosulate and perennial. 233 backcross plants were genotyped, and the data from 200 or 150 individuals used to generate the maps. The RAD reads were assembled using a de novo approach, aligned on the reference draft genome sequences of S. rexii using BWA or Stampy, and finally all markers were analyzed combined. In the latter, 599 mappable markers were generated which were assigned to 16 linkage groups, with an average inter-marker distance of 2.6 cM. This work provides a valuable genetic resource for further studies in the large and morphologically diverse genus Streptocarpus.

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential w-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization.

Key message: A begomovirus resistance gene pepy-1, which encodes the messenger RNA surveillance factor Pelota, was identified in pepper (C. annuum) through map-based cloning and functional characterization. Abstract: Pepper yellow leaf curl disease caused by begomoviruses seriously affects pepper (Capsicum spp.) production in a number of regions around the world. Ty genes of tomato, which confer resistance to the tomato yellow leaf curl virus, are the only begomovirus resistance genes cloned to date. In this study, we focused on the identification of begomovirus resistance genes in Capsicum annuum. BaPep-5 was identified as a novel source of resistance against pepper yellow leaf curl Indonesia virus (PepYLCIV) and pepper yellow leaf curl Aceh virus (PepYLCAV). A single recessive locus, which we named as pepper yellow leaf curl disease virus resistance 1 (pepy-1), responsible for PepYLCAV resistance in BaPep-5 was identified on chromosome 5 in an F2 population derived from a cross between BaPep-5 and the begomovirus susceptible accession BaPep-4. In the target region spanning 34 kb, a single candidate gene, the messenger RNA surveillance factor Pelota, was identified. Whole-genome resequencing of BaPep-4 and BaPep-5 and comparison of their genomic DNA sequences revealed a single nucleotide polymorphism (A to G) located at the splice site of the 9th intron of CaPelota in BaPep-5, which caused the insertion of the 9th intron into the transcript, resulting in the addition of 28 amino acids to CaPelota protein without causing a frameshift. Virus-induced gene silencing of CaPelota in the begomovirus susceptible pepper No.218 resulted in the gain of resistance against PepYLCIV, a phenotype consistent with BaPep-5. The DNA marker developed in this study will greatly facilitate marker-assisted breeding of begomovirus resistance in peppers.

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N50 = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.

Studies on plant viruses are biased towards crop diseases and little is known about viruses in natural vegetation. We conducted extensive surveys of plant viruses in wild Brassicaceae plants occurring in three local plant communities in central Japan. We applied RNA-Seq with selective depletion of rRNA, which allowed us to detect infections of all genome-reported viruses simultaneously. Infections of Turnip mosaic virus (TuMV), Cucumber mosaic virus (CMV), Brassica yellows virus, Pelargonium zonate spot virus, and Arabidopsis halleri partitivirus 1 were detected from the two perennial species, Arabidopsis halleri subsp. gemmifera and Rorippa indica. De novo assembly further detected partial sequences of a putative novel virus in Arabis fragellosa. Virus species composition and infection rate differed depending on site and plant species. Viruses were most frequently detected from the perennial clonal plant, A. halleri, in which a high clonal transmission rate of viruses across multiple years was confirmed. Phylogenetic analysis of TuMV and CMV showed that virus strains from wild Brassicaceae were included as a major clade of these viruses with other reported strains from crop plants, suggesting that viruses were shared among wild plants and crops. Our studies indicated that distribution of viruses in natural plant populations are determined by the combinations of life histories of viruses and hosts. Revealing viral distribution in the natural plant communities improves our knowledge on the ecology of plant viruses.

To obtain genetic information about the germplasm of tea (Camellia sinensis L.) in Japan, 167 accessions including 138 var. sinensis (96 Japanese var. sinensis and 42 exotic var. sinensis) and 29 Assam hybrids were analyzed using single nucleotide polymorphisms (SNPs) markers identified by double-digest restriction-site-associated DNA sequencing (ddRAD-seq) analysis. Approximately 10,000 SNPs were identified by ddRAD-seq and were mapped across the whole genome. The 167 tea accessions were classified into three genetic subgroups: (1) Japanese var. sinensis; (2) Japanese and exotic var. sinensis; (3) Assam hybrids and exotic var. sinensis. Leaf morphology varied widely within each genetic subgroups. The 96 Japanese var. sinensis were classified into four genetic subgroups as follows; two subgroups of Shizuoka (the largest tea production region) landraces, Uji (most ancient tea production region) landraces, and the pedigree of ‘Yabukita’, the leading green tea cultivar in Japan. These results indicated that the SNP markers obtained from ddRAD-seq are a useful tool to investigate the geographical background and breeding history of Japanese tea. This genetic information revealed the ancestral admixture situation of the ‘Yabukita’ pedigree, and showed that the genome structure of ‘Yabukita’ is clearly different from those of other Japanese accessions.

Orobanche foetida Poiret is a holoparasitic plant that lacks chlorophyll and totally depending on its host for its growth. Orobanche foetida parasitizes host plant roots and extract nutrient and water via a haustorium. Although O. foetida distributes in the Mediterranean region as a wild plant parasite, it parasitizes faba bean causing serious damages which may reach 90% yield losses in Tunisia. Analysis of genetic diversity of the parasite is important to better understand its evolution and spread, remained largely unknown. In this work, we present the first study on genetic diversity and population structure using the robust technique Restriction-site-Associated DNA sequencing (RADseq) for Orobanche spp. We collected 244 samples of O. foetida from 18 faba bean fields in the north of Tunisia including 17 populations from the north-west and one population form the north-east. To overcome the difficulty of SNP discovery in O. foetida genome as a non-model and tetraploid plant, we utilized three different informatics pipelines, namely UNEAK, pyRAD and Stacks. This study showed that genetic differentiation occurred in the Tunisian O. foetida emphasizing the isolation by distance effect. However, no strong population clustering was detected in this work basing on the three data sets and clustering methods used. The present study shed the light on the current distribution and the genetic variation situation of the fetid broomrape in Tunisia, highlighting the importance of understanding the evolution of this parasite and its genetic background. This will aid in developing efficient strategies to control this parasite and its expansion in Tunisia and worldwide.

As most organisms have evolved in seasonal environments, their environmental responses should be adapted to seasonal transitions. Here we show that the combination of temperature and day length shapes the seasonal dynamics of the transcriptome and adaptation to seasonal environments in a natural habitat of a perennial plant Arabidopsis halleri subsp. gemmifera. Weekly transcriptomes for two years and bihourly diurnal transcriptomes on the four equinoxes/solstices, identified 2,879 and 7,185 seasonally- and diurnally-oscillating genes, respectively. Dominance of annual temperature changes for defining seasonal oscillations of gene expressions was indicated by controlled environment experiments manipulating the natural 1.5-month lag of temperature behind day length. We found that plants have higher fitness in ‘natural’ chambers than in ‘unnatural’ chambers simulating in-phase and anti-phase oscillations between temperature and day length. Seasonal temperature responses were disturbed in unnatural chambers. Our results demonstrate how plants use multiple types of environmental information to adapt to seasonal environments.

Objective: Selective breeding for desirable traits is becoming popular in aquaculture. In Miyagi prefecture, Japan, a selectively bred population of Coho salmon (Oncorhynchus kisutch) has been established with the original, randomly breeding population maintained separately. Since they have been bred without family records, the genetic diversity within these populations remains unknown. In this study, we estimated the genetic diversity and key quantitative genetic parameters such as heritability and genomic breeding value for body size traits by means of genomic best linear unbiased prediction to assess the genetic health of these populations. Results: Ninety-nine and 83 females from the selective and random groups, respectively, were genotyped at 2350 putative SNPs by means of double digest restriction associated DNA sequencing. The genetic diversity in the selectively bred group was low, as were the estimated heritability and prediction accuracy for length and weight (h 2 = 0.26-0.28; accuracy = 0.34), compared to the randomly bred group (h 2 = 0.50-0.60; accuracy = 0.51-0.54). Although the tested sample size was small, these results suggest that further selection is difficult for the selectively bred population, while there is some potential for the randomly bred group, especially with the aid of genomic information.

Plants sense and respond to light via multiple photoreceptors including phytochrome. The decreased ratio of red to far-red light that occurs under a canopy triggers shade-avoidance responses, which allow plants to compete with neighboring plants. The leaf acts as a photoperceptive organ in this response. In this study, we investigated how the shade stimulus is spatially processed within the cotyledon. We performed transcriptome analysis on microtissue samples collected from vascular and nonvascular regions of Arabidopsis (Arabidopsis thaliana) cotyledons. In addition, we mechanically isolated and analyzed the vascular tissue. More genes were up-regulated by the shade stimulus in vascular tissues than in mesophyll and epidermal tissues. The genes up-regulated in the vasculature were functionally divergent and included many auxin-responsive genes, suggesting that various physiological/developmental processes might be controlled by shade stimulus in the vasculature. We then investigated the spatial regulation of these genes in the vascular tissues. A small vascular region within a cotyledon was irradiated with far-red light, and the response was compared with that when the whole seedling was irradiated with far-red light. Most of the auxin-responsive genes were not fully induced by the local irradiation, suggesting that perception of the shade stimulus requires that a wider area be exposed to far-red light or that a certain position in the mesophyll and epidermis of the cotyledon be irradiated. This result was consistent with a previous report that auxin synthesis genes are up-regulated in the periphery of the cotyledon. Hence, auxin acts as an important intraorgan signaling factor that controls the vascular shade response within the cotyledon.