A genomics approach reveals insights into the importance of gene losses for mammalian adaptations

Evolutionary repeatability of chromatin accessibility in the hippocampus between echolocating bats and soft-furred tree mice

Echolocation is a compelling example of behavioral convergence in mammals and has evolved independently in at least five mammalian lineages. Despite increasing evidence of convergent evolution in protein-coding regions among echolocating mammals, the evolution of gene regulatory activity underlying mammalian echolocation remains undetermined. Here, we systematically investigate chromatin accessibility and transcriptomic data from the hippocampi of echolocating bats and soft-furred tree mice, as the hippocampus plays a crucial role in processing echolocation signals. Our results showed significant evolutionary repeatability in accessible chromatin regions among echolocating mammals compared with their non-echolocating counterparts. An analysis that paired gene expression with chromatin accessibility revealed that echolocating mammals have more complex gene regulatory networks in the hippocampus than non-echolocating mammals. This complexity is primarily attributed to the observed evolutionary repeatability in accessible chromatin regions among echolocating mammals. The shared gene regulations among echolocating mammals are particularly enriched in functional terms related to synaptic function, which is supported by the higher synaptic and mitochondrial densities in soft-furred tree mice than in laboratory mice. These findings demonstrate significant evolutionary repeatability of gene regulatory activity in the hippocampus among echolocating mammalian species, suggesting crucial contributions of gene regulatory activities to convergent phenotypes in mammals.

Intragenomic mutational heterogeneity: structural and functional insights from gene evolution

Variation of mutation rates between species has been documented over decades, but the variation between different regions of a genome has been less often discussed. Recent studies using high-quality sequence data have revealed previously unknown levels of intragenomic heterogeneity of mutation rates and their association with other structural and functional features of DNA sequences. This article reviews accumulating evidence of this intragenomic heterogeneity and speculates its cause and influence on organismal phenotypes.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Genomic signatures associated with the evolutionary loss of egg yolk in parasitoid wasps

Trait regression and loss have occurred repeatedly in numerous lineages in response to environmental changes. In parasitoid wasps, a megadiverse group of hymenopteran insects, yolk protein reduction or loss has been observed in many species, likely linked to the transition from ectoparasitism to endoparasitism. However, the genetic basis of this trait and the impact of its loss on genome evolution remain poorly understood. Here, we performed a comparative genomic analysis of 64 hymenopteran insects. The conserved insect yolk protein gene vitellogenin (Vg) underwent five independent loss events in four families, involving 23 of the analyzed endoparasitoid species. Whole-genome alignment suggested that Vg loss occurred during genome rearrangement events. Analysis of Vg receptor gene (VgR) loss, selection, and structural variation in lineages lacking Vg demonstrated functional biases in the patterns of gene loss. The ectoparasitism to endoparasitism transition did not appear to be the primary driver of Vg loss or the subsequent VgR evolution. Together, these findings reveal the genomic changes underlying a unique trait loss in parasitoid wasps. More broadly, this study enhances our understanding of yolk protein loss evolution outside the class Mammalia, highlighting a potential evolutionary trend arising from the availability of an alternative nutrient source for embryonic development.

HEK-omics: The promise of omics to optimize HEK293 for recombinant adeno-associated virus (rAAV) gene therapy manufacturing

Gene therapy is poised to transition from niche to mainstream medicine, with recombinant adeno-associated virus (rAAV) as the vector of choice. However, robust, scalable, industrialized production is required to meet demand and provide affordable patient access, which has not yet materialized. Closing the chasm between demand and supply requires innovation in biomanufacturing to achieve the essential step change in rAAV product yield and quality. Omics provides a rich source of mechanistic knowledge that can be applied to HEK293, the most commonly used cell line for rAAV production. In this review, the findings from a growing number of diverse studies that apply genomics, epigenomics, transcriptomics, proteomics, and metabolomics to HEK293 bioproduction are explored. Learnings from CHO-omics, application of omics approaches to improve CHO bioproduction, provide a framework to explore the potential of "HEK-omics" as a multi-omics-informed approach providing actionable mechanistic insights for improved transient and stable production of rAAV and other recombinant products in HEK293.

The unexpected loss of the 'hunger hormone' ghrelin in true passerines: a game changer in migration physiology

Migratory birds must accumulate large amounts of fat prior to migration to sustain long flights. In passerines, the small body size limits the amount of energy stores that can be transported, and therefore birds undergo cycles of extreme fattening and rapid exhaustion of reserves. Research on these physiological adaptations was rattled by the discovery that birds have lost the main vertebrate regulator of fat deposition, leptin. Recent studies have thus focused on ghrelin, known as 'hunger hormone', a peptide secreted by the gastrointestinal tract to regulate, e.g. food intake and body mass in vertebrates. Studies on domestic species showed that, in birds, ghrelin has effects opposite to those described in mammals such as inhibiting instead of promoting food intake. Furthermore, recent studies have shown that ghrelin administration influences migratory behaviour in passerine birds. Using comparative genomics and immunoaffinity chromatography, we show that ghrelin has been lost in Eupasseres after the basic split from Acanthisitti about 50 Ma. We found that the ghrelin receptor is still conserved in passerines. The maintenance of a functional receptor system suggests that in Eupasserines, another ligand has replaced ghrelin, perhaps to bypass the feedback system that would hinder the large pre-migratory accumulation of subcutaneous fat.

Less, but More: New Insights From Appendicularians on Chordate Fgf Evolution and the Divergence of Tunicate Lifestyles

The impact of gene loss on the diversification of taxa and the emergence of evolutionary innovations remains poorly understood. Here, our investigation on the evolution of the Fibroblast Growth Factors (FGFs) in appendicularian tunicates as a case study reveals a scenario of "less, but more" characterized by massive losses of all Fgf gene subfamilies, except for the Fgf9/16/20 and Fgf11/12/13/14, which in turn underwent two bursts of duplications. Through phylogenetic analysis, synteny conservation, and gene and protein structure, we reconstruct the history of appendicularian Fgf genes, highlighting their paracrine and intracellular functions. An exhaustive analysis of developmental Fgf expression in Oikopleura dioica allows us to identify four associated evolutionary patterns characterizing the "less, but more" conceptual framework: conservation of ancestral functions; function shuffling between paralogs linked to gene losses; innovation of new functions after the duplication bursts; and function extinctions linked to gene losses. Our findings allow us to formulate novel hypotheses about the impact of Fgf losses and duplications on the transition from an ancestral ascidian-like biphasic lifestyle to the fully free-living appendicularians. These hypotheses include massive co-options of Fgfs for the development of the oikoblast and the tail fin; recruitment of Fgf11/12/13/14s into the evolution of a new mouth, and their role modulating neuronal excitability; the evolutionary innovation of an anterior tail FGF signaling source upon the loss of retinoic acid signaling; and the potential link between the loss of Fgf7/10/22 and Fgf8/17/18 and the loss of drastic metamorphosis and tail absorption in appendicularians, in contrast to ascidians.

Modulating vertebrate physiology by genomic fine-tuning of GPCR functions

G protein-coupled receptors (GPCRs) play a crucial role as membrane receptors, facilitating the communication of eukaryotic species with their environment and regulating cellular and organ interactions. Consequently, GPCRs hold immense potential in contributing to adaptation to ecological niches and responding to environmental shifts. Comparative analyses of vertebrate genomes reveal patterns of GPCR gene loss, expansion, and signatures of selection. Integrating these genomic data with insights from functional analyses of gene variants enables the interpretation of genotype-phenotype correlations. This review underscores the involvement of GPCRs in adaptive processes, presenting numerous examples of how alterations in GPCR functionality influence vertebrate physiology or, conversely, how environmental changes impact GPCR functions. The findings demonstrate that modifications in GPCR function contribute to adapting to aquatic, arid, and nocturnal habitats, influencing camouflage strategies, and specializing in particular dietary preferences. Furthermore, the adaptability of GPCR functions provides an effective mechanism in facilitating past, recent, or ongoing adaptations in animal domestication and human evolution and should be considered in therapeutic strategies and drug development.

Sperm RNA code in spermatogenesis and male infertility

Spermatozoa are traditionally thought to be transcriptionally inert, but recent studies have revealed the presence of sperm RNA, some of which is derived from the residues of spermatocyte transcription and some from epididymosomes. Paternal sperm RNA can be affected by external factors and further modified at the post-transcriptional level, for example N6-methyladenosine (m6A), thus shaping spermatogenesis and reproductive outcome. This review briefly introduces the origin of sperm RNA and, on this basis, summarizes the current knowledge on RNA modifications and their functional role in spermatogenesis and male infertility. The bottlenecks and knowledge gaps in the current research on RNA modification in male reproduction have also been indicated. Further investigations are needed to elucidate the functional consequences of these modifications, providing new therapeutic and preventive strategies for reproductive health and genetic inheritance.

Dissection of the synthesis of polyunsaturated fatty acids in nematodes and Collembola of the soil fauna

It is becoming increasingly clear that not only unicellular, photoautotrophic eukaryotes, plants, and fungi, but also invertebrates are capable of synthesizing ω3 long-chain polyunsaturated fatty acids (LC-PUFA) de novo. However, the distribution of this anabolic capacity among different invertebrate groups and its implementation at the gene and protein level are often still unknown. This study investigated the PUFA pathways in common soil fauna, i.e. two nematode and two Collembola species. Of these, one species each (Panagrellus redivivus, Folsomia candida) was assumed to produce ω3 LC-PUFA de novo, while the others (Acrobeloides bodenheimeri, Isotoma caerulea) were supposed to be unable to do so. A highly labeled oleic acid (99 % 13C) was supplemented and the isotopic signal was used to trace its metabolic path. All species followed the main pathway of lipid biosynthesis. However, in A. bodenheimeri this terminated at arachidonic acid (ω6 PUFA), whereas the other three species continued the pathway to eicosapentaenoic acid (ω3 PUFA), including I. caerulea. For the nematode P. redivivus the identification and functional characterization of four new fatty acid desaturase (FAD) genes was performed. These genes encode the FAD activities Δ9, Δ6, and Δ5, respectively. Additionally, the Δ12 desaturase was analyzed, yet the observed activity of an ω3 FAD could not be attributed to a coding gene. In the Collembola F. candida, 11 potential first desaturases (Δ9) and 13 front-end desaturases (Δ6 or Δ5 FADs) have been found. Further sequence analysis indicates the presence of omega FADs, specifically Δ12, which are likely derived from Δ9 FADs.

Progressive transcriptomic shifts in evolved yeast strains following gene knockout

Gene knockout disrupts cellular homeostasis, altering gene expression, and phenotypes. We investigated whether cells return to their pre-knockout transcriptomic state through adaptive evolution experiments on hap4Δ and ade1Δ yeast strains. Analysis revealed that genes with higher expression levels and more physical interaction partners in wild-type strains were more likely to be restored, suggesting that genes of significant functional importance have increased resilience to genetic perturbations. However, as the experiment progressed, most initially restored genes became unrestored. Over 60% of differentially expressed genes in knockout strains remained unrestored in evolved strains. Evolved strains exhibited distinct transcriptomic states, diverging from the original strain over time. Ribosome biogenesis components exhibited systematic sequential changes during the evolution. Our findings suggest the knockout strain transcriptomes struggle to return to the original state even after 28 days of culture. Instead, compensatory mechanisms lead to distinct suboptimal states, highlighting the complex transcriptomic dynamics following genetic perturbations.

Comparative Genomics of the World's Smallest Mammals Reveals Links to Echolocation, Metabolism, and Body Size Plasticity

Originating 30 million years ago, shrews (Soricidae) have diversified into around 400 species worldwide. Shrews display a wide array of adaptations, with some species having developed distinctive traits such as echolocation, underwater diving, and venomous saliva. Accordingly, these tiny insectivores are ideal to study the genomic mechanisms of evolution and adaptation. We conducted a comparative genomic analysis of four shrew species and 16 other mammals to identify genomic variations unique to shrews. Using two existing shrew genomes and two de novo assemblies for the maritime (Sorex maritimensis) and smoky (Sorex fumeus) shrews, we identified mutations in conserved regions of the genomes, also known as accelerated regions, gene families that underwent significant expansion, and positively selected genes. Our analyses unveiled shrew-specific genomic variants in genes associated with the nervous, metabolic, and auditory systems, which can be linked to unique traits in shrews. Notably, genes suggested to be under convergent evolution in echolocating mammals exhibited accelerated regions in shrews, and pathways linked to putative body size plasticity were detected. These findings provide insight into the evolutionary mechanisms shaping shrew species, shedding light on their adaptation and divergence over time.

Genome-Wide Comparisons Reveal Extensive Divergence Within the Lichen Photobiont Genus, Trebouxia

The green algal genus Trebouxia is the most frequently encountered photobiont of the lichen symbiosis. The single-celled symbionts have a worldwide distribution, including all continents and climate zones. The vast, largely undescribed, diversity of Trebouxia lineages is currently grouped into four phylogenetic clades (A, C, I, and S), based on a multilocus phylogeny. Genomes are still scarce, however, and it is unclear how the phylogenetic diversity, the broad ecological tolerances, and the ability to form symbioses with many different fungal host species are reflected in genome-wide differences. Here, we generated PacBio-based de novo genomes of six Trebouxia lineages belonging to the Clades A and S, isolated from lichen individuals of the genus Umbilicaria. Sequences belonging to Clade S have been reported in a previous study, but were reassembled and reanalyzed here. Genome sizes ranged between 63.08 and 73.88 Mb. Repeat content accounted for 9% to 16% of the genome sequences. Based on RNA evidence, we predicted 14,109 to 16,701 gene models per genome, of which 5,203 belonged to a core set of gene families shared by all 6 lineages. Between 121 and 454, gene families are specific to each lineage. About 53% of the genes could be functionally annotated. The presence of biosynthetic gene clusters (6 to 17 per genome) suggests that Trebouxia algae are able to synthesize alkaloids, saccharides, terpenes, NRPSs, and T3PKSs. Phylogenomic comparisons of the six strains indicate prevalent gene gain during Trebouxia evolution. Some of the gene families that exhibited significant evolutionary changes (i.e. gene expansion and contraction) are associated with metabolic processes linked to protein phosphorylation, which is known to have a role in photosynthesis regulation, particularly under changing light conditions. Overall, there is substantial genomic divergence within the algal genus Trebouxia, which may contribute to the genus' large ecological amplitude concerning fungal host diversity and climatic niches.

Distinct Genes with Similar Functions Underlie Convergent Evolution in Myotis Bat Ecomorphs

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes.

Diverse signatures of convergent evolution in cactus-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently approximately 17 times. Using a machine learning-based approach, we further found that cactophily can be predicted with 76% accuracy from both functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which we found to be likely associated with altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved independently through disparate molecular mechanisms. Notably, we found that multiple cactophilic species and their close relatives have been reported as emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-might preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high-throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Unraveling genomic features and phylogenomics through the analysis of three Mexican endemic Myotis genomes

Background

Genomic resource development for non-model organisms is rapidly progressing, seeking to uncover molecular mechanisms and evolutionary adaptations enabling thriving in diverse environments. Limited genomic data for bat species hinder insights into their evolutionary processes, particularly within the diverse Myotis genus of the Vespertilionidae family. In Mexico, 15 Myotis species exist, with three-M. vivesi, M. findleyi, and M. planiceps-being endemic and of conservation concern.

Methods

We obtained samples of Myotis vivesi, M. findleyi, and M. planiceps for genomic analysis. Each of three genomic DNA was extracted, sequenced, and assembled. The scaffolding was carried out utilizing the M. yumanensis genome via a genome-referenced approach within the ntJoin program. GapCloser was employed to fill gaps. Repeat elements were characterized, and gene prediction was done via ab initio and homology methods with MAKER pipeline. Functional annotation involved InterproScan, BLASTp, and KEGG. Non-coding RNAs were annotated with INFERNAL, and tRNAscan-SE. Orthologous genes were clustered using Orthofinder, and a phylogenomic tree was reconstructed using IQ-TREE.

Results

We present genome assemblies of these endemic species using Illumina NovaSeq 6000, each exceeding 2.0 Gb, with over 90% representing single-copy genes according to BUSCO analyses. Transposable elements, including LINEs and SINEs, constitute over 30% of each genome. Helitrons, consistent with Vespertilionids, were identified. Values around 20,000 genes from each of the three assemblies were derived from gene annotation and their correlation with specific functions. Comparative analysis of orthologs among eight Myotis species revealed 20,820 groups, with 4,789 being single copy orthogroups. Non-coding RNA elements were annotated. Phylogenomic tree analysis supported evolutionary chiropterans' relationships. These resources contribute significantly to understanding gene evolution, diversification patterns, and aiding conservation efforts for these endangered bat species.

New Insight into Nucleotide Changes on Irradiated Bactrocera dorsalis (Hendel), A Pest of Horticultural Importance

Bactrocera dorsalis (Hendel) is a major quarantine pest species infesting most of the tropical fruits. Its infestation had significantly reduced and disrupted the export market trade, thus, very crucial to be controlled during the preharvest and postharvest. One of the most sustainable control methods is by using the radiation technique to reduce the pest population, thus curbing the spread of this pest to new geographical areas. The objective of this study was to measure the nucleotide changes in B. dorsalis (larval, pupal and adult stages) which had been irradiated with 50 to 400 Gray, using Gamma Cell Biobeam GM8000 irradiator with Cesium-137 source at the Malaysian Nuclear Agency, Selangor, Malaysia. Data from the treated samples (with and without morphological changes) were analysed using cytochrome oxidase subunit I (COI). The alignment of 59 sequences resulted in 0.92% variables with only four characters that were parsimony informative, and six sites (30, 60, 234, 282, 483 and 589) which had nucleotide changes, but had not been translated to another protein. Low polymorphism was presented on the sample groups, with only four haplotypes, but with high diversity value (Hd) = 0.5885. The phylogeny trees formed soft polytomy in both trees [neighbour joining (NJ) and maximum parsimony (MP)] presenting a mixture of individuals but did not show any significant difference between treatments. This finding concluded that low mutation had occurred on the treated B. dorsalis and this information is very valuable in getting new insight on the survival of B. dorsalis in the horticulture industry.

Identification and characterization of STAT family in silver pomfret (Pampus argenteus) involved in different exogenous stresses

Members of the Signal Transducer and Activator of Transcription (STAT) family function pivotally as transcriptional activators integral to the modulation of inflammatory responses. The aquaculture of silver pomfret is frequently compromised by the imposition of exogenous stressors, which include thermal fluctuations, notably low-temperatures, diminished oxygen levels, and the onslaught of bacterial pathogens. Notwithstanding the critical impact of these stressors, the scientific literature presents a notable gap in our understanding of the STAT pathway's role in the silver pomfret's adaptive response mechanisms. To address this lacuna, we identified stat genes in the silver pomfret-denominated as Pastat1, Pastat2, Pastat3, Pastat4, and Pastat5-through a thorough and systematic bioinformatics analysis. Further scrutiny of the gene configurations and constituent motifs has elucidated that STAT proteins possess analogous structural frameworks and exhibit significant evolutionary preservation. Subsequently, the expression patterns of five stat genes were verified by RT-qPCR in twelve different tissues and four growth periods in healthy fish, showing that the expression of Pastat genes was temporally and spatially specific, with most of the stat genes expressed at higher levels in the spleen, following muscle, gill, and liver. Transcriptomic analysis of exposure to exogenous stressors, specifically formaldehyde and low-temperature conditions, elucidated that Pastat1 and Pastat2 genes exhibited a heightened sensitivity to these environmental challenges. RT-qPCR assays demonstrated a marked alteration in the expression profiles of jak1 and Pastat gene suites in PaS upon prolonged bacterial infection subsequent to these exogenous insults. Moreover, the gene expression of the downstream effectors involved in innate immunity and apoptosis displayed marked deviations. This study additionally elucidated the Pastat gene family's role in modulating the innate immune response and apoptotic regulation within the silver pomfret during exogenous stressors and subsequent pathogenic incursions.

The Evolution of Transglutaminases Underlies the Origin and Loss of Cornified Skin Appendages in Vertebrates

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes.

Adaptive Evolution of the Greater Horseshoe Bat AANAT: Insights into the Link between AANAT and Hibernation Rhythms

Arylalkylamine N-acetyltransferase (AANAT) is a crucial rate-limiting enzyme in the synthesis of melatonin. AANAT has been confirmed to be independently duplicated and inactivated in different animal taxa in order to adapt to the environment. However, the evolutionary forces associated with having a single copy of AANAT remain unclear. The greater horseshoe bat has a single copy of AANAT but exhibits different hibernation rhythms in various populations. We analyzed the adaptive evolution at the gene and protein levels of AANAT from three distinct genetic lineages in China: northeast (NE), central east (CE), and southwest (SW). The results revealed greater genetic diversity in the AANAT loci of the NE and CE lineage populations that have longer hibernation times, and there were two positive selection loci. The catalytic capacity of AANAT in the Liaoning population that underwent positive selection was significantly higher than that of the Yunnan population (p < 0.05). This difference may be related to the lower proportion of α helix and the variation in two interface residues. The adaptive evolution of AANAT was significantly correlated with climate and environment (p < 0.05). After controlling for geographical factors (latitude and altitude), the evolution of AANAT by the negative temperature factor was represented by the monthly mean temperature (r = -0.6, p < 0.05). The results identified the gene level variation, functional adaptation, and evolutionary driving factors of AANAT, provide an important foundation for further understanding the adaptive evolution of the single copy of AANAT in pteropods, and may offer evidence for adaptive hibernation rhythms in bats.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

Pervasive loss of regulated necrotic cell death genes in elephants, hyraxes, and sea cows (Paenungualta)

Gene loss can promote phenotypic differences between species, for example, if a gene constrains phenotypic variation in a trait, its loss allows for the evolution of a greater range of variation or even new phenotypes. Here, we explore the contribution of gene loss to the evolution of large bodies and augmented cancer resistance in elephants. We used genomes from 17 Afrotherian and Xenarthran species to identify lost genes, i.e., genes that have pseudogenized or been completely lost, and Dollo parsimony to reconstruct the evolutionary history of gene loss across species. We unexpectedly discovered a burst of gene losses in the Afrotherian stem lineage and found that the loss of genes with functions in regulated necrotic cell death modes was pervasive in elephants, hyraxes, and sea cows (Paenungulata). Among the lost genes are MLKL and RIPK3, which mediate necroptosis, and sensors that activate inflammasomes to induce pyroptosis, including AIM2, MEFV, NLRC4, NLRP1, and NLRP6. These data suggest that the mechanisms that regulate necrosis and pyroptosis are either extremely derived or potentially lost in these lineages, which may contribute to the repeated evolution of large bodies and cancer resistance in Paenungulates as well as susceptibility to pathogen infection.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Functional analysis of CYP71AV1 reveals the evolutionary landscape of artemisinin biosynthesis

Artemisinin biosynthesis, unique to Artemisia annua, is suggested to have evolved from the ancestral costunolide biosynthetic pathway commonly found in the Asteraceae family. However, the evolutionary landscape of this process is not fully understood. The first oxidase in artemisinin biosynthesis, CYP71AV1, also known as amorpha-4,11-diene oxidase (AMO), has specialized from ancestral germacrene A oxidases (GAOs). Unlike GAO, which exhibits catalytic promiscuity toward amorpha-4,11-diene, the natural substrate of AMO, AMO has lost its ancestral activity on germacrene A. Previous studies have suggested that the loss of the GAO copy in A. annua is responsible for the abolishment of the costunolide pathway. In the genome of A. annua, there are two copies of AMO, each of which has been reported to be responsible for the different product profiles of high- and low-artemisinin production chemotypes. Through analysis of their tissue-specific expression and comparison of their sequences with those of other GAOs, it was discovered that one copy of AMO (AMOHAP) exhibits a different transcript compared to the reported artemisinin biosynthetic genes and shows more sequence similarity to other GAOs in the catalytic regions. Furthermore, in a subsequent in vitro enzymatic assay, the recombinant protein of AMOHAP unequivocally demonstrated GAO activity. This result clearly indicates that AMOHAP is a GAO rather than an AMO and that its promiscuous activity on amorpha-4,11-diene has led to its misidentification as an AMO in previous studies. In addition, the divergent expression pattern of AMOHAP compared to that of the upstream germacrene A synthase may have contributed to the abolishment of costunolide biosynthesis in A. annua. Our findings reveal a complex evolutionary landscape in which the emergence of a new metabolic pathway replaces an ancestral one.

Adaptive strategies of aquatic mammals: Exploring the role of the HIF pathway and hypoxia tolerance

Aquatic mammals (marine and freshwater species) share significant and similar adaptations, enabling them to tolerate hypoxia during regular breath-hold diving. Despite the established importance of HIF1A, a master regulator in the molecular mechanism of hypoxia response, and other associated genes, their role in the evolutionary adaptation of aquatic mammals is not fully understood. In this study, we investigated this topic by employing a candidate gene approach to analyze 11 critical genes involved in the HIF1A signaling pathway in aquatic mammals. Our gene analyses included evaluating positive and negative selection, relaxation or constriction of selection, and molecular convergence compared to other terrestrial mammals, including subterranean mammals. Evidence of selection suggested a significant role of negative selection, as well as relaxation of the selective regime in cetaceans for most of these genes. We found that the glutamine 68 variant in the HIF3α protein is unique to cetaceans and initial evaluations indicated a destabilizing effect on protein structure. However, further analyses are necessary to evaluate its functional impact and adaptive relevance in this taxon.

Integrative single-cell characterization of a frugivorous and an insectivorous bat kidney and pancreas

Frugivory evolved multiple times in mammals, including bats. However, the cellular and molecular components driving it remain largely unknown. Here, we use integrative single-cell sequencing (scRNA-seq and scATAC-seq) on insectivorous (Eptesicus fuscus; big brown bat) and frugivorous (Artibeus jamaicensis; Jamaican fruit bat) bat kidneys and pancreases and identify key cell population, gene expression and regulatory differences associated with the Jamaican fruit bat that also relate to human disease, particularly diabetes. We find a decrease in loop of Henle and an increase in collecting duct cells, and differentially active genes and regulatory elements involved in fluid and electrolyte balance in the Jamaican fruit bat kidney. The Jamaican fruit bat pancreas shows an increase in endocrine and a decrease in exocrine cells, and differences in genes and regulatory elements involved in insulin regulation. We also find that these frugivorous bats share several molecular characteristics with human diabetes. Combined, our work provides insights from a frugivorous mammal that could be leveraged for therapeutic purposes.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Genome evolution in plants and the origins of innovation

Plant evolution has been characterised by a series of major novelties in their vegetative and reproductive traits that have led to greater complexity. Underpinning this diversification has been the evolution of the genome. When viewed at the scale of the plant kingdom, plant genome evolution has been punctuated by conspicuous instances of gene and whole-genome duplication, horizontal gene transfer and extensive gene loss. The periods of dynamic genome evolution often coincide with the evolution of key traits, demonstrating the coevolution of plant genomes and phenotypes at a macroevolutionary scale. Conventionally, plant complexity and diversity have been considered through the lens of gene duplication and the role of gene loss in plant evolution remains comparatively unexplored. However, in light of reductive evolution across multiple plant lineages, the association between gene loss and plant phenotypic diversity warrants greater attention.

Genome Assembly of Cordia subcordata, a Coastal Protection Species in Tropical Coral Islands

Cordia subcordata trees or shrubs, belonging to the Boraginaceae family, have strong resistance and have adapted to their habitat on a tropical coral island in China, but the lack of genome information regarding its genetic background is unclear. In this study, the genome was assembled using both short/long whole genome sequencing reads and Hi-C reads. The assembled genome was 475.3 Mb, with 468.7 Mb (99.22%) of the sequences assembled into 16 chromosomes. Repeat sequences accounted for 54.41% of the assembled genome. A total of 26,615 genes were predicted, and 25,730 genes were functionally annotated using different annotation databases. Based on its genome and the other 17 species, phylogenetic analysis using 336 single-copy genes obtained from ortholog analysis showed that C. subcordata was a sister to Coffea eugenioides, and the divergence time was estimated to be 77 MYA between the two species. Gene family evolution analysis indicated that the significantly expanded gene families were functionally related to chemical defenses against diseases. These results can provide a reference to a deeper understanding of the genetic background of C. subcordata and can be helpful in exploring its adaptation mechanism on tropical coral islands in the future.

Three Blind Moles: Molecular Evolutionary Insights on the Tempo and Mode of Convergent Eye Degeneration in Notoryctes typhlops (Southern Marsupial Mole) and Two Chrysochlorids (Golden Moles)

Golden moles (Chrysochloridae) and marsupial moles (Notoryctidae) are textbook examples of convergent evolution. Both taxa are highly adapted to subterranean lifestyles and have powerful limbs for digging through the soil/sand, ears that are adapted for low-frequency hearing, vestigial eyes that are covered by skin and fur, and the absence of optic nerve connections between the eyes and the brain. The eyes of marsupial moles also lack a lens as well as retinal rods and cones. Two hypotheses have been proposed to account for the greater degeneracy of the eyes of marsupial moles than golden moles. First, marsupial moles may have had more time to adapt to their underground habitat than other moles. Second, the eyes of marsupial moles may have been rapidly and recently vestigialized to (1) reduce the injurious effects of sand getting into the eyes and (2) accommodate the enlargement of lacrimal glands that keep the nasal cavity moist and prevent the entry of sand into the nasal passages during burrowing. Here, we employ molecular evolutionary methods on DNA sequences for 38 eye genes, most of which are eye-specific, to investigate the timing of relaxed selection (=neutral evolution) for different groups of eye-specific genes that serve as proxies for distinct functional components of the eye (rod phototransduction, cone phototransduction, lens/cornea). Our taxon sampling included 12 afrothere species, of which two are golden moles (Amblysomus hottentotus, Chrysochloris asiatica), and 28 marsupial species including two individuals of the southern marsupial mole (Notoryctes typhlops). Most of the sequences were mined from databases, but we also provide new genome data for A. hottentotus and one of the two N. typhlops individuals. Even though the eyes of golden moles are less degenerate than the eyes of marsupial moles, there are more inactivating mutations (e.g., frameshift indels, premature stop codons) in their cone phototransduction and lens/cornea genes than in orthologous genes of the marsupial mole. We estimate that cone phototransduction recovery genes were inactivated first in each group, followed by lens/cornea genes and then cone phototransduction activation genes. All three groups of genes were inactivated earlier in golden moles than in marsupial moles. For the latter, we estimate that lens/cornea genes were inactivated ~17.8 million years ago (MYA) when stem notoryctids were burrowing in the soft soils of Australian rainforests. Selection on phototransduction activation genes was relaxed much later (5.38 MYA), during the early stages of Australia's aridification that produced coastal sand plains and eventually sand dunes. Unlike cone phototransduction activation genes, rod phototransduction activation genes are intact in both golden moles and one of the two individuals of N. typhlops. A second marsupial mole individual has just a single inactivating mutation in one of the rod phototransduction activation genes (PDE6B). One explanation for this result is that some rod phototransduction activation genes are pleiotropic and are expressed in extraocular tissues, possibly in conjunction with sperm thermotaxis.

Identification and characterization of toll-like receptor genes in silver pomfret (Pampus argenteus) and their involvement in the host immune response to Photobacterium damselae subsp. Damselae and Nocardia seriolae infection

Toll-like receptors (TLRs) are vital pattern recognition receptors that play a critical role in the innate immune response against pathogenic attack. Among the bacteria commonly found in the culture process of silver pomfret, Photobacterium damselae subsp. Damselae (PDD, gram-negative) and Nocardia seriolae (NS, gram-positive), can cause large-scale mortality in this fish species. However, there is currently no research on the role of TLRs in mediating the immune response of silver pomfret to these two bacterial infections. Therefore, in this study, we identified nine PaTLRs family members, including several fish-specific TLRs (TLR14 and TLR21). Phylogenetic analysis revealed that these PaTLRs genes could be classified into five subfamilies, namely TLR1, TLR3, TLR5, TLR7, and TLR11, indicating their evolutionary conservation. To further explore the interactions of TLR genes with immune-related mediators, protein and protein interaction network (PPI) results were generated to explain the association of TLR genes with TNF receptor-associated factor 6 (TRAF6) and other relevant genes in the MyD88-dependent pathway and NF-κb signaling pathway. Subsequently, RT-qPCR was conducted to verify the expression patterns of the nine TLR genes in the gills, skin, kidney, liver, and spleen of healthy fish, with most of the TLRs showing high expression levels in the spleen. Following infection with PDD and NS, these PaTLRs exhibited different expression patterns in the spleen, with PaTLR2, PaTLR3, PaTLR5, PaTLR7, PaTLR9, and PaTLR14 being significantly up-regulated. Furthermore, when spleen cells were treated with bacterial compositions, the majority of PaTLRs expression was up-regulated in response to Lipopolysaccharide (LPS) and lipophosphorylcholic acid (LTA) treatment, except for PaTLR21. Finally, changes in the expression levels of TLR-interacting genes were also observed under the stimulation of bacteria and bacterial compositions. The results of this study provide a preliminary reference for further understanding the mechanism of the innate immune response of the TLR gene family in silver pomfret and offer theoretical support for addressing the disease problems encountered during large-scale fish breeding.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Pseudogenization of the Hair-Related Genes PADI3 and S100A3 in Cetaceans and Hippopotamus amphibius

Hair-related genes in mammals play important roles in the development and maintenance of hair and other keratinous structures in mammals. The peptidyl arginine deiminase 3 (PADI3) gene encodes an enzyme that catalyzes the conversion of arginine residues to citrulline. The S100 calcium binding protein A3 (S100A3) gene encodes a protein that is highly expressed in the hair cuticle and contains arginine residues that are converted to citrullines by PADI enzymes. In this study, we investigated the pseudogenization events of PADI3 and S100A3 in cetaceans and Hippopotamus amphibius. We found that PADI3 underwent three independent pseudogenization events during cetacean evolution, in baleen whales, toothed cetaceans other than Physeter catodon, and P. catodon. Notably, the entire PADI3 gene is absent in the baleen whales. Pseudogenization of S100A3 occurred independently in cetaceans and H. amphibius. Interestingly, we found that in cetaceans S100A3 underwent pseudogenization before PADI3, suggesting that differential selection pressures were acting on the two genes. Our findings provide valuable insights into the molecular evolution of these genes in cetaceans and hippopotamuses, highlighting their importance for understanding the evolution of hair-related genes.

Diverse signatures of convergent evolution in cacti-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently ~17 times. Using machine-learning, we further found that cactophily can be predicted with 76% accuracy from functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which is likely associated with duplication and altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved through disparate molecular mechanisms. Remarkably, multiple cactophilic lineages and their close relatives are emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-may preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Wide-scale identification of novel/eliminated genes responsible for evolutionary transformations

BACKGROUND

It is generally accepted that most evolutionary transformations at the phenotype level are associated either with rearrangements of genomic regulatory elements, which control the activity of gene networks, or with changes in the amino acid contents of proteins. Recently, evidence has accumulated that significant evolutionary transformations could also be associated with the loss/emergence of whole genes. The targeted identification of such genes is a challenging problem for both bioinformatics and evo-devo research.

RESULTS

To solve this problem we propose the WINEGRET method, named after the first letters of the title. Its main idea is to search for genes that satisfy two requirements: first, the desired genes were lost/emerged at the same evolutionary stage at which the phenotypic trait of interest was lost/emerged, and second, the expression of these genes changes significantly during the development of the trait of interest in the model organism. To verify the first requirement, we do not use existing databases of orthologs, but rely purely on gene homology and local synteny by using some novel quickly computable conditions. Genes satisfying the second requirement are found by deep RNA sequencing. As a proof of principle, we used our method to find genes absent in extant amniotes (reptiles, birds, mammals) but present in anamniotes (fish and amphibians), in which these genes are involved in the regeneration of large body appendages. As a result, 57 genes were identified. For three of them, c-c motif chemokine 4, eotaxin-like, and a previously unknown gene called here sod4, essential roles for tail regeneration were demonstrated. Noteworthy, we established that the latter gene belongs to a novel family of Cu/Zn-superoxide dismutases lost by amniotes, SOD4.

CONCLUSIONS

We present a method for targeted identification of genes whose loss/emergence in evolution could be associated with the loss/emergence of a phenotypic trait of interest. In a proof-of-principle study, we identified genes absent in amniotes that participate in body appendage regeneration in anamniotes. Our method provides a wide range of opportunities for studying the relationship between the loss/emergence of phenotypic traits and the loss/emergence of specific genes in evolution.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

Genome-wide analysis of cuticle protein family genes in rice stem borer Chilo suppressalis: Insights into their role in environmental adaptation and insecticidal stress response

Insect cuticle plays a key role in insect survival, adaptation and prosperity by serving as the exoskeleton and the first barrier against environmental stresses. As the major components of insect cuticle, the diverse structural cuticle proteins (CPs) contribute to variation in physical properties and functions of cuticle. However, the roles of CPs in cuticular versatility, especially in the stress response or adaption, remain incompletely understood. In this study, we performed a genome-wide analysis of CP superfamily in the rice-boring pest Chilo suppressalis. A total of 211 CP genes were identified and their encoding proteins were classified into eleven families and three subfamilies (RR1, RR2, and RR3). The comparative genomic analysis of CPs revealed that C. suppressalis had fewer CP genes compared to other lepidopteran species, which largely resulted from a less expansion of his-rich RR2 genes involved in cuticular sclerotization, suggesting long-term boring life of C. suppressalis inside rice hosts might evolutionarily prefer cuticular elasticity rather than cuticular sclerotization. We also investigated the response pattern of all CP genes under insecticidal stresses. >50 % CsCPs were upregulated at least 2-fold under insecticidal stresses. Notably, the majority of the highly upregulated CsCPs formed gene pairs or gene clusters on chromosomes, indicating the rapid response of adjacent CsCPs to insecticidal stress. Most high-response CsCPs encoded AAPA/V/L motifs that are related to cuticular elasticity and >50 % of the sclerotization-related his-rich RR2 genes were also upregulated. These results suggested the potential roles of CsCPs in balancing the elasticity and sclerotization of cuticles, which is essential for the survival and adaptation of plant borers including C. suppressalis. Our study provides valuable information for further developing cuticle-based strategies of both pest management and biomimetic applications.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Adaptive Selection in the Evolution of Aquaglyceroporins in Mammals

Aquaporins (AQPs) are integral membrane proteins responsible for water transport across cellular membranes in both prokaryotes and eukaryotes. A subfamily of AQPs, known as aquaglyceroporins (AQGPs), facilitate the transport of small solutes such as glycerol, water, and other solutes across cellular membranes. These proteins are involved in a variety of physiological processes, such as organogenesis, wound healing, and hydration. Although AQPs have been studied extensively in different species, their conservation patterns, phylogenetic relationships, and evolution in mammals remain unexplored. In the present study, 119 AQGP coding sequences from 31 mammalian species were analysed to identify conserved residues, gene organisation, and most importantly, the nature of AQGP gene selection. Repertoire analysis revealed the absence of AQP7, 9, and 10 genes in certain species of Primates, Rodentia, and Diprotodontia, although not all three genes were absent in a single species. Two Asparagine-Proline-Alanine (NPA) motifs located at the N- and C-terminal ends, aspartic acid (D) residues, and the ar/R region were conserved in AQP3, 9, and 10. Six exons encoding the functional MIP domain of AQGP genes were found to be conserved across mammalian species. Evolutionary analysis indicated signatures of positive selection in AQP7, 9, and 10 amongst different mammalian lineages. Furthermore, substitutions of certain amino acids located close to critical residues may alter AQGP functionality, which is crucial for substrate selectivity, pore formation, and transport efficiency required for the maintenance of homeostasis in different mammalian species.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Integrative single-cell characterization of frugivory adaptations in the bat kidney and pancreas

Frugivory evolved multiple times in mammals, including bats. However, the cellular and molecular components driving it remain largely unknown. Here, we used integrative single-cell sequencing on insectivorous and frugivorous bat kidneys and pancreases and identified key cell population, gene expression and regulatory element differences associated with frugivorous adaptation that also relate to human disease, particularly diabetes. We found an increase in collecting duct cells and differentially active genes and regulatory elements involved in fluid and electrolyte balance in the frugivore kidney. In the frugivorous pancreas, we observed an increase in endocrine and a decrease in exocrine cells and differences in genes and regulatory elements involved in insulin regulation. Combined, our work provides novel insights into frugivorous adaptation that also could be leveraged for therapeutic purposes.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Loss of a gluconeogenic muscle enzyme contributed to adaptive metabolic traits in hummingbirds

Hummingbirds possess distinct metabolic adaptations to fuel their energy-demanding hovering flight, but the underlying genomic changes are largely unknown. Here, we generated a chromosome-level genome assembly of the long-tailed hermit and screened for genes that have been specifically inactivated in the ancestral hummingbird lineage. We discovered that FBP2 (fructose-bisphosphatase 2), which encodes a gluconeogenic muscle enzyme, was lost during a time period when hovering flight evolved. We show that FBP2 knockdown in an avian muscle cell line up-regulates glycolysis and enhances mitochondrial respiration, coincident with an increased mitochondria number. Furthermore, genes involved in mitochondrial respiration and organization have up-regulated expression in hummingbird flight muscle. Together, these results suggest that FBP2 loss was likely a key step in the evolution of metabolic muscle adaptations required for true hovering flight.

Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

Molecular evolution and signatures of selective pressures on Bos, focusing on the Nelore breed (Bos indicus)

Evolutionary history leads to genome changes over time, especially for species that have experienced intense selective pressures over a short period. Here, we investigated the genomic evolution of Bos species by searching for potential selection signatures, focusing on Nelore, an economically relevant cattle breed in Brazil. We assessed the genomic processes determining the molecular evolution across Nelore and thirteen other related taxa by evaluating (i) amino acid sequence conservation, (ii) the dN/dS ratio, and (iii) gene families' turnover rate (λ). Low conserved regions potentially associated with fatty acid metabolism seem to reflect differences in meat fat content in taxa with different evolutionary histories. All Bos species presented genes under positive selection, especially B. indicus and Nelore, which include transport protein cobalamin, glycolipid metabolism, and hormone signaling. These findings could be explained by constant selective pressures to obtain higher immune resistance and efficient metabolism. The gene contraction rate across the Nelore + B. indicus branch was almost nine times higher than that in other lineages (λ = 0.01043 vs. 0.00121), indicating gene losses during the domestication process. Amino acid biosynthesis, reproductive and innate immune system-related pathways were associated with genes recognized within the most frequent rapidly evolving gene families and in genes under positive selection, supporting the substantial relevance of such traits from a domestication perspective. Our data provide new insights into how the genome may respond to intense artificial selection in distinct taxa, and reinforces the presence of selective pressures on traits potentially relevant for future animal breeding investments.

Endocardium in Hypoplastic Left Heart Syndrome: Implications from In Vitro Study

Endocardium lines the inner layer of the heart ventricle and serves as the source of valve endothelial cells and interstitial cells. Previously, endocardium-associated abnormalities in hypoplastic left heart syndrome (HLHS) have been reported, including endocardial fibroelastosis (EFE) and mitral and aortic valve malformation. However, few mechanistic studies have investigated the molecular pathological changes in endocardial cells. Recently, the emergence of a powerful in vitro system-induced pluripotent stem cells (iPSCs)-was applied to study various genetic diseases, including HLHS. This review summarized current in vitro studies in understanding the endocardial pathology in HLHS, emphasizing new findings of the cellular phenotypes and underlying molecular mechanisms. Lastly, a future perspective is provided regarding the better recapitulation of endocardial phenotypes in a dish.

Birth-and-death evolution of ribonuclease 9 genes in Cetartiodactyla

RNase9 plays a reproductive function and has been recognized as an important member of the ribonuclease (RNase) A superfamily, a gene family that is widely used as a model for molecular evolutionary studies. Here, we identified 178 RNase9 genes from 95 Cetartiodactyla species that represent all four lineages and 21 families of this clade. Unexpectedly, RNase9 experienced an evolutionary scenario of "birth and death" in Ruminantia, and expression analyses showed that duplicated RNase9A and RNase9B genes are expressed in reproductive tissues (epididymis, vas deferens or prostate). This expression pattern combined with the estimate that these genes duplicated during the middle Eocene, a time when Ruminantia become a successful lineage, suggests that the RNase9 gene duplication might have been advantageous for promoting sperm motility and male fertility as an adaptation to climate seasonality changes of this period. In contrast, all RNase9 genes were lost in the Cetacean lineage, which might be associated with their high levels of prostatic lesions and lower reproductive rates as adaptations to a fully aquatic environment and a balance to the demands of ocean resources. This study reveals a complex and intriguing evolutionary history and functional divergence for RNase9 in Cetartiodactyla, providing new insights into the evolution of the RNaseA superfamily and molecular mechanisms for organismal adaptations to the environment.