Published papers

Understanding how genetic variation within a foundation species determines the structure of associated communities and ecosystem processes has been an emerging frontier in ecology. Previous studies in common gardens identified close links between intraspecific variation and multispecies community structure, and these findings are now being evaluated directly in the complex natural ecosystem. In this study, we examined to what extent genomic variation in a foundation tree species explains the structure of associated arthropod communities in the field, comparing with spatial, temporal and environmental factors. In a continuous mixed forest, arthropods were surveyed on 85 mature alders (Alnus hirsuta) in 2 years. Moreover, we estimated Nei's genetic distance among the alders based on 1,077 single nucleotide polymorphisms obtained from restricted-site-associated DNA sequencing of the alders’ genome. In both years, we detected significant correlations between genetic distance and dissimilarity of arthropod communities. A generalized dissimilarity modelling indicated that the genetic distance of alder populations was the most important predictor to explain the variance of arthropod communities. Among arthropod functional groups, carnivores were consistently correlated with genetic distance of the foundation species in both years. Furthermore, the extent of year-to-year changes in arthropod communities was more similar between more genetically closed alder populations. This study demonstrates that the genetic similarity rule would be primarily prominent in community assembly of plant-associated arthropods under temporally and spatially variable environments in the field.

Japanese pepper, Zanthoxylum piperitum, is native to Japan and has four well-known lineages (Asakura, Takahara, Budou, and Arima), which are named after their production area or morphology. Restriction-site associated DNA sequencing (RAD-Seq) was used to analyse 93 accessions from various areas, including these four lineages. Single nucleotide variant analysis was used to classify the plants into eight groups: the Asakura and Arima lineages each had two groups, the Takahara and Budou lineages each had one group, and two additional groups were present. In one Asakura group and two Arima groups, the plants were present in agricultural fields and mountains, thus representing the early stage of domestication of the Japanese pepper. The second Asakura lineage group was closely related to plants present in various areas, and this represents the second stage of domestication of this plant because, after early domestication, genetically related lineages with desirable traits spread to the periphery. These results demonstrate that domestication of Japanese pepper is ongoing. In addition, this study shows that spineless plants are polyphyletic, despite the spineless lineage being considered a subspecies of Japanese pepper.

Brassicaceae plants have a dual-cell type of chemical defense against herbivory. Here, we show a novel single-cell defense involving endoplasmic reticulum (ER)-derived organelles (ER bodies) and the vacuoles. We identify various glucosinolates as endogenous substrates of the ER-body β-glucosidases BGLU23 and BGLU21. Woodlice strongly prefer to eat seedlings of bglu23 bglu21 or a glucosinolate-deficient mutant over wild-type seedlings, confirming that the β-glucosidases have a role in chemical defense: production of toxic compounds upon organellar damage. Deficiency of the Brassicaceae-specific protein NAI2 prevents ER-body formation, which results in a loss of BGLU23 and a loss of resistance to woodlice. Hence, NAI2 that interacts with BGLU23 is essential for sequestering BGLU23 in ER bodies and preventing its degradation. Artificial expression of NAI2 and BGLU23 in non-Brassicaceae plants results in the formation of ER bodies, indicating that acquisition of NAI2 by Brassicaceae plants is a key step in developing their single-cell defense system.

The phase of the circadian clock (body time) is related to various aspects of physiology and metabolism in plants. Therefore, it is important to evaluate body time to control plant quality and growth. The molecular timetable method, a statistical analysis of the transcriptome, is useful for plant body time estimation. However, it is necessary to verify some assumptions in this method for high resolution estimation in plants under various conditions. To investigate whether the assumptions are valid, we analyzed the period length and expression pattern of genes. We used Arabidopsis thaliana as a model plant, and both lettuce and tomato as model crops. A. thaliana and lettuce were cultivated under conditions of constant light in a closed plant factory, and tomato was cultivated in a sunlight-type plant factory. We found that estimation of the period length of the plant circadian rhythm was possible and that the gene expression pattern correlates well with test-fit cosine curves. Thus, we successfully demonstrated high-resolution estimation of body time in plants using a molecular timetable method.

Sympatric speciation has been demonstrated in few empirical case studies, despite intense searches, because of difficulties in testing the criteria for this mode of speciation. Here, we report a possible case of sympatric speciation in ricefishes of the genus Oryzias on Sulawesi, an island of Wallacea. Three species of Oryzias are known to be endemic to Lake Poso, an ancient tectonic lake in central Sulawesi. Phylogenetic analyses using RAD-seq-derived single nucleotide polymorphisms (SNPs) revealed that these species are monophyletic. We also found that the three species are morphologically distinguishable and clearly separated by population-structure analyses based on the SNPs, suggesting that they are reproductively isolated from each other. A mitochondrial DNA chronogram suggested that their speciation events occurred after formation of the tectonic lake, and existence of a historical allopatric phase was not supported by coalescent-based demographic inference. Demographic inference also suggested introgressive hybridization from an outgroup population. However, differential admixture among the sympatric species was not supported by any statistical tests. These results all concur with criteria necessary to demonstrate sympatric speciation. Ricefishes in this Wallacean lake provide a promising new model system for the study of sympatric speciation.

Geographic and environmental isolations of islands and the mainland offer excellent opportunity to investigate colonization and survival dynamics of island populations. We inferred and compared evolutionary processes and the demographic history of Rhododendron tsusiophyllum, in the Izu Islands and the much larger island Honshu, treated here as the mainland, using thousands of nuclear SNPs obtained by ddRAD-seq from eight populations of R. tsusiophyllum and three populations of R. tschonoskii as an outgroup. Phylogenetic relationships and their habitats suggest that R. tsusiophyllum had evolved and migrated from cold north to warm south regions. We detected clear genetic divergence among populations in three regions of Honshu and the Izu Islands, suggesting restricted migration between them due to isolated habitats on mountains even in the mainland. The three regions have different changes in effective population size, especially, genetic diversity and population size of the Izu Islands are small compared to the others. Further, habitats of populations in the Izu Islands are warmer than those in Honshu, suggesting that they have undergone adaptive evolution. Our study provides evidences of montane rather than insular isolation on genetic divergence, survival of populations and significance of adaptive evolution for island populations with small population size and low genetic diversity, despite close proximity to mainland populations.

The plant root-associated environments such as the rhizosphere, rhizoplane, and endosphere are different from the outer soil region (bulk soil). They establish characteristic conditions including microbiota, metabolites, and minerals, and they can directly affect plant growth and development. However, comprehensive insights into those characteristic environments, especially the rhizosphere, and molecular mechanisms of their formation are not well understood. In the present study, we investigated the spatiotemporal dynamics of the root-associated environment in actual field conditions by multi-omics analyses (mineral, microbiome, and transcriptome) of soybean plants. Mineral and microbiome analyses demonstrated a characteristic rhizosphere environment in which most of the minerals were highly accumulated and bacterial communities were distinct from those in the bulk soil. Mantel’s test and co-abundance network analysis revealed that characteristic community structures and dominant bacterial taxa in the rhizosphere significantly interact with mineral contents in the rhizosphere, but not in the bulk soil. Our field multi-omics analysis suggests a rhizosphere-specific close association between the microbiota and mineral environment.

Hydathodes are typically found at leaf teeth in vascular plants and are involved in water release to the outside. Although morphological and physiological analysis of hydathodes has been performed in various plants, little is known about the genes involved in hydathode function. In this study, we performed fluorescent protein-based imaging and tissue-specific RNA-seq analysis in Arabidopsis hydathodes. We used the enhancer trap line E325, which has been reported to express green fluorescent protein (GFP) at its hydathodes. We found that E325-GFP was expressed in small cells found inside the hydathodes (named E cells) that were distributed between the water pores and xylem ends. No fluorescence of the phloem markers pSUC2:GFP and pSEOR1:SEOR1-YFP was observed in the hydathodes. These observations indicate that Arabidopsis hydathodes are composed of three major components: Water pores, xylem ends, and E cells. In addition, we performed transcriptome analysis of the hydathode using the E325-GFP line. Microsamples were collected from GFP-positive or-negative regions of E325 leaf margins with a needle-based device (~130 μm in diameter). RNA-seq was performed with each single microsample using a high-throughput library preparation method called Lasy-Seq. We identified 72 differentially expressed genes. Among them, 68 genes showed significantly higher and four genes showed significantly lower expression in the hydathode. Our results provide new insights into the molecular basis for hydathode physiology and development.

A recent phylogenetic study based only on chloroplast DNA (cpDNA) variation revealed that populations of an Isodon species are frequently embedded paraphyletically among other Isodon species. This phylogenetic discrepancy between species taxonomy and molecular phylogeny was considered to have resulted from chloroplast DNA captures and/or incomplete lineage sorting. To elucidate which of these factors was mainly responsible for the observed phylogenetic pattern, we performed phylogenetic analyses of multiple populations of Isodon species in Japan using cpDNA variation, three single-copy nuclear genes, and double-digest restriction-site-associated DNA sequencing (ddRAD-seq). Although a species often shared chlorotypes with other species, our phylogenetical analyses based on variation in the three single-copy nuclear genes and the ddRAD-seq data showed that most populations belonging to the same species were monophyletic at the species level, suggesting that chloroplast capture may have frequently occurred between Isodon species. Some populations of an intraspecific taxon were embedded paraphyletically within the species, regardless of the large amount of phylogenetic information in nuclear DNA; this incongruity may have resulted from incomplete lineage sorting.

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was 10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/ KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Increase in the nitrogen (N)-use efficiency and optimization of N response in crop species are urgently needed. Although transcription factor-based genetic engineering is a promising approach for achieving these goals, transcription factors that play key roles in the response to N deficiency have not been studied extensively. Here, we performed RNA-seq analysis of root samples of 20 Asian rice (Oryza sativa) accessions with differential nutrient uptake. Data obtained from plants exposed to N-replete and N-deficient conditions were subjected to coexpression analysis and machine learning-based pathway inference to dissect the gene regulatory network required for the response to N deficiency. Four transcription factors, including members of the G2-like and bZIP families, were predicted to function as key regulators of gene transcription within the network in response to N deficiency. Cotransfection assays validated inferred novel regulatory pathways, and further analyses using genome-edited knockout lines suggested that these transcription factors are important for N-deficiency responses in planta. Many of the N deficiency-responsive genes, including those encoding key regulators within the network, were coordinately regulated by transcription factors belonging to different families. Transcription factors identified in this study could be valuable for the modification of N response and metabolism.

To study the genetic structure of clonal plant populations, genotyping and genet detection using genetic markers are necessary to assign ramets to corresponding genets. Assignment is difficult as it involves setting a robust threshold of genetic distance for genet distinction as neighbouring genets in a plant population are often genetically related. Here, we used restriction site-associated DNA sequencing (RAD-seq) for a rhizomatous clonal herb, Cardamine leucantha [Brassicaceae] to accurately determine genet structure in a natural population. We determined a draft genome sequence of this species for the first time, which resulted in 66 617 scaffolds with N50 = 6086 bp and an estimated genome size of approximately 253 Mbp. Using genetic distances based on the RAD-seq analysis, we successfully distinguished ramets that belonged to distinct genets even from a half-sib family. We applied these methods to 372 samples of C. leucantha collected at 1-m interval grids within a 20 ×20 m plot in a natural population in Hokkaido, Japan. From these samples, we identified 61 genets with high inequality in terms of genet size and patchy distribution. Spatial autocorrelation analyses indicated significant aggregation within 7 and 4 m at ramet and genet levels, respectively. An analysis of parallel DNA microsatellite loci (simple sequence repeats) suggested that RAD-seq can provide data that allows robust genet assignment. It remains unclear whether the large genets identified here became dominant stochastically or deterministically. Precise identification of genets will assist further study and characterization of dominant genets.

The genetic bases of growth and body weight are of economic and scientific interest, and teleost fish models have proven useful in such investigations. The Oryzias latipes species complex (medaka) is an abundant freshwater fish in Japan and suitable for genetic studies. We compared two wild medaka stocks originating from different latitudes. The Maizuru population from higher latitudes weighed more than the Ginoza population. We investigated the genetic basis of body weight, using quantitative trait locus (QTL) analysis of the F2 offspring of these populations. We detected one statistically significant QTL for body weight on medaka chromosome 4 and identified 12 candidate genes that might be associated with body weight or growth. Nine of these 12 genes had at least one single nucleotide polymorphism that caused amino acid substitutions in protein-coding regions, and we estimated the effects of these substitutions. The present findings might contribute to the marker-assisted selection of economically important aquaculture species.

Taishu horses are a native Japanese breed, of which only 41 individuals remained on Tsushima Island in 2018. Their genetic diversity is considered lower than that of other Japanese native horse breeds; thus, it needs to be investigated for sustainable conservation of this breed. Historical records revealed that several Taishu individuals were released areas off-Tsushima Island in mid-1980s. At present, Taishu horses living outside of Tsushima Island, hereafter referred to as Non-Tsushima Taishus (NTTs), are tagged. However, the genetic structure of the NTT individuals remains unclear, and such individuals are not included in the current mating plans for Taishu horses. Herein, we examined the genetic structure of 18 NTT individuals by comparing their genomic (SNP) information with that of individuals on Tsushima Island (TT), four other native Japanese breeds, and one Anglo- Arabian breed by using ddRAD-seq. We found that all individuals related to the Taishu can be grouped in one cluster, which was separated from other horse breeds. Patterns of specific and shared SNPs in NTT individuals closely resembled those of TT individuals, suggesting very minor genetic differences. Meanwhile, the heterozygosity of NTT individuals was slightly higher than that of TT individuals, and many NTT individuals were of fertile age, suggesting that the pedigree of NTT individuals would be useful in breed conservation plans for Taishu horses. Based on their genomic information, we also reconstructed the pedigree structures of four NTT individuals with no family information. The inclusion of NTT individuals in future mating plans on Tsushima Island may be an effective and feasible method for conserving the Taishu horse breed in Japan.

The genetic architecture of innate fear behavior in chickens is poorly understood. Here, we performed quantitative trait loci (QTL) analysis of innate responses to tonic immobility (TI) and open field (OF) fears in 242 newly hatched chicks of an F2 population between the native Japanese Nagoya breed and the White Leghorn breed using 881 single nucleotide polymorphism markers obtained by restriction site-associated DNA sequencing. At genome-wide 5% significance levels, four QTL for TI traits were revealed on chromosomes 1–3 and 24. Two of these loci had sex-specific effects on the traits. For OF traits, three QTL were revealed on chromosomes 2, 4 and 7. The TI and OF QTL identified showed no overlaps in genomic regions and different modes of inheritance. The three TI QTL and one OF QTL exerted antagonistic effects on the traits. The results demonstrated that context-dependent QTL underlie the variations in innate TI and OF behaviors.

Populations or species exploiting different habitats can differ in sensory perception as a result of divergent adaptation. In bony fish, the water current is perceived via neuromasts, the end organ of the lateral line system. Although fish in different habitats are known to vary in neuromasts, we know little about the genetic basis for such variation. Here, we investigate the genetic basis for variation in supraorbital neuromasts in a hybrid zone between the Japanese gobies Gymnogobius breunigii and Gymnogobius castaneus. The former has supraorbital canal neuromasts with six cephalic pores, whereas the latter has only superficial neuromasts with no canals or pores in the supraorbital region. Our genomic analysis showed that G. breunigii and G. castaneus occur mainly in the lower and mid/upper reaches, respectively. In a river in northern Japan, hybrids were found at the sites between the habitats of the two species. These hybrids exhibited anomalies of cephalic pores. Using this hybrid zone, we conducted genome-wide association studies and identified one locus significantly associated with the number of pores. Genomic cline analysis in the hybrid zone demonstrated that this locus exhibited a higher introgression rate compared with the genomic background, indicating the possibility of adaptive introgression.

Background and Aims: Contrasting life-history traits can evolve through generations of dwarf plant ecotypes, yet such phenotypic changes often involve decreased plant size and reproductive allocation, which can configure seed dispersal patterns and, subsequently, population demography. Therefore, evolutionary transitions to dwarfism can represent good study systems to test the roles of life-history traits in population demography by comparing genetic structure between related but phenotypically divergent ecotypes. Methods: In this study, we examined an ecotypic taxon pair of the world's smallest goldenrod (stem height 2.6 cm) in alpine habitats and its closely related lowland taxon (30-40 cm) found on Yakushima Island, Japan. Genetic variation in chloroplast DNA sequences, nuclear microsatellites and genome-wide single-nucleotide polymorphisms were used to investigate 197 samples from 16 populations, to infer the population genetic demography and compare local genetic structure of the ecotypes. Key Results: We found a pronounced level of genetic differentiation among alpine dwarf populations, which were much less geographically isolated than their lowland counterparts. In particular, several neighbouring dwarf populations (located ~500 m apart) harboured completely different sets of chloroplast haplotypes and nuclear genetic clusters. Demographic modelling revealed that the dwarf populations have not exchanged genes at significant levels after population divergence. Conclusions: These lines of evidence suggest that substantial effects of genetic drift have operated on these dwarf populations. The low-growing stature and reduced fecundity (only 3.1 heads per plant) of the dwarf plants may have reduced gene flow and rare long-distance seed dispersal among habitat patches, although the effects of life-history traits require further evaluation using ecological approaches.

Population genetic studies of marine species sometimes find no clear patterns of genetic differentiation, but signals of geographic gene flow typified by isolation by distance (IBD) due to high dispersal ability. A previous study of mitochondrial DNA of the yellowfin goby Acanthogobius flavimanus indicated that this species also does not have clearly differentiated population groups around the Japanese archipelago, in contrast with other coastal goby species that have highly differentiated groups between the Sea of Japan and Pacific Ocean populations. Here, to clearly reveal phylogeographic structures of the yellowfin goby around the Japanese archipelago, we conducted large-scale RAD-sequencing analysis of their nuclear genomic SNP loci. As a result, IBD patterns were confirmed across the yellowfin goby populations, and the larval dispersal distance of this species was revealed to be approximately 19 km on average. In addition, different larval dispersal distances were estimated for its Sea of Japan and Pacific Ocean populations, suggesting that different coastal environments and/or ecological features may affect its dispersal rates. These results would provide foundations for ecological and conservation studies of the yellowfin goby and coastal fishes having similar dispersal strategy to this species.

Transgressive segregation produces hybrid progeny phenotypes that exceed the parental phenotypes. Unlike heterosis, extreme phenotypes caused by transgressive segregation are heritably stable. We examined transgressive phenotypes of flowering time in rice, and revealed transgressive segregation in F2 populations derived from a cross between parents with similar (proximal) days to heading (DTH). The DTH phenotypes of the A58 · Kitaake F2 progenies were frequently more extreme than those of either parent. These transgressive phenotypes were maintained in the F3 and F4 populations. Both A58 and Kitaake are japonica rice cultivars adapted to Hokkaido, Japan, which is a high-latitude region, and have a short DTH. Among the four known loci required for a short DTH, three loci had common alleles in A58 and Kitaake, implying there is a similar genetic basis for DTH between the two varieties. A genome-wide single nucleotide polymorphism (SNP) analysis based on the F4 population identified five new quantitative trait loci (QTL) associated with transgressive DTH phenotypes. Each of these QTL had different degrees of additive effects on DTH, and two QTL had an epistatic effect on each other. Thus, a genome-wide SNP analysis facilitated the detection of genetic loci associated with extreme DTH phenotypes, and revealed that the transgressive phenotypes were produced by exchanging the complementary alleles of a few minor QTL in the similar parental phenotypes.

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.