1
|
Gable SM, Bushroe NA, Mendez JM, Wilson A, Pinto BJ, Gamble T, Tollis M. Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles. Genome Biol Evol 2024; 16:evae157. [PMID: 39031594 DOI: 10.1093/gbe/evae157] [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: 02/14/2024] [Revised: 07/09/2024] [Accepted: 07/11/2024] [Indexed: 07/22/2024] Open
Abstract
Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.
Collapse
Affiliation(s)
- Simone M Gable
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA
| | - Nicholas A Bushroe
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA
| | - Jasmine M Mendez
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA
| | - Adam Wilson
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA
| | - Brendan J Pinto
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, USA
| | - Tony Gamble
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, USA
- Department of Biological Sciences, Marquette University, Milwaukee, WI, USA
- Bell Museum of Natural History, University of Minnesota, St. Paul, MN, USA
| | - Marc Tollis
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA
| |
Collapse
|
2
|
Kuhl H, Euclide PT, Klopp C, Cabau C, Zahm M, Lopez-Roques C, Iampietro C, Kuchly C, Donnadieu C, Feron R, Parrinello H, Poncet C, Jaffrelo L, Confolent C, Wen M, Herpin A, Jouanno E, Bestin A, Haffray P, Morvezen R, de Almeida TR, Lecocq T, Schaerlinger B, Chardard D, Żarski D, Larson WA, Postlethwait JH, Timirkhanov S, Kloas W, Wuertz S, Stöck M, Guiguen Y. Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae. BMC Biol 2024; 22:141. [PMID: 38926709 PMCID: PMC11209984 DOI: 10.1186/s12915-024-01935-9] [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: 01/19/2024] [Accepted: 06/06/2024] [Indexed: 06/28/2024] Open
Abstract
BACKGROUND The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.
Collapse
Affiliation(s)
- Heiner Kuhl
- Leibniz-Institute of Freshwater Ecology and Inland Fisheries - IGB (Forschungsverbund Berlin), Müggelseedamm 301/310, D-12587, Berlin, Germany.
| | - Peter T Euclide
- Department of Forestry and Natural Resources | Illinois-Indiana Sea Grant, Purdue University, West Lafayette, USA
| | - Christophe Klopp
- Sigenae, Plateforme Bioinformatique, Genotoul, BioinfoMics, UR875 Biométrie et Intelligence Artificielle, INRAE, Castanet-Tolosan, France
| | - Cédric Cabau
- Sigenae, GenPhySE, Université de Toulouse, INRAE, ENVT, Castanet-Tolosan, France
| | - Margot Zahm
- Sigenae, Plateforme Bioinformatique, Genotoul, BioinfoMics, UR875 Biométrie et Intelligence Artificielle, INRAE, Castanet-Tolosan, France
| | | | | | - Claire Kuchly
- INRAE, US 1426, GeT-PlaGe, Genotoul, Castanet-Tolosan, France
| | | | - Romain Feron
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Hugues Parrinello
- Montpellier GenomiX (MGX), c/o Institut de Génomique Fonctionnelle, 141 rue de la Cardonille, 34094, Montpellier Cedex 05, France
| | - Charles Poncet
- GDEC Gentyane, INRAE, Université Clermont Auvergne, Clermont-Ferrand, France
| | - Lydia Jaffrelo
- GDEC Gentyane, INRAE, Université Clermont Auvergne, Clermont-Ferrand, France
| | - Carole Confolent
- GDEC Gentyane, INRAE, Université Clermont Auvergne, Clermont-Ferrand, France
| | - Ming Wen
- INRAE, LPGP, 35000, Rennes, France
- State Key Laboratory of Developmental Biology of Freshwater Fish, College of Life Science, Hunan Normal University, Changsha, China
| | | | | | - Anastasia Bestin
- SYSAAF, Station INRAE-LPGP, Campus de Beaulieu, 35042, Rennes Cedex, France
| | - Pierrick Haffray
- SYSAAF, Station INRAE-LPGP, Campus de Beaulieu, 35042, Rennes Cedex, France
| | - Romain Morvezen
- SYSAAF, Station INRAE-LPGP, Campus de Beaulieu, 35042, Rennes Cedex, France
| | | | | | | | | | - Daniel Żarski
- Department of Gamete and Embryo Biology, Institute of Animal Reproduction and Food Research, Polish Academy of Sciences, ul. Tuwima 10, 10-748, Olsztyn, Poland
| | - Wesley A Larson
- National Oceanographic and Atmospheric Administration, National Marine Fisheries Service, Alaska Fisheries Science Center, 17109 Point Lena Loop Road, Auke Bay LaboratoriesJuneau, AK, 99801, USA
| | | | | | - Werner Kloas
- Leibniz-Institute of Freshwater Ecology and Inland Fisheries - IGB (Forschungsverbund Berlin), Müggelseedamm 301/310, D-12587, Berlin, Germany
| | - Sven Wuertz
- Leibniz-Institute of Freshwater Ecology and Inland Fisheries - IGB (Forschungsverbund Berlin), Müggelseedamm 301/310, D-12587, Berlin, Germany
| | - Matthias Stöck
- Leibniz-Institute of Freshwater Ecology and Inland Fisheries - IGB (Forschungsverbund Berlin), Müggelseedamm 301/310, D-12587, Berlin, Germany
| | | |
Collapse
|
3
|
Mezzasalma M, Macirella R, Odierna G, Brunelli E. Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles. Genes (Basel) 2024; 15:371. [PMID: 38540430 PMCID: PMC10970613 DOI: 10.3390/genes15030371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Revised: 03/12/2024] [Accepted: 03/15/2024] [Indexed: 06/14/2024] Open
Abstract
Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.
Collapse
Affiliation(s)
- Marcello Mezzasalma
- Department of Biology, Ecology and Earth Science, University of Calabria, Via P. Bucci 4/B, 87036 Rende, Italy; (R.M.); (E.B.)
| | - Rachele Macirella
- Department of Biology, Ecology and Earth Science, University of Calabria, Via P. Bucci 4/B, 87036 Rende, Italy; (R.M.); (E.B.)
| | - Gaetano Odierna
- Independent Researcher, Via Michelangelo 123, 81031 Aversa, Italy;
| | - Elvira Brunelli
- Department of Biology, Ecology and Earth Science, University of Calabria, Via P. Bucci 4/B, 87036 Rende, Italy; (R.M.); (E.B.)
| |
Collapse
|
4
|
Webster TH, Vannan A, Pinto BJ, Denbrock G, Morales M, Dolby GA, Fiddes IT, DeNardo DF, Wilson MA. Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly. Genome Biol Evol 2024; 16:evae018. [PMID: 38319079 PMCID: PMC10950046 DOI: 10.1093/gbe/evae018] [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: 06/15/2023] [Revised: 01/23/2024] [Accepted: 01/23/2024] [Indexed: 02/07/2024] Open
Abstract
Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.
Collapse
Affiliation(s)
- Timothy H Webster
- Department of Anthropology, University of Utah, Salt Lake City, UT, USA
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Annika Vannan
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Brendan J Pinto
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, USA
| | - Grant Denbrock
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Matheo Morales
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Department of Genetics, Yale University, New Haven, CT, USA
| | - Greer A Dolby
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Department of Biology, University of Alabama at Birmingham, Birmingham, AL, USA
| | | | - Dale F DeNardo
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Melissa A Wilson
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- Center for Mechanisms of Evolution, Biodesign Institute, Tempe, AZ, USA
| |
Collapse
|
5
|
Pinto BJ, Nielsen SV, Sullivan KA, Behere A, Keating SE, van Schingen-Khan M, Nguyen TQ, Ziegler T, Pramuk J, Wilson MA, Gamble T. It's a trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex-determining gene in a lizard (Anguimorpha; Shinisauridae). Evolution 2024; 78:355-363. [PMID: 37952174 PMCID: PMC10834058 DOI: 10.1093/evolut/qpad205] [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: 07/07/2023] [Revised: 10/25/2023] [Accepted: 11/08/2023] [Indexed: 11/14/2023]
Abstract
Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally to genetically influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades are not conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 with a novel ZW system on LG3. The linkage group, LG3, corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1 Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2, critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex-determining gene for Shinisaurus.
Collapse
Affiliation(s)
- Brendan J Pinto
- School of Life Sciences, Arizona State University, Tempe, AZ, United States
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, United States
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, United States
| | - Stuart V Nielsen
- Department of Biological Sciences, Museum of Life Sciences, Louisiana State University-Shreveport, Shreveport, LA, United States
- Florida Museum of Natural History, University of Florida, Gainesville, FL, United States
| | - Kathryn A Sullivan
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, United States
- Department of Biological Sciences, Marquette University, Milwaukee, WI, United States
| | - Ashmika Behere
- Department of Biological Sciences, Marquette University, Milwaukee, WI, United States
| | - Shannon E Keating
- Department of Biological Sciences, Marquette University, Milwaukee, WI, United States
| | | | - Truong Q Nguyen
- Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Thomas Ziegler
- Cologne Zoo, Cologne, Germany
- Department of Biology, Institute of Zoology, University of Cologne, Cologne, Germany
| | - Jennifer Pramuk
- Former affiliation: Woodland Park Zoo, Seattle, WA, United States
| | - Melissa A Wilson
- School of Life Sciences, Arizona State University, Tempe, AZ, United States
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, United States
- Center for Mechanisms of Evolution, Biodesign Institute, Tempe, AZ, United States
| | - Tony Gamble
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, United States
- Department of Biological Sciences, Marquette University, Milwaukee, WI, United States
- Bell Museum of Natural History, University of Minnesota, St Paul, MN, United States
| |
Collapse
|
6
|
Pinto BJ, Gamble T, Smith CH, Keating SE, Havird JC, Chiari Y. The revised reference genome of the leopard gecko (Eublepharis macularius) provides insight into the considerations of genome phasing and assembly. J Hered 2023; 114:513-520. [PMID: 36869788 PMCID: PMC10445513 DOI: 10.1093/jhered/esad016] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Accepted: 02/28/2023] [Indexed: 03/05/2023] Open
Abstract
Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000.
Collapse
Affiliation(s)
- Brendan J Pinto
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, USA
| | - Tony Gamble
- Department of Zoology, Milwaukee Public Museum, Milwaukee, WI, USA
- Department of Biological Sciences, Marquette University, Milwaukee WI, USA
- Bell Museum of Natural History, University of Minnesota, St Paul, MN, USA
| | - Chase H Smith
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| | - Shannon E Keating
- Department of Biological Sciences, Marquette University, Milwaukee WI, USA
| | - Justin C Havird
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| | - Ylenia Chiari
- Department of Biology, George Mason University, Fairfax, VA, USA
| |
Collapse
|
7
|
Xie H, Chen Z, Pang S, Du W. Efficient and Highly Continuous Chromosome-Level Genome Assembly of the First Chameleon Genome. Genome Biol Evol 2023; 15:evad131. [PMID: 37481259 PMCID: PMC10410292 DOI: 10.1093/gbe/evad131] [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: 06/06/2022] [Revised: 07/06/2023] [Accepted: 07/10/2023] [Indexed: 07/24/2023] Open
Abstract
Most amniote genomes are diploid and moderate in size (∼1-6 Gbp) and contain a large proportion of repetitive sequences. The development of next-generation sequencing technology, especially the emergence of high-fidelity (HiFi) long-read data, has made it feasible to resolve high-quality genome assembly for nonmodel species efficiently. However, reference genomes for squamate reptiles have lagged behind other amniote lineages. Here, we de novo assembled the first genome from the Chameleonidae family, the panther chameleon (Furcifer pardalis). We obtained telomere-to-telomere contigs using only HiFi data, reaching a contig N50 of 158.72 Mbp. The final chromosome-level assembly is 1.61 Gbp in size, and 100% of primary contigs were placed to pseudochromosomes using Hi-C interaction data. We also found that sequencing depth > 30 folds can ensure both the integrity and accuracy of the genome, whereas insufficient depth led to false increase in genome size and proportion of duplicated genes. We provide a high-quality reference genome valuable for evolutionary and ecological studies in chameleons as well as provide comparative genomic resources for squamate reptiles.
Collapse
Affiliation(s)
- Hongxin Xie
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Zixuan Chen
- College of Biological Sciences, China Agricultural University, Beijing, China
| | | | - Weiguo Du
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|