1
|
Song H, Dong T, Wang W, Yan X, Geng C, Bai S, Hu H. GWAS Enhances Genomic Prediction Accuracy of Caviar Yield, Caviar Color and Body Weight Traits in Sturgeons Using Whole-Genome Sequencing Data. Int J Mol Sci 2024; 25:9756. [PMID: 39273703 PMCID: PMC11395957 DOI: 10.3390/ijms25179756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2024] [Revised: 09/05/2024] [Accepted: 09/07/2024] [Indexed: 09/15/2024] Open
Abstract
Caviar yield, caviar color, and body weight are crucial economic traits in sturgeon breeding. Understanding the molecular mechanisms behind these traits is essential for their genetic improvement. In this study, we performed whole-genome sequencing on 673 Russian sturgeons, renowned for their high-quality caviar. With an average sequencing depth of 13.69×, we obtained approximately 10.41 million high-quality single nucleotide polymorphisms (SNPs). Using a genome-wide association study (GWAS) with a single-marker regression model, we identified SNPs and genes associated with these traits. Our findings revealed several candidate genes for each trait: caviar yield: TFAP2A, RPS6KA3, CRB3, TUBB, H2AFX, morc3, BAG1, RANBP2, PLA2G1B, and NYAP1; caviar color: NFX1, OTULIN, SRFBP1, PLEK, INHBA, and NARS; body weight: ACVR1, HTR4, fmnl2, INSIG2, GPD2, ACVR1C, TANC1, KCNH7, SLC16A13, XKR4, GALR2, RPL39, ACVR2A, ADCY10, and ZEB2. Additionally, using the genomic feature BLUP (GFBLUP) method, which combines linkage disequilibrium (LD) pruning markers with GWAS prior information, we improved genomic prediction accuracy by 2%, 1.9%, and 3.1% for caviar yield, caviar color, and body weight traits, respectively, compared to the GBLUP method. In conclusion, this study enhances our understanding of the genetic mechanisms underlying caviar yield, caviar color, and body weight traits in sturgeons, providing opportunities for genetic improvement of these traits through genomic selection.
Collapse
Affiliation(s)
- Hailiang Song
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
- National Innovation Center for Digital Seed Industry, Beijing 100097, China
| | - Tian Dong
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
| | - Wei Wang
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
| | - Xiaoyu Yan
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
| | - Chenfan Geng
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
| | - Song Bai
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
| | - Hongxia Hu
- Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences & Beijing Key Laboratory of Fisheries Biotechnology, Beijing 100068, China
- Key Laboratory of Sturgeon Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Hangzhou 311799, China
- National Innovation Center for Digital Seed Industry, Beijing 100097, China
| |
Collapse
|
2
|
Haque MA, Alam MZ, Iqbal A, Lee YM, Dang CG, Kim JJ. Evaluation of accuracies of genomic predictions for body conformation traits in Korean Holstein. Anim Biosci 2024; 37:555-566. [PMID: 38271974 PMCID: PMC10915218 DOI: 10.5713/ab.23.0237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 08/31/2023] [Accepted: 11/22/2023] [Indexed: 01/27/2024] Open
Abstract
OBJECTIVE This study aimed to assess the genetic parameters and accuracy of genomic predictions for twenty-four linear body conformation traits and overall conformation scores in Korean Holstein dairy cows. METHODS A dataset of 2,206 Korean Holsteins was collected, and genotyping was performed using the Illumina Bovine 50K single nucleotide polymorphism (SNP) chip. The traits investigated included body traits (stature, height at front end, chest width, body depth, angularity, body condition score, and locomotion), rump traits (rump angle, rump width, and loin strength), feet and leg traits (rear leg set, rear leg rear view, foot angle, heel depth, and bone quality), udder traits (udder depth, udder texture, udder support, fore udder attachment, front teat placement, front teat length, rear udder height, rear udder width, and rear teat placement), and overall conformation score. Accuracy of genomic predictions was assessed using the single-trait animal model genomic best linear unbiased prediction method implemented in the ASReml-SA v4.2 software. RESULTS Heritability estimates ranged from 0.10 to 0.50 for body traits, 0.21 to 0.35 for rump traits, 0.13 to 0.29 for feet and leg traits, and 0.05 to 0.46 for udder traits. Rump traits exhibited the highest average heritability (0.29), while feet and leg traits had the lowest estimates (0.21). Accuracy of genomic predictions varied among the twenty-four linear body conformation traits, ranging from 0.26 to 0.49. The heritability and prediction accuracy of genomic estimated breeding value (GEBV) for the overall conformation score were 0.45 and 0.46, respectively. The GEBVs for body conformation traits in Korean Holstein cows had low accuracy, falling below the 50% threshold. CONCLUSION The limited response to selection for body conformation traits in Korean Holsteins may be attributed to both the low heritability of these traits and the lower accuracy estimates for GEBVs. Further research is needed to enhance the accuracy of GEBVs and improve the selection response for these traits.
Collapse
Affiliation(s)
- Md Azizul Haque
- Department of Biotechnology, Yeungnam University, Gyeongsan 38541,
Korea
| | | | - Asif Iqbal
- Department of Biotechnology, Yeungnam University, Gyeongsan 38541,
Korea
| | - Yun Mi Lee
- Department of Biotechnology, Yeungnam University, Gyeongsan 38541,
Korea
| | - Chang Gwon Dang
- Animal Breeding and Genetics Division, National Institute of Animal Science, Cheonan, 31000,
Korea
| | - Jong Joo Kim
- Department of Biotechnology, Yeungnam University, Gyeongsan 38541,
Korea
| |
Collapse
|
3
|
Zhu D, Zhao Y, Zhang R, Wu H, Cai G, Wu Z, Wang Y, Hu X. Genomic prediction based on selective linkage disequilibrium pruning of low-coverage whole-genome sequence variants in a pure Duroc population. Genet Sel Evol 2023; 55:72. [PMID: 37853325 PMCID: PMC10583454 DOI: 10.1186/s12711-023-00843-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 09/14/2023] [Indexed: 10/20/2023] Open
Abstract
BACKGROUND Although the accumulation of whole-genome sequencing (WGS) data has accelerated the identification of mutations underlying complex traits, its impact on the accuracy of genomic predictions is limited. Reliable genotyping data and pre-selected beneficial loci can be used to improve prediction accuracy. Previously, we reported a low-coverage sequencing genotyping method that yielded 11.3 million highly accurate single-nucleotide polymorphisms (SNPs) in pigs. Here, we introduce a method termed selective linkage disequilibrium pruning (SLDP), which refines the set of SNPs that show a large gain during prediction of complex traits using whole-genome SNP data. RESULTS We used the SLDP method to identify and select markers among millions of SNPs based on genome-wide association study (GWAS) prior information. We evaluated the performance of SLDP with respect to three real traits and six simulated traits with varying genetic architectures using two representative models (genomic best linear unbiased prediction and BayesR) on samples from 3579 Duroc boars. SLDP was determined by testing 180 combinations of two core parameters (GWAS P-value thresholds and linkage disequilibrium r2). The parameters for each trait were optimized in the training population by five fold cross-validation and then tested in the validation population. Similar to previous GWAS prior-based methods, the performance of SLDP was mainly affected by the genetic architecture of the traits analyzed. Specifically, SLDP performed better for traits controlled by major quantitative trait loci (QTL) or a small number of quantitative trait nucleotides (QTN). Compared with two commercial SNP chips, genotyping-by-sequencing data, and an unselected whole-genome SNP panel, the SLDP strategy led to significant improvements in prediction accuracy, which ranged from 0.84 to 3.22% for real traits controlled by major or moderate QTL and from 1.23 to 11.47% for simulated traits controlled by a small number of QTN. CONCLUSIONS The SLDP marker selection method can be incorporated into mainstream prediction models to yield accuracy improvements for traits with a relatively simple genetic architecture, however, it has no significant advantage for traits not controlled by major QTL. The main factors that affect its performance are the genetic architecture of traits and the reliability of GWAS prior information. Our findings can facilitate the application of WGS-based genomic selection.
Collapse
Affiliation(s)
- Di Zhu
- State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing, China
| | - Yiqiang Zhao
- State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing, China
| | - Ran Zhang
- State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing, China
| | - Hanyu Wu
- State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing, China
- National Research Facility for Phenotypic and Genotypic Analysis of Model Animals (Beijing), China Agricultural University, Beijing, China
| | - Gengyuan Cai
- National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangdong, China
| | - Zhenfang Wu
- National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangdong, China.
| | - Yuzhe Wang
- National Research Facility for Phenotypic and Genotypic Analysis of Model Animals (Beijing), China Agricultural University, Beijing, China.
| | - Xiaoxiang Hu
- State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing, China.
| |
Collapse
|
4
|
Gayathiri E, Prakash P, Pratheep T, Ramasubburayan R, Thirumalaivasan N, Gaur A, Govindasamy R, Rengasamy KRR. Bio surfactants from lactic acid bacteria: an in-depth analysis of therapeutic properties and food formulation. Crit Rev Food Sci Nutr 2023; 64:10925-10949. [PMID: 37401803 DOI: 10.1080/10408398.2023.2230491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Healthy humans and animals commonly harbor lactic acid bacteria (LAB) on their mucosal surfaces, which are often associated with food fermentation. These microorganisms can produce amphiphilic compounds, known as microbial surface-active agents, that exhibit remarkable emulsifying activity. However, the exact functions of these microbial surfactants within the producer cells remain unclear. Consequently, there is a growing urgency to develop biosurfactant production from nonpathogenic microbes, particularly those derived from LAB. This approach aims to harness the benefits of biosurfactants while ensuring their safety and applicability. This review encompasses a comprehensive analysis of native and genetically modified LAB biosurfactants, shedding light on microbial interactions, cell signaling, pathogenicity, and biofilm development. It aims to provide valuable insights into the applications of these active substances in therapeutic use and food formulation as well as their potential biological and other benefits. By synthesizing the latest knowledge and advancements, this review contributes to the understanding and utilization of LAB biosurfactants in the food and nutritional areas.
Collapse
Affiliation(s)
- Ekambaram Gayathiri
- Department of Plant Biology and Plant Biotechnology, Guru Nanak College (Autonomous), Chennai, Tamil Nadu, India
| | | | - Thangaraj Pratheep
- Department of Biotechnology, Rathinam College of Arts and Science, Coimbatore, Tamil Nadu, India
| | - Ramasamy Ramasubburayan
- Department of Prosthodontics, Saveetha Dental College and Hospitals, Saveetha Institute of Medical and Technical Sciences, Saveetha University, Chennai, Tamil Nadu, India
| | - Natesan Thirumalaivasan
- Department of Periodontics, Saveetha Dental College and Hospitals, Saveetha Institute of Medical and Technical Sciences, Saveetha University, Chennai, Tamil Nadu, India
| | - Arti Gaur
- Department of Applied Sciences, Parul University, Vadodara, Gujarat, India
| | - Rajakumar Govindasamy
- Department of Orthodontics, Saveetha Dental College and Hospitals, Saveetha Institute of Medical and Technical Sciences, Saveetha University, Chennai, Tamil Nadu, India
| | - Kannan R R Rengasamy
- Laboratory of Natural Products and Medicinal Chemistry (LNPMC), Department of Pharmacology, Saveetha Dental College, Saveetha Institute of Medical and Technical Sciences (SIMATS), Chennai, Tamil Nadu, India
| |
Collapse
|
5
|
Yan Z, Wang Z, Zhang Q, Yue S, Yin B, Jiang Y, Shi K. Identification of whole-genome significant single nucleotide polymorphisms in candidate genes associated with body conformation traits in Chinese Holstein cattle. Anim Genet 2019; 51:141-146. [PMID: 31633203 PMCID: PMC7003999 DOI: 10.1111/age.12865] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/16/2019] [Indexed: 11/29/2022]
Affiliation(s)
- Zhengui Yan
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Zhonghua Wang
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Qin Zhang
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Shujian Yue
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Bin Yin
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Yunliang Jiang
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| | - Kerong Shi
- College of Animal Science and Technology, Shandong Key Laboratory of Animal Bioengineering and Disease Prevention, Shandong Agricultural University, Taian, Shandong, 271018, China
| |
Collapse
|
6
|
Song H, Ye S, Jiang Y, Zhang Z, Zhang Q, Ding X. Using imputation-based whole-genome sequencing data to improve the accuracy of genomic prediction for combined populations in pigs. Genet Sel Evol 2019; 51:58. [PMID: 31638889 PMCID: PMC6805481 DOI: 10.1186/s12711-019-0500-8] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Accepted: 10/07/2019] [Indexed: 11/17/2022] Open
Abstract
BACKGROUND For genomic selection in populations with a small reference population, combining populations of the same breed or populations of related breeds is an effective way to increase the size of the reference population. However, genomic predictions based on single nucleotide polymorphism (SNP)-chip genotype data using combined populations with different genetic backgrounds or from different breeds have not shown a clear advantage over using within-population or within-breed predictions. The increasing availability of whole-genome sequencing (WGS) data provides new opportunities for combined population genomic prediction. Our objective was to investigate the accuracy of genomic prediction using imputation-based WGS data from combined populations in pigs. Using 80K SNP panel genotypes, WGS genotypes, or genotypes on WGS variants that were pruned based on linkage disequilibrium (LD), three methods [genomic best linear unbiased prediction (GBLUP), single-step (ss)GBLUP, and genomic feature (GF)BLUP] were implemented with different prior information to identify the best method to improve the accuracy of genomic prediction for combined populations in pigs. RESULTS In total, 2089 and 2043 individuals with production and reproduction phenotypes, respectively, from three Yorkshire populations with different genetic backgrounds were genotyped with the PorcineSNP80 panel. Imputation accuracy from 80K to WGS variants reached 92%. The results showed that use of the WGS data compared to the 80K SNP panel did not increase the accuracy of genomic prediction in a single population, but using WGS data with LD pruning and GFBLUP with prior information did yield higher accuracy than the 80K SNP panel. For the 80K SNP panel genotypes, using the combined population resulted in a slight improvement, no change, or even a slight decrease in accuracy in comparison with the single population for GBLUP and ssGBLUP, while accuracy increased by 1 to 2.4% when using WGS data. Notably, the GFBLUP method did not perform well for both the combined population and the single populations. CONCLUSIONS The use of WGS data was beneficial for combined population genomic prediction. Simply increasing the number of SNPs to the WGS level did not increase accuracy for a single population, while using pruned WGS data based on LD and GFBLUP with prior information could yield higher accuracy than the 80K SNP panel.
Collapse
Affiliation(s)
- Hailiang Song
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Shaopan Ye
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong China
| | - Yifan Jiang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Zhe Zhang
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong China
| | - Qin Zhang
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, Taian, China
| | - Xiangdong Ding
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
7
|
Iung LHS, Petrini J, Ramírez-Díaz J, Salvian M, Rovadoscki GA, Pilonetto F, Dauria BD, Machado PF, Coutinho LL, Wiggans GR, Mourão GB. Genome-wide association study for milk production traits in a Brazilian Holstein population. J Dairy Sci 2019; 102:5305-5314. [PMID: 30904307 DOI: 10.3168/jds.2018-14811] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Accepted: 10/19/2018] [Indexed: 12/19/2022]
Abstract
Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.
Collapse
Affiliation(s)
- L H S Iung
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - J Petrini
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - J Ramírez-Díaz
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - M Salvian
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - G A Rovadoscki
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - F Pilonetto
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - B D Dauria
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - P F Machado
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - L L Coutinho
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - G R Wiggans
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, Beltsville, MD 20705-2350
| | - G B Mourão
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil.
| |
Collapse
|