Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.).

KEY MESSAGE: Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

Extracellular oil production by Rhodotorula paludigena BS15 for biorefinery without complex downstream processes.

To realize biomass refinery without complex downstream processes, we extensively screened for microbial strains that efficiently produce extracellular oil from sugars. Rhodotorula paludigena (formerly Rhodosporidium paludigenum) BS15 was found to efficiently produce polyol esters of fatty acids (PEFAs), which mainly comprised of 3-acetoxypalmitic acid and partially acetylated mannitol/arabinitol. To evaluate the performance of this strain, fed-batch fermentation was demonstrated on a flask scale, and 110 g/L PEFA and 103 g/L dry cells were produced in 12 days. To the best of our knowledge, the strain BS15 exhibited the highest PEFA titer (g/L) ever to be reported so far. Because the PEFA precipitated at the bottom of the culture broth, it could be easily recovered by simply discarding the upper phase. Various carbon sources can be utilized for cell growth and/or PEFA production, which signifies the potential for converting diverse biomass sources. Two different types of next-generation sequencers, Illumina HiSeq and Oxford Nanopore PromethION, were used to analyze the whole-genome sequence of the strain BS15. The integrative data analysis generated a high-quality and reliable reference genome for PEFA-producing R. paludigena. The 22.5-M base genome sequence and the estimated genes were registered in Genbank (accession numbers BQKY01000001-BQKY01000019). KEY POINTS: • R. paludigena BS15 was isolated after an extensive screening of extracellular oil producers from natural sources. • Fed-batch fermentation of R. paludigena BS15 yielded 110 g/L of PEFA, which is the highest titer ever reported to date. • Combined analysis using Illumina and Oxford Nanopore sequencers produced the near-complete genome sequence.

Starch-dependent sodium accumulation in the leaves of Vigna riukiuensis.

This research provides insight into a unique salt tolerance mechanism of Vigna riukiuensis. V. riukiuensis is one of the salt-tolerant species identified from the genus Vigna. We have previously reported that V. riukiuensis accumulates a higher amount of sodium in the leaves, whereas V. nakashimae, a close relative of V. riukiuensis, suppresses sodium allocation to the leaves. We first suspected that V. riukiuensis would have developed vacuoles for sodium sequestration, but there were no differences compared to a salt-sensitive species V. angularis. However, many starch granules were observed in the chloroplasts of V. riukiuensis. In addition, forced degradation of leaf starch by shading treatment resulted in no radio-Na (22Na) accumulation in the leaves. We performed SEM-EDX to locate Na in leaf sections and detected Na in chloroplasts of V. riukiuensis, especially around the starch granules but not in the middle of. Our results could provide the second evidence of the Na-trapping system by starch granules, following the case of common reed that accumulates starch granule at the shoot base for binding Na.

Diversity of Na+ allocation in salt-tolerant species of the genus Vigna.

Wild species in the genus Vigna are a great resource of tolerance to various stresses including salinity. We have previously screened the genetic resources of the genus Vigna and identified several accessions that have independently evolved salt tolerance. However, many aspects of such tolerance have remained unknown. Thus, we used autoradiography with radioactive sodium (22Na+) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS) to visualize and compare Na+ allocation in Vigna angularis (Willd.) Ohwi & H.Ohashi (azuki bean), Vigna nakashimae (Ohwi) Ohwi & H.Ohashi, Vigna riukiuensis (Ohwi) Ohwi & H.Ohashi, Vigna luteola (Jacq.) Benth. and Vigna marina (Burm.) Merr.. The results indicated: 1) Tolerant accessions suppress Na+ accumulation compared to azuki bean. 2) V. nakashimae and V. marina does so by accumulating higher amount of K+, whereas V. riukiuensis and V. luteola does so by other mechanisms. 3) V. luteola avoids salt-shedding by allocating excess Na+ to newly expanded leaves. As the mechanisms of the tolerant species were different, they could be piled up in a single crop via classical breeding or by genetic engineering or genome editing.

Simultaneous site-directed mutagenesis of duplicated loci in soybean using a single guide RNA.

KEY MESSAGE: Using a gRNA and Agrobacterium-mediated transformation, we performed simultaneous site-directed mutagenesis of two GmPPD loci in soybean. Mutations in GmPPD loci were confirmed in at least 33% of T2 seeds. The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated endonuclease 9 (Cas9) system is a powerful tool for site-directed mutagenesis in crops. Using a single guide RNA (gRNA) and Agrobacterium-mediated transformation, we performed simultaneous site-directed mutagenesis of two homoeologous loci in soybean (Glycine max), GmPPD1 and GmPPD2, which encode the orthologs of Arabidopsis thaliana PEAPOD (PPD). Most of the T1 plants had heterozygous and/or chimeric mutations for the targeted loci. The sequencing analysis of T1 and T2 generations indicates that putative mutation induced in the T0 plant is transmitted to the T1 generation. The inheritable mutation induced in the T1 plant was also detected. This result indicates that continuous induction of mutations during T1 plant development increases the occurrence of mutations in germ cells, which ensures the transmission of mutations to the next generation. Simultaneous site-directed mutagenesis in both GmPPD loci was confirmed in at least 33% of T2 seeds examined. Approximately 19% of double mutants did not contain the Cas9/gRNA expression construct. Double mutants with frameshift mutations in both GmPPD1 and GmPPD2 had dome-shaped trifoliate leaves, extremely twisted pods, and produced few seeds. Taken together, our data indicate that continuous induction of mutations in the whole plant and advancing generations of transgenic plants enable efficient simultaneous site-directed mutagenesis in duplicated loci in soybean.

Characterization of a heat-activated retrotransposon in Vigna angularis.

In plants, several transposable elements are conserved across species. We found a homolog of ONSEN, which is a heat-activated retrotransposon originally isolated from Arabidopsis thaliana, in Vigna. The ONSEN-like elements (VaONS) were detected in all the analyzed Japanese accessions of Vigna angularis (adzuki bean) by Southern blot analysis. However, VaONS sequences were observed to be polymorphic in the different accessions. Interestingly, extrachromosomal DNA (ecDNA) was detected in some accessions of adzuki bean, indicating the conserved heat-activation of VaONS. Furthermore, we successfully induced retrotransposition of VaONS in adzuki plant regenerated through callus. Findings of our study should provide a new tool for molecular breeding of adzuki bean.

Multiple organ gigantism caused by mutation in VmPPD gene in blackgram (Vigna mungo).

Seed size is one of the most important traits in leguminous crops. We obtained a recessive mutant of blackgram that had greatly enlarged leaves, stems and seeds. The mutant produced 100% bigger leaves, 50% more biomass and 70% larger seeds though it produced 40% less number of seeds. We designated the mutant as multiple-organ-gigantism (mog) and found the mog phenotype was due to increase in cell numbers but not in cell size. We also found the mog mutant showed a rippled leaf (rl) phenotype, which was probably caused by a pleiotropic effect of the mutation. We performed a map-based cloning and successfully identified an 8 bp deletion in the coding sequence of VmPPD gene, an orthologue of Arabidopsis PEAPOD (PPD) that regulates arrest of cell divisions in meristematic cells. We found no other mutations in the neighboring genes between the mutant and the wild type. We also knocked down GmPPD genes and reproduced both the mog and rl phenotypes in soybean. Controlling PPD genes to produce the mog phenotype is highly valuable for breeding since larger seed size could directly increase the commercial values of grain legumes.

The Vigna Genome Server, 'VigGS': A Genomic Knowledge Base of the Genus Vigna Based on High-Quality, Annotated Genome Sequence of the Azuki Bean, Vigna angularis (Willd.) Ohwi & Ohashi.

The genus Vigna includes legume crops such as cowpea, mungbean and azuki bean, as well as >100 wild species. A number of the wild species are highly tolerant to severe environmental conditions including high-salinity, acid or alkaline soil; drought; flooding; and pests and diseases. These features of the genus Vigna make it a good target for investigation of genetic diversity in adaptation to stressful environments; however, a lack of genomic information has hindered such research in this genus. Here, we present a genome database of the genus Vigna, Vigna Genome Server ('VigGS', http://viggs.dna.affrc.go.jp), based on the recently sequenced azuki bean genome, which incorporates annotated exon-intron structures, along with evidence for transcripts and proteins, visualized in GBrowse. VigGS also facilitates user construction of multiple alignments between azuki bean genes and those of six related dicot species. In addition, the database displays sequence polymorphisms between azuki bean and its wild relatives and enables users to design primer sequences targeting any variant site. VigGS offers a simple keyword search in addition to sequence similarity searches using BLAST and BLAT. To incorporate up to date genomic information, VigGS automatically receives newly deposited mRNA sequences of pre-set species from the public database once a week. Users can refer to not only gene structures mapped on the azuki bean genome on GBrowse but also relevant literature of the genes. VigGS will contribute to genomic research into plant biotic and abiotic stresses and to the future development of new stress-tolerant crops.

Unexpected consequences of a sudden and massive transposon amplification on rice gene expression.

High-copy-number transposable elements comprise the majority of eukaryotic genomes where they are major contributors to gene and genome evolution. However, it remains unclear how a host genome can survive a rapid burst of hundreds or thousands of insertions because such bursts are exceedingly rare in nature and therefore difficult to observe in real time. In a previous study we reported that in a few rice strains the DNA transposon mPing was increasing its copy number by approximately 40 per plant per generation. Here we exploit the completely sequenced rice genome to determine 1,664 insertion sites using high-throughput sequencing of 24 individual rice plants and assess the impact of insertion on the expression of 710 genes by comparative microarray analysis. We find that the vast majority of transposable element insertions either upregulate or have no detectable effect on gene transcription. This modest impact reflects a surprising avoidance of exon insertions by mPing and a preference for insertion into 5' flanking sequences of genes. Furthermore, we document the generation of new regulatory networks by a subset of mPing insertions that render adjacent genes stress inducible. As such, this study provides evidence for models first proposed previously for the involvement of transposable elements and other repetitive sequences in genome restructuring and gene regulation.

Reciprocal translocation identified in Vigna angularis dominates the wild population in East Japan.

Using an F2 population derived from cultivated and wild azuki bean, we previously detected a reciprocal translocation and a seed size QTL near the translocation site. To test the hypothesis that the translocation in the cultivated variety contributed to the larger seed size, we performed further linkage analyses with several cross combinations between cultivated and wild azuki beans. In addition, we visually confirmed the translocation by cytogenetic approach using 25 wild and cultivated accessions. As a result, we found the translocation-type chromosomes in none of the cultivated accessions, but in a number of the wild accessions. Interestingly, all the wild accessions with the translocation were originally collected from East Japan, while all the accessions with normal chromosomes were from West Japan or the Sea of Japan-side region. Such biased geographical distribution could be explained by the glacial refugium hypothesis, and supported narrowing down the domestication origin of cultivated azuki bean.

QTL mapping for salt tolerance and domestication-related traits in Vigna marina subsp. oblonga, a halophytic species.

QTL mapping in F2 population [V. luteola × V. marina subsp. oblonga] revealed that the salt tolerance in V. marina subsp. oblonga is controlled by a single major QTL. The habitats of beach cowpea (Vigna marina) are sandy beaches in tropical and subtropical regions. As a species that grows closest to the sea, it has potential to be a gene source for breeding salt-tolerant crops. We reported here for the first time, quantitative trait loci (QTLs) mapping for salt tolerance in V. marina. A genetic linkage map was constructed from an F2 population of 120 plants derived from an interspecific cross between V. luteola and V. marina subsp. oblonga. The map comprised 150 SSR markers. The markers were clustered into 11 linkage groups spanning 777.6 cM in length with a mean distance between the adjacent markers of 5.59 cM. The F2:3 population was evaluated for salt tolerance under hydroponic conditions at the seedling and developmental stages. Segregation analysis indicated that salt tolerance in V. marina is controlled by a few genes. Multiple interval mapping consistently identified one major QTL which can explain about 50% of phenotypic variance. The flanking markers may facilitate transfer of the salt tolerance allele from V. marina subsp. oblonga into related Vigna crops. The QTL for domestication-related traits from V. marina are also discussed.

mPing: The bursting transposon.

Though transposable elements (TEs) have been considered as an efficient source of evolution, it has never been possible to test this hypothesis because most of TE insertions had occurred millions of years ago, or because currently active TEs have very few copies in a host genome. However, mPing, the first active DNA transposon in rice, was revealed to hold a key to answer this question. mPing has attained high copy numbers and still retained very high activity in a traditional rice strain, which enabled direct observation of behavior and impact of a bursting TE. A comprehensive analysis of mPing insertion sites has revealed it avoids exons but prefers promoter regions and thus moderately affects transcription of neighboring genes. Some of the mPing insertions have introduced possibly useful expression profile to adjacent genes that indicated TE's potential in de novo formation of gene regulatory network.

De novo assembly of the complete organelle genome sequences of azuki bean (Vigna angularis) using next-generation sequencers.

Since chloroplasts and mitochondria are maternally inherited and have unique features in evolution, DNA sequences of those organelle genomes have been broadly used in phylogenetic studies. Thanks to recent progress in next-generation sequencer (NGS) technology, whole-genome sequencing can be easily performed. Here, using NGS data generated by Roche GS Titanium and Illumina Hiseq 2000, we performed a hybrid assembly of organelle genome sequences of Vigna angularis (azuki bean). Both the mitochondrial genome (mtDNA) and the chloroplast genome (cpDNA) of V. angularis have very similar size and gene content to those of V. radiata (mungbean). However, in structure, mtDNA sequences have undergone many recombination events after divergence from the common ancestor of V. angularis and V. radiata, whereas cpDNAs are almost identical between the two. The stability of cpDNAs and the variability of mtDNAs was further confirmed by comparative analysis of Vigna organelles with model plants Lotus japonicus and Arabidopsis thaliana.

How to Sequence and Assemble Plant Genomes.

Although nanopore sequencer is a great tool, many plant scientists have suffered from bad sequencing results, even though they have exactly followed the official protocol in preparing a library. This is because the protocol is not optimized for plant genomic DNA. The protocol may be good for sequencing animal or bacterial genomes, but not for plants. However, if the protocol is properly modified, one can obtain lots of long reads and achieve a telomere-to-telomere assembly. Here I present a protocol to that end.

Domesticating Vigna stipulacea: Chromosome-Level genome assembly reveals VsPSAT1 as a candidate gene decreasing hard-seededness.

To increase food production under the challenges presented by global climate change, the concept of de novo domestication-utilizing stress-tolerant wild species as new crops-has recently gained considerable attention. We had previously identified mutants with desired domestication traits in a mutagenized population of the legume Vigna stipulacea Kuntze (minni payaru) as a pilot for de novo domestication. Given that there are multiple stress-tolerant wild legume species, it is important to establish efficient domestication processes using reverse genetics and identify the genes responsible for domestication traits. In this study, we identified VsPSAT1 as the candidate gene responsible for decreased hard-seededness, using a Vigna stipulacea isi2 mutant that takes up water from the lens groove. Scanning electron microscopy and computed tomography revealed that the isi2 mutant has lesser honeycomb-like wax sealing the lens groove than the wild-type, and takes up water from the lens groove. We also identified the pleiotropic effects of the isi2 mutant: accelerating leaf senescence, increasing seed size, and decreasing numbers of seeds per pod. While doing so, we produced a V. stipulacea whole-genome assembly of 441 Mbp in 11 chromosomes and 30,963 annotated protein-coding sequences. This study highlights the importance of wild legumes, especially those of the genus Vigna with pre-existing tolerance to biotic and abiotic stresses, for global food security during climate change.

Genetic structure of landraces in foxtail millet (Setaria italica (L.) P. Beauv.) revealed with transposon display and interpretation to crop evolution of foxtail millet.

Although the origin and domestication process of foxtail millet (Setaria italica subsp. italica (L.) P. Beauv.) has been studied by several groups, the issue is still ambiguous. It is essential to resolve this issue by studying a large number of accessions with sufficient markers covering the entire genome. Genetic structures were analyzed by transposon display (TD) using 425 accessions of foxtail millet and 12 of the wild ancestor green foxtail (Setaria italica subsp. viridis (L.) P. Beauv.). We used three recently active transposons (TSI-1, TSI-7, and TSI-10) as genome-wide markers and succeeded in demonstrating geographical structures of the foxtail millet. A neighbor-joining dendrogram based on TD grouped the foxtail millet accessions into eight major clusters, each of which consisted of accessions collected from adjacent geographical areas. Eleven out of 12 green foxtail accessions were grouped separately from the clusters of foxtail millet. These results indicated strong regional differentiations and a long history of cultivation in each region. Furthermore, we discuss the relationship between foxtail millet and green foxtail and suggest a monophyletic origin of foxtail millet domestication.

Whole genome sequencing data of Vigna nakashimae var. Ukushima and G418.

Vigna nakashimae is one of the closely related species of Vigna angularis (Adzuki bean). Two strain of 'Ukushima' and 'G418' were identified as salt tolerance strains in Vigna nakashimae from gene bank collection. F2 populations from an inter- or intra-specific cross between the sensitive and tolerant strains are useful for the detection of salt tolerance QTL in Vigna nakashimae. Although Vigna angularis reference genome is available and useful for genetic analysis by genotyping-by-sequencing/RADseq in closely related species, it is not enough for isolation of responsible genes. To reveal sequence variation in Vigna nakashimae "Ukushima" and "G418", the whole genome sequencing was performed using Illumina HiSeq X Ten system (411,174,986 and 478,116,282 read). NGS data was deposited in the DNA Data Bank of Japan (DDBJ) under accession number DRA009307.

Mapping patterns of abiotic and biotic stress resilience uncovers conservation gaps and breeding potential of Vigna wild relatives.

This study provides insights in patterns of distribution of abiotic and biotic stress resilience across Vigna gene pools to enhance the use and conservation of these genetic resources for legume breeding. Vigna is a pantropical genus with more than 88 taxa including important crops such as V. radiata (mung bean) and V. unguiculata (cowpea). Our results show that sources of pest and disease resistance occur in at least 75 percent of the Vigna taxa, which were part of screening assessments, while sources of abiotic stress resilience occur in less than 30 percent of screened taxa. This difference in levels of resilience suggests that Vigna taxa co-evolve with pests and diseases while taxa are more conservative to adapt to climatic changes and salinization. Twenty-two Vigna taxa are poorly conserved in genebanks or not at all. This germplasm is not available for legume breeding and requires urgent germplasm collecting before these taxa extirpate on farm and in the wild. Vigna taxa, which tolerate heat and drought stress are rare compared with taxa, which escape these stresses because of short growing seasons or with taxa, which tolerate salinity. We recommend prioritizing these rare Vigna taxa for conservation and screening for combined abiotic and biotic stress resilience resulting from stacked or multifunctional traits. The high presence of salinity tolerance compared with drought stress tolerance, suggests that Vigna taxa are good at developing salt-tolerant traits. Vigna taxa are therefore of high value for legume production in areas that will suffer from salinization under global climate change.