Selective Sweeps Research Articles

Population Genomics of Plasmodium falciparum and Malaria Control Implications in Abidjan (Cote d’ivoire)

Introduction: The onset of Plasmodium falciparum (P. falciparum) resistance to antimalarial drugs requires the careful surveillance of African parasite populations. Genomic tools are implemented to detect evolutionary changes that could impact malaria control and elimination strategies. Here, we evaluate the genome-wide pattern of selection and sequence variation in P. falciparum populations in Abidjan, Côte d’Ivoire. Methods: The study was conducted in three localities of Abidjan from 2013 to 2014. We collected 70 blood samples after receiving written informed consent from patients above two years of age. After extracting P. falciparum and human DNA from isolates, we performed whole-genome sequencing and used population genomics approaches to investigate the genetic diversity and complexity of infections and identify loci under positive directional selection. Results: We observed an excess of rare variants in the population, showing a clear mutation process in the isolates. Moderate Fst estimates (0.3) was detected for surfin, an immune invasion gene family. Seven iHS regions that had at least two SNPs with a score > 3.2 were identified. These regions code for genes that have been under strong directional selection. Two of these genes were the chloroquine resistance transporter (crt) on chromosome 7 and the dihydropteroate reductase (dhps) on chromosome 8. Our analyses showed that a recent selective sweep occurred for the erythrocyte membrane protein (Pfemp1). Conclusion: Our analyses identified genes under selective drug pressure and balancing selection on protective immune-specific genes. These findings demonstrate the effectiveness of genomics analyses to follow the evolution of malaria parasites and adopt appropriate strategies to eliminate malaria in Côte d’Ivoire.

Natural variations in the PbCPK28 promoter regulate sugar content through interaction with PbTST4 and PbVHA-A1 in pear.

Soluble sugars play an important role in plant growth, development and fruit quality. Pear fruits have demonstrated a considerable improvement in sugar quality during their long history of selection. However, little is known about the underlying molecular mechanisms accompanying the changes in fruit sugar content as a result of selection by horticulturists. Here, we identified a calcium-dependent protein kinase (PbCPK28), which is located on LG15 and is present within a selective sweep region, thus linked to the quantitative trait loci for soluble solids. Association analysis indicates that a single nucleotide polymorphism-13 variation (SNP13T/C ) in the PbCPK28 regulatory region led to fructose content diversity in pear. Elevated expression of PbCPK28 resulted in significantly increased fructose levels in pear fruits. Furthermore, PbCPK28 interacts with and phosphorylates PbTST4, a proton antiporter, thereby coupling the sugar import into the vacuole with proton export. We demonstrated that residues S277 and S314 of PbTST4 are crucial for its function. Additionally, PbCPK28 interacts with and phosphorylates the vacuolar hydrogen proton pump PbVHA-A1, which could provide proton motive forces for PbTST4. We also found that the T11 and Y120 phosphorylation sites in PbVHA-A1 are essential for its function. Evolution analysis and yeast-two-hybrid results support that the CPK-TST/CPK-VHA-A regulatory network is highly conserved in plants, especially the corresponding phosphorylation sites. Together, our work identifies an agriculturally important natural variation and an important regulatory network, allowing genetic improvement of fruit sugar contents in pears through modulation of PbCPK28 expression and phosphorylation of PbTST4 and PbVHA-A1.

A key gene for the climatic adaptation of Apis cerana populations in China according to selective sweep analysis

BackgroundApis cerana is widely distributed in China and, prior to the introduction of western honeybees, was the only bee species kept in China. During the long-term natural evolutionary process, many unique phenotypic variations have occurred among A. cerana populations in different geographical regions under varied climates. Understanding the molecular genetic basis and the effects of climate change on the adaptive evolution of A. cerana can promote A. cerana conservation in face of climate change and allow for the effective utilization of its genetic resources.ResultTo investigate the genetic basis of phenotypic variations and the impact of climate change on adaptive evolution, A. cerana workers from 100 colonies located at similar geographical latitudes or longitudes were analyzed. Our results revealed an important relationship between climate types and the genetic variation of A. cerana in China, and a greater influence of latitude compared with longitude was observed. Upon selection and morphometry analyses combination for populations under different climate types, we identified a key gene RAPTOR, which was deeply involved in developmental processes and influenced the body size.ConclusionThe selection of RAPTOR at the genomic level during adaptive evolution could allow A. cerana to actively regulate its metabolism, thereby fine-tuning body sizes in response to harsh conditions caused by climate change, such as food shortages and extreme temperatures, which may partially elucidate the size differences of A. cerana populations. This study provides crucial support for the molecular genetic basis of the expansion and evolution of naturally distributed honeybee populations.

Recombination Variation Shapes Phylogeny and Introgression in Wild Diploid Strawberries.

Introgressive hybridization is widespread in wild plants and has important consequences. However, frequent hybridization between species makes the estimation of the species' phylogeny challenging, and little is known about the genomic landscape of introgression as it results from complex interactions of multiple evolutionary processes. Here, we reconstructed the phylogeny of ten wild diploid strawberries with whole genome resequencing data and then investigated the influence of recombination rate variation on phylogeny and introgression. We found that genomic regions with low recombination showed reduced levels of incomplete lineage sorting and introgression, and concentrated phylogenetic signals, thus contributing to the most likely species tree of wild diploid strawberries. We revealed complex and widespread introgression across the genus Fragaria, with an average proportion of approximately 4.1% of the extant genome. Introgression tends to be retained in the regions with high recombination rates and low gene density. Furthermore, we identified four SLF genes under selective sweeps that may play potential roles in the possible regain of self-incompatibility by ancient introgression. Altogether, our study yielded novel insights into the evolutionary history and genomic characteristics of introgression in wild diploid strawberries and provides evidence for the role of introgression in plant mating system transitions.

Extraordinary selection on the human X chromosome associated with archaic admixture.

The X chromosome in non-African humans shows less diversity and less Neanderthal introgression than expected from neutral evolution. Analyzing 162 human male X chromosomes worldwide, we identified fourteen chromosomal regions where nearly identical haplotypes spanning several hundred kilobases are found at high frequencies in non-Africans. Genetic drift alone cannot explain the existence of these haplotypes, which must have been associated with strong positive selection in partial selective sweeps. Moreover, the swept haplotypes are entirely devoid of archaic ancestry as opposed to the non-swept haplotypes in the same genomic regions. The ancient Ust'-Ishim male dated at 45,000 before the present (BP)also carries the swept haplotypes, implying that selection on the haplotypes must have occurred between 45,000 and 55,000 years ago. Finally, we find that the chromosomal positions of sweeps overlap previously reported hotspots of selectivesweeps in great ape evolution, suggesting a mechanism of selection unique to X chromosomes.

Genomic diversity and signals of selection processes in wild and farm-reared red-legged partridges (Alectoris rufa)

The genetic dynamics of wild populations with releases of farm-reared reinforcements are very complex. These releases can endanger wild populations through genetic swamping or by displacing them. We assessed the genomic differences between wild and farm-reared red-legged partridges (Alectoris rufa) and described differential selection signals between both populations. We sequenced the whole genome of 30 wild and 30 farm-reared partridges. Both partridges had similar nucleotide diversity (π). Farm-reared partridges had a more negative Tajima's D and more and longer regions of extended haplotype homozygosity than wild partridges. We observed higher inbreeding coefficients (FIS and FROH) in wild partridges. Selective sweeps (Rsb) were enriched with genes that contribute to the reproductive, skin and feather colouring, and behavioural differences between wild and farm-reared partridges. The analysis of genomic diversity should inform future decisions for the preservation of wild populations.

Integrating Whole-Genome Resequencing and RNA Sequencing Data Reveals Selective Sweeps and Differentially Expressed Genes Related to Nervous System Changes in Luxi Gamecocks.

The Luxi gamecock developed very unique morphological and behavioral features under the special artificial selection of the most famous Chinese gamecocks. There are very few research studies on the genetics and selection of the Luxi gamecock. We used six methods (Fst, Tajima's D, hapFLK, iHS, XP-EHH, and Runs of homozygosity) to detect selective sweeps in whole-genome resequencing data of 19 Luxi gamecocks compared to other Chinese indigenous chickens. Eleven genes that were highly related to nervous system development (CDH18, SLITRK1, SLITRK6, NDST3, ATP23, LRIG3, IL1RAPL1, GADL1, C5orf22, UGT8, WISP1, and WNT9A) appeared in at least four methods and were regarded as the most significant genes under selection. Differentially expressed gene (DEG) analysis based on the RNA sequencing data of the cerebral cortex and middle brain between six Luxi gamecocks, Tibetan chickens, and white leghorns found that most differentially expressed genes were enriched in pathways with nervous system functions. Genes associated with aggressiveness-related neurotransmitters (SLC4A2, DRD1, DRD2, ADRA2A, and ADRA2B) showed differential expression rates in Luxi gamecocks as well. Combined results showed that most genes in selective sweep regions were also differentially expressed in Luxi gamecocks including the most significant genes (SLITRK6, IL1RAPL1, GADL1, WISP1, and LRIG3). This study provides more insight into molecular mechanisms of the aggressiveness of gamecocks and aims to promote further studies on animal and human aggression.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BackgroundThe genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.ResultsUsing the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) FST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).ConclusionBased on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Revealing Genetic Variations Associated with Chip-Processing Properties in Potato (Solanum tuberosum L.)

Although the potato chip industry is booming, and distinct chip-processing clones have been released over the past 60 years, the genetic architecture of their chip-processing characteristics remains largely unknown. Case-control genome-wide association studies (GWAS) with SolCAP SNP array data for chip-processing clones versus all other market classes in the 393-line potato diversity panel were performed using the GWASpoly R package, enabling detection of significant signals on chromosome 10. Our results were replicated using internal replication of a strata-corrected 190-line panel. Furthermore, the genomic scans employing selective sweep approaches such as the cross-population composite likelihood ratio method (XP-CLR) and PCAdapt redetected the same signals as those in our GWAS. Through applications of four selective sweep approaches, various genetic variants were found across the genome that had been differentially selected. These genomic regions under selection along with transcriptomic data analysis are involved in carbohydrate metabolism-related genes or loci and transcription factors, indicating to be associated with the improvement of chip-processing performance of potato cultivars. Kompetitive allele-specific PCR (KASP) assays were designed for the causal SNPs to use in validating the chip-processing clones. The results could have implications for genomics-assisted breeding of the promising chip-processing cultivars in potato.

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield.

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker-trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population.

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Sweepstakes reproductive success via pervasive and recurrent selective sweeps.

Highly fecund natural populations characterized by high early mortality abound, yet our knowledge about their recruitment dynamics is somewhat rudimentary. This knowledge gap has implications for our understanding of genetic variation, population connectivity, local adaptation, and the resilience of highly fecund populations. The concept of sweepstakes reproductive success, which posits a considerable variance and skew in individual reproductive output, is key to understanding the distribution of individual reproductive success. However, it still needs to be determined whether highly fecund organisms reproduce through sweepstakes and, if they do, the relative roles of neutral and selective sweepstakes. Here, we use coalescent-based statistical analysis of population genomic data to show that selective sweepstakes likely explain recruitment dynamics in the highly fecund Atlantic cod. We show that the Kingman coalescent (modelling no sweepstakes) and the Xi-Beta coalescent (modelling random sweepstakes), including complex demography and background selection, do not provide an adequate fit for the data. The Durrett-Schweinsberg coalescent, in which selective sweepstakes result from recurrent and pervasive selective sweeps of new mutations, offers greater explanatory power. Our results show that models of sweepstakes reproduction and multiple-merger coalescents are relevant and necessary for understanding genetic diversity in highly fecund natural populations. These findings have fundamental implications for understanding the recruitment variation of fish stocks and general evolutionary genomics of high-fecundity organisms.

Chromosome-level genome assembly and population genomics of Mongolian racerunner (Eremias argus) provide insights into high-altitude adaptation in lizards

BackgroundAlthough the extreme environmental adaptation of organisms is a hot topic in evolutionary biology, genetic adaptation to high-altitude environment remains poorly characterized in ectothermic animals. Squamates are among the most diverse terrestrial vertebrates, with tremendous ecological plasticity and karyotype diversity, and are a unique model system to investigate the genetic footprints of adaptation.ResultsWe report the first chromosome-level assembly of the Mongolian racerunner (Eremias argus) and our comparative genomics analyses found that multiple chromosome fissions/fusions events are unique to lizards. We further sequenced the genomes of 61 Mongolian racerunner individuals that were collected from altitudes ranging from ~ 80 to ~ 2600 m above sea level (m.a.s.l.). Population genomic analyses revealed many novel genomic regions under strong selective sweeps in populations endemic to high altitudes. Genes embedded in those genomic regions are mainly associated with energy metabolism and DNA damage repair pathways. Moreover, we identified and validated two substitutions of PHF14 that may enhance the lizards’ tolerance to hypoxia at high altitudes.ConclusionsOur study reveals the molecular mechanism of high-altitude adaptation in ectothermic animal using lizard as a research subject and provides a high-quality lizard genomic resource for future research.

Genome Polymorphism Analysis and Selected Sweep Regions Detection via the Genome Resequencing of 91 Cabbage (Brassica oleracea) Accessions

The completion of more and more high-quality cabbage genome sequencing attempts provides an important reference genome sequence for resequencing a large proportion of B. oleracea populations. This has laid a foundation for the study of the population diversity of B. oleracea and the excavation of genes related to important agronomic traits. Here, we performed genome resequencing of 91 B. oleracea accessions from 14 countries. We produced a total of 983.84 Gb of high-quality sequencing data, and the average sequencing depth of the genomes was over 15.71 among the 91 accessions. A total of 3,432,341 high-quality SNPs were detected in these B. oleracea accessions. A phylogenetic tree indicated that the leaf ball shapes of “Gaobian” accessions and most “Jianqiu” accessions converge with “Bianqiu” accessions, while “Yuanqiu” accessions were individually assigned to a branch. The principal component analysis (PCA) result was consistent with the phylogenetic relationships. The LD distance of “Yuanqiu” was 38.69 Kb, while that of “Bianqiu” was 30.16 Kb. This result indicated that the linkage degree of the “Yuanqiu” population is greater than that of the “Bianqiu” population, which may indicate that the “Yuanqiu” population has received more positive selection. Our analysis showed that, for the 91 accessions, the best population structure number was four. A total of 560 genes were identified across the 203 selected sweep regions identified in the “Yuanqiu” accessions. Similarly, 682 genes were identified across the 304 selective sweep regions in the “Bianqiu” accessions. Finally, several functional terms were identified via enrichment analysis of the genes in the selective sweep regions. In conclusion, this study will provide rich resources for studies on gene function related to leaf ball development and the population diversity of B. olearcea.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds.

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.

A heterophil/lymphocyte-selected population reveals the phosphatase PTPRJ is associated with immune defense in chickens

Quantification of leukocyte profiles is among the simplest measures of animal immune function. However, the relationship between H/L ratio and innate immunity and the measure’s utility as an index for heterophil function remains to be analyzed. Variants associated with H/L ratio were fine-mapped based on the resequencing of 249 chickens of different generations and an F2 segregating population generated by crossing selection and control lines. H/L ratio in the selection line was associated with a selective sweep of mutations in protein tyrosine phosphatase, receptor type J (PTPRJ), which affects proliferation and differentiation of heterophils through its downstream regulatory genes. The SNP downstream of PTPRJ (rs736799474) have a universal effect on H/L, with CC homozygotes exhibiting improved heterophil function because of downregulated PTPRJ expression. In short, we systematically elucidated the genetic basis of the change in heterophil function resulting from H/L selection by identifying the regulatory gene (PTPRJ) and causative SNP.

Modeling the genetic footprint of fluctuating balancing selection: from the local to the genomic scale.

Natural selection not only affects the actual loci under selection but also leaves "footprints" in patterns of genetic variation in linked genetic regions. This offers exciting opportunities for inferring selection and for understanding the processes shaping levels of genetic variation in natural populations. Here, we develop analytical approximations based on coalescent theory to characterize the genetic footprint of a complex, but potentially common type of natural selection: balancing selection with seasonally fluctuating allele frequencies. As we show analytically and confirm with stochastic simulations, seasonal allele frequency fluctuations can have important (and partly unexpected) consequences for the genetic footprint of balancing selection. Fluctuating balancing selection generally leads to an increase in genetic diversity close to the selected site, the effect of balancing selection, but reduces diversity further away from the selected site, which is a consequence of the allele-frequency fluctuations effectively producing recurrent bottlenecks of allelic backgrounds. This medium- and long-range reduction usually outweighs the short-range increase when averaging diversity levels across the entire chromosome. Strong fluctuating balancing selection even induces a loss of genetic variation in unlinked regions, e.g. on different chromosomes. If many loci in the genome are simultaneously under fluctuating balancing selection this can lead to substantial genome-wide reductions in genetic diversity, even when allele-frequency fluctuations are small and local footprints are difficult to detect. Thus, together with genetic drift, selective sweeps and background selection, fluctuating selection could be a major force shaping levels of genetic diversity in natural populations.

Nonparallel genome changes within subpopulations over time contributed to genetic diversity within the US Holstein population

Maintaining genetic variation in a population is important for long-term genetic gain. The existence of subpopulations within a breed helps maintain genetic variation and diversity. The 20,990 genotyped animals, representing the breeding animals in the year 2014, were identified as the sires of animals born after 2010 with at least 25 progenies, and females measured for type traits within the last 2 yr of data. K-means clustering with 5 clusters (C1, C2, C3, C4, and C5) was applied to the genomic relationship matrix based on 58,990 SNP markers to stratify the selected candidates into subpopulations. The general higher inbreeding resulting from within-cluster mating than across-cluster mating suggests the successful stratification into genetically different groups. The largest cluster (C4) contained animals that were less related to each animal within and across clusters. The average fixation index was 0.03, indicating that the populations were differentiated, and allele differences across the subpopulations were not due to drift alone. Starting with the selected candidates within each cluster, a family unit was identified by tracing back through the pedigree, identifying the genotyped ancestors, and assigning them to a pseudogeneration. Each of the 5 families (F1, F2, F3, F4, and F5) was traced back for 10 generations, allowing for changes in frequency of individual SNPs over time to be observed, which we call allele frequencies change. Alternative procedures were used to identify SNPs changing in a parallel or nonparallel way across families. For example, markers that have changed the most in the whole population, markers that have changed differently across families, and genes previously identified as those that have changed in allele frequency. The genomic trajectory taken by each family involves selective sweeps, polygenic changes, hitchhiking, and epistasis. The replicate frequency spectrum was used to measure the similarity of change across families and showed that populations have changed differently. The proportion of markers that reversed direction in allele frequency change varied from 0.00 to 0.02 if the rate of change was greater than 0.02 per generation, or from 0.14 to 0.24 if the rate of change was greater than 0.005 per generation within each family. Cluster-specific SNP effects for stature were estimated using only females and applied to obtain indirect genomic predictions for males. Reranking occurs depending on SNP effects used. Additive genetic correlations between clusters show possible differences in populations. Further research is required to determine how this knowledge can be applied to maintain diversity and optimize selection decisions in the future.

Genetic responses in sexual diploid and unisexual triploid goldfish (Carassius auratus) introduced into a high-altitude environment.

Anthropogenic biological invasions represent major concerns but enable us to investigate rapid evolutionary changes and adaptation to novel environments. The goldfish Carassius auratus with sexual diploids and unisexual triploids coexisting in natural waters is one of the most widespread invasive fishes in Tibet, providing an ideal model to study evolutionary processes during invasion in different reproductive forms from the same vertebrate. Here, using whole-genome resequencing data of 151 C.auratus individuals from invasive and native ranges, we found different patterns of genomic responses between diploid and triploid populations during their invasion of Tibet. For diploids, although invasive individuals derived from two different genetically distinct sources had a relative higher diversity (π) at the population level, their individual genetic diversity (genome-wide observed heterozygosity) was significantly lower (21.4%) than that of source individuals. Population structure analysis revealed that the invasive individuals formed a specific genetic cluster distinct from the source populations. Runs of homozygosity analysis showed low inbreeding only in invasive individuals, and only the invasive population experienced a recent decline in effective population size reflecting founder events. For triploids, however, invasive populations showed no loss of individual genetic diversity and no genetic differentiation relative to source populations. Regions of putative selective sweeps between invasive and source populations of diploids mainly involved genes associated with mannosidase activity and embryo development. Our results suggest that invasive diploids deriving from distinct sources still lost individual genetic diversity resulting from recent inbreeding and founder events and selective sweeps, and invasive triploids experienced no change in genetic diversity owing to their reproduction mode of gynogenesis that precludes inbreeding and founder effects and may make them more powerful invaders.

Whole-body adipose tissue multi-omic analyses in sheep reveal molecular mechanisms underlying local adaptation to extreme environments

The fat tail of sheep is an important organ that has evolved to adapt to extreme environments. However, the genetic mechanisms underlying the fat tail phenotype remain poorly understood. Here, we characterize transcriptome and lipidome profiles and morphological changes in 250 adipose tissues from two thin-tailed and three fat-tailed sheep populations in summer and winter. We implement whole-genome selective sweep tests to identify genetic variants related to fat-tails. We identify a set of functional genes that show differential expression in the tail fat of fat-tailed and thin-tailed sheep in summer and winter. These genes are significantly enriched in pathways, such as lipid metabolism, extracellular matrix (ECM) remodeling, molecular transport, and inflammatory response. In contrast to thin-tailed sheep, tail fat from fat-tailed sheep show slighter changes in adipocyte size, ECM remodeling, and lipid metabolism, and had less inflammation in response to seasonal changes, indicating improved homeostasis. Whole-genome selective sweep tests identify genes involved in preadipocyte commitment (e.g., BMP2, PDGFD) and terminal adipogenic differentiation (e.g., VEGFA), which could contribute to enhanced adipocyte hyperplasia. Altogether, we establish a model of regulatory networks regulating adipose homeostasis in sheep tails. These findings improve our understanding of how adipose homeostasis is maintained, in response to extreme environments in animals.