1
|
Manggabarani AM, Hashiguchi T, Hashiguchi M, Hayashi A, Kikuchi M, Mustamin Y, Bamba M, Kodama K, Tanabata T, Isobe S, Tanaka H, Akashi R, Nakaya A, Sato S. Construction of prediction models for growth traits of soybean cultivars based on phenotyping in diverse genotype and environment combinations. DNA Res 2022; 29:6653298. [PMID: 35916715 PMCID: PMC9358015 DOI: 10.1093/dnares/dsac024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Indexed: 11/14/2022] Open
Abstract
Abstract
As soybean cultivars are adapted to a relatively narrow range of latitude, the effects of climate changes are estimated to be severe. To address this issue, it is important to improve our understanding of the effects of climate change by applying the simulation model including both genetic and environmental factors with their interactions (G×E). To achieve this goal, we conducted the field experiments for soybean core collections using multiple sowing times in multi-latitudinal fields. Sowing time shifts altered the flowering time (FT) and growth phenotypes, and resulted in increasing the combinations of genotypes and environments. Genome-wide association studies for the obtained phenotypes revealed the effects of field and sowing time to the significance of detected alleles, indicating the presence of G×E. By using accumulated phenotypic and environmental data in 2018 and 2019, we constructed multiple regression models for FT and growth pattern. Applicability of the constructed models was evaluated by the field experiments in 2020 including a novel field, and high correlation between the predicted and measured values was observed, suggesting the robustness of the models. The models presented here would allow us to predict the phenotype of the core collections in a given environment.
Collapse
Affiliation(s)
| | - Takuyu Hashiguchi
- Faculty of Agriculture, University of Miyazaki , Miyazaki 889-2192, Japan
| | | | - Atsushi Hayashi
- Kazusa DNA Research Institute , Kisarazu, Chiba 292-0818, Japan
| | - Masataka Kikuchi
- Graduate School of Medicine, Osaka University , Suita, Osaka 565-0871, Japan
| | - Yusdar Mustamin
- Graduate School of Life Sciences, Tohoku University , Sendai, Miyagi 980-8577, Japan
| | - Masaru Bamba
- Graduate School of Life Sciences, Tohoku University , Sendai, Miyagi 980-8577, Japan
| | - Kunihiro Kodama
- Kazusa DNA Research Institute , Kisarazu, Chiba 292-0818, Japan
| | | | - Sachiko Isobe
- Kazusa DNA Research Institute , Kisarazu, Chiba 292-0818, Japan
| | - Hidenori Tanaka
- Faculty of Agriculture, University of Miyazaki , Miyazaki 889-2192, Japan
| | - Ryo Akashi
- Faculty of Agriculture, University of Miyazaki , Miyazaki 889-2192, Japan
| | - Akihiro Nakaya
- Graduate School of Medicine, Osaka University , Suita, Osaka 565-0871, Japan
- Graduate School of Frontier Sciences, The University of Tokyo , Kashiwa, Chiba 277-0882, Japan
| | - Shusei Sato
- Graduate School of Life Sciences, Tohoku University , Sendai, Miyagi 980-8577, Japan
| |
Collapse
|
2
|
Mendonça HC, Pereira LFP, Maldonado dos Santos JV, Meda AR, Sant’ Ana GC. Genetic Diversity and Selection Footprints in the Genome of Brazilian Soybean Cultivars. FRONTIERS IN PLANT SCIENCE 2022; 13:842571. [PMID: 35432410 PMCID: PMC9006619 DOI: 10.3389/fpls.2022.842571] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Accepted: 02/14/2022] [Indexed: 06/14/2023]
Abstract
Although Brazil is currently the largest soybean producer in the world, only a small number of studies have analyzed the genetic diversity of Brazilian soybean. These studies have shown the existence of a narrow genetic base. The objectives of this work were to analyze the population structure and genetic diversity, and to identify selection signatures in the genome of soybean germplasms from different companies in Brazil. A panel consisting of 343 soybean lines from Brazil, North America, and Asia was genotyped using genotyping by sequencing (GBS). Population structure was assessed by Bayesian and multivariate approaches. Genetic diversity was analyzed using metrics such as the fixation index, nucleotide diversity, genetic dissimilarity, and linkage disequilibrium. The software BayeScan was used to detect selection signatures between Brazilian and Asian accessions as well as among Brazilian germplasms. Region of origin, company of origin, and relative maturity group (RMG) all had a significant influence on population structure. Varieties belonging to the same company and especially to the same RMG exhibited a high level of genetic similarity. This result was exacerbated among early maturing accessions. Brazilian soybean showed significantly lower genetic diversity when compared to Asian accessions. This was expected, because the crop's region of origin is its main genetic diversity reserve. We identified 7 genomic regions under selection between the Brazilian and Asian accessions, and 27 among Brazilian varieties developed by different companies. Associated with these genomic regions, we found 96 quantitative trait loci (QTLs) for important soybean breeding traits such as flowering, maturity, plant architecture, productivity components, pathogen resistance, and seed composition. Some of the QTLs associated with the markers under selection have genes of great importance to soybean's regional adaptation. The results reported herein allowed to expand the knowledge about the organization of the genetic variability of the Brazilian soybean germplasm. Furthermore, it was possible to identify genomic regions under selection possibly associated with the adaptation of soybean to Brazilian environments.
Collapse
Affiliation(s)
| | - Luiz Filipe Protasio Pereira
- Centro de Ciências Biológicas, State University of Londrina, Londrina, Brazil
- Laboratório de Biotecnologia, Instituto de Desenvolvimento Rural do Paraná, Embrapa Café, Londrina, Brazil
| | | | | | | |
Collapse
|
3
|
Díaz P, Sarmiento F, Mathew B, Ballvora A, Mosquera Vásquez T. Genomic regions associated with physiological, biochemical and yield-related responses under water deficit in diploid potato at the tuber initiation stage revealed by GWAS. PLoS One 2021; 16:e0259690. [PMID: 34748612 PMCID: PMC8575265 DOI: 10.1371/journal.pone.0259690] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Accepted: 10/24/2021] [Indexed: 11/19/2022] Open
Abstract
Water deficit, which is increasing with climate change, is a serious threat to agricultural sustainability worldwide. Dissection of the genetic architecture of water deficit responses is highly desirable for developing water-deficit tolerant potato cultivars and enhancing the resilience of existing cultivars. This study examined genetic variation in response to water deficit in a panel of diploid potato and identified the QTL governing this trait via a genome-wide association study (GWAS). A panel of 104 diploid potato accessions were evaluated under both well-watered and water deficit treatments at tuber initiation stage. Drought stress index (DTI) was calculated to assess tolerance of the diploid potato genotypes to water deficit. The GWAS was conducted using a matrix of 47K single nucleotide polymorphisms (SNP), recently available for this population. We are reporting 38 QTL, seven for well-watered conditions, twenty-two for water deficit conditions and nine for DTI which explain between 12.6% and 44.1% of the phenotypic variance. A set of 6 QTL were found to be associated with more than one variable. Marker WDP-9.21 was found associated with tuber fresh weigh under WD and gene annotation analysis revealed co-localization with the Glucan/water dikinase (GWD) gene. Of the nine QTL detected from DTI on chromosomes 2,3,5,8,10 and 12, three candidate genes with a feasible role in water deficit response were identified. The findings of this study can be used in marker-assisted selection (MAS) for water- deficit tolerance breeding in potato.
Collapse
Affiliation(s)
- Paula Díaz
- Departamento de Agronomía, Facultad de Ciencias Agrarias, Universidad Nacional de Colombia-Sede Bogotá, Bogotá, Colombia
| | - Felipe Sarmiento
- Departamento de Biología, Facultad de Ciencias, Universidad Nacional de Colombia-Sede Bogotá, Bogotá, Colombia
| | - Boby Mathew
- Bayer CropScience, Monheim am Rhein, Germany
| | - Agim Ballvora
- Institute of Crop Science and Resource Conservation Plant Breeding, University of Bonn, Bonn, Germany
| | - Teresa Mosquera Vásquez
- Departamento de Agronomía, Facultad de Ciencias Agrarias, Universidad Nacional de Colombia-Sede Bogotá, Bogotá, Colombia
| |
Collapse
|
4
|
Kanapin AA, Sokolkova AB, Samsonova AA, Shchegolkov AV, Boldyrev SV, Aupova AF, Khaitovich PE, Nuzhdin SV, Samsonova MG. Genetic Variants Associated with Productivity and Contents of Protein and Oil in Soybeans. Biophysics (Nagoya-shi) 2020. [DOI: 10.1134/s0006350920020074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
5
|
Sun F, Xu M, Park C, Dwiyanti MS, Nagano AJ, Zhu J, Watanabe S, Kong F, Liu B, Yamada T, Abe J. Characterization and quantitative trait locus mapping of late-flowering from a Thai soybean cultivar introduced into a photoperiod-insensitive genetic background. PLoS One 2019; 14:e0226116. [PMID: 31805143 PMCID: PMC6894811 DOI: 10.1371/journal.pone.0226116] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 11/19/2019] [Indexed: 11/18/2022] Open
Abstract
The timing of both flowering and maturation determine crop adaptability and productivity. Soybean (Glycine max) is cultivated across a wide range of latitudes. The molecular-genetic mechanisms for flowering in soybean have been determined for photoperiodic responses to long days (LDs), but remain only partially determined for the delay of flowering under short-day conditions, an adaptive trait of cultivars grown in lower latitudes. Here, we characterized the late-flowering (LF) habit introduced from the Thai cultivar K3 into a photoperiod-insensitive genetic background under different photo-thermal conditions, and we analyzed the genetic basis using quantitative trait locus (QTL) mapping. The LF habit resulted from a basic difference in the floral induction activity and from the suppression of flowering, which was caused by red light-enriched LD lengths and higher temperatures, during which FLOWERING LOCUS T (FT) orthologs, FT2a and FT5a, were strongly down-regulated. QTL mapping using gene-specific markers for flowering genes E2, FT2a and FT5a and 829 single nucleotide polymorphisms obtained from restriction-site associated DNA sequencing detected three QTLs controlling the LF habit. Of these, a QTL harboring FT2a exhibited large and stable effects under all the conditions tested. A resequencing analysis detected a nonsynonymous substitution in exon 4 of FT2a from K3, which converted the glycine conserved in FT-like proteins to the aspartic acid conserved in TERMINAL FLOWER 1-like proteins (floral repressors), suggesting a functional depression in the FT2a protein from K3. The effects of the remaining two QTLs, likely corresponding to E2 and FT5a, were environment dependent. Thus, the LF habit from K3 may be caused by the functional depression of FT2a and the down-regulation of two FT genes by red light-enriched LD conditions and high temperatures.
Collapse
Affiliation(s)
- Fei Sun
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Meilan Xu
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
- Key Laboratory of Soybean Molecular Design Breeding, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Harbin, China
| | - Cheolwoo Park
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
| | | | | | - Jianghui Zhu
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
| | | | - Fanjiang Kong
- School of Life Sciences, Guangzhou University, Guangzhou, China
| | - Baohui Liu
- School of Life Sciences, Guangzhou University, Guangzhou, China
| | - Tetsuya Yamada
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Jun Abe
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Hokkaido, Japan
| |
Collapse
|
6
|
Wang J, Zhao X, Wang W, Qu Y, Teng W, Qiu L, Zheng H, Han Y, Li W. Genome-wide association study of inflorescence length of cultivated soybean based on the high-throughout single-nucleotide markers. Mol Genet Genomics 2019; 294:607-620. [PMID: 30739204 DOI: 10.1007/s00438-019-01533-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Accepted: 01/31/2019] [Indexed: 11/25/2022]
Abstract
As an important and complex trait, inflorescence length (IL) of soybean [Glycine max (L.) Merr.] significantly affected seed yields. Therefore, elucidating molecular basis of inflorescence architecture, especially for IL, was important for improving soybean yield potentials. Longer IL meaned to have more pod and seed in soybean. Hence, increasing IL and improving yield are targets for soybean breeding. In this study, a association panel, comprising 283 diverse samples, was used to dissect the genetic basis of IL based on genome-wide association analysis (GWAS) and haplotype analysis. GWAS and haplotype analysis were conducted through high-throughout single-nucleotide polymorphisms (SNP) developed by SLAF-seq methodology. A total of 39, 057 SNPs (minor allele frequency ≥ 0.2 and missing data ≤ 10%) were utilized to evaluate linkage disequilibrium (LD) level in the tested association panel. A total of 30 association signals were identified to be associated with IL via GWAS. Among them, 13 SNPs were novel, and another 17 SNPs were overlapped or located near the linked regions of known quantitative trait nucleotide (QTN) with soybean seed yield or yield component. The functional genes, located in the 200-kb genomic region of each peak SNP, were considered as candidate genes, such as the cell division/ elongation, specific enzymes, and signaling or transport of specific proteins. These genes have been reported to participant in the regulation of IL. Ten typical long-IL lines and ten typical short-IL lines were re-sequencing, and then, six SNPs from five genes were obtained based on candidate gene-based association. In addition, 42 haplotypes were defined based on haplotype analysis. Of them, 11 haplotypes were found to regulate long IL (> 14 mm) in soybean. The identified 30 QTN with beneficial alleles and their candidate genes might be valuable for dissecting the molecular mechanisms of IL and further improving the yield potential of soybean.
Collapse
Affiliation(s)
- Jinyang Wang
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Xue Zhao
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Wei Wang
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Yingfan Qu
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Weili Teng
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Lijuan Qiu
- Institute of Crop Science, National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI) Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Hongkun Zheng
- Bioinformatics Division, Biomarker Technologies Corporation, Beijing, 101300, China
| | - Yingpeng Han
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China.
| | - Wenbin Li
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China.
| |
Collapse
|
7
|
Zatybekov A, Abugalieva S, Didorenko S, Rsaliyev A, Turuspekov Y. GWAS of a soybean breeding collection from South East and South Kazakhstan for resistance to fungal diseases. Vavilovskii Zhurnal Genet Selektsii 2018. [DOI: 10.18699/vj18.392] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Soybean (Glycine max(L.) Merr) is an essential food, feed, and technical culture. In Kazakhstan the area under soybean is increasing every year, helping to solve the problem of protein deficiency in human nutrition and animal feeding. One of the main problems of soybean production is fungal diseases causing yields losses of up to 30 %. Modern genomic studies can be applied to facilitate efficient breeding research for improvement of soybean fungal disease tolerance. Therefore, the objective of this genome-wide association study (GWAS) was analysis of a soybean collection consisting of 182 accessions in relation to fungal diseases in the conditions of South East and South Kazakhstan. Field evaluation of the soybean collection suggested thatFusariumspp. andCercospora sojinaaffected plants in the South region (RIBSP), andSeptoria glycines– in the South East region (KRIAPP). The major objective of the study was identification of QTL associated with resistance to fusarium root rot (FUS), frogeye leaf spot (FLS), and brown spot (BS). GWAS using 4 442 SNP (single nucleotide polymorphism) markers of Illumina iSelect array allowed for identification of fifteen marker trait associations (MTA) resistant to the three diseases at two different stages of growth. Two QTL both for FUS (chromosomes 13 and 17) and BS (chromosomes 14 and 17) were genetically mapped, including one presumably novel QTL for BS (chromosome 17). Also, five presumably novel QTL for FLS were genetically mapped on chromosomes 2, 7, and 15. The results can be used for improvement of the local breeding projects based on marker-assisted selection approach.
Collapse
|
8
|
Shen Y, Liu J, Geng H, Zhang J, Liu Y, Zhang H, Xing S, Du J, Ma S, Tian Z. De novo assembly of a Chinese soybean genome. SCIENCE CHINA. LIFE SCIENCES 2018; 61:871-884. [PMID: 30062469 DOI: 10.1007/s11427-018-9360-0] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Accepted: 07/05/2018] [Indexed: 10/28/2022]
Abstract
Soybean was domesticated in China and has become one of the most important oilseed crops. Due to bottlenecks in their introduction and dissemination, soybeans from different geographic areas exhibit extensive genetic diversity. Asia is the largest soybean market; therefore, a high-quality soybean reference genome from this area is critical for soybean research and breeding. Here, we report the de novo assembly and sequence analysis of a Chinese soybean genome for "Zhonghuang 13" by a combination of SMRT, Hi-C and optical mapping data. The assembled genome size is 1.025 Gb with a contig N50 of 3.46 Mb and a scaffold N50 of 51.87 Mb. Comparisons between this genome and the previously reported reference genome (cv. Williams 82) uncovered more than 250,000 structure variations. A total of 52,051 protein coding genes and 36,429 transposable elements were annotated for this genome, and a gene co-expression network including 39,967 genes was also established. This high quality Chinese soybean genome and its sequence analysis will provide valuable information for soybean improvement in the future.
Collapse
Affiliation(s)
- Yanting Shen
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, 100101, China
- University of Chinese Academy of Sciences, Beijing, 100039, China
| | - Jing Liu
- Provincial Key Laboratory of Agrobiology, Institute of Crop Germplasm and Biotechnology, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
| | - Haiying Geng
- School of Life Sciences, University of Science and Technology of China, Hefei, 230027, China
| | - Jixiang Zhang
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, 100101, China
- University of Chinese Academy of Sciences, Beijing, 100039, China
| | - Yucheng Liu
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, 100101, China
- University of Chinese Academy of Sciences, Beijing, 100039, China
| | | | - Shilai Xing
- Berry Genomics Corporation, Beijing, 100015, China
| | - Jianchang Du
- Provincial Key Laboratory of Agrobiology, Institute of Crop Germplasm and Biotechnology, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China.
| | - Shisong Ma
- School of Life Sciences, University of Science and Technology of China, Hefei, 230027, China.
| | - Zhixi Tian
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, 100101, China.
- University of Chinese Academy of Sciences, Beijing, 100039, China.
| |
Collapse
|
9
|
Arriagada O, do Amaral Junior AT, Mora F. Thirteen years under arid conditions: exploring marker-trait associations in Eucalyptus cladocalyx for complex traits related to flowering, stem form and growth. BREEDING SCIENCE 2018; 68:367-374. [PMID: 30100804 PMCID: PMC6081299 DOI: 10.1270/jsbbs.17131] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/27/2017] [Accepted: 04/01/2018] [Indexed: 06/08/2023]
Abstract
We present an association analysis for seven key traits related to flowering, stem form and growth in Eucalyptus cladocalyx, a tree species suitable for low rainfall sites, using a long-term progeny trial with 49 open-pollinated maternal families in the southern Atacama Desert, Chile. The progeny trial was carried out in an arid environment with a mean annual rainfall of 152 mm. Simple sequence repeats (SSR) from a full consensus map of Eucalyptus were used for genotyping 245 individual trees. Twenty-three significant marker-trait associations were identified, explaining between 5.9 and 23.7% of the phenotypic variance. The marker EMBRA101 located on LG10 at 56.5 cM was concomitantly associated with diameter at breast height and tree height. Nine SSR were significantly associated with stem forking and stem straightness, explaining between 5.9 and 14.8% of the phenotypic variation. To our knowledge, this is the first study reporting a SSR-based association mapping analysis for stem form traits in Eucalyptus. These results provide novel and valuable information for understanding the genetic base of key traits in E. cladocalyx for breeding purposes under arid conditions.
Collapse
Affiliation(s)
- Osvin Arriagada
- Institute of Biological Sciences, University of Talca,
2 Norte 685, 3460000 Talca,
Chile
| | - Antonio Teixeira do Amaral Junior
- Laboratório de Melhoramento Genético Vegetal, Universidade Estadual do Norte Fluminense Darcy Ribeiro,
Av. Alberto Lamego 2000, 28013-602 Campos dos Goytacazes,
Brazil
| | - Freddy Mora
- Institute of Biological Sciences, University of Talca,
2 Norte 685, 3460000 Talca,
Chile
| |
Collapse
|