Evolvability predicts macroevolution under fluctuating selection.

Heritable variation is a prerequisite for evolutionary change, but the relevance of genetic constraints on macroevolutionary timescales is debated. By using two datasets on fossil and contemporary taxa, we show that evolutionary divergence among populations, and to a lesser extent among species, increases with microevolutionary evolvability. We evaluate and reject several hypotheses to explain this relationship and propose that an effect of evolvability on population and species divergence can be explained by the influence of genetic constraints on the ability of populations to track rapid, stationary environmental fluctuations.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution.

BACKGROUND: Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected. METHODS: We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes. RESULTS: Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host. CONCLUSIONS: Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure.

BACKGROUND: Shallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon's true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India. RESULTS: Shallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia. CONCLUSION: This study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Evolutionary and functional analyses of LRP5 in archaic and extant modern humans.

BACKGROUND: The human lineage has undergone a postcranial skeleton gracilization (i.e. lower bone mass and strength relative to body size) compared to other primates and archaic populations such as the Neanderthals. This gracilization has been traditionally explained by differences in the mechanical load that our ancestors exercised. However, there is growing evidence that gracilization could also be genetically influenced. RESULTS: We have analyzed the LRP5 gene, which is known to be associated with high bone mineral density conditions, from an evolutionary and functional point of view. Taking advantage of the published genomes of archaic Homo populations, our results suggest that this gene has a complex evolutionary history both between archaic and living humans and within living human populations. In particular, we identified the presence of different selective pressures in archaics and extant modern humans, as well as evidence of positive selection in the African and South East Asian populations from the 1000 Genomes Project. Furthermore, we observed a very limited evidence of archaic introgression in this gene (only at three haplotypes of East Asian ancestry out of the 1000 Genomes), compatible with a general erasing of the fingerprint of archaic introgression due to functional differences in archaics compared to extant modern humans. In agreement with this hypothesis, we observed private mutations in the archaic genomes that we experimentally validated as putatively increasing bone mineral density. In particular, four of five archaic missense mutations affecting the first ß-propeller of LRP5 displayed enhanced Wnt pathway activation, of which two also displayed reduced negative regulation. CONCLUSIONS: In summary, these data suggest a genetic component contributing to the understanding of skeletal differences between extant modern humans and archaic Homo populations.

Multiscale selection in spatially structured populations.

The spatial structure of populations is key to many (eco-)evolutionary processes. In such cases, the strength and sign of selection on a trait may depend on the spatial scale considered. An example is the evolution of altruism: selection in local environments often favours cheaters over altruists, but this can be outweighed by selection at larger scales, favouring clusters of altruists over clusters of cheaters. For populations subdivided into distinct groups, this effect is described formally by multilevel selection theory. However, many populations do not consist of non-overlapping groups but rather (self-)organize into other ecological patterns. We therefore present a mathematical framework for multiscale selection. This framework decomposes natural selection into two parts: local selection, acting within environments of a certain size, and interlocal selection, acting among them. Varying the size of the local environments subsequently allows one to measure the contribution to selection of each spatial scale. To illustrate the use of this framework, we apply it to models of the evolution of altruism and pathogen transmissibility. The analysis identifies how and to what extent ecological processes at different spatial scales contribute to selection and compete, thus providing a rigorous underpinning to eco-evolutionary intuitions.

Rapid vertebrate speciation via isolation, bottlenecks, and drift.

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Comparative transcriptome analysis reveals lineage- and environment-specific adaptations in cacti from the Brazilian Atlantic Forest.

MAIN CONCLUSION: Natural selection influenced adaptive divergence between Cereus fernambucensis and Cereus insularis, revealing key genes governing abiotic stress responses and supporting neoteny in C. insularis. Uncovering the molecular mechanisms driving adaptive divergence in traits related to habitat adaptation remains a central challenge. In this study, we focused on the cactus clade, which includes Cereus sericifer F.Ritter, Cereus fernambucensis Lem., and Cereus insularis Hemsley. These allopatric species inhabit distinct relatively drier regions within the Brazilian Atlantic Forest, each facing unique abiotic conditions. We leveraged whole transcriptome data and abiotic variables datasets to explore lineage-specific and environment-specific adaptations in these species. Employing comparative phylogenetic methods, we identified genes under positive selection (PSG) and examined their association with non-synonymous genetic variants and abiotic attributes through a PhyloGWAS approach. Our analysis unveiled signatures of selection in all studied lineages, with C. fernambucensis northern populations and C. insularis showing the most PSGs. These PSGs predominantly govern abiotic stress regulation, encompassing heat tolerance, UV stress response, and soil salinity adaptation. Our exclusive observation of gene expression tied to early developmental stages in C. insularis supports the hypothesis of neoteny in this species. We also identified genes associated with abiotic variables in independent lineages, suggesting their role as environmental filters on genetic diversity. Overall, our findings suggest that natural selection played a pivotal role in the geographic range of these species in response to environmental and biogeographic transitions.

The molecular evolution of cancer associated genes in mammals.

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Improving selection decisions with mating information by accounting for Mendelian sampling variances looking two generations ahead.

BACKGROUND: Breeding programs are judged by the genetic level of animals that are used to disseminate genetic progress. These animals are typically the best ones of the population. To maximise the genetic level of very good animals in the next generation, parents that are more likely to produce top performing offspring need to be selected. The ability of individuals to produce high-performing progeny differs because of differences in their breeding values and gametic variances. Differences in gametic variances among individuals are caused by differences in heterozygosity and linkage. The use of the gametic Mendelian sampling variance has been proposed before, for use in the usefulness criterion or Index5, and in this work, we extend existing approaches by not only considering the gametic Mendelian sampling variance of individuals, but also of their potential offspring. Thus, the criteria developed in this study plan one additional generation ahead. For simplicity, we assumed that the true quantitative trait loci (QTL) effects, genetic map and the haplotypes of all animals are known. RESULTS: In this study, we propose a new selection criterion, ExpBVSelGrOff, which describes the genetic level of selected grand-offspring that are produced by selected offspring of a particular mating. We compare our criterion with other published criteria in a stochastic simulation of an ongoing breeding program for 21 generations for proof of concept. ExpBVSelGrOff performed better than all other tested criteria, like the usefulness criterion or Index5 which have been proposed in the literature, without compromising short-term gains. After only five generations, when selection is strong (1%), selection based on ExpBVSelGrOff achieved 5.8% more commercial genetic gain and retained 25% more genetic variance without compromising inbreeding rate compared to selection based only on breeding values. CONCLUSIONS: Our proposed selection criterion offers a new tool to accelerate genetic progress for contemporary genomic breeding programs. It retains more genetic variance than previously published criteria that plan less far ahead. Considering future gametic Mendelian sampling variances in the selection process also seems promising for maintaining more genetic variance.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits.

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.

High-throughput DNA sequence analysis elucidates novel insight into the genetic basis of adaptation in local sheep.

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

Experimental Evaluation of a Direct Fitness Effect of the De Novo Evolved Mouse Gene Pldi.

De novo evolved genes emerge from random parts of noncoding sequences and have, therefore, no homologs from which a function could be inferred. While expression analysis and knockout experiments can provide insights into the function, they do not directly test whether the gene is beneficial for its carrier. Here, we have used a seminatural environment experiment to test the fitness of the previously identified de novo evolved mouse gene Pldi, which has been implicated to have a role in sperm differentiation. We used a knockout mouse strain for this gene and competed it against its parental wildtype strain for several generations of free reproduction. We found that the knockout (ko) allele frequency decreased consistently across three replicates of the experiment. Using an approximate Bayesian computation framework that simulated the data under a demographic scenario mimicking the experiment's demography, we could estimate a selection coefficient ranging between 0.21 and 0.61 for the wildtype allele compared to the ko allele in males, under various models. This implies a relatively strong selective advantage, which would fix the new gene in less than hundred generations after its emergence.

Genome-wide scans for selection signatures in indigenous chickens reveal candidate genes associated with local adaptation.

Population growth and climate change pose challenges to the sustainability of poultry farming. The emphasis on high-yield traits in commercialized breeds has led to a decline in their adaptability. Chicken varieties adapted to the local environment, possessing traits that facilitate adaptation to climate change, such as disease resistance and tolerance to extreme weather conditions, can improve hybridization outcomes. In this study, we conducted an analysis of the population structure and genetic diversity of 110 chickens representing indigenous breeds from southern China and two different commercial breeds. Further, we performed comparative population genomics, utilizing nucleotide diversity and fixation statistics, to characterize genomic features of natural selection and to identify unique genetic traits and potential selection markers developed by indigenous breeds after adapting to the local environment. Results based on genetic diversity and population structure analyses showed that indigenous varieties exhibited high levels of genetic diversity. Commercial breeds that have been indigenously bred demonstrated higher levels of genetic diversity than those that have not, and breeds with different selection practices displayed significant differences in genetic structure. Additionally, we further searched for potential genomic regions in native chicken ecotypes, uncovering several candidate genes related to ecological adaptations affecting local breeds, such as IKBKB, S1PR1, TSHR, IL1RAPL1 and AMY2A, which are involved in disease resistance, heat tolerance, immune regulation and behavioral traits. This work provides important insights into the genomic characterization of ecotypes of native chicken in southern China. The identification of candidate genes associated with traits such as disease resistance, heat tolerance, immunomodulation, and behavioral traits is a significant outcome. These candidate genes may contribute to the understanding of the molecular basis of the adaptation of the southern native chicken to the local environment. It is recommended that these genes be integrated into chicken breeding programs to enhance sustainable agriculture and promote effective conservation and utilization strategies.

Unraveling inbreeding patterns and selection signals in Alpine Grey cattle.

Inbreeding plays a crucial role in livestock breeding, influencing genetic diversity and phenotypic traits. Genomic data have helped address limitations posed by incomplete pedigrees, providing deeper insights into breed genetic diversity. This study assesses inbreeding levels via pedigree and genomic approaches and analyzes old and recent inbreeding using runs of homozygosity (ROH), and selection signals in Alpine Grey cattle. Pedigree data from 165 575 individuals, analyzed with INBUPGF90 software, computed inbreeding coefficients. Genomic-based coefficients derived from PLINK v1.9. or DetectRUNS R package analyses of 1 180 individuals' genotypes. Common single nucleotide polymorphisms within ROH pinpointed genomic regions, aggregating into "ROH islands" indicative of selection pressure. Overlaps with USCS Genome Browser unveiled gene presence. Moderate correlations (0.20-0.54) existed between pedigree and genomic coefficients, with most genomic estimators having higher (>0.8) correlation values. Inbreeding averaged 0.04 in < 8 Mb ROH segments, and 0.03 in > 16 Mb segments; > 90% of ROHs were < 8 Mb, indicating ancient inbreeding prevalence. Recent inbreeding proved less detrimental than in cosmopolitan breeds. Two major ROH islands on chromosomes 6 and 7 harbored genes linked to immune response, disease resistance (PYURF, HERC3), and fertility (EIF4EBP3, SRA1). This study underscores the need for detailed inbreeding analyses to understand genetic characteristics and historical changes in local breeds like Alpine Grey cattle. Genomic insights, especially from ROH, facilitated overcoming pedigree limitations, illuminating breed genetic diversity. Our findings reveal ancient inbreeding's enduring genetic impact and ROH islands potential for selective sweeps, elucidating traits in Alpine Grey cattle.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Group-selection via aggregative propagule-formation enables cooperative multicellularity in an individual based, spatial model.

The emergence of multicellularity is one of the major transitions in evolution that happened multiple times independently. During aggregative multicellularity, genetically potentially unrelated lineages cooperate to form transient multicellular groups. Unlike clonal multicellularity, aggregative multicellular organisms do not rely on kin selection instead other mechanisms maintain cooperation against cheater phenotypes that benefit from cooperators but do not contribute to groups. Spatiality with limited diffusion can facilitate group selection, as interactions among individuals are restricted to local neighbourhoods only. Selection for larger size (e.g. avoiding predation) may facilitate the emergence of aggregation, though it is unknown, whether and how much role such selection played during the evolution of aggregative multicellularity. We have investigated the effect of spatiality and the necessity of predation on the stability of aggregative multicellularity via individual-based modelling on the ecological timescale. We have examined whether aggregation facilitates the survival of cooperators in a temporally heterogeneous environment against cheaters, where only a subset of the population is allowed to periodically colonize a new, resource-rich habitat. Cooperators constitutively produce adhesive molecules to promote aggregation and propagule-formation while cheaters spare this expense to grow faster but cannot aggregate on their own, hence depending on cooperators for long-term survival. We have compared different population-level reproduction modes with and without individual selection (predation) to evaluate the different hypotheses. In a temporally homogeneous environment without propagule-based colonization, cheaters always win. Predation can benefit cooperators, but it is not enough to maintain the necessary cooperator amount in successive dispersals, either randomly or by fragmentation. Aggregation-based propagation however can ensure the adequate ratio of cooperators-to-cheaters in the propagule and is sufficient to do so even without predation. Spatiality combined with temporal heterogeneity helps cooperators via group selection, thus facilitating aggregative multicellularity. External stress selecting for larger size (e.g. predation) may facilitate aggregation, however, according to our results, it is neither necessary nor sufficient for aggregative multicellularity to be maintained when there is effective group-selection.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality.

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens.

BACKGROUND: Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS: To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION: Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds.

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.