Search

Postmating isolation is thought to be an important driver of the late stages of speciation. However, relatively little is empirically known about the process compared with other isolating mechanisms that drive the early stages of speciation, especially in non-model organisms. We characterized the genetic architecture of postmating isolation between 2 rockfishes, Sebastes schlegelii and S. trivittatus, whose reproductive isolation is complete. We examined transmission ratio distortion (TRD) patterns of genetic markers in 2 reciprocal backcross populations. Markers showing either of the 2 types of TRD was widespread across the genome, with some of the distorted markers forming extensive clusters around the recombination coldspots. These suggest that the postmating isolation effectively prevents gene flow across the genome and the recombination landscape contributes to the genetic architecture. Comparisons between 2 backcross families and 2 developmental stages showed little similarity in the distorted markers, suggesting asymmetry and stage specificity of the isolation. This may be due to hybrid incompatibility involving maternal factors or extrinsic selection. The lack of sex-ratio distortion in the mapping families suggested that Haldane’s rule in terms of hybrid inviability does not hold. Additionally, quantitative trait loci (QTL) mapping detected significant QTLs for sex and the morphological traits relevant to speciation and convergence of rockfishes, including body coloration. Genes in the melanocortin system, including agouti-signaling protein 1 (asip1) and melanocortin 1 receptor (mc1r), might underlie the horizontal and vertical color patterns on the body, respectively. These findings constitute an essential step toward a comprehensive understanding of speciation and morphological diversification of rockfishes.

Background and Aims Plant propagules often possess specialized morphologies that facilitate dispersal across specific landscapes. In the fruit dimorphism of a coastal shrub, Scaevola taccada, individual plants produce either cork-morph or pulp-morph fruits. The former is buoyant and common on sandy beaches, whereas the latter does not float, is bird-dispersed, and is common on elevated sites such as slopes on sea cliffs and behind rocky shores. We hypothesized that beach populations bridge the heterogeneous landscapes by serving as a source of both fruit types, while dispersal is biased for the pulp morph on elevated sites within the islands and for the cork morph between beaches of different islands. Based on this hypothesis, we predicted that populations in elevated sites would diverge genetically over time due to isolation by distance, whereas beach populations would maintain high genetic similarity via current gene flow. Methods The genetic structure and gene flow in S. taccada were evaluated by investigating genome-wide single nucleotide polymorphisms in plants from 17 sampling sites on six islands (belonging to the Ryukyu, Daito and Ogasawara Islands) in Japan. Key Results Geographical isolation was detected among the three distant island groups. Analyses within the Ryukyu Islands suggested that sandy beach populations were characterized by genetic admixture, whereas populations in elevated sites were relatively isolated between the islands. Pairwise FST values between islands were lowest between sandy beaches, intermediate between sandy beaches and elevated sites, and highest between elevated sites. Conclusions Dispersal across the ocean by cork morphs is sufficiently frequent to prevent genetic divergence between beaches of different islands. Stronger genetic isolation of elevated sites between islands suggests that bird dispersal by pulp morphs is restricted mainly within islands. These contrasting patterns of gene flow realized by fruit dimorphism provide evidence that fruit characteristics can strongly mediate genetic structure.

The glucosinolate–myrosinase defense system (GMDS), characteristic of Brassicales, is involved in plant defense. Previous single-cell transcriptomic analyses have reported the expression profiles of multiple GMDS-related cell types (i.e. myrosinase-rich myrosin idioblasts and multiple types of potential glucosinolate synthetic cells as well as a candidate S-cell for glucosinolate accumulation). However, differences in plant stages and cell-type annotation methods have hindered comparisons among studies. Here, we used the single-cell transcriptome profiles of extended Arabidopsis leaves and verified the distribution of previously used markers to refine the expression profiles of GMDS-associated cell types. Moreover, we performed beta-glucuronidase promoter assays to confirm the histological expression patterns of newly obtained markers for GMDS-associated candidates. As a result, we found a set of new specific reporters for myrosin cells and potential glucosinolate-producing cells.

The long-term dynamics of the transcriptome under natural field conditions remain unclear. We conducted comprehensive gene expression analyses of rice leaves and roots grown under natural field conditions for a long period, from the tillering stage to the ripening stage. In this experiment, changes in the transcriptome were captured in relation to microclimatic parameters, particularly potential evaporation (Ep), which is a multiple meteorological factor and acts as an indicator of transpirational demand. The results indicated that many genes were regulated by changes in temperature and Ep in both leaves and roots. Furthermore, the correlation between gene expression and meteorological factors differed significantly between the vegetative and reproductive stages. Since Ep triggers transpiration, we analyzed aquaporin gene expression, which is responsible for water transport, and found that many aquaporin genes in leaves were positively correlated with Ep throughout the growth period, whereas in roots, two plasma membrane intrinsic aquaporins, PIP2;4 and PIP2;5 were strongly correlated with Ep during reproductive growth. Other genes closely related to productivity, such as those involved in nutrient absorption and photosynthesis, exhibited different responses to meteorological factors at different growth stages. The stage-dependent shift in the microclimate response provides an important perspective on crop physiology in light of climate change.

Water scarcity is a serious agricultural problem causing significant losses to crop yield and product quality. The development of technologies to mitigate the damage caused by drought stress is essential for ensuring a sustainable food supply for the increasing global population. We herein report that the exogenous application of ethanol, an inexpensive and environmentally friendly chemical, significantly enhances drought tolerance in Arabidopsis thaliana, rice and wheat. The transcriptomic analyses of ethanol-treated plants revealed the upregulation of genes related to sucrose and starch metabolism, phenylpropanoids and glucosinolate biosynthesis, while metabolomic analysis showed an increased accumulation of sugars, glucosinolates and drought-tolerance-related amino acids. The phenotyping analysis indicated that drought-induced water loss was delayed in the ethanol-treated plants. Furthermore, ethanol treatment induced stomatal closure, resulting in decreased transpiration rate and increased leaf water contents under drought stress conditions. The ethanol treatment did not enhance drought tolerance in the mutant of ABI1, a negative regulator of abscisic acid (ABA) signaling in Arabidopsis, indicating that ABA signaling contributes to ethanol-mediated drought tolerance. The nuclear magnetic resonance analysis using 13C-labeled ethanol indicated that gluconeogenesis is involved in the accumulation of sugars. The ethanol treatment did not enhance the drought tolerance in the aldehyde dehydrogenase (aldh) triple mutant (aldh2b4/aldh2b7/aldh2c4). These results show that ABA signaling and acetic acid biosynthesis are involved in ethanol-mediated drought tolerance and that chemical priming through ethanol application regulates sugar accumulation and gluconeogenesis, leading to enhanced drought tolerance and sustained plant growth. These findings highlight a new survival strategy for increasing crop production under water-limited conditions.

Environmental DNA (eDNA) metabarcoding, a method that applies high-throughput sequencing and universal primer sets to eDNA analysis, has been a promising approach for efficient, comprehensive biodiversity monitoring. However, significant money-, labor-, and time-costs are still required for performing eDNA metabarcoding. In this study, we assessed the performance of an “early-pooling” protocol (a protocol based on 1st PCR tagging) to reduce the experimental costs of library preparation for eDNA metabarcoding. Specifically, we performed three experiments to investigate the effects of 1st PCR-tagging and 2nd PCR-indexing protocols on the community composition revealed by eDNA metabarcoding, the effects of post-1st PCR exonuclease purification on tag jumping (corresponds to index hopping in 2nd PCR indexing), and the effects of the number of PCR replicates and the eDNA template volume on the number of detected OTUs. Analyses of 204 eDNA libraries from three natural aquatic ecosystems and one mock eDNA sample showed that (i) 1st PCR tagging does not cause clear biases in the outcomes of eDNA metabarcoding, (ii) post-1st PCR exonuclease purification reduces the risk of tag jumping, and (iii) increasing the eDNA template volume may increase the number of detected OTUs and reduce variations in the detected community compositions, similar to increasing the number of 1st PCR replicates. Our results show that an early-pooling protocol with post-1st PCR exonuclease purification and an increased amount of the DNA template reduces the risk of tag jumping, the costs for consumables and reagents (except for many tagged 1st PCR primers), and the handling time in library preparation, and produces similar results to a 2nd PCR-indexing protocol. Therefore, once a target metabarcoding region is selected and a set of tagged-1st PCR primers is prepared, the early-pooling protocol provides a cost, labor, and time-efficient approach for processing a large number of samples.

SummaryGenetically diverse populations can increase plant resistance to natural enemies. Yet, beneficial genotype pairs remain elusive due to the occurrence of both positive and negative effects of mixed planting on plant resistance, called associational resistance and susceptibility. We used genome-wide polymorphisms of the plant speciesArabidopsis thalianato identify genotype pairs that enhance associational resistance to herbivory. By quantifying neighbor interactions among 199 genotypes grown in a randomized block design, we predicted that 823 of the 19,701 candidate pairs could reduce herbivory through associational resistance. We planted such pairs with predicted associational resistance in mixtures and monocultures and found a significant reduction in herbivore damage in the mixtures. Our study highlights the potential application to assemble genotype mixtures with positive biodiversity effects.

Population sizes of the Japanese red coral Corallium japonicum have been severely affected by poaching and overfishing. Although genetic structure and connectivity patterns are considered important parameters for conservation strategies, there are few studies focusing on the population genetics of C. japonicum in the Northwest Pacific. We examined the genetic population structure of C. japonicum, in the Northwest Pacific. We used restriction-site-associated DNA sequencing (RAD-seq), which can be used to identify genome-wide single-nucleotide polymorphism (SNPs), to reveal detailed within-species genetic variations. Using the variable SNP loci identified from this analysis, we successfully evaluated the population-level genetic diversity and patterns of gene flow among multiple populations of C. japonicum around Japan. The results of genetic analysis basically showed that gene flow is widely maintained in the geographic range examined in this study, but the analysis in combination with larval dispersal simulations revealed several populations that were genetically distinct from the other populations, suggesting geographically limited gene flows. The information obtained from this study will be useful for the design of effective management schemes for C. japonicum, which is under threat from overfishing.

The importance of hybridization in morphological diversification is a fundamental topic in evolutionary biology. However, despite accumulated knowledge concerning adult hybrid variations, how hybridization affects ontogenetic allometry is less well understood. Herein, we assessed the effects of hybridization on the postnatal ontogenetic allometry in skulls of a putative hybrid population of introduced Taiwanese macaques (Macaca cyclopis) and native Japanese macaques (M. fuscata). Genomic analyses indicated that the population consisted of individuals with various degrees of admixture, formed by male migration of Japanese to Taiwanese macaques. Hybridization parallelly shifted the ontogenetic trajectories of the overall skull shape in an almost additive manner. Transgressive variations in the overall skull shape were observed mainly in non-allometric components. Hybridization lifted the restriction found in Japanese macaques on the growth rate of the maxillary sinus (the hollow space in the face), producing hybrids with a mosaic pattern, i.e., the maxillary sinus is as large as that in Taiwanese macaques, although the overall skull shape is intermediate. Our findings suggest that the transgressive variations can be caused by prenatal shape modifications and the lifting of the genetic restriction on the regional growth rate, highlighting the complex genetic and ontogenetic bases underlying hybridization-induced morphological diversification.

The differences between plants grown in field and in controlled environments have long been recognized. However, few studies have addressed the underlying molecular mechanisms. To evaluate plant responses to fluctuating environments using laboratory equipment, we developed SmartGC, a high-performance growth chamber that reproduces the fluctuating irradiance, temperature and humidity of field environments. We analysed massive transcriptome data of rice plants grown under field and SmartGC conditions to clarify the differences in plant responses to field and controlled environments. Rice transcriptome dynamics in SmartGC mimicked those in the field, particularly during the morning and evening but those in conventional growth chamber conditions did not. Further analysis revealed that fluctuation of irradiance affects transcriptome dynamics in the morning and evening, while fluctuation of temperature affects transcriptome dynamics only in the morning. We found upregulation of genes related to biotic and abiotic stress, and their expression was affected by environmental factors that cannot be mimicked by SmartGC. Our results reveal fillable and unfillable gaps in the transcriptomes of rice grown in field and controlled environments and can accelerate the understanding of plant responses to field environments for both basic biology and agricultural applications.

Modes of reproduction in animals are diverse, with different modes having evolved independently in multiple lineages across a variety of taxa. However, an understanding of the genomic change driving the transition between different modes of reproduction is limited. Several ricefishes (Adrianichthyidae) on the island of Sulawesi have a unique mode of reproduction called "pelvic-fin brooding," wherein females carry externally fertilized eggs until hatching using their pelvic fins. Phylogenomic analysis demonstrated pelvic-fin brooders to have evolved at least twice in two distant clades of the Adrianichthyidae. We investigated the genetic architecture of the evolution of this unique mode of reproduction. Morphological analyses and laboratory observations revealed that females of pelvic-fin brooders have longer pelvic fins and a deeper abdominal concavity, and that they can carry an egg clutch for longer than nonbrooding adrianichthyids, suggesting that these traits play important roles in this reproductive mode. Quantitative trait locus mapping using a cross between a pelvic-fin brooder Oryzias eversi and a nonbrooding O. dopingdopingensis reveals different traits involved in pelvic-fin brooding to be controlled by different loci on different chromosomes. Genomic analyses of admixture detected no signatures of introgression between two lineages with pelvic-fin brooders, indicating that introgression is unlikely to be responsible for repeated evolution of pelvic-fin brooding. These findings suggest that multiple independent mutations may have contributed to the convergent evolution of this novel mode of reproduction.

Analyses of phylogeographic patterns and genetic diversity provide fundamental information for the management and conservation of species. However, little is published about these patterns in Japanese kelp species. In this study, we conducted phylogeographic analyses of a canopy-forming kelp, Eisenia bicyclis, based on genome-wide SNPs identified by ddRAD-seq. We obtained 1,299 SNPs for 76 samples from nine localities across the distribution. STRUCTURE, NeighborNet, and discriminant analysis of principal components consistently showed high genetic differentiation among the Eastern Pacific, Central Pacific, and Sea of Japan coastal regions. Relatively strong gene flow was detected only within populations in the Eastern Pacific and in the Sea of Japan. Genetic diversity and genetic uniqueness were high in the Central Pacific and low in the Sea of Japan. These results suggest that there were at least three independent refugia corresponding to the three regions during the Last Glacial Maximum (LGM). Furthermore, relatively larger populations in the Central Pacific and smaller populations in the Sea of Japan have been maintained in the demographic history from before the LGM to the present. These phylogeographic histories were supported by an Approximate Bayesian Computation analysis. From a conservation genetics perspective, the loss of southern populations in the Central Pacific would greatly reduce the total genetic diversity of the species. Southern populations in the Sea of Japan, which have relatively low genetic diversity, may be highly vulnerable to environmental change, such as heat waves and increased feeding. Therefore, careful monitoring and conservation are needed in the two regions.

Severe genetic redundancy is particularly clear in gene families encoding plant hormone receptors, each subtype sharing redundant and specific functions. Genetic redundancy of receptor family members represents a major challenge for the functional dissection of each receptor subtype. A paradigmatic example is the perception of the hormone (+)-7-iso-jasmonoyl-L-isoleucine, perceived by several COI1-JAZ complexes; the specific role of each receptor subtype still remains elusive. Subtype-selective agonists of the receptor are valuable tools for analyzing the responses regulated by individual receptor subtypes. We constructed a stereoisomer library consisting of all stereochemical isomers of coronatine (COR), a mimic of the plant hormone (+)-7-iso-jasmonoyl-L-isoleucine, to identify subtype-selective agonists for COI1-JAZ co-receptors in Arabidopsis thaliana and Solanum lycopersicum. An agonist selective for the Arabidopsis COI1-JAZ9 co-receptor efficiently revealed that JAZ9 is not involved in most of the gene downregulation caused by COR, and the degradation of JAZ9-induced defense without inhibiting growth.

Increasing the lodging resistance of rice through genetic improvement has been an important target in breeding. To further enhance the lodging resistance of high-yielding rice varieties amidst climate change, it is necessary to not only shorten culms but strengthen them as well. A landrace rice variety, Omachi, which was established more than 100 years ago, has the largest culm diameter and bending moment at breaking in the basal internodes among 135 temperate japonica accessions. Using unused alleles in such a landrace is an effective way to strengthen the culm. In this study, we performed quantitative trait locus (QTL) analysis to identify the genetic factors of culm strength of Omachi using recombinant inbred lines (RILs) derived from a cross between Omachi and Koshihikari, a standard variety in Japan. We identified three QTLs for the culm diameter of the 5th internode on chromosomes 3 (qCD3) and 7 (qCD7-1, qCD7-2). Among them, qCD7-2 was verified by QTL analysis using the F2 population derived from a cross between one of the RILs and Koshihikari. RNA-seq analysis of shoot apex raised 10 candidate genes underlying the region of qCD7-2. The increase in culm strength by accumulating Omachi alleles of qCD3, qCD7-1 and qCD7-2 was 25.0% in 2020. These QTLs for culm diameter pleiotropically increased spikelet number per panicle but did not affect days to heading or culm length. These results suggest that the Omachi alleles of qCD3, qCD7-1 and qCD7-2 are useful for breeding to increase lodging resistance and yield.

As biodiversity loss continues, there is an urgent need to develop efficient conservation measures to protect diversity with limited conservation resources. Conservation targets have generally been selected based on their population size, but more detailed assessments clarifying the phylogenetic genetic status, history, and phylogenetic uniqueness of rare species is crucial to set more appropriate and effective conservation measures. In Japan, the Act on Conservation of Endangered Species of Wild Fauna and Flora designated endangered plants with high conservation priority, but >40% of these species also grow overseas. We conducted comparative analyses based on ddRADseq and MIG-seq to evaluate the population conservation status and value of Vaccinium emarginatum and Elatostema platyphyllum which are growing across national borders at the eastern edge of their species distribution range. The analyses revealed contrasting conservation status between the two species; the Japanese population of V. emarginatum had lower genetic diversity at the individual level and phylogenetically differentiated from Taiwanese populations, while that of E. platyphyllum had higher diversity at the individual level and is a relatively recent migrant with little phylogenetical differentiation from Taiwanese populations. The two species, which share the common feature of being critically rare in Japan, showed contrasting genetic/phylogenetic characteristics. This study provided useful information for appropriate conservation measures based on species’ phylogenetic traits and genetic diversity.

Plant specialized metabolites (PSMs) are often stored as glycosides within cells and released from the roots with some chemical modifications. While isoflavones are known to function as symbiotic signals with rhizobia and to modulate the soybean rhizosphere microbiome, the underlying mechanisms of root-to-soil delivery are poorly understood. In addition to transporter-mediated secretion, the hydrolysis of isoflavone glycosides in the apoplast by an isoflavone conjugate-hydrolyzing β-glucosidase (ICHG) has been proposed but not yet verified. To clarify the role of ICHG in isoflavone supply to the rhizosphere, we have isolated two independent mutants defective in ICHG activity from a soybean high-density mutant library. In the root apoplastic fraction of ichg mutants, the isoflavone glycosides contents were significantly increased while isoflavone aglycone contents were decreased, indicating that ICHG hydrolyzes isoflavone glycosides into aglycones in the root apoplast. When grown in a field, the lack of ICHG activity considerably reduced isoflavone aglycone contents in roots and the rhizosphere soil, although the transcriptomes showed no distinct differences between the ichg mutants and WTs. Despite the change in isoflavone contents and composition of the root and rhizosphere of the mutants, root and rhizosphere bacterial communities were not distinctive from those of the WTs. Root bacterial communities and nodulation capacities of the ichg mutants did not differ from the WTs under nitrogen-deficient conditions, either. Taken together, these results indicate that ICHG elevates the accumulation of isoflavones in the soybean rhizosphere but is not essential in isoflavone-mediated plant-microbe interactions., ダイズ根圏へのイソフラボン供給量を増やす酵素を発見 --植物が機能性成分を根から土壌へ分泌するメカニズムの理解に貢献--. 京都大学プレスリリース. 2023-02-10.

This study aimed to identify quantitative trait loci (QTLs) for growth-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers in Japanese quail. A QTL mapping population of 277 F2 birds was obtained from an intercross between a male of a large-sized strain and three females of a normal-sized strain. Body weight (BW) was measured weekly from hatching to 16 weeks of age. Non-linear regression growth models of Weibull, Logistic, Gompertz, Richards, and Brody were analyzed, and growth curve parameters of Richards was selected as the best model to describe the quail growth curve of the F2 birds. Restriction-site associated DNA sequencing developed 125 SNP markers that were informative between their parental strains. The SNP markers were distributed on 16 linkage groups that spanned 795.9 centiMorgan (cM) with an average marker interval of 7.3 cM. QTL analysis of phenotypic traits revealed four main-effect QTLs. Detected QTLs were located on chromosomes 1 and 3 and were associated with BW from 4 to 16 weeks of age and asymptotic weight of Richards model at genome-wide significant at 1% or 5% level. No QTL was detected for BW from 0 to 3 weeks of age. This is the first report identified QTLs for asymptotic weight of the Richards parameter in Japanese quail. These results highlight that the combination of QTL studies and the RAD-seq method will aid future breeding programs identify genes underlying the QTL and the application of marker-assisted selection in the poultry industry, particularly the Japanese quail.

Pacific crown-of-thorns starfish (Acanthaster cf. solaris) are known to be widely distributed, and previous studies (showing limited genetic structure) suggest that there are high levels of connectivity at regional scales. It is possible, however, that methods used in previous studies (e.g., microsatellites) did not effectively resolve genetic differentiation and structure for crown-of-thorns starfish. This study examined patterns of genetic structure for populations of the A. cf. solaris in the Ryukyu Archipelago and at Miyazaki, Japan, using restriction-site-associated DNA sequencing (RADseq) analysis, to identify genome-wide single nucleotide polymorphism (SNP) loci to examine genetic diversity and gene flow of CoTS populations. We obtained 186 biallelic SNPs completely shared among 33 individuals by utilizing the published genome sequence of CoTS. Population genomic analysis with biallelic SNPs confirmed that there is limited genetic differentiation and seemingly strong connectivity and gene flow among CoTS populations in the Ryukyu Archipelago and at Miyazaki, Japan. Our findings likely reflect the rapid range expansion of CoTS through southern Japan or panmixia facilitated by the northward flowing Kuroshio Current.

How genetic variations affect gene expression dynamics of field-grown plants remains unclear. Expression quantitative trait loci (eQTL) analysis is frequently used to find genomic regions underlying gene expression polymorphisms. This approach requires transcriptome data for the complete set of the QTL mapping population under the given conditions. Therefore, only a limited range of environmental conditions is covered by a conventional eQTL analysis. We sampled sparse time series of field-grown rice from chromosome segment substitution lines (CSSLs) and conducted RNA sequencing (RNA-Seq). Then, by using statistical analysis integrating meteorological data and the RNA-Seq data, we identified 1,675 eQTLs leading to polymorphisms in expression dynamics under field conditions. A genomic region on chromosome 11 influences the expression of several defense-related genes in a time-of-day- and scaled-age-dependent manner. This includes the eQTLs that possibly influence the time-of-day- and scaled-age-dependent differences in the innate immunity between Koshihikari and Takanari. Based on the eQTL and meteorological data, we successfully predicted gene expression under environments different from training environments and in rice cultivars with more complex genotypes than the CSSLs. Our novel approach of eQTL identification facilitated the understanding of the genetic architecture of expression dynamics under field conditions, which is difficult to assess by conventional eQTL studies. The prediction of expression based on eQTLs and environmental information could contribute to the understanding of plant traits under diverse field conditions.

Latitude is correlated with environmental components that determine the distribution of biodiversity. In combination with geographic factors, latitude‐associated environmental variables are expected to influence speciation, but empirical evidence on how those factors interplay is scarce. We evaluated the genetic and environmental variation among populations in the pair of sister species Dioon sonorense–D. vovidesii, two cycads distributed along a latitudinal environmental gradient in northwestern Mexico, to reveal their demographic histories and the environmental factors involved in their divergence. Using genome‐wide loci data, we determined the species delimitation, estimated the gene flow, and compared multiple demographic scenarios of divergence. Also, we estimated the variation of climatic variables among populations and used ecological niche models to test niche overlap between species. The effect of geographic and environmental variables on the genetic variation among populations was evaluated using linear models. Our results suggest the existence of a widespread ancestral population that split into the two species ~829 ky ago. The geographic delimitation along the environmental gradient occurs in the absence of major geographic barriers, near the 28th parallel north, where a zonation of environmental seasonality exists. The northern species, D. vovidesii, occurs in more seasonal environments but retains the same niche of the southern species, D. sonorense. The genetic variation throughout populations cannot be solely explained by stochastic processes; the latitudinal‐associated seasonality has been an additive factor that strengthened the species divergence. This study represents an example of how speciation can be achieved by the effect of the latitude‐associated factors on the genetic divergence among populations.

The least absolute shrinkage and selection operator (lasso) and principal component regression (PCR) are popular methods of estimating traits from high-dimensional omics data, such as transcriptomes. The prediction accuracy of these estimation methods is highly dependent on the covariance structure, which is characterized by gene regulation networks. However, the manner in which the structure of a gene regulation network together with the sample size affects prediction accuracy has not yet been sufficiently investigated. In this study, Monte Carlo simulations are conducted to investigate the prediction accuracy for several network structures under various sample sizes. When the gene regulation network is a random graph, a sufficiently large number of observations are required to ensure good prediction accuracy with the lasso. The PCR provided poor prediction accuracy regardless of the sample size. However, a real gene regulation network is likely to exhibit a scale-free structure. In such cases, the simulation indicates that a relatively small number of observations, such as $$N=300$$ N = 300 , is sufficient to allow the accurate prediction of traits from a transcriptome with the lasso.

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.

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.

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.

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).