1
|
Sukhija N, Malik AA, Devadasan JM, Dash A, Bidyalaxmi K, Ravi Kumar D, Kousalaya Devi M, Choudhary A, Kanaka KK, Sharma R, Tripathi SB, Niranjan SK, Sivalingam J, Verma A. Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India. Anim Biotechnol 2024; 35:2290521. [PMID: 38088885 DOI: 10.1080/10495398.2023.2290521] [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: 02/22/2024]
Abstract
The peculiarity of Indian cattle lies in milk quality, resistance to diseases and stressors as well as adaptability. The investigation addressed selection signatures in Gir and Tharparkar cattle, belonging to arid ecotypes of India. Double digest restriction-site associated DNA sequencing (ddRAD-seq) yielded nearly 26 million high-quality reads from unrelated seven Gir and seven Tharparkar cows. In all, 19,127 high-quality SNPs were processed for selection signature analysis. An approach involving within-population composite likelihood ratio (CLR) statistics and between-population FST statistics was used to capture selection signatures within and between the breeds, respectively. A total of 191 selection signatures were addressed using CLR and FST approaches. Selection signatures overlapping 86 and 73 genes were detected as Gir- and Tharparkar-specific, respectively. Notably, genes related to production (CACNA1D, GHRHR), reproduction (ESR1, RBMS3), immunity (NOSTRIN, IL12B) and adaptation (ADAM22, ASL) were annotated to selection signatures. Gene pathway analysis revealed genes in insulin/IGF pathway for milk production, gonadotropin releasing hormone pathway for reproduction, Wnt signalling pathway and chemokine and cytokine signalling pathway for adaptation. This is the first study where selection signatures are identified using ddRAD-seq in indicine cattle breeds. The study shall help in conservation and leveraging genetic improvements in Gir and Tharparkar cattle.
Collapse
Affiliation(s)
- Nidhi Sukhija
- ICAR-National Dairy Research Institute, Karnal, India
| | - Anoop Anand Malik
- TERI School of Advanced Studies, Delhi, India
- The Energy and Resources Institute, North Eastern Regional Centre, Guwahati, India
| | | | | | - Kangabam Bidyalaxmi
- ICAR-National Dairy Research Institute, Karnal, India
- ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - D Ravi Kumar
- ICAR-National Dairy Research Institute, Karnal, India
| | | | | | - K K Kanaka
- ICAR-National Dairy Research Institute, Karnal, India
- ICAR- Indian Institute of Agricultural Biotechnology, Ranchi, India
| | - Rekha Sharma
- ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | | | | | | | - Archana Verma
- ICAR-National Dairy Research Institute, Karnal, India
| |
Collapse
|
2
|
Moore T. X centromeric drive may explain the prevalence of polycystic ovary syndrome and other conditions: Genomic structure of the human X chromosome pericentromeric region is consistent with meiotic drive associated with PCOS and other conditions. Bioessays 2024; 46:e2400056. [PMID: 39072829 DOI: 10.1002/bies.202400056] [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: 03/12/2024] [Revised: 07/02/2024] [Accepted: 07/04/2024] [Indexed: 07/30/2024]
Abstract
X chromosome centromeric drive may explain the prevalence of polycystic ovary syndrome and contribute to oocyte aneuploidy, menopause, and other conditions. The mammalian X chromosome may be vulnerable to meiotic drive because of X inactivation in the female germline. The human X pericentromeric region contains genes potentially involved in meiotic mechanisms, including multiple SPIN1 and ZXDC paralogs. This is consistent with a multigenic drive system comprising differential modification of the active and inactive X chromosome centromeres in female primordial germ cells and preferential segregation of the previously inactivated X chromosome centromere to the polar body at meiosis I. The drive mechanism may explain differences in X chromosome regulation in the female germlines of the human and mouse and, based on the functions encoded by the genes in the region, the transmission of X pericentromeric genetic or epigenetic variants to progeny could contribute to preeclampsia, autism, and differences in sexual differentiation.
Collapse
Affiliation(s)
- Tom Moore
- School of Biochemistry and Cell Biology, University College Cork, Cork, Ireland
| |
Collapse
|
3
|
Whitehouse LS, Ray D, Schrider DR. Tree sequences as a general-purpose tool for population genetic inference. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.20.581288. [PMID: 39185244 PMCID: PMC11343121 DOI: 10.1101/2024.02.20.581288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/27/2024]
Abstract
As population genetics data increases in size new methods have been developed to store genetic information in efficient ways, such as tree sequences. These data structures are computationally and storage efficient, but are not interchangeable with existing data structures used for many population genetic inference methodologies such as the use of convolutional neural networks (CNNs) applied to population genetic alignments. To better utilize these new data structures we propose and implement a graph convolutional network (GCN) to directly learn from tree sequence topology and node data, allowing for the use of neural network applications without an intermediate step of converting tree sequences to population genetic alignment format. We then compare our approach to standard CNN approaches on a set of previously defined benchmarking tasks including recombination rate estimation, positive selection detection, introgression detection, and demographic model parameter inference. We show that tree sequences can be directly learned from using a GCN approach and can be used to perform well on these common population genetics inference tasks with accuracies roughly matching or even exceeding that of a CNN-based method. As tree sequences become more widely used in population genetics research we foresee developments and optimizations of this work to provide a foundation for population genetics inference moving forward.
Collapse
Affiliation(s)
- Logan S. Whitehouse
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina, USA
| | - Dylan Ray
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina, USA
| | - Daniel R. Schrider
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina, USA
| |
Collapse
|
4
|
Götsch H, Bürger R. Polygenic dynamics underlying the response of quantitative traits to directional selection. Theor Popul Biol 2024; 158:21-59. [PMID: 38677378 DOI: 10.1016/j.tpb.2024.04.006] [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: 02/22/2023] [Revised: 04/14/2024] [Accepted: 04/19/2024] [Indexed: 04/29/2024]
Abstract
We study the response of a quantitative trait to exponential directional selection in a finite haploid population, both at the genetic and the phenotypic level. We assume an infinite sites model, in which the number of new mutations per generation in the population follows a Poisson distribution (with mean Θ) and each mutation occurs at a new, previously monomorphic site. Mutation effects are beneficial and drawn from a distribution. Sites are unlinked and contribute additively to the trait. Assuming that selection is stronger than random genetic drift, we model the initial phase of the dynamics by a supercritical Galton-Watson process. This enables us to obtain time-dependent results. We show that the copy-number distribution of the mutant in generation n, conditioned on non-extinction until n, is described accurately by the deterministic increase from an initial distribution with mean 1. This distribution is related to the absolutely continuous part W+ of the random variable, typically denoted W, that characterizes the stochasticity accumulating during the mutant's sweep. A suitable transformation yields the approximate dynamics of the mutant frequency distribution in a Wright-Fisher population of size N. Our expression provides a very accurate approximation except when mutant frequencies are close to 1. On this basis, we derive explicitly the (approximate) time dependence of the expected mean and variance of the trait and of the expected number of segregating sites. Unexpectedly, we obtain highly accurate approximations for all times, even for the quasi-stationary phase when the expected per-generation response and the trait variance have equilibrated. The latter refine classical results. In addition, we find that Θ is the main determinant of the pattern of adaptation at the genetic level, i.e., whether the initial allele-frequency dynamics are best described by sweep-like patterns at few loci or small allele-frequency shifts at many. The number of segregating sites is an appropriate indicator for these patterns. The selection strength determines primarily the rate of adaptation. The accuracy of our results is tested by comprehensive simulations in a Wright-Fisher framework. We argue that our results apply to more complex forms of directional selection.
Collapse
Affiliation(s)
- Hannah Götsch
- Faculty of Mathematics, University of Vienna, 1090 Vienna, Austria; Vienna Graduate School of Population Genetics, Austria.
| | - Reinhard Bürger
- Faculty of Mathematics, University of Vienna, 1090 Vienna, Austria
| |
Collapse
|
5
|
Samano A, Kumar N, Liao Y, Ishtiaq F, Chakraborty M. Genome structural variants shape adaptive success of an invasive urban malaria vector Anopheles stephensi. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.29.605641. [PMID: 39211149 PMCID: PMC11360885 DOI: 10.1101/2024.07.29.605641] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Global changes are associated with the emergence of several invasive species. However, the genomic determinants of the adaptive success of an invasive species in a new environment remain poorly understood. Genomic structural variants (SVs), consisting of copy number variants, play an important role in adaptation. SVs often cause large adaptive shifts in ecologically important traits, which makes SVs compelling candidates for driving rapid adaptations to environmental changes, which is critical to invasive success. To address this problem, we investigated the role SVs play in the adaptive success of Anopheles stephensi , a primary vector of urban malaria in South Asia and an invasive malaria vector in several South Asian islands and Africa. We collected whole genome sequencing data from 115 mosquitoes from invasive island populations and four locations from mainland India, an ancestral range for the species. We identified 2,988 duplication copy number variants and 16,038 deletions in these strains, with ∼50% overlapping genes. SVs are enriched in genomic regions with signatures of selective sweeps in the mainland and invasive island populations, implying a putative adaptive role of SVs. Nearly all high-frequency SVs, including the candidate adaptive variants, in the invasive island populations are present on the mainland, suggesting a major contribution of existing variation to the success of the island populations. Among the candidate adaptive SVs, three duplications involving toxin-resistance genes evolved, likely due to the widespread application of insecticides in India since the 1950s. We also identify two SVs associated with the adaptation of An. stephensi larvae to brackish water in the island and two coastal mainland populations, where the mutations likely originated. Our results suggest that existing SVs play a vital role in the evolutionary success of An. stephensi in new environmental conditions.
Collapse
|
6
|
Poulin R, Salloum PM, Bennett J. Evolution of parasites in the Anthropocene: new pressures, new adaptive directions. Biol Rev Camb Philos Soc 2024. [PMID: 38984760 DOI: 10.1111/brv.13118] [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: 01/23/2024] [Revised: 06/27/2024] [Accepted: 07/02/2024] [Indexed: 07/11/2024]
Abstract
The Anthropocene is seeing the human footprint rapidly spreading to all of Earth's ecosystems. The fast-changing biotic and abiotic conditions experienced by all organisms are exerting new and strong selective pressures, and there is a growing list of examples of human-induced evolution in response to anthropogenic impacts. No organism is exempt from these novel selective pressures. Here, we synthesise current knowledge on human-induced evolution in eukaryotic parasites of animals, and present a multidisciplinary framework for its study and monitoring. Parasites generally have short generation times and huge fecundity, features that predispose them for rapid evolution. We begin by reviewing evidence that parasites often have substantial standing genetic variation, and examples of their rapid evolution both under conditions of livestock production and in serial passage experiments. We then present a two-step conceptual overview of the causal chain linking anthropogenic impacts to parasite evolution. First, we review the major anthropogenic factors impacting parasites, and identify the selective pressures they exert on parasites through increased mortality of either infective stages or adult parasites, or through changes in host density, quality or immunity. Second, we discuss what new phenotypic traits are likely to be favoured by the new selective pressures resulting from altered parasite mortality or host changes; we focus mostly on parasite virulence and basic life-history traits, as these most directly influence the transmission success of parasites and the pathology they induce. To illustrate the kinds of evolutionary changes in parasites anticipated in the Anthropocene, we present a few scenarios, either already documented or hypothetical but plausible, involving parasite taxa in livestock, aquaculture and natural systems. Finally, we offer several approaches for investigations and real-time monitoring of rapid, human-induced evolution in parasites, ranging from controlled experiments to the use of state-of-the-art genomic tools. The implications of fast-evolving parasites in the Anthropocene for disease emergence and the dynamics of infections in domestic animals and wildlife are concerning. Broader recognition that it is not only the conditions for parasite transmission that are changing, but the parasites themselves, is needed to meet better the challenges ahead.
Collapse
Affiliation(s)
- Robert Poulin
- Department of Zoology, University of Otago, P.O. Box 56, Dunedin, New Zealand
| | - Priscila M Salloum
- Department of Zoology, University of Otago, P.O. Box 56, Dunedin, New Zealand
| | - Jerusha Bennett
- Department of Zoology, University of Otago, P.O. Box 56, Dunedin, New Zealand
| |
Collapse
|
7
|
Ndiaye YD, Wong W, Thwing J, Schaffner SF, Brenneman KV, Tine A, Diallo MA, Deme AB, Sy M, Bei AK, Thiaw AB, Daniels R, Ndiaye T, Gaye A, Ndiaye IM, Toure M, Gadiaga N, Sene A, Sow D, Garba MN, Yade MS, Dieye B, Diongue K, Zoumarou D, Ndiaye A, Gomis JF, Fall FB, Ndiop M, Diallo I, Sene D, Macinnis B, Seck MC, Ndiaye M, Ngom B, Diedhiou Y, Mbaye AM, Ndiaye L, Sy N, Badiane AS, Hartl DL, Wirth DF, Volkman SK, Ndiaye D. Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time. Malar J 2024; 23:205. [PMID: 38982475 PMCID: PMC11234717 DOI: 10.1186/s12936-024-05024-8] [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: 04/17/2024] [Accepted: 06/25/2024] [Indexed: 07/11/2024] Open
Abstract
BACKGROUND Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance. METHODS Data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics. RESULTS This dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014. DISCUSSION (CONCLUSION) The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.
Collapse
Affiliation(s)
- Yaye D Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Wesley Wong
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
| | - Julie Thwing
- Malaria Branch, Division of Parasitic Diseases and Malaria, Global Health Center, Centers for Disease Control and Prevention, Atlanta, GA, USA
| | - Stephen F Schaffner
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Katelyn Vendrely Brenneman
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
| | - Abdoulaye Tine
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mamadou A Diallo
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Awa B Deme
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mouhamad Sy
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Amy K Bei
- Yale School of Public Health, 60 College St, New Haven, CT, 06510, USA
| | - Alphonse B Thiaw
- Department of Biochemistry and Functional Genomics, Sherbrooke University, 2500 Bd de L'Universite, Sherbrooke, QC, J1K 2R1, Canada
| | - Rachel Daniels
- RNA Therapeutics Institute, UMass Chan Medical School, 368 Plantation Street, Worcester, MA, 01605, USA
| | - Tolla Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Amy Gaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Ibrahima M Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mariama Toure
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Nogaye Gadiaga
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Aita Sene
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Djiby Sow
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mamane N Garba
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mamadou S Yade
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Baba Dieye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Khadim Diongue
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Daba Zoumarou
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Aliou Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Jules F Gomis
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Fatou B Fall
- National Malaria Control Programme (NMCP), 25270, Dakar, Senegal
| | - Medoune Ndiop
- National Malaria Control Programme (NMCP), 25270, Dakar, Senegal
| | - Ibrahima Diallo
- National Malaria Control Programme (NMCP), 25270, Dakar, Senegal
| | - Doudou Sene
- National Malaria Control Programme (NMCP), 25270, Dakar, Senegal
| | - Bronwyn Macinnis
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Mame C Seck
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Mouhamadou Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Bassirou Ngom
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Younouss Diedhiou
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Amadou M Mbaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Lamine Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Ngayo Sy
- Service de Lutte Antiparasitaire (SLAP), Thiès, Senegal
| | - Aida S Badiane
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
| | - Daniel L Hartl
- Department of Organismic and Evolutionary Biology, Harvard University, 16 Divinity Avenue, Cambridge, MA, 02138, USA
| | - Dyann F Wirth
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Sarah K Volkman
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA.
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA.
- Simmons University, 300 The Fenway, Boston, MA, 02115, USA.
| | - Daouda Ndiaye
- International Research Training Center On Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, 16477, Dakar, Senegal
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
| |
Collapse
|
8
|
Li T, Wang Y, Zhang Z, Ji C, Zheng N, Huang Y. A comparative analysis reveals the genomic diversity among 8 Muscovy duck populations. G3 (BETHESDA, MD.) 2024; 14:jkae112. [PMID: 38789099 PMCID: PMC11228869 DOI: 10.1093/g3journal/jkae112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Revised: 03/05/2024] [Accepted: 05/10/2024] [Indexed: 05/26/2024]
Abstract
The Muscovy duck (Cairina moschata) is a waterfowl indigenous to the neotropical regions of Central and South America. It has low demand for concentrated feed and strong adaptability to different rearing conditions. After introduced to China through Eurasian commercial trade, Muscovy ducks have a domestication history of around 300 years in the Fujian Province of China. In the 1990s, the commodity Muscovy duck breed "Crimo," cultivated in Europe, entered the Chinese market for consumption and breeding purposes. Due to the different selective breeding processes, Muscovy ducks have various populational traits and lack transparency of their genetic background. To remove this burden in the Muscovy duck breeding process, we analyzed genomic data from 8 populations totaling 83 individuals. We identify 11.24 million single nucleotide polymorphisms (SNPs) and categorized these individuals into the Fujian-bred and the Crimo populations according to phylogenetic analyses. We then delved deeper into their evolutionary relationships through assessing population structure, calculating fixation index (FST) values, and measuring genetic distances. Our exploration of runs of homozygosity (ROHs) and homozygous-by-descent (HBD) uncovered genomic regions enriched for genes implicated in fatty acid metabolism, development, and immunity pathways. Selective sweep analyses further indicated strong selective pressures exerted on genes including TECR, STAT2, and TRAF5. These findings provide insights into genetic variations of Muscovy ducks, thus offering valuable information regarding genetic diversity, population conservation, and genome associated with the breeding of Muscovy ducks.
Collapse
Affiliation(s)
- Te Li
- State Key Laboratory of Farm Animal Biotech Breeding, College of Biology Sciences, China Agricultural University, No.2 Yuan Ming Yuan West Road, Hai Dian District, Beijing 100193, China
| | - Yiming Wang
- State Key Laboratory of Farm Animal Biotech Breeding, College of Biology Sciences, China Agricultural University, No.2 Yuan Ming Yuan West Road, Hai Dian District, Beijing 100193, China
| | - Zhou Zhang
- National Key Laboratory for Swine Genetic Improvement and Production Technology, Jiangxi Agricultural University, Nanchang 330045, China
| | - Congliang Ji
- Technology Department (Research Institute) Livestock and Poultry Breeding Research Office, Wens Foodstuff Group Co. Ltd, Huineng North Road, Xincheng Town, Xinxing County, Yunfu City, Guangdong Province 527400, China
| | - Nengzhu Zheng
- Institute of Animal Husbandry and Veterinary Medicine, Fujian Academy of Agricultural Sciences, Fuzhou 350013, China
| | - Yinhua Huang
- State Key Laboratory of Farm Animal Biotech Breeding, College of Biology Sciences, China Agricultural University, No.2 Yuan Ming Yuan West Road, Hai Dian District, Beijing 100193, China
| |
Collapse
|
9
|
She H, Liu Z, Xu Z, Zhang H, Wu J, Wang X, Cheng F, Charlesworth D, Qian W. Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra. THE NEW PHYTOLOGIST 2024; 243:477-494. [PMID: 38715078 DOI: 10.1111/nph.19799] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Accepted: 04/18/2024] [Indexed: 06/07/2024]
Abstract
Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.
Collapse
Affiliation(s)
- Hongbing She
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Zhiyuan Liu
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Zhaosheng Xu
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Helong Zhang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Jian Wu
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Xiaowu Wang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Feng Cheng
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Deborah Charlesworth
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Charlotte Auerbach Road, Edinburgh, EH9 3FL, UK
| | - Wei Qian
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| |
Collapse
|
10
|
Zhang X, Blaxter M, Wood JMD, Tracey A, McCarthy S, Thorpe P, Rayner JG, Zhang S, Sikkink KL, Balenger SL, Bailey NW. Temporal genomics in Hawaiian crickets reveals compensatory intragenomic coadaptation during adaptive evolution. Nat Commun 2024; 15:5001. [PMID: 38866741 PMCID: PMC11169259 DOI: 10.1038/s41467-024-49344-4] [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/20/2023] [Accepted: 05/24/2024] [Indexed: 06/14/2024] Open
Abstract
Theory predicts that compensatory genetic changes reduce negative indirect effects of selected variants during adaptive evolution, but evidence is scarce. Here, we test this in a wild population of Hawaiian crickets using temporal genomics and a high-quality chromosome-level cricket genome. In this population, a mutation, flatwing, silences males and rapidly spread due to an acoustically-orienting parasitoid. Our sampling spanned a social transition during which flatwing fixed and the population went silent. We find long-range linkage disequilibrium around the putative flatwing locus was maintained over time, and hitchhiking genes had functions related to negative flatwing-associated effects. We develop a combinatorial enrichment approach using transcriptome data to test for compensatory, intragenomic coevolution. Temporal changes in genomic selection were distributed genome-wide and functionally associated with the population's transition to silence, particularly behavioural responses to silent environments. Our results demonstrate how 'adaptation begets adaptation'; changes to the sociogenetic environment accompanying rapid trait evolution can generate selection provoking further, compensatory adaptation.
Collapse
Affiliation(s)
- Xiao Zhang
- Tianjin Key Laboratory of Conservation and Utilization of Animal Diversity, College of Life Sciences, Tianjin Normal University, Tianjin, China.
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, Fife, UK.
| | - Mark Blaxter
- Tree of Life, Wellcome Sanger Institute, Cambridge, UK
| | | | - Alan Tracey
- Tree of Life, Wellcome Sanger Institute, Cambridge, UK
| | | | - Peter Thorpe
- School of Medicine, University of St Andrews, St Andrews, Fife, UK
- Data Analysis Group, Division of Computational Biology, School of Life Sciences, University of Dundee, Dundee, UK
| | - Jack G Rayner
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, Fife, UK
| | - Shangzhe Zhang
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, Fife, UK
| | | | - Susan L Balenger
- College of Biological Sciences, University of Minnesota, Saint Paul, MN, USA
| | - Nathan W Bailey
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, Fife, UK.
| |
Collapse
|
11
|
Zhang X, Yang L, Xie Z, Gan J, Zhu P, Song J, Kang H, Zhang Z, Liu L, Xiang H, Li H. Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity. BMC Genomics 2024; 25:577. [PMID: 38858651 PMCID: PMC11163754 DOI: 10.1186/s12864-024-10492-y] [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: 10/17/2023] [Accepted: 06/03/2024] [Indexed: 06/12/2024] Open
Abstract
BACKGROUND Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH. RESULTS The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands. CONCLUSIONS These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.
Collapse
Affiliation(s)
- Xing Zhang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Liu Yang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
- Key Laboratory of National Forestry and Grassland Administration on Conservation Biology of Rare Animals in the Giant, China Conservation and Research Centre for the Giant Panda, Panda National Park, Chengdu, 611830, China
| | - Zhuojun Xie
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Jiankang Gan
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China
| | - Piao Zhu
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Jiani Song
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Huimin Kang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Zhengfen Zhang
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China
| | - Lingbin Liu
- College of Animal Science and Technology, Southwest University, Chongqing, 400715, China
| | - Hai Xiang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China.
| | - Hua Li
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China.
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China.
| |
Collapse
|
12
|
Soni V, Jensen JD. Temporal challenges in detecting balancing selection from population genomic data. G3 (BETHESDA, MD.) 2024; 14:jkae069. [PMID: 38551137 DOI: 10.1093/g3journal/jkae069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 12/21/2023] [Accepted: 03/19/2024] [Indexed: 04/28/2024]
Abstract
The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ 85281, USA
| | - Jeffrey D Jensen
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ 85281, USA
| |
Collapse
|
13
|
Yu H, Zhang K, Cheng G, Mei C, Wang H, Zan L. Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle. BMC Genomics 2024; 25:558. [PMID: 38834950 DOI: 10.1186/s12864-024-10482-0] [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/26/2023] [Accepted: 05/30/2024] [Indexed: 06/06/2024] Open
Abstract
BACKGROUND Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19). RESULTS A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima's D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle. CONCLUSIONS We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.
Collapse
Affiliation(s)
- Hengwei Yu
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, 712100, China
| | - Ke Zhang
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, 712100, China
| | - Gong Cheng
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, 712100, China
| | - Chugang Mei
- College of Grassland Agriculture, Northwest A&F University, No.22 Xinong Road, Yangling, 712100, China
- National Beef Cattle Improvement Center, Yangling, 712100, China
| | - Hongbao Wang
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, 712100, China
| | - Linsen Zan
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, 712100, China.
- National Beef Cattle Improvement Center, Yangling, 712100, China.
| |
Collapse
|
14
|
Yang X, Su Y, Huang S, Hou Q, Wei P, Hao Y, Huang J, Xiao H, Ma Z, Xu X, Wang X, Cao S, Cao X, Zhang M, Wen X, Ma Y, Peng Y, Zhou Y, Cao K, Qiao G. Comparative population genomics reveals convergent and divergent selection in the apricot-peach-plum-mei complex. HORTICULTURE RESEARCH 2024; 11:uhae109. [PMID: 38883333 PMCID: PMC11179850 DOI: 10.1093/hr/uhae109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 04/06/2024] [Indexed: 06/18/2024]
Abstract
The economically significant genus Prunus includes fruit and nut crops that have been domesticated for shared and specific agronomic traits; however, the genomic signals of convergent and divergent selection have not been elucidated. In this study, we aimed to detect genomic signatures of convergent and divergent selection by conducting comparative population genomic analyses of the apricot-peach-plum-mei (APPM) complex, utilizing a haplotype-resolved telomere-to-telomere (T2T) genome assembly and population resequencing data. The haplotype-resolved T2T reference genome for the plum cultivar was assembled through HiFi and Hi-C reads, resulting in two haplotypes 251.25 and 251.29 Mb in size, respectively. Comparative genomics reveals a chromosomal translocation of ~1.17 Mb in the apricot genomes compared with peach, plum, and mei. Notably, the translocation involves the D locus, significantly impacting titratable acidity (TA), pH, and sugar content. Population genetic analysis detected substantial gene flow between plum and apricot, with introgression regions enriched in post-embryonic development and pollen germination processes. Comparative population genetic analyses revealed convergent selection for stress tolerance, flower development, and fruit ripening, along with divergent selection shaping specific crop, such as somatic embryogenesis in plum, pollen germination in mei, and hormone regulation in peach. Notably, selective sweeps on chromosome 7 coincide with a chromosomal collinearity from the comparative genomics, impacting key fruit-softening genes such as PG, regulated by ERF and RMA1H1. Overall, this study provides insights into the genetic diversity, evolutionary history, and domestication of the APPM complex, offering valuable implications for genetic studies and breeding programs of Prunus crops.
Collapse
Affiliation(s)
- Xuanwen Yang
- Key Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), Institute of Agro-bioengineering/College of Life Sciences, Guizhou University, Guiyang 550025, China
- National Key Laboratory for Germplasm Innovation & Utilization of Horticultural Crops, Zhengzhou Fruit Research Institute, Chinese Academy of Agricultural Science, Zhengzhou 450009, China
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- College of Horticulture and Forestry Sciences, Huazhong Agricultural University, Wuhan 430070, China
| | - Ying Su
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Xinjiang, Urumqi 830046, China
| | - Siyang Huang
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Qiandong Hou
- Key Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), Institute of Agro-bioengineering/College of Life Sciences, Guizhou University, Guiyang 550025, China
| | - Pengcheng Wei
- National Key Laboratory for Germplasm Innovation & Utilization of Horticultural Crops, Zhengzhou Fruit Research Institute, Chinese Academy of Agricultural Science, Zhengzhou 450009, China
| | - Yani Hao
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- Department of Bioinformatics, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou 215123, China
| | - Jiaqi Huang
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- College of Life Sciences, Wuhan University, Wuhan 430072, China
| | - Hua Xiao
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Zhiyao Ma
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Xiaodong Xu
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Xu Wang
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Shuo Cao
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- College of Horticulture and Forestry Sciences, Huazhong Agricultural University, Wuhan 430070, China
| | - Xuejing Cao
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Mengyan Zhang
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Xiaopeng Wen
- Key Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), Institute of Agro-bioengineering/College of Life Sciences, Guizhou University, Guiyang 550025, China
| | - Yuhua Ma
- Institute of Pomology Science, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China
| | - Yanling Peng
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Yongfeng Zhou
- National Key Laboratory of Tropical Crop Breeding, Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
- National Key Laboratory of Tropical Crop Breeding, Tropical Crops Genetic Resources Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 570100, China
| | - Ke Cao
- National Key Laboratory for Germplasm Innovation & Utilization of Horticultural Crops, Zhengzhou Fruit Research Institute, Chinese Academy of Agricultural Science, Zhengzhou 450009, China
| | - Guang Qiao
- Key Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), Institute of Agro-bioengineering/College of Life Sciences, Guizhou University, Guiyang 550025, China
| |
Collapse
|
15
|
Fonseca PAS, Suárez-Vega A, Arranz JJ, Gutiérrez-Gil B. Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits. Genet Sel Evol 2024; 56:40. [PMID: 38773423 PMCID: PMC11106937 DOI: 10.1186/s12711-024-00910-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: 12/08/2023] [Accepted: 05/07/2024] [Indexed: 05/23/2024] Open
Abstract
BACKGROUND Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS). RESULTS In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits. CONCLUSIONS The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.
Collapse
Affiliation(s)
- Pablo A S Fonseca
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana S/N, 24071, León, Spain
| | - Aroa Suárez-Vega
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana S/N, 24071, León, Spain
| | - Juan J Arranz
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana S/N, 24071, León, Spain
| | - Beatriz Gutiérrez-Gil
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana S/N, 24071, León, Spain.
| |
Collapse
|
16
|
Korfmann K, Temple-Boyer M, Sellinger T, Tellier A. Determinants of rapid adaptation in species with large variance in offspring production. Mol Ecol 2024; 33:e16982. [PMID: 37199145 DOI: 10.1111/mec.16982] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Revised: 04/26/2023] [Accepted: 05/02/2023] [Indexed: 05/19/2023]
Abstract
The speed of population adaptation to changing biotic and abiotic environments is determined by the interaction between genetic drift, positive selection and linkage effects. Many marine species (fish, crustaceans), invertebrates and pathogens of humans and crops, exhibit sweepstakes reproduction characterized by the production of a very large amount of offspring (fecundity phase) from which only a small fraction may survive to the next generation (viability phase). Using stochastic simulations, we investigate whether the occurrence of sweepstakes reproduction affects the efficiency of a positively selected unlinked locus, and thus, the speed of adaptation since fecundity and/or viability have distinguishable consequences on mutation rate, probability and fixation time of advantageous alleles. We observe that the mean number of mutations at the next generation is always the function of the population size, but the variance increases with stronger sweepstakes reproduction when mutations occur in the parents. On the one hand, stronger sweepstakes reproduction magnifies the effect of genetic drift thus increasing the probability of fixation of neutral allele and decreasing that of selected alleles. On the other hand, the time to fixation of advantageous (as well as neutral) alleles is shortened by stronger sweepstakes reproduction. Importantly, fecundity and viability selection exhibit different probabilities and times to fixation of advantageous alleles under intermediate and weak sweepstakes reproduction. Finally, alleles under both strong fecundity and viability selection display a synergistic efficiency of selection. We conclude that measuring and modelling accurately fecundity and/or viability selection are crucial to predict the adaptive potential of species with sweepstakes reproduction.
Collapse
Affiliation(s)
- Kevin Korfmann
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| | - Marie Temple-Boyer
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| | - Thibaut Sellinger
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
- Department of Environment and Biodiversity, Paris Lodron University of Salzburg, Salzburg, Austria
| | - Aurélien Tellier
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| |
Collapse
|
17
|
Le Pennec G, Retel C, Kowallik V, Becks L, Feulner PGD. Demographic fluctuations and selection during host-parasite co-evolution interactively increase genetic diversity. Mol Ecol 2024; 33:e16939. [PMID: 36997280 DOI: 10.1111/mec.16939] [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: 10/28/2022] [Revised: 03/17/2023] [Accepted: 03/22/2023] [Indexed: 04/01/2023]
Abstract
Host-parasite interactions can cause strong demographic fluctuations accompanied by selective sweeps of resistance/infectivity alleles. Both demographic bottlenecks and frequent sweeps are expected to reduce the amount of segregating genetic variation and therefore might constrain adaptation during co-evolution. Recent studies, however, suggest that the interaction of demographic and selective processes is a key component of co-evolutionary dynamics and may rather positively affect levels of genetic diversity available for adaptation. Here, we provide direct experimental testing of this hypothesis by disentangling the effects of demography, selection and their interaction in an experimental host-parasite system. We grew 12 populations of a unicellular, asexually reproducing algae (Chlorella variabilis) that experienced either growth followed by constant population sizes (three populations), demographic fluctuations (three populations), selection induced by exposure to a virus (three populations), or demographic fluctuations together with virus-induced selection (three populations). After 50 days (~50 generations), we conducted whole-genome sequencing of each algal host population. We observed more genetic diversity in populations that jointly experienced selection and demographic fluctuations than in populations where these processes were experimentally separated. In addition, in those three populations that jointly experienced selection and demographic fluctuations, experimentally measured diversity exceeds expected values of diversity that account for the cultures' population sizes. Our results suggest that eco-evolutionary feedbacks can positively affect genetic diversity and provide the necessary empirical measures to guide further improvements of theoretical models of adaptation during host-parasite co-evolution.
Collapse
Affiliation(s)
- Guénolé Le Pennec
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Cas Retel
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Vienna Kowallik
- Community Dynamics Group, Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
- Albert-Ludwigs University Freiburg, Faculty of Environment and Natural Resources, Professorship of Forest Entomology and Protection, Stegen-Wittental, Germany
| | - Lutz Becks
- Community Dynamics Group, Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
- Aquatic Ecology and Evolution, Limnological Institute University of Konstanz, Konstanz, Germany
| | - Philine G D Feulner
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
- Division of Aquatic Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| |
Collapse
|
18
|
Eldon B, Stephan W. Sweepstakes reproduction facilitates rapid adaptation in highly fecund populations. Mol Ecol 2024; 33:e16903. [PMID: 36896794 DOI: 10.1111/mec.16903] [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: 05/23/2022] [Revised: 02/21/2023] [Accepted: 02/23/2023] [Indexed: 03/11/2023]
Abstract
Adaptation enables natural populations to survive in a changing environment. Understanding the mechanics of adaptation is therefore crucial for learning about the evolution and ecology of natural populations. We focus on the impact of random sweepstakes on selection in highly fecund haploid and diploid populations partitioned into two genetic types, with one type conferring selective advantage. For the diploid populations, we incorporate various dominance mechanisms. We assume that the populations may experience recurrent bottlenecks. In random sweepstakes, the distribution of individual recruitment success is highly skewed, resulting in a huge variance in the number of offspring contributed by the individuals present in any given generation. Using computer simulations, we investigate the joint effects of random sweepstakes, recurrent bottlenecks and dominance mechanisms on selection. In our framework, bottlenecks allow random sweepstakes to have an effect on the time to fixation, and in diploid populations, the effect of random sweepstakes depends on the dominance mechanism. We describe selective sweepstakes that are approximated by recurrent sweeps of strongly beneficial allelic types arising by mutation. We demonstrate that both types of sweepstakes reproduction may facilitate rapid adaptation (as defined based on the average time to fixation of a type conferring selective advantage conditioned on fixation of the type). However, whether random sweepstakes cause rapid adaptation depends also on their interactions with bottlenecks and dominance mechanisms. Finally, we review a case study in which a model of recurrent sweeps is shown to essentially explain population genomic data from Atlantic cod.
Collapse
Affiliation(s)
- Bjarki Eldon
- Institute of Evolution and Biodiversity Science, Natural History Museum Berlin, Berlin, Germany
| | | |
Collapse
|
19
|
Otalora K, Riera JL, Tavecchia G, Rotger A, Igual JM, Trotta JP, Baldo L. Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data. Ecol Evol 2024; 14:e11407. [PMID: 38799398 PMCID: PMC11116764 DOI: 10.1002/ece3.11407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 04/12/2024] [Accepted: 04/30/2024] [Indexed: 05/29/2024] Open
Abstract
Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.
Collapse
Affiliation(s)
- Katherin Otalora
- Department of Evolutionary Biology, Ecology and Environmental Sciences, Biology FacultyUniversity of Barcelona (UB)BarcelonaSpain
- Fundación FUNMAJO, EBA, RAIEC, Biodiversity BranchTunjaBoyacáColombia
| | - Joan Lluís Riera
- Department of Evolutionary Biology, Ecology and Environmental Sciences, Biology FacultyUniversity of Barcelona (UB)BarcelonaSpain
| | - Giacomo Tavecchia
- Animal Demography and Ecology Unit (GEDA‐IMEDEA, CSIC‐UIB)EsporlesSpain
| | - Andreu Rotger
- Animal Demography and Ecology Unit (GEDA‐IMEDEA, CSIC‐UIB)EsporlesSpain
| | - José Manuel Igual
- Animal Demography and Ecology Unit (GEDA‐IMEDEA, CSIC‐UIB)EsporlesSpain
| | - Jean‐Remi Paul Trotta
- CNAG‐CRG, Centre for Genomic Regulation (CRG)Barcelona Institute of Science and Technology (BIST)BarcelonaSpain
| | - Laura Baldo
- Department of Evolutionary Biology, Ecology and Environmental Sciences, Biology FacultyUniversity of Barcelona (UB)BarcelonaSpain
- Institute for Research on Biodiversity (IRBio)University of Barcelona (UB)BarcelonaSpain
| |
Collapse
|
20
|
Soni V, Terbot JW, Jensen JD. Population genetic considerations regarding the interpretation of within-patient SARS-CoV-2 polymorphism data. Nat Commun 2024; 15:3240. [PMID: 38627371 PMCID: PMC11021480 DOI: 10.1038/s41467-024-46261-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Accepted: 01/29/2024] [Indexed: 04/19/2024] Open
Affiliation(s)
- Vivak Soni
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA
| | - John W Terbot
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA
- Division of Biological Sciences, University of Montana, Missoula, MT, USA
| | - Jeffrey D Jensen
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA.
| |
Collapse
|
21
|
Rajawat D, Ghildiyal K, Sonejita Nayak S, Sharma A, Parida S, Kumar S, Ghosh AK, Singh U, Sivalingam J, Bhushan B, Dutt T, Panigrahi M. Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle. Gene 2024; 901:148178. [PMID: 38242377 DOI: 10.1016/j.gene.2024.148178] [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: 09/26/2023] [Revised: 01/10/2024] [Accepted: 01/16/2024] [Indexed: 01/21/2024]
Abstract
The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.
Collapse
Affiliation(s)
- Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Kanika Ghildiyal
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Anurodh Sharma
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Subhashree Parida
- Pharmacology & Toxicology Division, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Shive Kumar
- Department of Animal Genetics and Breeding, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, India
| | - A K Ghosh
- Department of Animal Genetics and Breeding, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, India
| | - Umesh Singh
- ICAR Central Institute for Research on Cattle, Meerut, UP, India
| | | | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India.
| |
Collapse
|
22
|
Rajawat D, Panigrahi M, Nayak SS, Bhushan B, Mishra BP, Dutt T. Dissecting the genomic regions of selection on the X chromosome in different cattle breeds. 3 Biotech 2024; 14:50. [PMID: 38268984 PMCID: PMC10803714 DOI: 10.1007/s13205-023-03905-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Accepted: 12/18/2023] [Indexed: 01/26/2024] Open
Abstract
Mammalian X and Y chromosomes independently evolved from various autosomes approximately 300 million years ago (MYA). To fully understand the relationship between genomic composition and phenotypic diversity arising due to the course of evolution, we have scanned regions of selection signatures on the X chromosome in different cattle breeds. In this study, we have prepared the datasets of 184 individuals of different cattle breeds and explored the complete X chromosome by utilizing four within-population and two between-population methods. There were 23, 25, 30, 17, 17, and 12 outlier regions identified in Tajima's D, CLR, iHS, ROH, FST, and XP-EHH. Bioinformatics analysis showed that these regions harbor important candidate genes like AKAP4 for reproduction in Brown Swiss, MBTS2 for production traits in Brown Swiss and Guernsey, CXCR3 and CITED1 for health traits in Jersey and Nelore, and BMX and CD40LG for regulation of X chromosome inactivation in Nelore and Gir. We identified genes shared among multiple methods, such as TRNAC-GCA and IL1RAPL1, which appeared in Tajima's D, ROH, and iHS analyses. The gene TRNAW-CCA was found in ROH, CLR and iHS analyses. The X chromosome exhibits a distinctive interaction between demographic factors and genetic variations, and these findings may provide new insight into the X-linked selection in different cattle breeds.
Collapse
Affiliation(s)
- Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, UP 243122 India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, UP 243122 India
| | - Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, UP 243122 India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, UP 243122 India
| | - B. P. Mishra
- ICAR-National Bureau of Animal Genetic Resources, Karnal, Karnal, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, Bareilly, UP 243122 India
| |
Collapse
|
23
|
Sadler DE, van Dijk S, Karjalainen J, Watts PC, Uusi‐Heikkilä S. Does size-selective harvesting erode adaptive potential to thermal stress? Ecol Evol 2024; 14:e11007. [PMID: 38333098 PMCID: PMC10850808 DOI: 10.1002/ece3.11007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Revised: 01/12/2024] [Accepted: 01/30/2024] [Indexed: 02/10/2024] Open
Abstract
Overharvesting is a serious threat to many fish populations. High mortality and directional selection on body size can cause evolutionary change in exploited populations via selection for a specific phenotype and a potential reduction in phenotypic diversity. Whether the loss of phenotypic diversity that accompanies directional selection impairs response to environmental stress is not known. To address this question, we exposed three zebrafish selection lines to thermal stress. Two lines had experienced directional selection for (1) large and (2) small body size, and one was (3) subject to random removal of individuals with respect to body size (i.e. line with no directional selection). Selection lines were exposed to three temperatures (elevated, 34°C; ambient, 28°C; low, 22°C) to determine the response to an environmental stressor (thermal stress). We assessed differences among selection lines in their life history (growth and reproduction), physiological traits (metabolic rate and critical thermal max) and behaviour (activity and feeding behaviour) when reared at different temperatures. Lines experiencing directional selection (i.e. size selected) showed reduced growth rate and a shift in average phenotype in response to lower or elevated thermal stress compared with fish from the random-selected line. Our data indicate that populations exposed to directional selection can have a more limited capacity to respond to thermal stress compared with fish that experience a comparable reduction in population size (but without directional selection). Future studies should aim to understand the impacts of environmental stressors on natural fish stocks.
Collapse
Affiliation(s)
- Daniel E. Sadler
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskyläFinland
| | - Stephan van Dijk
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskyläFinland
| | - Juha Karjalainen
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskyläFinland
| | - Phillip C. Watts
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskyläFinland
| | - Silva Uusi‐Heikkilä
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskyläFinland
| |
Collapse
|
24
|
Chase MA, Vilcot M, Mugal CF. The role of recombination dynamics in shaping signatures of direct and indirect selection across the Ficedula flycatcher genome †. Proc Biol Sci 2024; 291:20232382. [PMID: 38228173 DOI: 10.1098/rspb.2023.2382] [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/20/2022] [Accepted: 12/14/2023] [Indexed: 01/18/2024] Open
Abstract
Recombination is a central evolutionary process that reshuffles combinations of alleles along chromosomes, and consequently is expected to influence the efficacy of direct selection via Hill-Robertson interference. Additionally, the indirect effects of selection on neutral genetic diversity are expected to show a negative relationship with recombination rate, as background selection and genetic hitchhiking are stronger when recombination rate is low. However, owing to the limited availability of recombination rate estimates across divergent species, the impact of evolutionary changes in recombination rate on genomic signatures of selection remains largely unexplored. To address this question, we estimate recombination rate in two Ficedula flycatcher species, the taiga flycatcher (Ficedula albicilla) and collared flycatcher (Ficedula albicollis). We show that recombination rate is strongly correlated with signatures of indirect selection, and that evolutionary changes in recombination rate between species have observable impacts on this relationship. Conversely, signatures of direct selection on coding sequences show little to no relationship with recombination rate, even when restricted to genes where recombination rate is conserved between species. Thus, using measures of indirect and direct selection that bridge micro- and macro-evolutionary timescales, we demonstrate that the role of recombination rate and its dynamics varies for different signatures of selection.
Collapse
Affiliation(s)
- Madeline A Chase
- Department of Ecology and Genetics, Uppsala University, 75236 Uppsala, Sweden
- Swiss Ornithological Institute, 6204 Sempach, Switzerland
| | - Maurine Vilcot
- Department of Ecology and Genetics, Uppsala University, 75236 Uppsala, Sweden
- CEFE, University of Montpellier, CNRS, EPHE, IRD, 34293 Montpellier 5, France
| | - Carina F Mugal
- Department of Ecology and Genetics, Uppsala University, 75236 Uppsala, Sweden
- Laboratory of Biometry and Evolutionary Biology, University of Lyon 1, CNRS UMR 5558, 69622 Villeurbanne cedex, France
| |
Collapse
|
25
|
Diaz‐Martin Z, De Vitis M, Havens K, Kramer AT, MacKechnie LM, Fant J. Species-specific effects of production practices on genetic diversity in plant reintroduction programs. Evol Appl 2023; 16:1956-1968. [PMID: 38143906 PMCID: PMC10739063 DOI: 10.1111/eva.13614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 10/17/2023] [Accepted: 10/18/2023] [Indexed: 12/26/2023] Open
Abstract
Plant production practices can influence the genetic diversity of cultivated plant materials and, ultimately, their potential to adapt to a reintroduction site. A common step in the plant production process is the application of seed pretreatment to alleviate physiological seed dormancy and successfully germinate seeds. In production settings, the seeds that germinate more rapidly may be favored in order to fill plant quotas. In this study, we investigated how the application of cold-moist stratification treatments with different durations can lead to differences in the genetic diversity of the propagated plant materials. Specifically, we exposed seeds of three Viola species to two different cold stratification durations, and then we analyzed the genetic diversity of the resulting subpopulations through double-digestion restriction site-associated sequencing (ddRADseq). Our results show that, in two out of three species, utilizing a short stratification period will decrease the genetic diversity of neutral and expressed loci, likely due to the imposition of a genetic bottleneck and artificial selection. We conclude that, in some species, the use of minimal stratification practices in production may jeopardize the adaptive potential and long-term persistence of reintroduced populations and suggest that practitioners carefully consider the evolutionary implications of their production protocols. We highlight the need to consider the germination ecology of target species when selecting the length of dormancy-breaking pretreatments.
Collapse
Affiliation(s)
- Zoe Diaz‐Martin
- Department of BiologySpelman CollegeAtlantaGeorgiaUSA
- Chicago Botanic GardenNegaunee Institute for Plant Conservation Science and ActionGlencoeIllinoisUSA
- Plant Biology and ConservationNorthwestern UniversityEvanstonIllinoisUSA
| | - Marcello De Vitis
- Chicago Botanic GardenNegaunee Institute for Plant Conservation Science and ActionGlencoeIllinoisUSA
- Plant Biology and ConservationNorthwestern UniversityEvanstonIllinoisUSA
- Southeastern Grasslands InstituteAustin Peay State UniversityClarksvilleTennesseeUSA
| | - Kayri Havens
- Chicago Botanic GardenNegaunee Institute for Plant Conservation Science and ActionGlencoeIllinoisUSA
- Plant Biology and ConservationNorthwestern UniversityEvanstonIllinoisUSA
| | - Andrea T. Kramer
- Chicago Botanic GardenNegaunee Institute for Plant Conservation Science and ActionGlencoeIllinoisUSA
- Plant Biology and ConservationNorthwestern UniversityEvanstonIllinoisUSA
| | | | - Jeremie Fant
- Chicago Botanic GardenNegaunee Institute for Plant Conservation Science and ActionGlencoeIllinoisUSA
- Plant Biology and ConservationNorthwestern UniversityEvanstonIllinoisUSA
| |
Collapse
|
26
|
Höllinger I, Wölfl B, Hermisson J. A theory of oligogenic adaptation of a quantitative trait. Genetics 2023; 225:iyad139. [PMID: 37550847 PMCID: PMC10550320 DOI: 10.1093/genetics/iyad139] [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: 04/20/2023] [Revised: 04/20/2023] [Accepted: 07/13/2023] [Indexed: 08/09/2023] Open
Abstract
Rapid phenotypic adaptation is widespread in nature, but the underlying genetic dynamics remain controversial. Whereas population genetics envisages sequential beneficial substitutions, quantitative genetics assumes a collective response through subtle shifts in allele frequencies. This dichotomy of a monogenic and a highly polygenic view of adaptation raises the question of a middle ground, as well as the factors controlling the transition. Here, we consider an additive quantitative trait with equal locus effects under Gaussian stabilizing selection that adapts to a new trait optimum after an environmental change. We present an analytical framework based on Yule branching processes to describe how phenotypic adaptation is achieved by collective changes in allele frequencies at the underlying loci. In particular, we derive an approximation for the joint allele-frequency distribution conditioned on the trait mean as a comprehensive descriptor of the adaptive architecture. Depending on the model parameters, this architecture reproduces the well-known patterns of sequential, monogenic sweeps, or of subtle, polygenic frequency shifts. Between these endpoints, we observe oligogenic architecture types that exhibit characteristic patterns of partial sweeps. We find that a single compound parameter, the population-scaled background mutation rate Θbg, is the most important predictor of the type of adaptation, while selection strength, the number of loci in the genetic basis, and linkage only play a minor role.
Collapse
Affiliation(s)
- Ilse Höllinger
- Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria
| | - Benjamin Wölfl
- Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria
- Vienna Graduate School of Population Genetics, University of Vienna and Veterinary Medical University of Vienna, Vienna, Austria
- Vienna Doctoral School of Ecology and Evolution, University of Vienna, Vienna, Austria
| | - Joachim Hermisson
- Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria
- Max Perutz Labs, Vienna Biocenter Campus (VBC), Dr.-Bohr-Gasse 9, 1030 Vienna, Austria
| |
Collapse
|
27
|
Soni V, Johri P, Jensen JD. Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models. Evolution 2023; 77:2113-2127. [PMID: 37395482 PMCID: PMC10547124 DOI: 10.1093/evolut/qpad120] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Revised: 06/15/2023] [Accepted: 06/30/2023] [Indexed: 07/04/2023]
Abstract
The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, Arizona State University, Tempe, AZ, United States
| | - Parul Johri
- School of Life Sciences, Arizona State University, Tempe, AZ, United States
| | - Jeffrey D Jensen
- School of Life Sciences, Arizona State University, Tempe, AZ, United States
| |
Collapse
|
28
|
Abstract
AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.
Collapse
|
29
|
Wang X, Ingvarsson PK. Quantifying adaptive evolution and the effects of natural selection across the Norway spruce genome. Mol Ecol 2023; 32:5288-5304. [PMID: 37622583 DOI: 10.1111/mec.17106] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Revised: 08/07/2023] [Accepted: 08/09/2023] [Indexed: 08/26/2023]
Abstract
Detecting natural selection is one of the major goals of evolutionary genomics. Here, we sequenced the whole genome of 25 Picea abies individuals and quantified the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we showed that both negative selection and the rate of positively selected substitutions are very limited in coding regions. We found a positive correlation between the rate of adaptive substitutions and recombination rate and a negative correlation between the rate of adaptive substitutions and gene density, suggesting a widespread influence from Hill-Robertson interference on the efficiency of protein adaptation in P. abies. Finally, the distinct population statistics between genomic regions under either positive or balancing selection with that under neutral regions indicated the impact of natural selection on the genomic architecture of Norway spruce. Further gene ontology enrichment analysis for genes located in regions identified as undergoing either positive or long-term balancing selection also highlighted the specific molecular functions and biological processes that appear to be targets of selection in Norway spruce.
Collapse
Affiliation(s)
- Xi Wang
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
30
|
Korfmann K, Abu Awad D, Tellier A. Weak seed banks influence the signature and detectability of selective sweeps. J Evol Biol 2023; 36:1282-1294. [PMID: 37551039 DOI: 10.1111/jeb.14204] [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/02/2023] [Revised: 06/20/2023] [Accepted: 06/27/2023] [Indexed: 08/09/2023]
Abstract
Seed banking (or dormancy) is a widespread bet-hedging strategy, generating a form of population overlap, which decreases the magnitude of genetic drift. The methodological complexity of integrating this trait implies it is ignored when developing tools to detect selective sweeps. But, as dormancy lengthens the ancestral recombination graph (ARG), increasing times to fixation, it can change the genomic signatures of selection. To detect genes under positive selection in seed banking species it is important to (1) determine whether the efficacy of selection is affected, and (2) predict the patterns of nucleotide diversity at and around positively selected alleles. We present the first tree sequence-based simulation program integrating a weak seed bank to examine the dynamics and genomic footprints of beneficial alleles in a finite population. We find that seed banking does not affect the probability of fixation and confirm expectations of increased times to fixation. We also confirm earlier findings that, for strong selection, the times to fixation are not scaled by the inbreeding effective population size in the presence of seed banks, but are shorter than would be expected. As seed banking increases the effective recombination rate, footprints of sweeps appear narrower around the selected sites and due to the scaling of the ARG are detectable for longer periods of time. The developed simulation tool can be used to predict the footprints of selection and draw statistical inference of past evolutionary events in plants, invertebrates, or fungi with seed banks.
Collapse
Affiliation(s)
- Kevin Korfmann
- Department of Life Science Systems, School of Life Sciences, Technical University of Munich, München, Germany
| | - Diala Abu Awad
- Department of Life Science Systems, School of Life Sciences, Technical University of Munich, München, Germany
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, GQE-Le Moulon, Gif-sur-Yvette, France
| | - Aurélien Tellier
- Department of Life Science Systems, School of Life Sciences, Technical University of Munich, München, Germany
| |
Collapse
|
31
|
Feng J, Zhu W, Shi H, Peng D, Zang L, Wang Y, ZhaXi L, BaiMa J, Amevor FK, Wang X, Ma X, Zhao X. Analysis of the Selection Signal of the Tibetan Black Chicken Genome Based on Whole-Genome Sequencing. Genes (Basel) 2023; 14:1672. [PMID: 37761812 PMCID: PMC10531317 DOI: 10.3390/genes14091672] [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: 07/25/2023] [Revised: 08/18/2023] [Accepted: 08/22/2023] [Indexed: 09/29/2023] Open
Abstract
BACKGROUND The Tibetan chicken has adapted well to high altitudes genetically after its long-term habitation in the plateau. In this study, we analyzed the selection signal of Tibetan black chickens (TBCs) and discovered genes associated with the characteristics of germplasm. METHODS Whole-genome sequencing (WGS) was used to identify the single-nucleotide polymorphism (SNP) markers and genetic structures in the genome of Tibetan black chickens. Further, we performed a comparative population genomics analysis between the genomic data obtained in this present study and the genomic data for five wild red jungle fowls (RJFs) accessed from the NCBI database (GenBank accession number PRJNA241474). Thereafter, the Fst and Pi selections were used to identify genes under positive selection in the Tibetan black chicken genome. RESULTS A total of 9,490,690 SNPs were identified in the Tibetan black chickens. In addition, the results from the gene ontology (GO) analysis showed that 732 genes of TBCs were enriched in a total of 210 GO terms with specific molecular functions such as regulation of cellular catabolic process, the MAPK signaling pathway, regulation of ion transport, growth, morphogenesis and lung alveolus development which may provide a better mechanism to facilitate oxygen transport and utilization in TBCs. Moreover, the results from the KEGG analysis showed that 732 genes of the TBCs were significantly enriched in the calcium signaling pathway, circadian entrainment (ADCY1, GNG7 and PER3), oxytocin signaling pathway and pathways of multiple neurodegeneration diseases. In addition, the CD86 antigen (CD86) was identified as a gene associated with the immune response in chickens. It was also revealed that genes such as TRIT1, HPCAL4, NT5C1A and HEYL were discovered under selection in Tibetan black chickens on chromosome 23. These genes may be related to the local adaptive characteristics of Tibetan black chickens, for instance, NT5C1A and HEYL may be involved in the high-altitude adaption of oxygen delivery in Tibetan black chickens. CONCLUSIONS In summary, we found that selection mainly affects the disease resistance and cold acclimatization of Tibetan black chickens. Hence, these results may provide important genetic information for the evolution and breeding of Tibetan black chickens.
Collapse
Affiliation(s)
- Jing Feng
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lhasa 850009, China
| | - Wei Zhu
- State Key Laboratory of Swine and Poultry Breeding Industry, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (W.Z.); (F.K.A.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| | - Hairen Shi
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
- College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China
| | - Da Peng
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
| | - Lei Zang
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
| | - Yan Wang
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
| | - Luobu ZhaXi
- Shannan Longzi County Agriculture and Animal Husbandry Comprehensive Service Center, Shannan 856600, China (J.B.)
| | - Jiancai BaiMa
- Shannan Longzi County Agriculture and Animal Husbandry Comprehensive Service Center, Shannan 856600, China (J.B.)
| | - Felix Kwame Amevor
- State Key Laboratory of Swine and Poultry Breeding Industry, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (W.Z.); (F.K.A.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| | - Xiaoqi Wang
- Agriculture and Animal Husbandry Comprehensive Service Center of Lazi County, Shigatse 858100, China;
| | - Xueying Ma
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agricultural and Animal Husbandry Science, Lhasa 850009, China; (H.S.); (D.P.); (Y.W.); (X.M.)
| | - Xiaoling Zhao
- State Key Laboratory of Swine and Poultry Breeding Industry, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (W.Z.); (F.K.A.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| |
Collapse
|
32
|
Shishir TA, Saha O, Rajia S, Mondol SM, Masum MHU, Rahaman MM, Hossen F, Bahadur NM, Ahmed F, Naser IB, Amin MR. Genome-wide study of globally distributed respiratory syncytial virus (RSV) strains implicates diversification utilizing phylodynamics and mutational analysis. Sci Rep 2023; 13:13531. [PMID: 37598270 PMCID: PMC10439963 DOI: 10.1038/s41598-023-40760-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Accepted: 08/16/2023] [Indexed: 08/21/2023] Open
Abstract
Respiratory syncytial virus (RSV) is a common respiratory pathogen that causes mild cold-like symptoms and severe lower respiratory tract infections, causing hospitalizations in children, the elderly and immunocompromised individuals. Due to genetic variability, this virus causes life-threatening pneumonia and bronchiolitis in young infants. Thus, we examined 3600 whole genome sequences submitted to GISAID by 31 December 2022 to examine the genetic variability of RSV. While RSVA and RSVB coexist throughout RSV seasons, RSVA is more prevalent, fatal, and epidemic-prone in several countries, including the United States, the United Kingdom, Australia, and China. Additionally, the virus's attachment glycoprotein and fusion protein were highly mutated, with RSVA having higher Shannon entropy than RSVB. The genetic makeup of these viruses contributes significantly to their prevalence and epidemic potential. Several strain-specific SNPs co-occurred with specific haplotypes of RSVA and RSVB, followed by different haplotypes of the viruses. RSVA and RSVB have the highest linkage probability at loci T12844A/T3483C and G13959T/C2198T, respectively. The results indicate that specific haplotypes and SNPs may significantly affect their spread. Overall, this analysis presents a promising strategy for tracking the evolving epidemic situation and genetic variants of RSV, which could aid in developing effective control, prophylactic, and treatment strategies.
Collapse
Affiliation(s)
- Tushar Ahmed Shishir
- Department of Mathematics and Natural Sciences, BRAC University, Dhaka, Bangladesh
| | - Otun Saha
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh.
| | - Sultana Rajia
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | | | - Md Habib Ullah Masum
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | | | - Foysal Hossen
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | | | - Firoz Ahmed
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | - Iftekhar Bin Naser
- Department of Mathematics and Natural Sciences, BRAC University, Dhaka, Bangladesh
| | - Mohammad Ruhul Amin
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh.
| |
Collapse
|
33
|
Nickles GR, Oestereicher B, Keller NP, Drott M. Mining for a new class of fungal natural products: the evolution, diversity, and distribution of isocyanide synthase biosynthetic gene clusters. Nucleic Acids Res 2023; 51:7220-7235. [PMID: 37427794 PMCID: PMC10415135 DOI: 10.1093/nar/gkad573] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Revised: 06/16/2023] [Accepted: 07/06/2023] [Indexed: 07/11/2023] Open
Abstract
The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We amalgamated a pipeline of tools to predict BGCs based on shared promoter motifs and located 3800 ICS BGCs in 3300 genomes, making ICS BGCs the fifth largest class of specialized metabolites compared to canonical classes found by antiSMASH. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1/2 gene cluster family (GCF), which was prior only studied in yeast, is present in ∼30% of all Ascomycetes. The dit variety ICS exhibits greater similarity to bacterial ICS than other fungal ICS, suggesting a potential convergence of the ICS backbone domain. The evolutionary origins of the dit GCF in Ascomycota are ancient and these genes are diversifying in some lineages. Our results create a roadmap for future research into ICS BGCs. We developed a website (https://isocyanides.fungi.wisc.edu/) that facilitates the exploration and downloading of all identified fungal ICS BGCs and GCFs.
Collapse
Affiliation(s)
- Grant R Nickles
- Department of Medical Microbiology and Immunology, University of Wisconsin—Madison, Madison, WI 53706, USA
| | | | - Nancy P Keller
- Department of Medical Microbiology and Immunology, University of Wisconsin—Madison, Madison, WI 53706, USA
- Department of Plant Pathology, University of Wisconsin—Madison, Madison, WI 53706, USA
| | - Milton T Drott
- USDA-ARS Cereal Disease Lab (CDL), St. Paul, MN 55108, USA
| |
Collapse
|
34
|
Soni V, Johri P, Jensen JD. Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.15.545166. [PMID: 37398347 PMCID: PMC10312679 DOI: 10.1101/2023.06.15.545166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/04/2023]
Abstract
The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modelled by a realistic mutation rate and as part of a realistic distribution of fitness effects (DFE), as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modelled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false positive rates are in excess of true positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong. Teaser Text Outlier-based genomic scans have proven a popular approach for identifying loci that have potentially experienced recent positive selection. However, it has previously been shown that an evolutionarily appropriate baseline model that incorporates non-equilibrium population histories, purifying and background selection, and variation in mutation and recombination rates is necessary to reduce often extreme false positive rates when performing genomic scans. Here we evaluate the power to detect recurrent selective sweeps using common SFS-based and haplotype-based methods under these increasingly realistic models. We find that while these appropriate evolutionary baselines are essential to reduce false positive rates, the power to accurately detect recurrent selective sweeps is generally low across much of the biologically relevant parameter space.
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Parul Johri
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
- Present address: Department of Biology, Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | | |
Collapse
|
35
|
Ndiaye YD, Wong W, Thwing J, Schaffner SS, Tine A, Diallo MA, Deme A, Sy M, Bei AK, Thiaw AB, Daniels R, Ndiaye T, Gaye A, Ndiaye IM, Toure M, Gadiaga N, Sene A, Sow D, Garba MN, Yade MS, Dieye B, Diongue K, Zoumarou D, Ndiaye A, Gomis J, Fall FB, Ndiop M, Diallo I, Sene D, Macinnis B, Seck MC, Ndiaye M, Badiane AS, Hartl DL, Volkman SK, Wirth DF, Ndiaye D. Two decades of molecular surveillance in Senegal reveal changes in known drug resistance mutations associated with historical drug use and seasonal malaria chemoprevention. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.04.24.23288820. [PMID: 37163114 PMCID: PMC10168519 DOI: 10.1101/2023.04.24.23288820] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. We analyzed data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasite strains in Senegal to determine how historical changes in drug administration policy may have affected parasite evolution. We profiled several known drug resistance markers and their surrounding haplotypes using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole-genome sequence (WGS) based population genomics. We observed rapid changes in drug resistance markers associated with the withdrawal of chloroquine and introduction of sulfadoxine-pyrimethamine in 2003. We also observed a rapid increase in Pfcrt K76T and decline in Pfdhps A437G starting in 2014, which we hypothesize may reflect changes in resistance or fitness caused by seasonal malaria chemoprevention (SMC). Parasite populations evolve rapidly in response to drug use, and SMC preventive efficacy should be closely monitored.
Collapse
Affiliation(s)
- Yaye Die Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Wesley Wong
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
| | - Julie Thwing
- Centers for Disease Control and Prevention, 1600 Clifton Road, Atlanta, GA ,30329, USA
| | - Stephen S Schaffner
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Abdoulaye Tine
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mamadou Alpha Diallo
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Awa Deme
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mouhammad Sy
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Amy K Bei
- Yale School of Public Health, 60 College St, New Haven, CT 06510
| | - Alphonse B Thiaw
- Department of biochemistry and Functional Genomics, Sherbrooke University, 2500 Bd de l'Universite, Sherbrooke, QC J1K 2R1, Canada
| | - Rachel Daniels
- RNA Therapeutics Institute, UMass Chan Medical School, 368 Plantation Street, Worcester MA 01605
| | - Tolla Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Amy Gaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Ibrahima Mbaye Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mariama Toure
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Nogaye Gadiaga
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Aita Sene
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Djiby Sow
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mamane N Garba
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mamadou Samba Yade
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Baba Dieye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Khadim Diongue
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Daba Zoumarou
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Aliou Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Jules Gomis
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Fatou Ba Fall
- Yale School of Public Health, 60 College St, New Haven, CT 06510
| | - Medoune Ndiop
- National Malaria Control Program (NMCP), Rue FN 20, Dakar 25270, Senegal
| | - Ibrahima Diallo
- National Malaria Control Program (NMCP), Rue FN 20, Dakar 25270, Senegal
| | - Doudou Sene
- National Malaria Control Program (NMCP), Rue FN 20, Dakar 25270, Senegal
| | - Bronwyn Macinnis
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Mame Cheikh Seck
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Mouhamadou Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Aida S Badiane
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
| | - Daniel L Hartl
- Department of Organismic and Evolutionary Biology, Harvard University, 16 Divinity Avenue, Cambridge, MA, 02138 USA
| | - Sarah K Volkman
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
- Simmons University, 300 The Fenway, Boston, MA, 02115, USA
| | - Dyann F Wirth
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA, 02142, USA
| | - Daouda Ndiaye
- International Research Training Center on Genomics and Health Surveillance (CIGASS), Cheikh Anta Diop University, Dakar, 16477, Senegal
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, 665 Huntington Ave, Boston, MA, 02115, USA
| |
Collapse
|
36
|
Nickles GR, Oestereicher B, Keller NP, Drott MT. Mining for a New Class of Fungal Natural Products: The Evolution, Diversity, and Distribution of Isocyanide Synthase Biosynthetic Gene Clusters. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.17.537281. [PMID: 37131656 PMCID: PMC10153163 DOI: 10.1101/2023.04.17.537281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) have notable bioactivities that mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We developed the first genome-mining pipeline to identify ICS BGCs, locating 3,800 ICS BGCs in 3,300 genomes. Genes in these clusters share promoter motifs and are maintained in contiguous groupings by natural selection. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1 / 2 gene cluster family (GCF), which was thought to only exist in yeast, is present in ∼30% of all Ascomycetes, including many filamentous fungi. The evolutionary history of the dit GCF is marked by deep divergences and phylogenetic incompatibilities that raise questions about convergent evolution and suggest selection or horizontal gene transfers have shaped the evolution of this cluster in some yeast and dimorphic fungi. Our results create a roadmap for future research into ICS BGCs. We developed a website ( www.isocyanides.fungi.wisc.edu ) that facilitates the exploration, filtering, and downloading of all identified fungal ICS BGCs and GCFs.
Collapse
Affiliation(s)
- Grant R. Nickles
- Department of Medical Microbiology and Immunology, University of Wisconsin—Madison, Madison, WI 53706, USA
| | | | - Nancy P. Keller
- Department of Medical Microbiology and Immunology, University of Wisconsin—Madison, Madison, WI 53706, USA
- Department of Plant Pathology, University of Wisconsin—Madison, Madison, WI 53706, USA
| | | |
Collapse
|
37
|
Kaushik S. Effect of beneficial sweeps and background selection on genetic diversity in changing environments. J Theor Biol 2023; 562:111431. [PMID: 36754344 DOI: 10.1016/j.jtbi.2023.111431] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 12/13/2022] [Accepted: 01/27/2023] [Indexed: 02/08/2023]
Abstract
Neutral theory predicts that the genetic diversity within a population is proportional to the effective population size. In contrast, the observed genetic diversity for various species vary in a narrow range for several orders of magnitude change in the population sizes (Lewontin's paradox). The selective sweeps and background selection, reduce the genetic variation at the linked neutral sites and have been studied considering the environment to be selectively constant. However, in a natural population, the selective environment varies with time. Here, we investigate the impact of selective sweeps and background selection on neutral genetic diversity when the selection coefficient changes periodically over time. The reduction in genetic variation due to selective sweeps is known to depend on the conditional fixation time. Here, we find that the effect of changing environment on conditional mean fixation time is most substantial for the randomly mating population than the inbreeding population with arbitrary inbreeding coefficient. We also study the effect of background selection on neutral sites when the selection coefficient of linked deleterious mutation change periodically in time. In the slowly changing environment, we find that neutral heterozygosity is significantly different, and the site frequency spectrum has a different shape than that in the static environment.
Collapse
Affiliation(s)
- Sachin Kaushik
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560064, India.
| |
Collapse
|
38
|
Thorn CS, Maness RW, Hulke JM, Delmore KE, Criscione CD. Population genomics of helminth parasites. J Helminthol 2023; 97:e29. [PMID: 36927601 DOI: 10.1017/s0022149x23000123] [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] [Indexed: 03/18/2023]
Abstract
Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.
Collapse
Affiliation(s)
- C S Thorn
- Department of Biology, Texas A&M University, 3258 TAMU, College Station, TX, 77843, USA
| | - R W Maness
- Department of Biology, Texas A&M University, 3258 TAMU, College Station, TX, 77843, USA
| | - J M Hulke
- Department of Biology, Texas A&M University, 3258 TAMU, College Station, TX, 77843, USA
| | - K E Delmore
- Department of Biology, Texas A&M University, 3258 TAMU, College Station, TX, 77843, USA
| | - C D Criscione
- Department of Biology, Texas A&M University, 3258 TAMU, College Station, TX, 77843, USA
| |
Collapse
|
39
|
Fluctuating selection and the determinants of genetic variation. Trends Genet 2023; 39:491-504. [PMID: 36890036 DOI: 10.1016/j.tig.2023.02.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Revised: 02/01/2023] [Accepted: 02/07/2023] [Indexed: 03/08/2023]
Abstract
Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.
Collapse
|
40
|
Cutter AD. Synthetic gene drives as an anthropogenic evolutionary force. Trends Genet 2023; 39:347-357. [PMID: 36997427 DOI: 10.1016/j.tig.2023.02.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Revised: 02/13/2023] [Accepted: 02/14/2023] [Indexed: 03/30/2023]
Abstract
Genetic drive represents a fundamental evolutionary force that can exact profound change to the genetic composition of populations by biasing allele transmission. Herein I propose that the use of synthetic homing gene drives, the human-mediated analog of endogenous genetic drives, warrants the designation of 'genetic welding' as an anthropogenic evolutionary force. Conceptually, this distinction parallels that of artificial and natural selection. Genetic welding is capable of imposing complex and rapid heritable phenotypic change on entire populations, whether motivated by biodiversity conservation or public health. Unanticipated possible long-term evolutionary outcomes, however, demand further investigation and bioethical consideration. The emerging importance of genetic welding also compels our explicit recognition of genetic drive as an addition to the other four fundamental forces of evolution.
Collapse
|
41
|
Chung MY, Merilä J, Kim Y, Mao K, López‐Pujol J, Chung MG. A review on
Q
ST
–
F
ST
comparisons of seed plants: Insights for conservation. Ecol Evol 2023; 13:e9926. [PMID: 37006890 PMCID: PMC10049885 DOI: 10.1002/ece3.9926] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2022] [Revised: 02/14/2023] [Accepted: 03/02/2023] [Indexed: 03/30/2023] Open
Abstract
Increased access to genome-wide data provides new opportunities for plant conservation. However, information on neutral genetic diversity in a small number of marker loci can still be valuable because genomic data are not available to most rare plant species. In the hope of bridging the gap between conservation science and practice, we outline how conservation practitioners can more efficiently employ population genetic information in plant conservation. We first review the current knowledge about neutral genetic variation (NGV) and adaptive genetic variation (AGV) in seed plants, regarding both within-population and among-population components. We then introduce the estimates of among-population genetic differentiation in quantitative traits (Q ST) and neutral markers (F ST) to plant biology and summarize conservation applications derived from Q ST-F ST comparisons, particularly on how to capture most AGV and NGV on both in-situ and ex-situ programs. Based on a review of published studies, we found that, on average, two and four populations would be needed for woody perennials (n = 18) to capture 99% of NGV and AGV, respectively, whereas four populations would be needed in case of herbaceous perennials (n = 14). On average, Q ST is about 3.6, 1.5, and 1.1 times greater than F ST in woody plants, annuals, and herbaceous perennials, respectively. Hence, conservation and management policies or suggestions based solely on inference on F ST could be misleading, particularly in woody species. To maximize the preservation of the maximum levels of both AGV and NGV, we suggest using maximum Q ST rather than average Q ST. We recommend conservation managers and practitioners consider this when formulating further conservation and restoration plans for plant species, particularly woody species.
Collapse
Affiliation(s)
- Mi Yoon Chung
- Department of Biological SciencesChungnam National UniversityDaejeon34134South Korea
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Program, Faculty of Biological and Environmental SciencesUniversity of HelsinkiHelsinkiFI‐00014Finland
- Area of Ecology & BiodiversitySchool of Biological SciencesThe University of Hong KongHong Kong SARChina
| | - Yuseob Kim
- Division of EcoScienceEwha Womans UniversitySeoul03760South Korea
- Department of Life ScienceEwha Womans UniversitySeoul03760South Korea
| | - Kangshan Mao
- Key Laboratory for Bio‐resources and Eco‐environment of Ministry of Education, College of Life Science, State Key Laboratory of Hydraulics and Mountain River EngineeringSichuan UniversityChengdu610065China
| | - Jordi López‐Pujol
- Botanic Institute of Barcelona (IBB), CSIC‐Ajuntament de BarcelonaBarcelona08038CataloniaSpain
- Universidad Espíritu Santo (UEES)Samborondón091650Ecuador
| | - Myong Gi Chung
- Division of Life Science and RINSGyeongsang National UniversityJinju52828South Korea
| |
Collapse
|
42
|
Jiranek J, Gibson A. Diet can alter the cost of resistance to a natural parasite in Caenorhabditis elegans. Ecol Evol 2023; 13:e9793. [PMID: 36789344 PMCID: PMC9911625 DOI: 10.1002/ece3.9793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Revised: 01/07/2023] [Accepted: 01/16/2023] [Indexed: 02/12/2023] Open
Abstract
Resistance to parasites confers a fitness advantage, yet hosts show substantial variation in resistance in natural populations. Evolutionary theory indicates that resistant and susceptible genotypes can coexist if resistance is costly, but there is mixed evidence that resistant individuals have lower fitness in the absence of parasites. One explanation for this discrepancy is that the cost of resistance varies with environmental context. We tested this hypothesis using Caenorhabditis elegans and its natural microsporidian parasite, Nematocida ironsii. We used multiple metrics to compare the fitness of two near-isogenic host genotypes differing at regions associated with resistance to N. ironsii. To quantify the effect of the environment on the cost associated with these known resistance regions, we measured fitness on three microbial diets. We found that the cost of resistance varied with both diet and the measure of fitness. We detected no cost to resistance, irrespective of diet, when fitness was measured as fecundity. However, we detected a cost when fitness was measured in terms of population growth, and the magnitude of this cost varied with diet. These results provide a proof of concept that, by mediating the cost of resistance, environmental context may govern the rate and nature of resistance evolution in heterogeneous environments.
Collapse
Affiliation(s)
- Juliana Jiranek
- Department of BiologyUniversity of VirginiaCharlottesvilleVirginiaUSA
| | - Amanda Gibson
- Department of BiologyUniversity of VirginiaCharlottesvilleVirginiaUSA
| |
Collapse
|
43
|
Bohutínská M, Vlček J, Monnahan P, Kolář F. Population Genomic Analysis of Diploid-Autopolyploid Species. Methods Mol Biol 2023; 2545:297-324. [PMID: 36720820 DOI: 10.1007/978-1-0716-2561-3_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
This chapter outlines an empirical analysis of genome-wide single-nucleotide polymorphism (SNP) variation and its underlying drivers among multiple natural populations within a diploid-autopolyploid species. The aim is to reconstruct the genetic structure among natural populations of varying ploidy and infer footprints of selection in these populations, framed around specific questions that are typically encountered when analyzing a mixed-ploidy data set,e.g., addressing the relevance of natural whole-genome duplication for speciation and adaptation. We briefly review the options for the analysis of polyploid population genomic data involving variant calling, population structure, demographic history inference, and selection scanning approaches. Further, we provide suggestions for methods and associated software, possible caveats, and examples of their application to mixed-ploidy and autopolyploid data sets.
Collapse
Affiliation(s)
- Magdalena Bohutínská
- Department of Botany, Faculty of Science, Charles University, Prague, Czech Republic.,Institute of Botany of the Czech Academy of Sciences, Průhonice, Czech Republic
| | - Jakub Vlček
- Department of Botany, Faculty of Science, Charles University, Prague, Czech Republic
| | - Patrick Monnahan
- Department of Pediatrics, University of Minnesota, Minneapolis, MN, USA
| | - Filip Kolář
- Department of Botany, Faculty of Science, Charles University, Prague, Czech Republic. .,Institute of Botany of the Czech Academy of Sciences, Průhonice, Czech Republic.
| |
Collapse
|
44
|
Frankham R. Effects of genomic homozygosity on total fitness in an invertebrate: lethal equivalent estimates for Drosophila melanogaster. CONSERV GENET 2022. [DOI: 10.1007/s10592-022-01493-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
45
|
Souilmi Y, Tobler R, Johar A, Williams M, Grey ST, Schmidt J, Teixeira JC, Rohrlach A, Tuke J, Johnson O, Gower G, Turney C, Cox M, Cooper A, Huber CD. Admixture has obscured signals of historical hard sweeps in humans. Nat Ecol Evol 2022; 6:2003-2015. [PMID: 36316412 PMCID: PMC9715430 DOI: 10.1038/s41559-022-01914-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2021] [Accepted: 09/16/2022] [Indexed: 11/06/2022]
Abstract
The role of natural selection in shaping biological diversity is an area of intense interest in modern biology. To date, studies of positive selection have primarily relied on genomic datasets from contemporary populations, which are susceptible to confounding factors associated with complex and often unknown aspects of population history. In particular, admixture between diverged populations can distort or hide prior selection events in modern genomes, though this process is not explicitly accounted for in most selection studies despite its apparent ubiquity in humans and other species. Through analyses of ancient and modern human genomes, we show that previously reported Holocene-era admixture has masked more than 50 historic hard sweeps in modern European genomes. Our results imply that this canonical mode of selection has probably been underappreciated in the evolutionary history of humans and suggest that our current understanding of the tempo and mode of selection in natural populations may be inaccurate.
Collapse
Affiliation(s)
- Yassine Souilmi
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia.
| | - Raymond Tobler
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia.
- Evolution of Cultural Diversity Initiative, Australian National University, Canberra, Australian Capital Territory, Australia.
| | - Angad Johar
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia.
- Department of Cardiovascular Diseases, Mayo Clinic, Rochester, MN, USA.
| | - Matthew Williams
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - Shane T Grey
- Transplantation Immunology Group, Immunology Division, Garvan Institute of Medical Research, Darlinghurst, New South Wales, Australia
- St Vincent's Clinical School, Faculty of Medicine, UNSW, Darlinghurst, New South Wales, Australia
| | - Joshua Schmidt
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - João C Teixeira
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - Adam Rohrlach
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, The University of Adelaide, Adelaide, South Australia, Australia
- Department of Archaeogenetics, Max Planck Institute for the Science of Human History, Jena, Germany
| | - Jonathan Tuke
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, The University of Adelaide, Adelaide, South Australia, Australia
- School of Mathematical Sciences, The University of Adelaide, Adelaide, South Australia, Australia
| | - Olivia Johnson
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - Graham Gower
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - Chris Turney
- Chronos 14Carbon-Cycle Facility and Earth and Sustainability Science Research Centre, University of New South Wales, Sydney, New South Wales, Australia
| | - Murray Cox
- Statistics and Bioinformatics Group, School of Fundamental Sciences, Massey University, Palmerston North, New Zealand
| | - Alan Cooper
- South Australian Museum, Adelaide, South Australia, Australia.
- BlueSky Genetics, Ashton, South Australia, Australia.
| | - Christian D Huber
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia.
- Department of Biology, Penn State University, University Park, PA, USA.
| |
Collapse
|
46
|
Schlichta F, Moinet A, Peischl S, Excoffier L. The Impact of Genetic Surfing on Neutral Genomic Diversity. Mol Biol Evol 2022; 39:msac249. [PMID: 36403964 PMCID: PMC9703594 DOI: 10.1093/molbev/msac249] [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] [Indexed: 11/22/2022] Open
Abstract
Range expansions have been common in the history of most species. Serial founder effects and subsequent population growth at expansion fronts typically lead to a loss of genomic diversity along the expansion axis. A frequent consequence is the phenomenon of "gene surfing," where variants located near the expanding front can reach high frequencies or even fix in newly colonized territories. Although gene surfing events have been characterized thoroughly for a specific locus, their effects on linked genomic regions and the overall patterns of genomic diversity have been little investigated. In this study, we simulated the evolution of whole genomes during several types of 1D and 2D range expansions differing by the extent of migration, founder events, and recombination rates. We focused on the characterization of local dips of diversity, or "troughs," taken as a proxy for surfing events. We find that, for a given recombination rate, once we consider the amount of diversity lost since the beginning of the expansion, it is possible to predict the initial evolution of trough density and their average width irrespective of the expansion condition. Furthermore, when recombination rates vary across the genome, we find that troughs are over-represented in regions of low recombination. Therefore, range expansions can leave local and global genomic signatures often interpreted as evidence of past selective events. Given the generality of our results, they could be used as a null model for species having gone through recent expansions, and thus be helpful to correctly interpret many evolutionary biology studies.
Collapse
Affiliation(s)
- Flávia Schlichta
- Computational and Molecular Population Genetics lab, Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Antoine Moinet
- Computational and Molecular Population Genetics lab, Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
- Interfaculty Bioinformatics Unit, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
| | - Stephan Peischl
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
- Interfaculty Bioinformatics Unit, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
| | - Laurent Excoffier
- Computational and Molecular Population Genetics lab, Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
47
|
Min J, Gupta M, Desai MM, Weissman DB. Spatial structure alters the site frequency spectrum produced by hitchhiking. Genetics 2022; 222:iyac139. [PMID: 36094352 PMCID: PMC9630975 DOI: 10.1093/genetics/iyac139] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Accepted: 08/24/2022] [Indexed: 11/13/2022] Open
Abstract
The reduction of genetic diversity due to genetic hitchhiking is widely used to find past selective sweeps from sequencing data, but very little is known about how spatial structure affects hitchhiking. We use mathematical modeling and simulations to find the unfolded site frequency spectrum left by hitchhiking in the genomic region of a sweep in a population occupying a 1D range. For such populations, sweeps spread as Fisher waves, rather than logistically. We find that this leaves a characteristic 3-part site frequency spectrum at loci very close to the swept locus. Very low frequencies are dominated by recent mutations that occurred after the sweep and are unaffected by hitchhiking. At moderately low frequencies, there is a transition zone primarily composed of alleles that briefly "surfed" on the wave of the sweep before falling out of the wavefront, leaving a spectrum close to that expected in well-mixed populations. However, for moderate-to-high frequencies, there is a distinctive scaling regime of the site frequency spectrum produced by alleles that drifted to fixation in the wavefront and then were carried throughout the population. For loci slightly farther away from the swept locus on the genome, recombination is much more effective at restoring diversity in 1D populations than it is in well-mixed ones. We find that these signatures of space can be strong even in apparently well-mixed populations with negligible spatial genetic differentiation, suggesting that spatial structure may frequently distort the signatures of hitchhiking in natural populations.
Collapse
Affiliation(s)
- Jiseon Min
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA
- NSF-Simons Center for Mathematical and Statistical Analysis of Biology, Harvard University, Cambridge, MA 02138, USA
- Quantitative Biology Initiative, Harvard University, Cambridge, MA 02138, USA
| | - Misha Gupta
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - Michael M Desai
- NSF-Simons Center for Mathematical and Statistical Analysis of Biology, Harvard University, Cambridge, MA 02138, USA
- Quantitative Biology Initiative, Harvard University, Cambridge, MA 02138, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
- Department of Physics, Harvard University, Cambridge, MA 02138, USA
| | | |
Collapse
|
48
|
Rushworth CA, Wagner MR, Mitchell-Olds T, Anderson JT. The Boechera model system for evolutionary ecology. AMERICAN JOURNAL OF BOTANY 2022; 109:1939-1961. [PMID: 36371714 DOI: 10.1002/ajb2.16090] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 08/27/2022] [Accepted: 08/30/2022] [Indexed: 06/16/2023]
Abstract
Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.
Collapse
Affiliation(s)
| | - Maggie R Wagner
- Department of Ecology and Evolutionary Biology, Kansas Biological Survey and Center for Ecological Research, University of Kansas, Lawrence, KS, 66045, USA
| | | | - Jill T Anderson
- Department of Genetics and Odum School of Ecology, University of Georgia, Athens, GA, 30602, USA
| |
Collapse
|
49
|
de Greef E, Einfeldt AL, Miller PJO, Ferguson SH, Garroway CJ, Lefort KJ, Paterson IG, Bentzen P, Feyrer LJ. Genomics reveal population structure, evolutionary history, and signatures of selection in the northern bottlenose whale, Hyperoodon ampullatus. Mol Ecol 2022; 31:4919-4931. [PMID: 35947506 PMCID: PMC9804413 DOI: 10.1111/mec.16643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 07/18/2022] [Accepted: 08/03/2022] [Indexed: 01/05/2023]
Abstract
Information on wildlife population structure, demographic history, and adaptations are fundamental to understanding species evolution and informing conservation strategies. To study this ecological context for a cetacean of conservation concern, we conducted the first genomic assessment of the northern bottlenose whale, Hyperoodon ampullatus, using whole-genome resequencing data (n = 37) from five regions across the North Atlantic Ocean. We found a range-wide pattern of isolation-by-distance with a genetic subdivision distinguishing three subgroups: the Scotian Shelf, western North Atlantic, and Jan Mayen regions. Signals of elevated levels of inbreeding in the Endangered Scotian Shelf population indicate this population may be more vulnerable than the other two subgroups. In addition to signatures of inbreeding, evidence of local adaptation in the Scotian Shelf was detected across the genome. We found a long-term decline in effective population size for the species, which poses risks to their genetic diversity and may be exacerbated by the isolating effects of population subdivision. Protecting important habitat and migratory corridors should be prioritized to rebuild population sizes that were diminished by commercial whaling, strengthen gene flow, and ensure animals can move across regions in response to environmental changes.
Collapse
Affiliation(s)
- Evelien de Greef
- Department of BiologyDalhousie UniversityHalifaxNova ScotiaCanada,Department of Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| | | | | | | | - Colin J. Garroway
- Department of Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| | - Kyle J. Lefort
- Department of Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| | - Ian G. Paterson
- Department of BiologyDalhousie UniversityHalifaxNova ScotiaCanada
| | - Paul Bentzen
- Department of BiologyDalhousie UniversityHalifaxNova ScotiaCanada
| | - Laura J. Feyrer
- Department of BiologyDalhousie UniversityHalifaxNova ScotiaCanada,Fisheries and Oceans CanadaBedford Institute of OceanographyDartmouthNova ScotiaCanada
| |
Collapse
|
50
|
Zhang Z, Sui Z, Zhang J, Li Q, Zhang Y, Wang C, Li X, Xing F. Identification of Signatures of Selection for Litter Size and Pubertal Initiation in Two Sheep Populations. Animals (Basel) 2022; 12:ani12192520. [PMID: 36230262 PMCID: PMC9559472 DOI: 10.3390/ani12192520] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Revised: 09/13/2022] [Accepted: 09/19/2022] [Indexed: 11/23/2022] Open
Abstract
Fecundity is an important economic trait in sheep that directly affects their economic and productive efficiency. Our study aimed to identify SNP loci associated with sheep puberty or litter size which could be used in future breeding programs to improve fertility. Genomic DNA was obtained from Hetian and Cele Black sheep breeds and used for reduced-representation genome sequencing to identify SNP loci associated with pubertal initiation and litter size. Selective signatures analysis was performed based on the fixation index and nucleotide diversity, followed by pathway analysis of the genes contained in the selected regions. The selected SNP loci in the genes associated with pubertal initiation and litter size were validated using both sheep breeds. In total, 384,718 high quality SNPs were obtained and 376 genes were selected. Functional annotation of genes and enrichment analysis identified 12 genes associated with pubertal initiation and 11 genes associated with litter size. SNP locus validation showed that two SNP on PAK1 and four on ADCY1 may be associated with pubertal initiation, and one SNP on GNAQ gene (NC_040253.1: g.62677376G > A) was associated with litter size in Cele Black sheep. Our results provide new theoretical support for sheep breeding.
Collapse
|