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

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.

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.

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.

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.

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.

Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants

The origin of plants and their colonization of land fundamentally transformed the terrestrial environment. Here we elucidate the basis of this formative episode in Earth history through patterns of lineage, gene and genome evolution. We use new fossil calibrations, a relative clade age calibration (informed by horizontal gene transfer) and new phylogenomic methods for mapping gene family origins. Distinct rooting strategies resolve tracheophytes (vascular plants) and bryophytes (non-vascular plants) as monophyletic sister groups that diverged during the Cambrian, 515-494 million years ago. The embryophyte stem is characterized by a burst of gene innovation, while bryophytes subsequently experienced an equally dramatic episode of reductive genome evolution in which they lost genes associated with the elaboration of vasculature and the stomatal complex. Overall, our analyses reveal that extant tracheophytes and bryophytes are both highly derived from a more complex ancestral land plant. Understanding the origin of land plants requires tracing character evolution across a diversity of modern lineages.

Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis

Land plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is a continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning spermatogenesis remain to be elucidated. We successfully identified candidate genes involved in spermatogenesis, deeply divergent BLD10s, by computational analyses combining multiple methods and omics data. We then examined the functions of BLD10s in the liverwort Marchantia polymorpha and the moss Physcomitrium patens. MpBLD10 and PpBLD10 are required for normal basal body (BB) and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, and thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. We identified orthologs of MpBLD10 and PpBLD10 in diverse Streptophyta and found that MpBLD10 and PpBLD10 are orthologous to BLD10/CEP135 family proteins, which function in BB assembly. However, BLD10s evolved especially quickly in land plants and MpBLD10 might have acquired additional functions in spermatozoid formation through rapid molecular evolution.

Glycan-Adhering Lectins and Experimental Evaluation of a Lectin FimH Inhibitor in Enterohemorrhagic Escherichia coli (EHEC) O157:H7 Strain EDL933

In this study, we tried to develop a FimH inhibitor that inhibits adhesion of enterohemorrhagic Escherichia coli (EHEC) on the epithelium of human intestine during the initial stage of infections. Using a T7 phage display method with a reference strain, EHEC EDL933, FimH was selected as an adherent lectin to GM1a and Gb3 glycans. In order to detect the ligand binding domain (LBD) of FimH, we used a docking simulation and found three binding site sequences of FimH, i.e., P1, P2, and P3. Among Gb3 mimic peptides, P2 was found to have the strongest binding strength. Moreover, in vitro treatment with peptide P2 inhibited binding activity in a concentration-dependent manner. Furthermore, we conducted confirmation experiments through several strains isolated from patients in Korea, EHEC NCCP15736, NCCP15737, and NCCP15739. In addition, we analyzed the evolutionary characteristics of the predicted FimH lectin-like adhesins to construct a lectin-glycan interaction (LGI). We selected 70 recently differentiated strains from the phylogenetic tree of 2240 strains with Shiga toxin in their genome. We can infer EHEC strains dynamically evolved but FimH was conserved during the evolution time according to the phylogenetic tree. Furthermore, FimH could be a reliable candidate of drug target in terms of evolution. We examined how pathogen lectins interact with host glycans early in infection in EDL933 as well as several field strains and confirmed that glycan-like peptides worked as an initial infection inhibitor.

Molecular evolution of vision-related genes may contribute to marsupial photic niche adaptations

Vision plays an essential role in the life of many animals. While most mammals are night-active (nocturnal), many have adapted to novel light environments. This includes diurnal (day-active) and crepuscular (twilight-active) species. Here, we used integrative approaches to investigate the molecular evolution of 112 vision-related genes across 19 genomes representing most marsupial orders. We found that four genes (GUCA1B, GUCY2F, RGR, and SWS2) involved in retinal phototransduction likely became functionally redundant in the ancestor of marsupials, a group of largely obligate nocturnal mammals. We also show evidence of rapid evolution and positive selection of bright-light vision genes in the common ancestor of Macropus (kangaroos, wallaroos, and wallabies). Macropus-specific amino acid substitutions in opsin genes (LWS and SWS1), in particular, may be an adaptation for crepuscular vision in this genus via opsin spectral sensitivity tuning. Our study set the stage for functional genetics studies and provides a stepping stone to future research efforts that fully capture the visual repertoire of marsupials.

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.

Vision-related convergent gene losses reveal SERPINE3's unknown role in the eye

Despite decades of research, knowledge about the genes that are important for development and function of the mammalian eye and are involved in human eye disorders remains incomplete. During mammalian evolution, mammals that naturally exhibit poor vision or regressive eye phenotypes have independently lost many eye-related genes. This provides an opportunity to predict novel eye-related genes based on specific evolutionary gene loss signatures. Building on these observations, we performed a genome-wide screen across 49 mammals for functionally uncharacterized genes that are preferentially lost in species exhibiting lower visual acuity values. The screen uncovered several genes, including SERPINE3, a putative serine proteinase inhibitor. A detailed investigation of 381 additional mammals revealed that SERPINE3 is independently lost in 18 lineages that typically do not primarily rely on vision, predicting a vision-related function for this gene. To test this, we show that SERPINE3 has the highest expression in eyes of zebrafish and mouse. In the zebrafish retina, serpine3 is expressed in Müller glia cells, a cell type essential for survival and maintenance of the retina. A CRISPR-mediated knockout of serpine3 in zebrafish resulted in alterations in eye shape and defects in retinal layering. Furthermore, two human polymorphisms that are in linkage with SERPINE3 are associated with eye-related traits. Together, these results suggest that SERPINE3 has a role in vertebrate eyes. More generally, by integrating comparative genomics with experiments in model organisms, we show that screens for specific phenotype-associated gene signatures can predict functions of uncharacterized genes.

Gene loss and compensatory evolution promotes the emergence of morphological novelties in budding yeast

Deleterious mutations are generally considered to be irrelevant for morphological evolution. However, they could be compensated by conditionally beneficial mutations, thereby providing access to new adaptive paths. Here we use high-dimensional phenotyping of laboratory-evolved budding yeast lineages to demonstrate that new cellular morphologies emerge exceptionally rapidly as a by-product of gene loss and subsequent compensatory evolution. Unexpectedly, the capacities for invasive growth, multicellular aggregation and biofilm formation also spontaneously evolve in response to gene loss. These multicellular phenotypes can be achieved by diverse mutational routes and without reactivating the canonical regulatory pathways. These ecologically and clinically relevant traits originate as pleiotropic side effects of compensatory evolution and have no obvious utility in the laboratory environment. The extent of morphological diversity in the evolved lineages is comparable to that of natural yeast isolates with diverse genetic backgrounds and lifestyles. Finally, we show that both the initial gene loss and subsequent compensatory mutations contribute to new morphologies, with their synergistic effects underlying specific morphological changes. We conclude that compensatory evolution is a previously unrecognized source of morphological diversity and phenotypic novelties.

Molecular evolutionary analyses of tooth genes support sequential loss of enamel and teeth in baleen whales (Mysticeti)

The loss of teeth and evolution of baleen racks in Mysticeti was a profound transformation that permitted baleen whales to radiate and diversify into a previously underutilized ecological niche of bulk filter-feeding on zooplankton and other small prey. Ancestral state reconstructions suggest that postnatal teeth were lost in the common ancestor of crown Mysticeti. Genomic studies provide some support for this hypothesis and suggest that the genetic toolkit for enamel production was inactivated in the common ancestor of living baleen whales. However, molecular studies to date have not provided direct evidence for the complete loss of teeth, including their dentin component, on the stem mysticete branch. Given these results, several questions remain unanswered: (1) Were teeth lost in a single step or did enamel loss precede dentin loss? (2) Was enamel lost early or late on the stem mysticete branch? (3) If enamel and dentin/tooth loss were decoupled in the ancestry of baleen whales, did dentin loss occur on the stem mysticete branch or independently in different crown mysticete lineages? To address these outstanding questions, we compiled and analyzed complete protein-coding sequences for nine tooth-related genes from cetaceans with available genome data. Seven of these genes are associated with enamel formation (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) whereas two other genes are either dentin-specific (DSPP) or tooth-specific (ODAPH) but not enamel-specific. Molecular evolutionary analyses indicate that all seven enamel-specific genes have inactivating mutations that are scattered across branches of the mysticete tree. Three of the enamel genes (ACP4, KLK4, MMP20) have inactivating mutations that are shared by all mysticetes. The two genes that are dentin-specific (DSPP) or tooth-specific (ODAPH) do not have any inactivating mutations that are shared by all mysticetes, but there are shared mutations in Balaenidae as well as in Plicogulae (Neobalaenidae + Balaenopteroidea). These shared mutations suggest that teeth were lost at most two times. Shared inactivating mutations and dN/dS analyses, in combination with cetacean divergence times, were used to estimate inactivation times of genes and by proxy enamel and tooth phenotypes at ancestral nodes. The results of these analyses are most compatible with a two-step model for the loss of teeth in the ancestry of living baleen whales: enamel was lost very early on the stem Mysticeti branch followed by the independent loss of dentin (and teeth) in the common ancestors of Balaenidae and Plicogulae, respectively. These results imply that some stem mysticetes, and even early crown mysticetes, may have had vestigial teeth comprised of dentin with no enamel. Our results also demonstrate that all odontocete species (in our study) with absent or degenerative enamel have inactivating mutations in one or more of their enamel genes.

Deciphering the RNA universe in sperm in its role as a vertical information carrier

The inheritance of neurophysiologic and neuropsychologic complex diseases can only partly be explained by the Mendelian concept of genetic inheritance. Previous research showed that both psychological disorders like post-traumatic stress disorder and metabolic diseases are more prevalent in the progeny of affected parents. This could suggest an epigenetic mode of transmission. Human studies give first insight into the scope of intergenerational influence of stressors but are limited in exploring the underlying mechanisms. Animal models have elucidated the mechanistic underpinnings of epigenetic transmission. In this review, we summarize progress on the mechanisms of paternal intergenerational transmission by means of sperm RNA in mouse models. We discuss relevant details for the modelling of RNA-mediated transmission, point towards currently unanswered questions and propose experimental considerations for tackling these questions.

Comparison of Reptilian Genomes Reveals Deletions Associated with the Natural Loss of γδ T Cells in Squamates

T lymphocytes or T cells are key components of the vertebrate response to pathogens and cancer. There are two T cell classes based on their TCRs, αβ T cells and γδ T cells, and each plays a critical role in immune responses. The squamate reptiles may be unique among the vertebrate lineages by lacking an entire class of T cells, the γδ T cells. In this study, we investigated the basis of the loss of the γδ T cells in squamates. The genome and transcriptome of a sleepy lizard, the skink Tiliqua rugosa, were compared with those of tuatara, Sphenodon punctatus, the last living member of the Rhynchocephalian reptiles. We demonstrate that the lack of TCRγ and TCRδ transcripts in the skink are due to large deletions in the T. rugosa genome. We also show that tuataras are on a growing list of species, including sharks, frogs, birds, alligators, and platypus, that can use an atypical TCRδ that appears to be a chimera of a TCR chain with an Ab-like Ag-binding domain. Tuatara represents the nearest living relative to squamates that retain γδ T cells. The loss of γδTCR in the skink is due to genomic deletions that appear to be conserved in other squamates. The genes encoding the αβTCR chains in the skink do not appear to have increased in complexity to compensate for the loss of γδ T cells.