Signatures Of Positive Selection Research Articles

Complete organelle genomes of the threatened aquatic species Scheuchzeria palustris (Scheuchzeriaceae): Insights into adaptation and phylogenomic placement.

Scheuchzeria palustris, the only species in the Scheuchzeriaceae family, plays a crucial role in methane production and transportation, influencing the global carbon cycle and maintaining ecosystem stability. However, it is now threatened by human activities and global warming. In this study, we generated new organelle genomes for S. palustris, with the plastome (pt) measuring 158,573 bp and the mitogenome (mt) measuring 420,724 bp. We predicted 296 RNA editing sites in mt protein-coding genes (PCGs) and 142 in pt-PCGs. Notably, abundant RNA editing sites in pt-PCGs likely originated from horizontal gene transfer between the plastome and mitogenome. Additionally, we identified positive selection signals in four mt-PCGs (atp4, ccmB, nad3, and sdh4) and one pt-PCG (rps7), which may contribute to the adaptation of S. palustris to low-temperature and high-altitude environments. Furthermore, we identified 35 mitochondrial plastid DNA (MTPT) segments totaling 58,479 bp, attributed to dispersed repeats near most MTPT. Phylogenetic trees reconstructed from mt- and pt-PCGs showed topologies consistent with the APG IV system. However, the conflicting position of S. palustris can be explained by significant differences in the substitution rates of its mt- and pt-PCGs (p < .001). In conclusion, our study provides vital genomic resources to support future conservation efforts and explores the adaptation mechanisms of S. palustris.

A Genetic Study of Spillovers in the Bean Common Mosaic Subgroup of Potyviruses.

Nine viruses of the bean common mosaic virus subgroup of potyviruses are major international crop pathogens, but their phylogenetically closest relatives from non-crop plants have mostly been found only in SE Asia and Oceania, which is thus likely to be their "centre of emergence". We have compared over 700 of the complete genomic ORFs of the crop pandemic and the non-crop viruses in various ways. Only one-third of crop virus genomes are non-recombinant, but more than half the non-crop virus genomes are. Four of the viruses were from crops domesticated in the Old World (Africa to SE Asia), and the other five were from New World crops. There was a temporal signal in only three of the crop virus datasets, but it confirmed that the most recent common ancestors of all the crop viruses were before inter-continental marine trade started after 1492 CE, whereas all the crown clusters of the phylogenies are from after that date. The non-crop virus datasets are genetically more diverse than those of the crop viruses, and Tajima's D analyses showed that their populations were contracting, and only one of the crop viruses had a significantly expanding population. dN/dS analyses showed that most of the genes and codons in all the viruses were under significant negative selection, and the few that were under significant positive selection were mostly in the PIPO-encoding region of the P3 protein, or the PIPO protein itself. Interestingly, more positively selected codons were found in non-crop than in crop viruses, and, as the hosts of the former were taxonomically more diverse than the latter, this may indicate that the positively selected codons are involved in host range determination; AlphaFold3 modelling was used to investigate this possibility.

Association between thermogenic brown fat and genes under positive natural selection in circumpolar populations

BackgroundAdaptation to cold was essential for human migration across Eurasia. Non-shivering thermogenesis through brown adipose tissue (BAT) participates in cold adaptation because some genes involved in the differentiation and function of BAT exhibit signatures of positive natural selection in populations at high latitudes. Whether these genes are associated with the inter-individual variability in BAT thermogenesis remains unclear. In this study, we evaluated the potential associations between BAT activity and single nucleotide polymorphisms (SNPs) in candidate gene regions in East Asian populations.MethodsBAT activity induced by mild cold exposure was measured in 399 healthy Japanese men and women using fluorodeoxyglucose-positron emission tomography and computed tomography (FDG-PET/CT). The capacity for cold-induced thermogenesis and fat oxidation was measured in 56 men. Association analyses with physiological traits were performed for 11 SNPs at six loci (LEPR, ANGPTL8, PLA2G2A, PLIN1, TBX15-WARS2, and FADS1) reported to be under positive natural selection. Associations found in the FDG-PET/CT population were further validated in 84 healthy East Asian men and women, in whom BAT activity was measured using infrared thermography. Associations between the SNP genotypes and BAT activity or other related traits were tested using multiple logistic and linear regression models.ResultsOf the 11 putative adaptive alleles of the six genes, two intronic SNPs in LEPR (rs1022981 and rs12405556) tended to be associated with higher BAT activity. However, these did not survive multiple test comparisons. Associations with lower body fat percentage, plasma triglyceride, insulin, and HOMA-IR levels were observed in the FDG-PET/CT population (P < 0.05). Other loci, including TBX15-WARS2, which is speculated to mediate cold adaptation in Greenland Inuits, did not show significant differences in BAT thermogenesis.ConclusionsOur results suggest a marginal but significant association between LEPR SNPs. However, robust supporting evidence was not established for the involvement of other loci under positive natural selection in cold adaptation through BAT thermogenesis in East Asian adults. Given the pleiotropic function of these genes, factors other than cold adaptation through BAT thermogenesis, such as diet adaptation, may contribute to positive natural selection at these loci.

Comparative Genomics Supports Ecologically Induced Selection as a Putative Driver of Banded Penguin Diversification.

The relative importance of genetic drift and local adaptation in facilitating speciation remains unclear. This is particularly true for seabirds, which can disperse over large geographic distances, providing opportunities for intermittent gene flow among distant colonies that span the temperature and salinity gradients of the oceans. Here, we delve into the genomic basis of adaptation and speciation of banded penguins, Galápagos (Spheniscus mendiculus), Humboldt (Spheniscus humboldti), Magellanic (Spheniscus magellanicus), and African penguins (Spheniscus demersus), by analyzing 114 genomes from the main 16 breeding colonies. We aim to identify the molecular mechanism and genomic adaptive traits that have facilitated their diversifications. Through positive selection and gene family expansion analyses, we identified candidate genes that may be related to reproductive isolation processes mediated by ecological thermal niche divergence. We recover signals of positive selection on key loci associated with spermatogenesis, especially during the recent peripatric divergence of the Galápagos penguin from the Humboldt penguin. High temperatures in tropical habitats may have favored selection on loci associated with spermatogenesis to maintain sperm viability, leading to reproductive isolation among young species. Our results suggest that genome-wide selection on loci associated with molecular pathways that underpin thermoregulation, osmoregulation, hypoxia, and social behavior appears to have been crucial in local adaptation of banded penguins. Overall, these results contribute to our understanding of how the complexity of biotic, but especially abiotic, factors, along with the high dispersal capabilities of these marine species, may promote both neutral and adaptive lineage divergence even in the presence of gene flow.

Positive selection in cilia-related genes may facilitate deep-sea adaptation of Thermocollonia jamsteci

Deep-sea hydrothermal vents are characterized by high hydrostatic pressure, hypoxia, darkness and toxic substances. However, how organisms adapt to such extreme marine ecosystems remain poorly understood. We hypothesize that adaptive evolution plays an essential role in generating novelty for evolutionary adaptation to the deep-sea environment because adaptive evolution has been found to be critical for species origin and evolution. In this project, the chromosome-level genome of the deep-sea hydrothermal vent gastropod T. jamsteci was constructed for the first time to examine molecular mechanisms of its adaptation to the deep-sea environment. The genome size was large (2.54 Gb), ranking at the top of all species in the Vetigastropoda subclass, driven primarily by the bursts of transposable elements (TEs). The transposition of TEs may also trigger chromosomal changes including both inter-chromosomal fusions and intra-chromosomal activities involving chromosome inversions, rearrangements and fusions, as revealed by comparing the genomes of T. jamsteci and its closely related shallow-sea species Gibbula magus. Innovative changes including the expansion of the ABC transporter gene family that may facilitate detoxification, duplication of genes related to endocytosis, immunity, apoptosis, and anti-apoptotic domains that may help T. jamsteci fight against microbial pathogens, were identified. Furthermore, comparative analysis identified positive selection signals in a large number of genes including the hypoxia up-regulated protein 1, which is a chaperone that may promote adaptation of the T. jamsteci to hypoxic deepsea environments, hox2, Rx2, Pax6 and cilia-related genes BBS1, BBS2, BBS9 and RFX4. Notably, because of the critical importance of cilia and IFT in development, positive selection in cilia-related genes may play a critical role in facilitating T. jamsteci to adapt to the high-pressure deep-sea ecosystem. Results from this study thus revealed important molecular clues that may facilitate further research on the adaptation of molluscs to deep-sea hydrothermal vents.

Genetic Signatures of Positive Selection in Human Populations Adapted to High Altitude in Papua New Guinea.

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

The complete chloroplast genome of Rhododendronambiguum and comparative genomics of related species.

Rhododendron Linnaeus, 1753, the largest genus of woody plants in the Northern Hemisphere, includes some of the most significant species in horticulture. Rhododendronambiguum Hemsl, 1911, a member of subsection Triflora Sleumer 1947, exemplifies typical alpine Rhododendron species. The analysis of the complete chloroplast genome of R.ambiguum offers new insights into the evolution of Rhododendron species and enhances the resolution of phylogenetic relationships. This genome is composed of 207,478 base pairs, including a pair of inverted repeats (IRs) of 47,249 bp each, separated by a large single-copy (LSC) region of 110,367 bp and a small single-copy (SSC) region of 2,613 bp. It contains 110 genes: 77 protein-coding genes, 29 tRNAs, four unique rRNAs (4.5S, 5S, 16S, and 23S), with 16 genes duplicated in the IRs. Comparative analyses reveal substantial diversity in the Rhododendron chloroplast genome structures, identifying a fourth variant pattern. Specifically, four highly divergent regions (trnI-rpoB, ndhE-psaC, rpl32-ndhF, rrn16S-trnI) were noted in the intergenic spacers. Additionally, 76 simple sequence repeats were identified. Positive selection signals were detected in four genes (cemA, rps4, rpl16, and rpl14), evidenced by high Ka/Ks ratios. Phylogenetic reconstruction based on two datasets (shared protein-coding genes and complete chloroplast genomes) suggests that R.ambiguum is closely related to R.concinnum Hemsley, 1889. However, the phylogenetic positions of subsection Triflora Pojarkova, 1952 species remain unresolved, indicating that the use of complete chloroplast genomes for phylogenetic research in Rhododendron requires careful consideration. Overall, our findings provide valuable genetic information that will enhance understanding of the evolution, molecular biology, and genetic improvement of Rhododendron spieces.

Moderate Genetic Diversity of MHC Genes in an Isolated Small Population of Black-and-White Snub-Nosed Monkeys (Rhinopithecus bieti).

Genetic diversity is an essential indicator that echoes the natural selection and environmental adaptation of a species. Isolated small populations are vulnerable to genetic drift, inbreeding, and limited gene flow; thus, assessing their genetic diversity is critical in conservation. In this study, we studied the genetic diversity of black-and-white snub-nosed monkeys (Rhinopithecus bieti) using neutral microsatellites and five adaptive major histocompatibility complex (MHC) genes. Two DQA1 alleles, two DQB1 alleles, two DRB1 alleles, two DRB5 alleles, and three DPB1 alleles were isolated from a population. The results indicate that neutral microsatellites demonstrate a high degree of heterozygosity and polymorphism, while adaptive MHC genes display a high degree of heterozygosity and moderate polymorphism. The results also show that balancing selection has prominently influenced the MHC diversity of the species during evolution: (1) significant positive selection is identified at several amino acid sites (primarily at and near antigen-binding sites) of the DRB1, DRB5, and DQB1 genes; (2) phylogenetic analyses display the patterns of trans-species evolution for all MHC loci. This study provides valuable genetic diversity insights into black-and-white snub-nosed monkeys, which dwell at the highest altitude and have experienced the harshest environmental selection of all primates globally since the Pleistocene. Such results provide valuable scientific evidence and a reference for making or amending conservation strategies for this endangered primate species.

Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

Genetic characterization of Tibetan pigs adapted to high altitude under natural selection based on a large whole-genome dataset

The Qinghai-Tibet Plateau is a valuable genetic resource pool, and the high-altitude adaptation of Tibetan pigs is a classic example of the adaptive evolution of domestic animals. Here, we report the presence of Darwinian positive selection signatures in Tibetan pigs (TBPs) using 348 genome-wide datasets (127 whole-genome sequence datasets (WGSs) and 221 whole-genome single-nucleotide polymorphism (SNP) chip datasets). We characterized a high-confidence list of genetic signatures related response to high-altitude adaptation in Tibetan pigs, including 4,598 candidate SNPs and 131 candidate genes. Functional annotation and enrichment analysis revealed that 131 candidate genes are related to multiple systems and organs in Tibetan pigs. Notably, eight of the top ten novel genes, RALB, NBEA, LIFR, CLEC17A, PRIM2, CDH7, GK5 and FAM83B, were highlighted and associated with improved adaptive heart functions in Tibetan pigs high-altitude adaptation. Moreover, genome-wide association analysis revealed that 29 SNPs were involved in 13 candidate genes associated with at least one adaptive trait. In particular, among the top ten candidate genes, CLEC17A is related to a reduction in hemoglobin (HGB) in Tibetan pigs. Overall, our study provides a robust SNP/gene list involving genetic adaptation for Tibetan pig high-altitude adaptation, and it will be a valuable resource for future Tibetan pig studies.

Evolution of ion channels in cetaceans: a natural experiment in the tree of life

Cetaceans represent a natural experiment within the tree of life in which a lineage changed from terrestrial to aquatic habitats. This shift involved phenotypic modifications, representing an opportunity to explore the genetic bases of phenotypic diversity. Among the different molecular systems that maintain cellular homeostasis, ion channels are crucial for the proper physiological functioning of all living species. This study aims to explore the evolution of ion channels during the evolutionary history of cetaceans. To do so, we created a bioinformatic pipeline to annotate the repertoire of ion channels in the genome of the species included in our sampling. Our main results show that cetaceans have, on average, fewer protein-coding genes and a higher percentage of annotated ion channels than non-cetacean mammals. Signals of positive selection were detected in ion channels related to the heart, locomotion, visual and neurological phenotypes. Interestingly, we predict that the NaV1.5 ion channel of most toothed whales (odontocetes) is sensitive to tetrodotoxin, similar to NaV1.7, given the presence of tyrosine instead of cysteine, in a specific position of the ion channel. Finally, the gene turnover rate of the cetacean crown group is more than three times faster than that of non-cetacean mammals.

Copy number losses of oncogenes and gains of tumor suppressor genes generate common driver mutations

Cancer driver genes can undergo positive selection for various types of genetic alterations, including gain-of-function or loss-of-function mutations and copy number alterations (CNA). We investigated the landscape of different types of alterations affecting driver genes in 17,644 cancer exomes and genomes. We find that oncogenes may simultaneously exhibit signatures of positive selection and also negative selection in different gene segments, suggesting a method to identify additional tumor types where an oncogene is a driver or a vulnerability. Next, we characterize the landscape of CNA-dependent selection effects, revealing a general trend of increased positive selection on oncogene mutations not only upon CNA gains but also upon CNA deletions. Similarly, we observe a positive interaction between mutations and CNA gains in tumor suppressor genes. Thus, two-hit events involving point mutations and CNA are universally observed regardless of the type of CNA and may signal new therapeutic opportunities. An analysis with focus on the somatic CNA two-hit events can help identify additional driver genes relevant to a tumor type. By a global inference of point mutation and CNA selection signatures and interactions thereof across genes and tissues, we identify 9 evolutionary archetypes of driver genes, representing different mechanisms of (in)activation by genetic alterations.

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining Plasmodium falciparum Genomes.

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is increasingly being used in Plasmodium population genomics to estimate genetic relatedness, effective population size (N e ), population structure, and signals of positive selection. Despite its potential, a thorough evaluation of IBD segment detection tools for species with high recombination rates, such as P. falciparum, remains absent. Here, we perform comprehensive benchmarking of IBD callers - probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD) - using population genetic simulations tailored for high recombination, and IBD quality metrics at both the IBD segment level and the IBD-based downstream inference level. Our results demonstrate that low marker density per genetic unit, related to high recombination relative to mutation, significantly compromises the accuracy of detected IBD segments. In genomes with high recombination rates resembling P. falciparum, most IBD callers exhibit high false negative rates for shorter IBD segments, which can be partially mitigated through optimization of IBD caller parameters, especially those related to marker density. Notably, IBD detected with optimized parameters allows for more accurate capture of selection signals and population structure; IBD-based N e inference is very sensitive to IBD detection errors, with IBD called from hmmIBD uniquely providing less biased estimates of N e in this context. Validation with empirical data from the MalariaGEN Pf 7 database, representing different transmission settings, corroborates these findings. We conclude that context-specific evaluation and parameter optimization are essential for accurate IBD detection in high-recombining species and recommend hmmIBD for quality-sensitive analysis, such as estimation of N e in these species. Our optimization and high-level benchmarking methods not only improve IBD segment detection in high-recombining genomes but also enhance overall genomic analysis, paving the way for more accurate genomic surveillance and targeted intervention strategies for malaria.

Structural and genetic diversity in the secreted mucins MUC5AC and MUC5B

The secreted mucins MUC5AC and MUC5B are large glycoproteins that play critical defensive roles in pathogen entrapment and mucociliary clearance. Their respective genes contain polymorphic and degenerate protein-coding variable number tandem repeats (VNTRs) that make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5,761-5,762 amino acids [aa]); however, seven haplotypes have expanded VNTRs (6,291-7,019 aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5,249-6,325 aa) with cysteine-rich domain and VNTR copy-number variation. We group MUC5AC alleles into three phylogenetic clades: H1 (46%, ∼5,654 aa), H2 (33%, ∼5,742 aa), and H3 (7%, ∼6,325 aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium and Tajima's D analyses reveal that East Asians carry exceptionally large blocks with an excess of rare variation (p<0.05) at MUC5AC. To validate this result, we use Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observe a signature of positive selection in H1 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium (p<0.05), consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein-coding VNTRs for improved disease associations.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Long-range linkage disequilibrium events on the genome of dromedary camels as a signal of epistatic and directional positive selection

The genome of dromedary camels has been subjected to various evolutionary forces, such as genetic admixture, natural positive selection, and epistatic selection. These forces are considered as main factors associated with the formation of long-range linkage disequilibrium (LRLD) events. We have analyzed whole-genome data of 56 dromedary camel samples from different geographical regions across the Arabian Peninsula for two main purposes: first, to assess the level of linkage disequilibrium, and second, to identify autosomal LRLD events. The analysis revealed a mean r2 value of 0.25 (±0.028) over the dromedary autosomes, with a continuous decay until reaching a plateau at inter-variant distances >400 kb. A total of 1847 LRLD events were identified within the dromedary autosomes, which harbor 36 prevalent haplotypes. A level of genetic admixture was observed among the dromedary populations analyzed, which might be a source for the observed LRLD events. Four functional interactions were revealed among the genes found within the LRLD events, with some genes overlapping with prevalent haplotypes, indicative of potential epistatic selection. Genes related to renal function, fertility, thermal regulation, bone structure, and insulin regulation were found among the LRLD genes. These genes, along with the defined prevalent haplotypes, can be considered as hotspots for natural positive selection associated with the LRLD distribution on dromedary genomes. In this study, we have for the first time analyzed the genome of dromedary camels for LRLD events possibly influenced by forces including genetic admixture, epistatic and positive selection. The revealed LRLD elements and prevalent haplotypes should be accounted for when designing breeding programmes to conserve the genetic stock of this well-adapted domestic species.

Exploring selection signatures in the divergence and evolution of lipid droplet (LD) associated genes in major oilseed crops

BackgroundOil bodies or lipid droplets (LDs) in the cytosol are the subcellular storage compartments of seeds and the sites of lipid metabolism providing energy to the germinating seeds. Major LD-associated proteins are lipoxygenases, phospholipaseD, oleosins, TAG-lipases, steroleosins, caleosins and SEIPINs; involved in facilitating germination and enhancing peroxidation resulting in off-flavours. However, how natural selection is balancing contradictory processes in lipid-rich seeds remains evasive. The present study was aimed at the prediction of selection signatures among orthologous clades in major oilseeds and the correlation of selection effect with gene expression.ResultsThe LD-associated genes from the major oil-bearing crops were analyzed to predict natural selection signatures in phylogenetically close-knit ortholog clusters to understand adaptive evolution. Positive selection was the major force driving the evolution and diversification of orthologs in a lineage-specific manner. Significant positive selection effects were found in 94 genes particularly in oleosin and TAG-lipases, purifying with excess of non-synonymous substitution in 44 genes while 35 genes were neutral to selection effects. No significant selection impact was noticed in Brassicaceae as against LOX genes of oil palm. A heavy load of deleterious mutations affecting selection signatures was detected in T-lineage oleosins and LOX genes of Arachis hypogaea. The T-lineage oleosin genes were involved in mainly anther, tapetum and anther wall morphogenesis. In Ricinus communis and Sesamum indicum > 85% of PLD genes were under selection whereas selection pressures were low in Brassica juncea and Helianthus annuus. Steroleosin, caleosin and SEIPINs with large roles in lipid droplet organization expressed mostly in seeds and were under considerable positive selection pressures. Expression divergence was evident among paralogs and homeologs with one gene attaining functional superiority compared to the other. The LOX gene Glyma.13g347500 associated with off-flavor was not expressed during germination, rather its paralog Glyma.13g347600 showed expression in Glycine max. PLD-α genes were expressed on all the tissues except the seed,δ genes in seed and meristem while β and γ genes expressed in the leaf.ConclusionsThe genes involved in seed germination and lipid metabolism were under strong positive selection, although species differences were discernable. The present study identifies suitable candidate genes enhancing seed oil content and germination wherein directional selection can become more fruitful.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets.

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

Genetic and developmental divergence in the neural crest programme between cichlid fish species.

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.