Complete Genome Sequence of Annamia dubia, filamentous colony-making Chroococcales with the analysis of FraC gene influencing filament integrity.

The complete genome of Annamia dubia was sequenced. The genome size is 4.02 Mbp, including 3886286 bp circular chromosome and four circular plasmids (31516, 42453, 38085 and 24903 bp). It included 3718 protein-coding sequences, 45 tRNA genes, three sets of rRNA genes, a microcystin biosynthesis gene cluster and six CRISPR (clustered regularly interspaced short palindromic repeat). Annamia is the only one genus in the Chroococcales that makes filamentous colonies. FraC and FraG were identified in the genome. These genes are required for the integrity of cell junctions and influencing filament integrity and are thought to be related to colony formation. These genes are first reported from Chroococcales, and may play a significant role in the colony formation of this species. In the phylogenetic tree of the FraC gene, A. dubia was located in the basal position of Oscillatoriales. The GC ratio of FraC gene of A. dubia is very low from the genome and the FraC gene of Microcoleaceae. The presence of these genes in the basal region and the low GC ratio suggests that the FraC gene in this species was introduced by horizontal gene transfer. Since the filamentous colony is a fundamental and important taxonomic feature for the classification of cyanobacteria, the possibility of horizontal transmission of genes involved in filamentous cyanobacterial colonies is an important discovery for the classification of cyanobacteria.

DNA metabarcoding focusing on the plankton community: an effective approach to reconstruct the paleo-environment.

DNA metabarcoding (DNA-MB) targeting the whole plankton community is a promising approach in studies of sediment samples from water bodies, but its effectiveness in ancient material is not well demonstrated. We applied DNA-MB of plankton in a sediment core to reconstruct the paleo-environment of Lake Shinji, Japan, through a marine lagoon/freshwater lake transition during the past 2300 years. We interpreted core-sample plankton taxonomy and habitat by reference to the modern plankton community in water samples. OTUs (operational taxonomic units) belonging to Dictyochophyceae were 81.05% of the total reads in sediments. However, Ciliophora, Copepoda and Labyrinthulea formed the majority of plankton taxa in the water samples, suggesting that they are under-represented in sediment. A drastic change in plankton composition correlated with a large decrease in sediment sulfur concentration, implying the change of aquatic environment from marine lagoon to freshwater lake. This event took place ca. 1200 CE in Lake Shinji. A 250 year-long transitional period followed, during which the total DNA sequence reads were very low. This suggests that salinity fluctuations created a hostile environment for both marine and freshwater plankton species. Our results show that DNA-MB of the whole plankton community is effective in reconstructing paleo-environments.

Development of 12 sets of chromosome segment substitution lines that enhance allele mining in Asian cultivated rice.

Many agronomic traits that are important in rice breeding are controlled by multiple genes. The extensive time and effort devoted so far to identifying and selecting such genes are still not enough to target multiple agronomic traits in practical breeding in Japan because of a lack of suitable plant materials in which to efficiently detect and validate beneficial alleles from diverse genetic resources. To facilitate the comprehensive analysis of genetic variation in agronomic traits among Asian cultivated rice, we developed 12 sets of chromosome segment substitution lines (CSSLs) with the japonica background, 11 of them in the same genetic background, using donors representing the genetic diversity of Asian cultivated rice. Using these materials, we overviewed the chromosomal locations of 1079 putative QTLs for seven agronomic traits and their allelic distribution in Asian cultivated rice through multiple linear regression analysis. The CSSLs will allow the effects of putative QTLs in the highly homogeneous japonica background to be validated.

DNA metabarcoding focused on difficult-to-culture protists: An effective approach to clarify biological interactions.

DNA metabarcoding on a single organism is a promising approach to clarify the biological interactions (e.g., predator-prey relationships and symbiosis, including parasitism) of difficult-to-culture protists. To evaluate the effectiveness of this method, Radiolaria and Phaeodaria, which are ecologically important protistan groups, were chosen as target taxa. DNA metabarcoding on a single organism focused on the V9 region of the 18S rRNA gene revealed potential symbionts, parasites and food sources of Radiolaria and Phaeodaria. Previously reported hosts and symbionts (parasites) were detected, and newly recognized combinations were also identified. The contained organisms largely differed between Radiolaria and Phaeodaria. In Radiolaria, members of the same order tended to contain similar organisms, and the taxonomic composition of possible symbionts, parasites, and food sources was fixed at the species level. Members of the same phaeodarian family, however, did not contain similar organisms, and body part (i.e., the central capsule or the phaeodium) was the most important factor that divided the taxonomic composition of detected organisms, implying that the selection of appropriate body part is important when trying to ascertain contained organisms, even for unicellular zooplankton. Our results show that DNA metabarcoding on a single organism is effective in revealing the biological interactions of difficult-to-culture protists.

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat.

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

Development of a High-Quality/Yield Long-Read Sequencing-Adaptable DNA Extraction Method for Crop Seeds.

Genome sequencing is important for discovering critical genes in crops and improving crop breeding efficiency. Generally, fresh, young leaves are used for DNA extraction from plants. However, seeds, the storage form, are more efficient because they do not require cultivation and can be ground at room temperature. Yet, only a few DNA extraction kits or methods suitable for seeds have been developed to date. In this study, we introduced an improved (IMP) Boom method that is relatively low-cost, simple to operate, and yields high-quality DNA that can withstand long-read sequencing. The method successfully extracted approximately 8 µg of DNA per gram of seed weight from soybean seeds at an average concentration of 48.3 ng/µL, approximately 40-fold higher than that extracted from seeds using a common extraction method kit. The A260/280 and A260/230 values of the DNA were 1.90 and 2.43, respectively, which exceeded the respective quality thresholds of 1.8 and 2.0. The DNA also had a DNA integrity number value (indicating the degree of DNA degradation) of 8.1, higher than that obtained using the kit and cetyltrimethylammonium bromide methods. Furthermore, the DNA showed a read length N50 of 20.96 kbp and a maximum read length of 127.8 kbp upon long-read sequencing using the Oxford Nanopore sequencer, with both values being higher than those obtained using the other methods. DNA extracted from seeds using the IMP Boom method showed an increase in the percentage of the nuclear genome with a decrease in the relative ratio of chloroplast DNA. These results suggested that the proposed IMP Boom method can extract high-quality and high-concentration DNA that can be used for long-read sequencing, which cannot be achieved from plant seeds using other conventional DNA extraction methods. The IMP Boom method could also be adapted to crop seeds other than soybeans, such as pea, okra, maize, and sunflower. This improved method is expected to improve the efficiency of various crop-breeding operations, including seed variety determination, testing of genetically modified seeds, and marker-assisted selection.

Phosphate starvation response precedes abscisic acid response under progressive mild drought in plants.

Drought severely damages crop production, even under conditions so mild that the leaves show no signs of wilting. However, it is unclear how field-grown plants respond to mild drought. Here, we show through six years of field trials that ridges are a useful experimental tool to mimic mild drought stress in the field. Mild drought reduces inorganic phosphate levels in the leaves to activate the phosphate starvation response (PSR) in soybean plants in the field. Using Arabidopsis thaliana and its mutant plants grown in pots under controlled environments, we demonstrate that PSR occurs before abscisic acid response under progressive mild drought and that PSR plays a crucial role in plant growth under mild drought. Our observations in the field and laboratory using model crop and experimental plants provide insight into the molecular response to mild drought in field-grown plants and the relationship between nutrition and drought stress response.

Unique Salt-Tolerance-Related QTLs, Evolved in Vigna riukiuensis (Na+ Includer) and V. nakashimae (Na+ Excluder), Shed Light on the Development of Super-Salt-Tolerant Azuki Bean (V. angularis) Cultivars.

Wild relatives of crops have the potential to improve food crops, especially in terms of improving abiotic stress tolerance. Two closely related wild species of the traditional East Asian legume crops, Azuki bean (Vigna angularis), V. riukiuensis "Tojinbaka" and V. nakashimae "Ukushima" were shown to have much higher levels of salt tolerance than azuki beans. To identify the genomic regions responsible for salt tolerance in "Tojinbaka" and "Ukushima", three interspecific hybrids were developed: (A) azuki bean cultivar "Kyoto Dainagon" × "Tojinbaka", (B) "Kyoto Dainagon" × "Ukushima" and (C) "Ukushima" × "Tojinbaka". Linkage maps were developed using SSR or restriction-site-associated DNA markers. There were three QTLs for "percentage of wilt leaves" in populations A, B and C, while populations A and B had three QTLs and population C had two QTLs for "days to wilt". In population C, four QTLs were detected for Na+ concentration in the primary leaf. Among the F2 individuals in population C, 24% showed higher salt tolerance than both wild parents, suggesting that the salt tolerance of azuki beans can be further improved by combining the QTL alleles of the two wild relatives. The marker information would facilitate the transfer of salt tolerance alleles from "Tojinbaka" and "Ukushima" to azuki beans.

Complete chloroplast and mitochondrial genomes of Ditrichum rhynchostegium Kindb. (Ditrichaceae, Bryophyta).

The moss family Pottiaceae is one of the most diverse lineages of the subclass Dicranidae (haplolepideous mosses). Nevertheless, the phylogenetic relationships of Pottiaceae with other Dicranidae families remain unclear. To better understand the ancestral genomic structure and evolution of the Pottiaceae, herein, we present the chloroplast and mitochondrial genomes of Ditrichum rhynchostegium (Ditrichaceae, Bryophyta). The chloroplast genome is 124,628 bp long and displayed a circular structure composed of a large single-copy region, a small single-copy region, and a pair of inverted repeats. It has 118 genes, including 82 protein-coding genes, 32 tRNA genes, and four rRNA genes. The mitochondrial genome is 106,246 bp long and has a circular structure. It contains 67 genes, including 40 protein-coding genes, 24 tRNA genes, and three rRNA genes. Phylogenetic trees based on the coding sequences strongly support the sister relationship of D. rhynchostegium with all Pottiaceous accessions, and the dextrosely arranged operculum cells suggest its affinity for Pottiaceae. This study also demonstrates that long-read sequencing employing the Nanopore platform facilitates the repair of unassembled or misassembled organellar genomic regions.

High source-sink ratio at and after sink capacity formation promotes green stem disorder in soybean.

Green stem disorder (GSD) of soybean is characterized by delayed leaf and stem maturation despite normal pod maturation. Previous studies have suggested that GSD occurrence is promoted by a high source-sink ratio, which is produced by thinning or shade removal at the R5 growth stage (the beginning of seed filling). Here the effects of different times and durations of shade removal after the R5 stage on GSD severity were analyzed. First, shade removal for more than 28 days after R5 increased GSD severity by more than 0.4 point in GSD score. Thinning treatment at R5 increased specific leaf weight by 23%, suppressed stem dry weight reduction, and upregulated 19 genes including those encoding vegetative storage proteins at R5 + 28d, indicating excess source ability relative to sink size. On the contrary, shade removal for 14 days after R5 decreased GSD severity by 0.5 point in GSD score. In this treatment, seed size was smaller, while seed number was significantly larger than control, suggesting that shortage of source ability relative to sink size. These results implied that soybean plants regulate GSD occurrences either positively or negatively according to a source-sink ratio during the R5 to R5 + 28d growth stages.

De Novo Genome Assembly of Stinkhorn Mushroom Clathrus columnatus (Basidiomycota, Fungi) Using Illumina and Nanopore Sequencing Data.

We report the reference genome of Clathrus columnatus isolate MO-923, which was isolated from Chichijima Island, the Ogasawara (Bonin) Islands, Japan. Oxford Nanopore Technologies MinION and Illumina sequence reads were assembled using NECAT and polished using Pilon to yield a 36.51-Mb genome with 10,625 predicted protein-coding genes.

Targeted amplicon sequencing + next-generation sequencing-based bulked segregant analysis identified genetic loci associated with preharvest sprouting tolerance in common buckwheat (Fagopyrum esculentum).

BACKGROUND: Common buckwheat (2n = 2x = 16) is an outcrossing pseudocereal whose seeds contain abundant nutrients and potential antioxidants. As these beneficial compounds are damaged by preharvest sprouting (PHS) and PHS is likely to increase with global warming, it is important to find efficient ways to develop new PHS-tolerant lines. However, genetic loci and selection markers associated with PHS in buckwheat have not been reported. RESULTS: By next-generation sequencing (NGS) of whole-genome of parental lines, we developed a genome-wide set of 300 markers. By NGS- based bulked segregant analysis (NGS-BSA), we developed 100 markers linked to PHS tolerance. To confirm the effectiveness of marker development from NGS-BSA data, we developed 100 markers linked to the self-compatibility (SC) trait from previous NGS-BSA data. Using these markers, we developed genetic maps with AmpliSeq technology, which can quickly detect polymorphisms by amplicon-based multiplex targeted NGS, and performed quantitative trait locus (QTL) analysis for PHS tolerance in combination with NGS-BSA. QTL analysis detected two major and two minor QTLs for PHS tolerance in a segregating population developed from a cross between the PHS-tolerant 'Kyukei 29' and the self-compatible susceptible 'Kyukei SC7'. We found different major and minor QTLs in other segregating populations developed from the PHS-tolerant lines 'Kyukei 28' and 'NARO-FE-1'. Candidate markers linked to PHS developed by NGS-BSA were located near these QTL regions. We also investigated the effectiveness of markers linked to these QTLs for selection of PHS-tolerant lines among other segregating populations. CONCLUSIONS: We efficiently developed genetic maps using a method combined with AmpliSeq technology and NGS-BSA, and detected QTLs associated with preharvest sprouting tolerance in common buckwheat. This is the first report to identify QTLs for PHS tolerance in buckwheat. Our marker development system will accelerate genetic research and breeding in common buckwheat.

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections.

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.

Characterization of the genomic sequence data around common cutworm resistance genes in soybean (Glycine max) using short- and long-read sequencing methods.

The common cutworm (CCW, Spodopteraab litura Fabricius) is one of the pests that most severely infect soybean (Glycine max L. Merr.). In a previous report, quantitative trait loci (QTL) analysis of CCW resistance using a recombinant inbred line derived from a cross between a susceptible cultivar 'Fukuyutaka' and a resistant cultivar 'Himeshirazu', identified two antixenosis resistance QTLs, CCW-1 and CCW-2. To reveal sequence variation between the aforementioned two cultivars, whole genome resequencing was performed using Illumina HiSeq2000 (75,632,747 and 91,540,849 reads). The generated datasets can be used for fine mapping and gene isolation of CCW-1 and CCW-2 as well as for revealing more detailed genetic differences between 'Fukuyutaka' and 'Himeshirazu' .

IonBreeders: bioinformatics plugins toward genomics-assisted breeding.

Polymorphism information generated by next-generation sequencing (NGS) technologies has enabled applications of genome-wide markers assisted breeding. However, handling such large-scale data remains a challenge for experimental researchers and breeders, calling for the urgent development of a flexible and straightforward analysis tool for NGS data. We developed "IonBreeders" as bioinformatics plugins that implement general analysis steps from genotyping to genomic prediction. IonBreeders comprises three plugins, "ABH", "IMPUTATION", and "GENOMIC PREDICTION", for format conversion of genotyping data, preprocessing and imputation of genotyping data, and genomic prediction, respectively. "ABH" converts genotyping data derived from NGS into the ABH format, which is acceptable for our further plugins and with other breeding software tools, R/qtl, MapMaker, and AntMap. "IMPUTATION" filters out non-informative markers and imputes missing marker genotypes. In "GENOMIC PREDICTION", users can use four statistical methods based on their target trait, quantitative trait locus effect, and number of markers, and construct a prediction model for genomic selection. IonBreeders is operated in Torrent Suite, but can also handle genotype data in standard formats, e.g., Variant Call Format (VCF), by format conversion using free software or our provided scripts.

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.

Detection of heading date QTLs in advanced-backcross populations of an elite indica rice cultivar, IR64.

IR64 is one of the world's most popular rice cultivars. To collect genetic factors involved in controlling its heading date, we developed 70 reciprocal advanced-backcross populations with a total of 6284 individuals at the BC4F2 generation from crosses between Koshihikari and IR64. We detected 29 QTLs associated with heading date on chromosomes 3, 5-8, 10, and 12. Twenty QTLs were located in the same chromosome regions as previously isolated heading date genes (Hd1, Hd6, Hd16, Ghd7, DTH8, Hd17, and Hd18). The rest were located in other chromosome regions. We found more number of QTLs than previous studies using mapping populations of IR64. Fine mapping in additional advanced-backcross populations clearly revealed that QTLs on the long arm of chromosome 7 are overlapping and seem to be a novel genetic factor for heading date because of their different locations from OsPRR37. Our results suggest that the difference in heading date between IR64 and Koshihikari is genetically controlled by many factors, and that a non-functional allele of Hd1 contributes to early heading of IR64 in the genetic background of functional alleles of other heading date QTLs and genes such as Hd6 and Hd16.

Highly multiplexed AmpliSeq technology identifies novel variation of flowering time-related genes in soybean (Glycine max).

Whole-genome re-sequencing is a powerful approach to detect gene variants, but it is expensive to analyse only the target genes. To circumvent this problem, we attempted to detect novel variants of flowering time-related genes and their homologues in soybean mini-core collection by target re-sequencing using AmpliSeq technology. The average depth of 382 amplicons targeting 29 genes was 1,237 with 99.85% of the sequence data mapped to the reference genome. Totally, 461 variants were detected, of which 150 sites were novel and not registered in dbSNP. Known and novel variants were detected in the classical maturity loci-E1, E2, E3, and E4. Additionally, large indel alleles, E1-nl and E3-tr, were successfully identified. Novel loss-of-function and missense variants were found in FT2a, MADS-box, WDR61, phytochromes, and two-component response regulators. The multiple regression analysis showed that four genes-E2, E3, Dt1, and two-component response regulator-can explain 51.1-52.3% of the variation in flowering time of the mini-core collection. Among them, the two-component response regulator with a premature stop codon is a novel gene that has not been reported as a soybean flowering time-related gene. These data suggest that the AmpliSeq technology is a powerful tool to identify novel alleles.

Development of Genome-Wide SNP Markers for Barley via Reference- Based RNA-Seq Analysis.

Marker-assisted selection of crop plants requires DNA markers that can distinguish between the closely related strains often used in breeding. The availability of reference genome sequence facilitates the generation of markers, by elucidating the genomic positions of new markers as well as of their neighboring sequences. In 2017, a high quality genome sequence was released for the six-row barley (Hordeum vulgare) cultivar Morex. Here, we developed a de novo RNA-Seq-based genotyping procedure for barley strains used in Japanese breeding programs. Using RNA samples from the seedling shoot, seedling root, and immature flower spike, we mapped next-generation sequencing reads onto the transcribed regions, which correspond to ∼590 Mb of the whole ∼4.8-Gbp reference genome sequence. Using 150 samples from 108 strains, we detected 181,567 SNPs and 45,135 indels located in the 28,939 transcribed regions distributed throughout the Morex genome. We evaluated the quality of this polymorphism detection approach by analyzing 387 RNA-Seq-derived SNPs using amplicon sequencing. More than 85% of the RNA-Seq SNPs were validated using the highly redundant reads from the amplicon sequencing, although half of the indels and multiple-allele loci showed different polymorphisms between the platforms. These results demonstrated that our RNA-Seq-based de novo polymorphism detection system generates genome-wide markers, even in the closely related barley genotypes used in breeding programs.

Soybean leaves transcriptomic data dissects the phenylpropanoid pathway genes as a defence response against Phakopsora pachyrhizi.

Asian soybean rust (ASR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is responsible for severe yield losses of up to 90% in all soybean producing countries. Till today, eight resistance to Phakopsora pachyrhizi (Rpp) loci have been mapped in soybean. Their resistance mechanism is race specific but largely unknown. The transcriptomes of susceptible BRS184 and Rpp3 with ASR isolates T1-2 at 24 h after inoculation (hai) and without ASR inoculation (mock) were annotated by similarity searching with different databases. A total of 4518 differentially expressed genes were identified. We found 70.89%, 56.61%, 32.13%, and 56.04% genes in the protein family databases (PFAM), Gene Ontology (GO), Eukaryotic clusters of Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG), respectively. KEGG disclosed that 52% of the phenylpropanoid pathway related genes were up-regulated. The relative gene expression study for selected genes of that pathway was conducted by RT-qPCR using Rpp1-Rpp4 carrying lines with T1-2 infection. The RT-qPCR results revealed that the Rpp lines utilized these genes in a rate limiting manner as a defence response. With the exception of glycinol 4-dimethylallyltransferase (G4DT) and chalcone reductase (CHR), all the genes showed the greatest expression at 12 hai, but the gene expressions which occur between 24 and 96 hai make these Rpp lines unique to their respective ASR isolates. Moreover, functional coordination of arogenate dehydratase 6 (ADT6) and 4-hydroxy-3-methylbut-2-enyl diphosphate synthase (ispG), chalcone synthase (CHS) and CHR, and G4DT and phytyltransferase 3 (PT3) may have a great impact on soybean resistance against ASR.