1
|
Delattre M, Toda Y, Tressou J, Iwata H. Modeling soybean growth: A mixed model approach. PLoS Comput Biol 2024; 20:e1011258. [PMID: 38990979 DOI: 10.1371/journal.pcbi.1011258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Accepted: 06/17/2024] [Indexed: 07/13/2024] Open
Abstract
The evaluation of plant and animal growth, separately for genetic and environmental effects, is necessary for genetic understanding and genetic improvement of environmental responses of plants and animals. We propose to extend an existing approach that combines nonlinear mixed-effects model (NLMEM) and the stochastic approximation of the Expectation-Maximization algorithm (SAEM) to analyze genetic and environmental effects on plant growth. These tools are widely used in many fields but very rarely in plant biology. During model formulation, a nonlinear function describes the shape of growth, and random effects describe genetic and environmental effects and their variability. Genetic relationships among the varieties were also integrated into the model using a genetic relationship matrix. The SAEM algorithm was chosen as an efficient alternative to MCMC methods, which are more commonly used in the domain. It was implemented to infer the expected growth patterns in the analyzed population and the expected curves for each variety through a maximum-likelihood and a maximum-a-posteriori approaches, respectively. The obtained estimates can be used to predict the growth curves for each variety. We illustrate the strengths of the proposed approach using simulated data and soybean plant growth data obtained from a soybean cultivation experiment conducted at the Arid Land Research Center, Tottori University. In this experiment, plant height was measured daily using drones, and the growth was monitored for approximately 200 soybean cultivars for which whole-genome sequence data were available. The NLMEM approach improved our understanding of the determinants of soybean growth and can be successfully used for the genomic prediction of growth pattern characteristics.
Collapse
Affiliation(s)
- Maud Delattre
- Université Paris-Saclay, INRAE, MaIAGE, Jouy-en-Josas, France
| | - Yusuke Toda
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Jessica Tressou
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
- Paris-Saclay University-AgroParisTech-INRAE, UMR MIA-Paris-Saclay, Palaiseau, France
| | - Hiroyoshi Iwata
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
2
|
Sabag I, Bi Y, Sahoo MM, Herrmann I, Morota G, Peleg Z. Leveraging genomics and temporal high-throughput phenotyping to enhance association mapping and yield prediction in sesame. THE PLANT GENOME 2024:e20481. [PMID: 38926134 DOI: 10.1002/tpg2.20481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Accepted: 05/16/2024] [Indexed: 06/28/2024]
Abstract
Sesame (Sesamum indicum) is an important oilseed crop with rising demand owing to its nutritional and health benefits. There is an urgent need to develop and integrate new genomic-based breeding strategies to meet these future demands. While genomic resources have advanced genetic research in sesame, the implementation of high-throughput phenotyping and genetic analysis of longitudinal traits remains limited. Here, we combined high-throughput phenotyping and random regression models to investigate the dynamics of plant height, leaf area index, and five spectral vegetation indices throughout the sesame growing seasons in a diversity panel. Modeling the temporal phenotypic and additive genetic trajectories revealed distinct patterns corresponding to the sesame growth cycle. We also conducted longitudinal genomic prediction and association mapping of plant height using various models and cross-validation schemes. Moderate prediction accuracy was obtained when predicting new genotypes at each time point, and moderate to high values were obtained when forecasting future phenotypes. Association mapping revealed three genomic regions in linkage groups 6, 8, and 11, conferring trait variation over time and growth rate. Furthermore, we leveraged correlations between the temporal trait and seed-yield and applied multi-trait genomic prediction. We obtained an improvement over single-trait analysis, especially when phenotypes from earlier time points were used, highlighting the potential of using a high-throughput phenotyping platform as a selection tool. Our results shed light on the genetic control of longitudinal traits in sesame and underscore the potential of high-throughput phenotyping to detect a wide range of traits and genotypes that can inform sesame breeding efforts to enhance yield.
Collapse
Affiliation(s)
- Idan Sabag
- The Robert H. Smith Institute of Plant Sciences and Genetics in Agriculture, The Hebrew University of Jerusalem, Rehovot, Israel
- School of Animal Sciences, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, USA
| | - Ye Bi
- School of Animal Sciences, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, USA
| | - Maitreya Mohan Sahoo
- The Robert H. Smith Institute of Plant Sciences and Genetics in Agriculture, The Hebrew University of Jerusalem, Rehovot, Israel
| | - Ittai Herrmann
- The Robert H. Smith Institute of Plant Sciences and Genetics in Agriculture, The Hebrew University of Jerusalem, Rehovot, Israel
| | - Gota Morota
- School of Animal Sciences, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, USA
- Center for Advanced Innovation in Agriculture, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, USA
| | - Zvi Peleg
- The Robert H. Smith Institute of Plant Sciences and Genetics in Agriculture, The Hebrew University of Jerusalem, Rehovot, Israel
| |
Collapse
|
3
|
Toda Y, Sasaki G, Ohmori Y, Yamasaki Y, Takahashi H, Takanashi H, Tsuda M, Kajiya-Kanegae H, Tsujimoto H, Kaga A, Hirai M, Nakazono M, Fujiwara T, Iwata H. Reaction norm for genomic prediction of plant growth: modeling drought stress response in soybean. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2024; 137:77. [PMID: 38460027 PMCID: PMC10924738 DOI: 10.1007/s00122-024-04565-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Accepted: 01/30/2024] [Indexed: 03/11/2024]
Abstract
KEY MESSAGE We proposed models to predict the effects of genomic and environmental factors on daily soybean growth and applied them to soybean growth data obtained with unmanned aerial vehicles. Advances in high-throughput phenotyping technology have made it possible to obtain time-series plant growth data in field trials, enabling genotype-by-environment interaction (G × E) modeling of plant growth. Although the reaction norm is an effective method for quantitatively evaluating G × E and has been implemented in genomic prediction models, no reaction norm models have been applied to plant growth data. Here, we propose a novel reaction norm model for plant growth using spline and random forest models, in which daily growth is explained by environmental factors one day prior. The proposed model was applied to soybean canopy area and height to evaluate the influence of drought stress levels. Changes in the canopy area and height of 198 cultivars were measured by remote sensing using unmanned aerial vehicles. Multiple drought stress levels were set as treatments, and their time-series soil moisture was measured. The models were evaluated using three cross-validation schemes. Although accuracy of the proposed models did not surpass that of single-trait genomic prediction, the results suggest that our model can capture G × E, especially the latter growth period for the random forest model. Also, significant variations in the G × E of the canopy height during the early growth period were visualized using the spline model. This result indicates the effectiveness of the proposed models on plant growth data and the possibility of revealing G × E in various growth stages in plant breeding by applying statistical or machine learning models to time-series phenotype data.
Collapse
Affiliation(s)
- Yusuke Toda
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Goshi Sasaki
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Yoshihiro Ohmori
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Yuji Yamasaki
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
- Arid Land Research Center, Tottori University, Tottori, Japan
| | - Hirokazu Takahashi
- Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya, Japan
| | - Hideki Takanashi
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Mai Tsuda
- Tsukuba-Plant Innovation Research Center (T-PIRC), University of Tsukuba, Tsukuba, Japan
| | | | | | - Akito Kaga
- Institute of Crop Science, NARO, Tsukuba, Japan
| | - Masami Hirai
- RIKEN Center for Sustainable Resource Science, Tsukuba, Japan
| | - Mikio Nakazono
- Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya, Japan
| | - Toru Fujiwara
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Hiroyoshi Iwata
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan.
| |
Collapse
|
4
|
Taniguchi S, Sakamoto T, Imase R, Nonoue Y, Tsunematsu H, Goto A, Matsushita K, Ohmori S, Maeda H, Takeuchi Y, Ishii T, Yonemaru JI, Ogawa D. Prediction of heading date, culm length, and biomass from canopy-height-related parameters derived from time-series UAV observations of rice. FRONTIERS IN PLANT SCIENCE 2022; 13:998803. [PMID: 36582650 PMCID: PMC9792801 DOI: 10.3389/fpls.2022.998803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Accepted: 11/28/2022] [Indexed: 06/17/2023]
Abstract
Unmanned aerial vehicles (UAVs) are powerful tools for monitoring crops for high-throughput phenotyping. Time-series aerial photography of fields can record the whole process of crop growth. Canopy height (CH), which is vertical plant growth, has been used as an indicator for the evaluation of lodging tolerance and the prediction of biomass and yield. However, there have been few attempts to use UAV-derived time-series CH data for field testing of crop lines. Here we provide a novel framework for trait prediction using CH data in rice. We generated UAV-based digital surface models of crops to extract CH data of 30 Japanese rice cultivars in 2019, 2020, and 2021. CH-related parameters were calculated in a non-linear time-series model as an S-shaped plant growth curve. The maximum saturation CH value was the most important predictor for culm length. The time point at the maximum CH contributed to the prediction of days to heading, and was able to predict stem and leaf weight and aboveground weight, possibly reflecting the association of biomass with duration of vegetative growth. These results indicate that the CH-related parameters acquired by UAV can be useful as predictors of traits typically measured by hand.
Collapse
Affiliation(s)
- Shoji Taniguchi
- Research Center for Agricultural Information Technology, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Toshihiro Sakamoto
- Institute for Agro-Environmental Sciences, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Ryoji Imase
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Yasunori Nonoue
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Hiroshi Tsunematsu
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Akitoshi Goto
- Research Center for Agricultural Information Technology, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Kei Matsushita
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Sinnosuke Ohmori
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Hideo Maeda
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Yoshinobu Takeuchi
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Takuro Ishii
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Jun-ichi Yonemaru
- Research Center for Agricultural Information Technology, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| | - Daisuke Ogawa
- Institute of Crop Science, National Agricultural and Food Research Organization (NARO), Tsukuba, Japan
| |
Collapse
|
5
|
Zheng C, Abd-Elrahman A, Whitaker V, Dalid C. Prediction of Strawberry Dry Biomass from UAV Multispectral Imagery Using Multiple Machine Learning Methods. REMOTE SENSING 2022; 14:4511. [DOI: 10.3390/rs14184511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Biomass is a key biophysical parameter for precision agriculture and plant breeding. Fast, accurate and non-destructive monitoring of biomass enables various applications related to crop growth. In this paper, strawberry dry biomass weight was modeled using 4 canopy geometric parameters (area, average height, volume, standard deviation of height) and 25 spectral variables (5 band original reflectance values and 20 vegetation indices (VIs)) extracted from the Unmanned Aerial Vehicle (UAV) multispectral imagery. Six regression techniques—multiple linear regression (MLR), random forest (RF), support vector machine (SVM), multivariate adaptive regression splines (MARS), eXtreme Gradient Boosting (XGBoost) and artificial neural network (ANN)—were employed and evaluated for biomass prediction. The ANN had the highest accuracy in a five-fold cross-validation, with R2 of 0.89~0.93, RMSE of 7.16~8.98 g and MAE of 5.06~6.29 g. As for the other five models, the addition of VIs increased the R2 from 0.77~0.80 to 0.83~0.86, and reduced the RMSE from 8.89~9.58 to 7.35~8.09 g and the MAE from 6.30~6.70 to 5.25~5.47 g, respectively. Red-edge-related VIs, including the normalized difference red-edge index (NDRE), simple ratio vegetation index red-edge (SRRedEdge), modified simple ratio red-edge (MSRRedEdge) and chlorophyll index red and red-edge (CIred&RE), were the most influential VIs for biomass modeling. In conclusion, the combination of canopy geometric parameters and VIs obtained from the UAV imagery was effective for strawberry dry biomass estimation using machine learning models.
Collapse
|
6
|
Optimization of UAV-Based Imaging and Image Processing Orthomosaic and Point Cloud Approaches for Estimating Biomass in a Forage Crop. REMOTE SENSING 2022. [DOI: 10.3390/rs14102396] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Forage and field peas provide essential nutrients for livestock diets, and high-quality field peas can influence livestock health and reduce greenhouse gas emissions. Above-ground biomass (AGBM) is one of the vital traits and the primary component of yield in forage pea breeding programs. However, a standard method of AGBM measurement is a destructive and labor-intensive process. This study utilized an unmanned aerial vehicle (UAV) equipped with a true-color RGB and a five-band multispectral camera to estimate the AGBM of winter pea in three breeding trials (two seed yields and one cover crop). Three processing techniques—vegetation index (VI), digital surface model (DSM), and 3D reconstruction model from point clouds—were used to extract the digital traits (height and volume) associated with AGBM. The digital traits were compared with the ground reference data (measured plant height and harvested AGBM). The results showed that the canopy volume estimated from the 3D model (alpha shape, α = 1.5) developed from UAV-based RGB imagery’s point clouds provided consistent and high correlation with fresh AGBM (r = 0.78–0.81, p < 0.001) and dry AGBM (r = 0.70–0.81, p < 0.001), compared with other techniques across the three trials. The DSM-based approach (height at 95th percentile) had consistent and high correlation (r = 0.71–0.95, p < 0.001) with canopy height estimation. Using the UAV imagery, the proposed approaches demonstrated the potential for estimating the crop AGBM across winter pea breeding trials.
Collapse
|
7
|
Cao P, Zhao Y, Wu F, Xin D, Liu C, Wu X, Lv J, Chen Q, Qi Z. Multi-Omics Techniques for Soybean Molecular Breeding. Int J Mol Sci 2022; 23:4994. [PMID: 35563386 PMCID: PMC9099442 DOI: 10.3390/ijms23094994] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 04/22/2022] [Accepted: 04/28/2022] [Indexed: 02/04/2023] Open
Abstract
Soybean is a major crop that provides essential protein and oil for food and feed. Since its origin in China over 5000 years ago, soybean has spread throughout the world, becoming the second most important vegetable oil crop and the primary source of plant protein for global consumption. From early domestication and artificial selection through hybridization and ultimately molecular breeding, the history of soybean breeding parallels major advances in plant science throughout the centuries. Now, rapid progress in plant omics is ushering in a new era of precision design breeding, exemplified by the engineering of elite soybean varieties with specific oil compositions to meet various end-use targets. The assembly of soybean reference genomes, made possible by the development of genome sequencing technology and bioinformatics over the past 20 years, was a great step forward in soybean research. It facilitated advances in soybean transcriptomics, proteomics, metabolomics, and phenomics, all of which paved the way for an integrated approach to molecular breeding in soybean. In this review, we summarize the latest progress in omics research, highlight novel findings made possible by omics techniques, note current drawbacks and areas for further research, and suggest that an efficient multi-omics approach may accelerate soybean breeding in the future. This review will be of interest not only to soybean breeders but also to researchers interested in the use of cutting-edge omics technologies for crop research and improvement.
Collapse
Affiliation(s)
- Pan Cao
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Ying Zhao
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Fengjiao Wu
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Dawei Xin
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Chunyan Liu
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Xiaoxia Wu
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Jian Lv
- Department of Innovation, Syngenta Biotechnology China, Beijing 102206, China
| | - Qingshan Chen
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| | - Zhaoming Qi
- College of Agriculture, Northeast Agricultural University, Harbin 150030, China; (P.C.); (Y.Z.); (F.W.); (D.X.); (C.L.)
| |
Collapse
|
8
|
Toda Y, Sasaki G, Ohmori Y, Yamasaki Y, Takahashi H, Takanashi H, Tsuda M, Kajiya-Kanegae H, Lopez-Lozano R, Tsujimoto H, Kaga A, Nakazono M, Fujiwara T, Baret F, Iwata H. Genomic Prediction of Green Fraction Dynamics in Soybean Using Unmanned Aerial Vehicles Observations. FRONTIERS IN PLANT SCIENCE 2022; 13:828864. [PMID: 35371133 PMCID: PMC8966771 DOI: 10.3389/fpls.2022.828864] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 02/21/2022] [Indexed: 05/25/2023]
Abstract
With the widespread use of high-throughput phenotyping systems, growth process data are expected to become more easily available. By applying genomic prediction to growth data, it will be possible to predict the growth of untested genotypes. Predicting the growth process will be useful for crop breeding, as variability in the growth process has a significant impact on the management of plant cultivation. However, the integration of growth modeling and genomic prediction has yet to be studied in depth. In this study, we implemented new prediction models to propose a novel growth prediction scheme. Phenotype data of 198 soybean germplasm genotypes were acquired for 3 years in experimental fields in Tottori, Japan. The longitudinal changes in the green fractions were measured using UAV remote sensing. Then, a dynamic model was fitted to the green fraction to extract the dynamic characteristics of the green fraction as five parameters. Using the estimated growth parameters, we developed models for genomic prediction of the growth process and tested whether the inclusion of the dynamic model contributed to better prediction of growth. Our proposed models consist of two steps: first, predicting the parameters of the dynamics model with genomic prediction, and then substituting the predicted values for the parameters of the dynamics model. By evaluating the heritability of the growth parameters, the dynamic model was able to effectively extract genetic diversity in the growth characteristics of the green fraction. In addition, the proposed prediction model showed higher prediction accuracy than conventional genomic prediction models, especially when the future growth of the test population is a prediction target given the observed values in the first half of growth as training data. This indicates that our model was able to successfully combine information from the early growth period with phenotypic data from the training population for prediction. This prediction method could be applied to selection at an early growth stage in crop breeding, and could reduce the cost and time of field trials.
Collapse
Affiliation(s)
- Yusuke Toda
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Goshi Sasaki
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Yoshihiro Ohmori
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Yuji Yamasaki
- Arid Land Research Center, Tottori University, Tottori, Japan
| | - Hirokazu Takahashi
- Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya, Japan
| | - Hideki Takanashi
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Mai Tsuda
- Tsukuba-Plant Innovation Research Center (T-PIRC), University of Tsukuba, Tsukuba, Japan
| | - Hiromi Kajiya-Kanegae
- Research Center for Agricultural Information Technology, National Agriculture and Food Research Organization (NARO), Tokyo, Japan
| | - Raul Lopez-Lozano
- Joint Research Unit of Mediterranean Environment and Modelling of Agroecosystems, National Research Institute for Agriculture, Food and Environment (INRAE), Avignon, France
| | | | - Akito Kaga
- Institute of Crop Science, National Agriculture and Food Research Organization (NARO), Tsukuba, Japan
| | - Mikio Nakazono
- Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya, Japan
| | - Toru Fujiwara
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Frederic Baret
- Joint Research Unit of Mediterranean Environment and Modelling of Agroecosystems, National Research Institute for Agriculture, Food and Environment (INRAE), Avignon, France
| | - Hiroyoshi Iwata
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| |
Collapse
|