Genomic selection for tolerance to aluminum toxicity in a synthetic population of upland rice.

Over half of the world's arable land is acidic, which constrains cereal production. In South America, different rice-growing regions (Cerrado in Brazil and Llanos in Colombia and Venezuela) are particularly affected due to high aluminum toxicity levels. For this reason, efforts have been made to breed for tolerance to aluminum toxicity using synthetic populations. The breeding program of CIAT-CIRAD is a good example of the use of recurrent selection to increase productivity for the Llanos in Colombia. In this study, we evaluated the performance of genomic prediction models to optimize the breeding scheme by hastening the development of an improved synthetic population and elite lines. We characterized 334 families at the S0:4 generation in two conditions. One condition was the control, managed with liming, while the other had high aluminum toxicity. Four traits were considered: days to flowering (FL), plant height (PH), grain yield (YLD), and zinc concentration in the polished grain (ZN). The population presented a high tolerance to aluminum toxicity, with more than 72% of the families showing a higher yield under aluminum conditions. The performance of the families under the aluminum toxicity condition was predicted using four different models: a single-environment model and three multi-environment models. The multi-environment models differed in the way they integrated genotype-by-environment interactions. The best predictive abilities were achieved using multi-environment models: 0.67 for FL, 0.60 for PH, 0.53 for YLD, and 0.65 for ZN. The gain of multi-environment over single-environment models ranged from 71% for YLD to 430% for FL. The selection of the best-performing families based on multi-trait indices, including the four traits mentioned above, facilitated the identification of suitable families for recombination. This information will be used to develop a new cycle of recurrent selection through genomic selection.

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.

Genomic stasis over millions of years in subseafloor sediment.

One of the significant challenges in microbiology is to understand the extent and mechanisms of evolution within life beneath the surface of the Earth. The population bottleneck that microbes in deep marine sediment experience implies that mutational and population genetic forces could lead to higher levels of relaxed selection and an increase in pseudogenes. To investigate this hypothesis, a group of Thalassospira strains were isolated from subseafloor sediment that is 3 to 6 million years old, as reported by Orsi and colleagues in 2021. These isolates, representing lineages that have been buried for millions of years, offer an excellent opportunity to study the evolution of life beneath the seafloor over a long period. The existence of closely related strains from environments on the surface of the Earth enabled us to examine the impact of selection within each group. We discovered that isolates from beneath the seafloor show lineage-specific similarities to Thalassospira from the surface world, both in the overall intensity of selection on the genome and in the specific genes affected by mutation. We found no signs of increased relaxed selection or other notable genomic changes in the genomes of the Thalassospira isolates from beneath the seafloor, suggesting that these subseafloor isolates were awakened from a million-year near-stasis. The unique genomic characteristics of each Thalassospira lineage from beneath the seafloor must then reflect genetic changes that surface-inhabiting decendants acquired in the past 3-6 million years. Remarkably, Thalassospira lineages beneath the surface appear to have stably maintained their genomes in the midst of metabolic dormancy and extremely long generation times.

RRmorph-a new R package to map phenotypic evolutionary rates and patterns on 3D meshes.

The study of evolutionary rates and patterns is the key to understand how natural selection shaped the current and past diversity of phenotypes. Phylogenetic comparative methods offer an array of solutions to undertake this challenging task, and help understanding phenotypic variation in full in most circumstances. However, complex, three-dimensional structures such as the skull and the brain serve disparate goals, and different portions of these phenotypes often fulfil different functions, making it hard to understand which parts truly were recruited by natural selection. In the recent past, we developed tools apt to chart evolutionary rate and patterns directly on three-dimensional shapes, according to their magnitude and direction. Here, we present further developments of these tools, which now allow to restitute the mapping of rates and patterns with full biological realism. The tools are condensed in a new R software package.

Mitochondrial, nuclear and morphological differentiation in the swimming crab Liocarcinus depurator along the Atlantic-Mediterranean transition.

Environmental gradients in the sea may coincide with phenotypic or genetic gradients resulting from an evolutionary balance between selection and dispersal. The population differentiation of the swimming crab, Liocarcinus depurator, an important by-catch species in the Mediterranean Sea and North-East Atlantic, was assessed using both genetic and morphometric approaches. A total of 472 specimens were collected along its distribution area, and 17 morphometric landmarks, one mitochondrial gene (COI) and 11 polymorphic microsatellite markers were scored in 350, 287 and 280 individuals, respectively. Morphometric data lacked significant differences, but genetic analyses showed significant genetic differentiation between Atlantic and Mediterranean populations, with a steeper gradient in COI compared to microsatellite markers. Interestingly, nuclear differentiation was due to an outlier locus with a gradient in the Atlantic-Mediterranean transition area overlapping with the mtDNA gradient. Such overlapping clines are likely to be maintained by natural selection. Our results suggest a scenario of past isolation with local adaptation and secondary contact between the two basins. Local adaptation during the process of vicariance may reinforce genetic differentiation at loci maintained by environmental selection even after secondary contact.

Diversifying selection identified in immune epitopes of bovine coronavirus isolates from Irish cattle.

Bovine betacoronavirus (BoCoV) is a pneumoenteric pathogen of cattle that is closely related to human coronavirus OC43. Vaccines are administered to protect against diseases caused by BoCoV, but knowledge gaps exist with regard to correlates of protection and the effect of immune evasion on driving evolution. In this study, immune epitopes were mapped onto BoCoV structural proteins, including spike and haemagglutinin esterase (HE), and then supported with targeted gene sequencing of Irish clinical isolates and selective pressure analysis. Increased prevalence of diversifying selection and amino acid changes in some mapped immune epitopes suggests that immune escape is selecting for non-synonymous mutations arising in these regions. Selection analysis and sequencing provided increased support for neutralising antibody (nAb) epitopes compared to others, suggesting that nAbs are an important arm of the immune response to BoCoV. Phylogenetic analysis of spike and HE sequences showed that Irish isolates from this study were in the European clade, except for one HE sequence that sat in the Asian/American clade, while the spike gene of this sample was in the European clade. Recombination between a European and an Asian/American isolate would give rise to such a sequence. This study has gathered evidence suggesting that pressure to evade the nAb response is contributing to BoCoV evolution.

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

BACKGROUND: Adaptation 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. METHODS: BAT 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. RESULTS: Of 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. CONCLUSIONS: Our 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.

Genetic divergence and truncation and simultaneous selection in inbred families (S1) of elephant grass for bioenergetic purposes via mixed models.

The State University of North Fluminense Darcy Ribeiro (UENF) has been developing for fifteen years a breeding program that aims at the development of new cultivars of elephant grass due to its high potential and the low availability of cultivars developed by genetic breeding programs that meet the needs of producers in the State of Rio de Janeiro. In this sense, inbred families were also obtained as a way of fixing potential alleles for traits related to production, as the inbreeding process apparently does not strongly affect elephant grass in aspects related to inbreeding depression. This study aimed to estimate genetic diversity, variance components and prediction of genotypic values in 11 (S1) elephant grass families, and perform the truncation and simultaneous selection of traits using the selection index, by mixed models. The experimental design consisted of randomized blocks with 11 (S1) families, three replications, and six plants per plot. For variables dry matter production, percentage of dry matter, plant height, stem diameter, number of tillers and leaf blade width, was performed the estimation of genetic parameters and selection of the best genotypes based selection index using mixed model. The descriptors were subjected to correlation analysis, distance matrices were generated by the Mahalanobis method, and individuals were grouped by the UPGMA method. In the selection via mixed models (REML/BLUP), families 6, 11, 8, 1, 3, 7, and 9 contributed most of the genotypes selected for the evaluated traits, indicating their high potential to generate superior genotype. The selection indices via mixed models indicated that the multiplicative index presented a greater selection gain.The phenotypic correlations showed the possibility of performing an indirect selection from six evaluated traits.The genotypes were separated into 18 groups by the Mahalanobis distance, allowing the observation of a wide genetic diversity. The most divergent and productive genotypes were self-fertilized to obtain the second generation (S2), continuing the development program.

Selection for stress tolerance and longevity in Drosophila melanogaster have strong impacts on microbiome profiles.

There is experimental evidence that microbiomes have a strong influence on a range of host traits. Understanding the basis and importance of symbiosis between host and associated microorganisms is a rapidly developing research field, and we still lack a mechanistic understanding of ecological and genetic pressures affecting host-microbiome associations. Here Drosophila melanogaster lines from a large-scale artificial selection experiment were used to investigate whether the microbiota differ in lines selected for different stress resistance traits and longevity. Following multiple generations of artificial selection all selection regimes and corresponding controls had their microbiomes assessed. The microbiome was interrogated based on 16S rRNA sequencing. We found that the microbiome of flies from the different selection regimes differed markedly from that of the unselected control regime, and microbial diversity was consistently higher in selected relative to control regimes. Several common Drosophila bacterial species showed differentially abundance in the different selection regimes despite flies being exposed to similar environmental conditions for two generations prior to assessment. Our findings provide strong evidence for symbiosis between host and microbiomes but we cannot reveal whether the interactions are adaptive, nor whether they are caused by genetic or ecological factors.

The Ruminant Telomere-to-Telomere (RT2T) Consortium.

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

Impact of truncating historical data on prediction ability of dairy sheep selection candidates.

Along the last decades, the genetic evaluation methodology has evolved, improving breeding value estimates. Many breeding programmes have historical phenotypic records and large number of generations, but to make use of them could result in more inconveniences than benefits. In this study, the prediction ability of genotyped young animals was assessed by simultaneously evaluating the removal of historical data, two pedigree deepness and two methodologies (traditional BLUP and single-step genomic BLUP or ssGBLUP), using milk yield records of 40 years of three Latxa dairy sheep populations. The linear regression method was used to compare predictions of young rams before and after progeny testing, with six cut-off points, by intervals of 4 years (from 1992 to 2012), and statistics of ratio of accuracies, bias, and dispersion were calculated. The prediction accuracy of selection candidates, when genomic information was included, was the highest in all Latxa populations (between 0.54 and 0.69 with full data set). Nevertheless, the deletion of historical phenotypic data resulted on moderate accuracy gain in the bigger data size populations (mean gain 2.5%), and the smaller population took advantage of a moderate data deletion (2.7% gain by removing data until 2004), reducing accuracy when more records were removed. The bias of validation individuals was lower when the breeding value was predicted based on genomic information (between 2.1 and 13.9), being lower when the biggest amount of data was deleted in the bigger data size populations (5.2% reduction), and the smaller population was benefited from data deletion between 1996 and 2008 (3.8% bias reduction). Meanwhile, the slope of estimated genetic trend was lower when less data were included, and an overestimation of the unknown parent group estimates was observed. The results indicated that ssGBLUP evaluations were outstanding, compared with traditional BLUP evaluations, while the depth of pedigree had a very small influence, and deletion of historical phenotypic data was beneficial. Thus, Latxa routine genetic evaluations would benefit from truncating phenotypic records between 2000 and 2004, the use of two pedigree generations and the implementation of ssGBLUP methodology.

Whole-genome scan for selection signature associated with temperature adaptation in Iranian sheep breeds.

The present study aimed to identify the selection signature associated with temperature adaptation in Iranian sheep breeds raised in cold and hot environments. The Illumina HD ovine SNP600K BeadChip genomic arrays were utilized to analyze 114 animals from eight Iranian sheep breeds, namely Ghezel, Afshari, Shall, Sanjabi, Lori-Bakhtiari, Karakul, Kermani, and Balochi. All animals were classified into two groups: cold-weather breeds and hot-weather breeds, based on the environments to which they are adapted and the regions where they have been raised for many years. The unbiased FST (Theta) and hapFLK tests were used to identify the selection signatures. The results revealed five genomic regions on chromosomes 2, 10, 11, 13, and 14 using the FST test, and three genomic regions on chromosomes 10, 14, and 15 using the hapFLK test to be under selection in cold and hot groups. Further exploration of these genomic regions revealed that most of these regions overlapped with genes previously identified to affect cold and heat stress, nervous system function, cell division and gene expression, skin growth and development, embryo and skeletal development, adaptation to hypoxia conditions, and the immune system. These regions overlapped with QTLs that had previously been identified as being associated with various important economic traits, such as body weight, skin color, and horn characteristics. The gene ontology and gene network analyses revealed significant pathways and networks that distinguished Iranian cold and hot climates sheep breeds from each other. We identified positively selected genomic regions in Iranian sheep associated with pathways related to cell division, biological processes, cellular responses to calcium ions, metal ions and inorganic substances. This study represents the initial effort to identify selective sweeps linked to temperature adaptation in Iranian indigenous sheep breeds. It may provide valuable insights into the genomic regions involved in climate adaptation in sheep.

Varying Selection Pressure for a Na+ Sensing Site in Epithelial Na+ Channel Subunits Reflect Divergent Roles in Na+ Homeostasis.

The epithelial Na+ channel (ENaC) emerged early in vertebrates and has played a role in Na+ and fluid homeostasis throughout vertebrate evolution. We previously showed that proteolytic activation of the channel evolved at the water-to-land transition of vertebrates. Sensitivity to extracellular Na+, known as Na+ self-inhibition, reduces ENaC function when Na+ concentrations are high and is a distinctive feature of the channel. A fourth ENaC subunit, δ, emerged in jawed fishes from an α subunit gene duplication. Here, we analyzed 849 α and δ subunit sequences and found that a key Asp in a postulated Na+ binding site was nearly always present in the α subunit, but frequently lost in the δ subunit (e.g. human). Analysis of site evolution and codon substitution rates provide evidence that the ancestral α subunit had the site and that purifying selection for the site relaxed in the δ subunit after its divergence from the α subunit, coinciding with a loss of δ subunit expression in renal tissues. We also show that the proposed Na+ binding site in the α subunit is a bona fide site by conferring novel function to channels comprising human δ subunits. Together, our findings provide evidence that ENaC Na+ self-inhibition improves fitness through its role in Na+ homeostasis in vertebrates.

Using viral diversity to identify HIV-1 variants under HLA-dependent selection in a systematic viral genome-wide screen.

The pathogenesis of HIV-1 infection is governed by a highly dynamic, time-dependent interaction between the host and the viral genome. In this study, we developed a novel systematic approach to assess the host-virus interaction, using average pairwise viral diversity as a proxy for time since infection, and applied this method to nearly whole viral genome sequences (n = 4,464), human leukocyte antigen (HLA) genotyping data (n = 1,044), and viral RNA load (VL) measurements during the untreated chronic phase (n = 829) of Swiss HIV Cohort Study participants. Our systematic genome-wide screen revealed for 98 HLA/viral-variant pairs a signature of immune-driven selection in the form of an HLA-dependent effect of infection time on the presence of HIV amino acid variants. Of these pairs, 12 were found to have an effect on VL. Furthermore, 28/58 pairs were validated by time-to-event analyses and 48/92 by computational HLA-epitope predictions. Our diversity-based approach allows a powerful and systematic investigation of the interaction between the virus and cellular immunity, revealing a notable subset of such interaction effects. From an evolutionary perspective, these observations underscore the complexity of HLA-mediated selection pressures on the virus that shape viral evolution and pathogenesis.

Into the Blue: Exploring genetic mechanisms behind the evolution of baleen whales.

Marine ecosystems are ideal for studying evolutionary adaptations involved in lineage diversification due to few physical barriers and reduced opportunities for strict allopatry compared to terrestrial ecosystems. Cetaceans (whales, dolphins, and porpoises) are a diverse group of mammals that successfully adapted to various habitats within the aquatic environment around 50 million years ago. While the overall adaptive transition from terrestrial to fully aquatic species is relatively well understood, the radiation of modern whales is still unclear. Here high-quality genomes derived from previously published data were used to identify genomic regions that potentially underpinned the diversification of baleen whales (Balaenopteridae). A robust molecular phylogeny was reconstructed based on 10,159 single copy and complete genes for eight mysticetes, seven odontocetes and two cetacean outgroups. Analysis of positive selection across 3,150 genes revealed that balaenopterids have undergone numerous idiosyncratic and convergent genomic variations that may explain their diversification. Genes associated with aging, survival and homeostasis were enriched in all species. Additionally, positive selection on genes involved in the immune system were disclosed for the two largest species, blue and fin whales. Such genes can potentially be ascribed to their morphological evolution, allowing them to attain greater length and increased cell number. Further evidence is presented about gene regions that might have contributed to the extensive anatomical changes shown by cetaceans, including adaptation to distinct environments and diets. This study contributes to our understanding of the genomic basis of diversification in baleen whales and the molecular changes linked to their adaptive radiation, thereby enhancing our understanding of cetacean evolution.