1
|
An ZX, Shi LG, Hou GY, Zhou HL, Xun WJ. Genetic diversity and selection signatures in Hainan black goats revealed by whole-genome sequencing data. Animal 2024; 18:101147. [PMID: 38843669 DOI: 10.1016/j.animal.2024.101147] [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] [Received: 10/30/2023] [Revised: 03/19/2024] [Accepted: 03/22/2024] [Indexed: 06/22/2024] Open
Abstract
Understanding the genetic characteristics of indigenous goat breeds is crucial for their conservation and breeding efforts. Hainan black goats, as a native breed of south China's tropical island province of Hainan, possess distinctive traits such as black hair, a moderate growth rate, good meat quality, and small body size. However, they exhibit exceptional resilience to rough feeding conditions, possess high-quality meat, and show remarkable resistance to stress and heat. In this study, we resequenced the whole genome of Hainan black goats to study the economic traits and genetic basis of these goats, we leveraged whole-genome sequencing data from 33 Hainan black goats to analyze single nucleotide polymorphism (SNP) density, Runs of homozygosity (ROH), Integrated Haplotype Score (iHS), effective population size (Ne), Nucleotide diversity Analysis (Pi) and selection characteristics. Our findings revealed that Hainan black goats harbor a substantial degree of genetic variation, with a total of 23 608 983 SNPs identified. Analysis of ROHs identified 53 710 segments, predominantly composed of short fragments, with inbreeding events mainly occurring in ancient ancestors, the estimates of inbreeding based on ROH in Hainan black goats typically exhibit moderate values ranging from 0.107 to 0.186. This is primarily attributed to significant declines in the effective population size over recent generations. Moreover, we identified 921 candidate genes within the intersection candidate region of ROH and iHS. Several of these genes are associated with crucial traits such as immunity (PTPRC, HYAL1, HYAL2, HYAL3, CENPE and PKN1), heat tolerance (GNG2, MAPK8, CAPN2, SLC1A1 and LEPR), meat quality (ACOX1, SSTR1, CAMK2B, PPP2CA and PGM1), cashmere production (AKT4, CHRM2, OXTR, AKT3, HMCN1 and CDK19), and stress resistance (TLR2, IFI44, ENPP1, STK3 and NFATC1). The presence of these genes may be attributed to the genetic adaptation of Hainan black goats to local climate conditions. The insights gained from this study provide valuable references and a solid foundation for the preservation, breeding, and utilization of Hainan black goats and their valuable genetic resources.
Collapse
Affiliation(s)
- Z X An
- School of Tropical Agriculture and Forestry, Hainan University, Haikou 570228, China; Tropical Crops Genetic Resources Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571100, China
| | - L G Shi
- Tropical Crops Genetic Resources Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571100, China
| | - G Y Hou
- Tropical Crops Genetic Resources Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571100, China
| | - H L Zhou
- Zhanjiang Experimental Station, Chinese Academy of Tropical Agricultural Sciences, Zhanjiang 524000, China
| | - W J Xun
- School of Tropical Agriculture and Forestry, Hainan University, Haikou 570228, China.
| |
Collapse
|
2
|
Tolone M, Sardina MT, Senczuk G, Chessari G, Criscione A, Moscarelli A, Riggio S, Rizzuto I, Di Gerlando R, Portolano B, Mastrangelo S. Genomic Tools for the Characterization of Local Animal Genetic Resources: Application in Mascaruna Goat. Animals (Basel) 2022; 12:2840. [PMID: 36290231 PMCID: PMC9597745 DOI: 10.3390/ani12202840] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Revised: 10/13/2022] [Accepted: 10/18/2022] [Indexed: 02/06/2024] Open
Abstract
Italy contains a large number of local goat populations, some of which do not have a recognized genetic structure. The "Mascaruna" is a goat population reared for milk production in Sicily. In this study, a total of 72 individuals were genotyped with the Illumina Goat_IGGC_65K_v2 BeadChip with the aim to characterize the genetic diversity, population structure and relatedness with another 31 Italian goat populations. The results displayed a moderate level of genetic variability for Mascaruna, in concordance with the estimated values for Italian goats. Runs of homozygosity islands are linked to genes involved in milk production, immune response and local adaptation. Population structure analyses separated Mascaruna from the other goat populations, indicating a clear genetic differentiation. Although they are not conclusive, our current results represent a starting point for the creation of monitoring and conservation plans. Additional analyses and a wider sampling would contribute to refine and validate these results. Finally, our study describing the diversity and structure of Mascaruna confirms the usefulness of applied genomic analyses as valid tools for the study of the local uncharacterized genetic resources.
Collapse
Affiliation(s)
- Marco Tolone
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Maria Teresa Sardina
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Gabriele Senczuk
- Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, 86100 Campobasso, Italy
| | - Giorgio Chessari
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, 95131 Catania, Italy
| | - Andrea Criscione
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, 95131 Catania, Italy
| | - Angelo Moscarelli
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Silvia Riggio
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Ilaria Rizzuto
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Rosalia Di Gerlando
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Baldassare Portolano
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
| |
Collapse
|
3
|
Torres-Hernández G, Maldonado-Jáquez JA, Granados-Rivera LD, Salinas-González H, Castillo-Hernández G. Status quo of genetic improvement in local goats: a review. Arch Anim Breed 2022; 65:207-221. [PMID: 35693297 PMCID: PMC9176210 DOI: 10.5194/aab-65-207-2022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Accepted: 04/14/2022] [Indexed: 11/26/2022] Open
Abstract
This review aims to summarize and synthesize the
fragmented information available on the genetic improvement of local goats
(criollo, indigenous, native) on the American and other continents, where
populations with these goats have an important role in food security and the
economy of rural communities, as well as in conservation of biodiversity and
productivity improvement. Topics such as the current state of goat
production globally, conservation programs, resistance to parasites and
diseases, use of phenotypical characteristics and genomic information, and
molecular markers for genetic improvement are addressed. The main
challenges, opportunities, and limitations described in recent literature
concerning local goats in the immediate future are discussed.
Collapse
Affiliation(s)
| | - Jorge Alonso Maldonado-Jáquez
- Colegio de Postgraduados-Campus Montecillo, 56230 Montecillo, Estado
de México, México
- Instituto Nacional de Investigaciones Forestales, Agrícolas y
Pecuarias, Centro de Investigación Regional Norte Centro, Campo
Experimental La Laguna, 27440 Matamoros, Coahuila, México
| | - Lorenzo Danilo Granados-Rivera
- Instituto Nacional de Investigaciones Forestales, Agrícolas y
Pecuarias, Centro de Investigación Regional Noreste, Campo Experimental
General Terán, 67400 General Terán, Nuevo León, México
| | | | - Gabriela Castillo-Hernández
- Colegio de Postgraduados-Campus Montecillo, 56230 Montecillo, Estado
de México, México
- Facultad de Estudios
Superiores Cuautitlán, Universidad Nacional Autónoma de México, 54714 Cuautitlán Izcalli, Estado de
México, México
| |
Collapse
|
4
|
Li G, Tang J, Huang J, Jiang Y, Fan Y, Wang X, Ren J. Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds. Front Genet 2022; 13:774196. [PMID: 35559012 PMCID: PMC9086400 DOI: 10.3389/fgene.2022.774196] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 03/14/2022] [Indexed: 11/13/2022] Open
Abstract
Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.
Collapse
Affiliation(s)
- Guixin Li
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Jianhong Tang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China.,Laboratory Animal Engineering Research Center of Ganzhou, Gannan Medical University, Ganzhou, China
| | - Jinyan Huang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Yongchuang Jiang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Yin Fan
- Department of Animal Science, Jiangxi Biotech Vocational College, Nanchang, China
| | - Xiaopeng Wang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Jun Ren
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| |
Collapse
|
5
|
Dadousis C, Ablondi M, Cipolat-Gotet C, van Kaam JT, Marusi M, Cassandro M, Sabbioni A, Summer A. Genomic inbreeding coefficients using imputed genotypes: Assessing different estimators in Holstein-Friesian dairy cows. J Dairy Sci 2022; 105:5926-5945. [DOI: 10.3168/jds.2021-21125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2021] [Accepted: 03/08/2022] [Indexed: 11/19/2022]
|
6
|
Wang X, Li G, Ruan D, Zhuang Z, Ding R, Quan J, Wang S, Jiang Y, Huang J, Gu T, Hong L, Zheng E, Li Z, Cai G, Wu Z, Yang J. Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines. Front Vet Sci 2022; 9:832633. [PMID: 35350434 PMCID: PMC8957889 DOI: 10.3389/fvets.2022.832633] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 01/24/2022] [Indexed: 12/29/2022] Open
Abstract
Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.
Collapse
Affiliation(s)
- Xiaopeng Wang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Guixin Li
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Donglin Ruan
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Zhanwei Zhuang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Rongrong Ding
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
- Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, China
| | - Jianping Quan
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Shiyuan Wang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Yongchuang Jiang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Jinyan Huang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Ting Gu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Linjun Hong
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Enqin Zheng
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Zicong Li
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
| | - Gengyuan Cai
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
- Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, China
| | - Zhenfang Wu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
- Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, China
- Guangdong Provincial Laboratory of Lingnan Modern Agricultural Science and Technology, Guangzhou, China
- *Correspondence: Zhenfang Wu
| | - Jie Yang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, China
- Guangdong Provincial Laboratory of Lingnan Modern Agricultural Science and Technology, Guangzhou, China
- Jie Yang
| |
Collapse
|
7
|
Ablondi M, Sabbioni A, Stocco G, Cipolat-Gotet C, Dadousis C, van Kaam JT, Finocchiaro R, Summer A. Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection. Front Vet Sci 2022; 8:773985. [PMID: 35097040 PMCID: PMC8792952 DOI: 10.3389/fvets.2021.773985] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 11/30/2021] [Indexed: 01/09/2023] Open
Abstract
Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (FROH). The average pedigree inbreeding was 0.07, while the average FROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.
Collapse
Affiliation(s)
- Michela Ablondi
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Alberto Sabbioni
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Giorgia Stocco
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Claudio Cipolat-Gotet
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
- *Correspondence: Claudio Cipolat-Gotet
| | - Christos Dadousis
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Raffaella Finocchiaro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Andrea Summer
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| |
Collapse
|
8
|
Fabbri MC, Dadousis C, Tiezzi F, Maltecca C, Lozada-Soto E, Biffani S, Bozzi R. Genetic diversity and population history of eight Italian beef cattle breeds using measures of autozygosity. PLoS One 2021; 16:e0248087. [PMID: 34695128 PMCID: PMC8544844 DOI: 10.1371/journal.pone.0248087] [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: 02/17/2021] [Accepted: 10/06/2021] [Indexed: 11/18/2022] Open
Abstract
In the present study, GeneSeek GGP-LDv4 33k single nucleotide polymorphism chip was used to detect runs of homozygosity (ROH) in eight Italian beef cattle breeds, six breeds with distribution limited to Tuscany (Calvana, Mucca Pisana, Pontremolese) or Sardinia (Sarda, Sardo Bruna and Sardo Modicana) and two cosmopolitan breeds (Charolais and Limousine). ROH detection analyses were used to estimate autozygosity and inbreeding and to identify genomic regions with high frequency of ROH, which might reflect selection signatures. Comparative analysis among breeds revealed differences in length and distribution of ROH and inbreeding levels. The Charolais, Limousine, Sarda, and Sardo Bruna breeds were found to have a high frequency of short ROH (~ 15.000); Calvana and Mucca Pisana presented also runs longer than 16 Mbp. The highest level of average genomic inbreeding was observed in Tuscan breeds, around 0.3, while Sardinian and cosmopolitan breeds showed values around 0.2. The population structure and genetic distances were analyzed through principal component and multidimensional scaling analyses, and resulted in a clear separation among the breeds, with clusters related to productive purposes. The frequency of ROH occurrence revealed eight breed-specific genomic regions where genes of potential selective and conservative interest are located (e.g. MYOG, CHI3L1, CHIT1 (BTA16), TIMELESS, APOF, OR10P1, OR6C4, OR2AP1, OR6C2, OR6C68, CACNG2 (BTA5), COL5A2 and COL3A1 (BTA2)). In all breeds, we found the largest proportion of homozygous by descent segments to be those that represent inbreeding events that occurred around 32 generations ago, with Tuscan breeds also having a significant proportion of segments relating to more recent inbreeding.
Collapse
Affiliation(s)
- Maria Chiara Fabbri
- Dipartimento di Scienze e Tecnologie Agrarie, Alimentari, Ambientali e Forestali, Università di Firenze, Firenze, Italy
- * E-mail:
| | - Christos Dadousis
- Dipartimento di Scienze Medico‐Veterinarie, Università di Parma, Parma, Italy
| | - Francesco Tiezzi
- Department of Animal Science, North Carolina State University, Raleigh, NC, United States of America
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh, NC, United States of America
| | - Emmanuel Lozada-Soto
- Department of Animal Science, North Carolina State University, Raleigh, NC, United States of America
| | - Stefano Biffani
- Institute of Agricultural Biology and Biotechnology (CNR), Milano, Italy
| | - Riccardo Bozzi
- Dipartimento di Scienze e Tecnologie Agrarie, Alimentari, Ambientali e Forestali, Università di Firenze, Firenze, Italy
| |
Collapse
|